perkin-elmer · glass bulk parameters. the agreement between theory and experiment is greatly...

50

Upload: others

Post on 23-Jan-2020

5 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for
Page 2: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER

I I 1 1 1 1 1 I

THE PERKIN-ELMER CORPORATION

REPORT NO. 8995

ACOUSTIC SELF-TRAPPING OF LASER BEAMS

First Semiannual Technical Summary

Through 31 October 1967

Submitted by

Edwin L. Kerr, Project Scientist (Phone 203-76 2-4650)

to the

OFFICE OF NAVAL RESEARCH

ARPA Order Nunber 306 Project Code Numbe Contract Number N00014-67-C-0468

Contract Duration 1 May 1967 through 29 February 1968

"Reproduction in Whole or in Part of this Report-1»- Permitted for any Purpcse of the United States Government."

This Research is part of Project DEF2NDER under the joint sponsorship of the Advanced Research Projects Agency, the Office of Naval Research, and the Depart- ment of Defense.

19$7

A

-

Page 3: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-E'LMER Report No. 8995

TABLE OF CONTENTS

faction

I

II

III

IV

VI

REFERENCES

Titlfl Page

ABSTRACT iv

PREFACE 1

INTRODUCTION 2

2.1 A New Power Limit 2.2 Threshold PreJlctable 2.3 Reconcilable with Steady State

Trapping Theory 2.4 Scaling Laws 2.5 Typical Results for Glass

2 2 3

3 5

DERIVATION OF THE ACOUSTIC BEAM TRAPPING THRESHOLD

3.1 The Driving Force 3.2 The Photoslastlc Effect 3.3 The Lcser Beam 3.4 The Acoustic Wave 3.5 The Numerical Integration 3.6 The Trapping Condition 3.'/ Th« Trapping Threshold 3.8 Numerical Calculations

AN ANALYTICAL SOLUTION OF THE LASER FEAM RADIUS EQUATION FOR KERR EFFECT TRAPPING

4.1 The Beam Radius Equation 4.2 Beam Propagation Without Self-Focu&lng 4.3 The Laser Beam Intensity Equation 4.4 Variation of Refractive Index with

Intensity

BEAM TRACING IN TF» ACOUSTIC TRAP

5.1 Other Features of the Glass Damage Phenomenon

5.2 Definition of the Problem 5.3 The Coupled Equations 5.4 The Step-By-Step Computer Solution

SUMMARY OF RESULTS AND FUTURE PLANS

11

8 8 9 9 14 16 16 20

22

22 24 25 27

32

32

34 34 35

5.5 Typical Results for Acoustic Beam Trapping 37

41

43

Page 4: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-£LMER Report No. SJ95

LIST OF ILLUSTRATIONS

i

I I I

Figure

1

2

3

4

5

6

7

8

9

10

11

1 12

■ 13

1 14

1 15

1 16

I 17

18

Title Page

Scaling Laws for Kerr Effect Trapping, Breakdown, 4 and Acoustic Trapping

Laser Damage Thresholds for Thr<se Optical Glasses 7

Laser Pulse Temporal Intensity Distribution, and the 10 Two-Impulse Approximation

Gaussian Spatial Intensity Distribution 11

Initial Force Distribution 11

Initial Rate of Change of Ccmpretision Distribution 11

Hankel Transfc m of Initial Ret« of Change of Compression 11

Acoustic Compression at T = 0.596. Proportional to Re- 11 fractive Index Change

Beam Center Compression as a Function of Time IS

Sound Wave Compression at Various Times After An Initial 17 Impulse

Sound Wave Compression at Various Times Near the First 18 Peak

Equivalent Beam Guides 23

Unit Quartic Pulse 26

Kerr Effect Trapping at 99 Percent Threshold 29

Kerr Effect Trapping at 95 Percent Threshold 30

Kerr Effect Trapping at 85 Percent Threshold 31

Geometry of the Beam Trapping Analysis 36

Acoustic Beam Trapping at 300 Percent Threshold 38

ill

Page 5: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

I I I I i I

ABSTRACT

Internal filamentary glass damage caused by high power Q-switched

pulse lasers has been reported in the literature. Also, filamentary trapping

in liquids has been observed. This report presents a theory of acoustic beam

trapping. In the theory, electrostriction is the sound wave driving force and

sound wave compressions cause focusing of the light wave fields. The theo-

retical trapping thresholds are calculated from the laser wavelength and the

density, refractive index. Young's modulus, and Poisson's raLi.o of a solid

material, or the density, refractive index, and speed of sound of a liquid

medium. The acoustic trapping thresholds agree with experimental glass

damage thresholds to within experimental accuracy, and they vary the same way

with initial beam size. Trapping movies showing the formation of strongly

focused regions are presented. An explanation is given for most of the salient

features observed in the damage phenomenon.

Ifl

I

An analysis of Kerr effect trapping is also presented for purposes

of comparison.

Iv

BBHjKSiSäSE^-iÄf«^ MHUJMIUWJ

Page 6: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

I 1

SECTION I

PREFACE

The goal of the research reported here is to continue the develop-

ment of a theory of acoustic be?m trapping. The theory is a possible explana-

tion of the mechanism for internal, filamentary glass damage by lasers. It

may also be an explanation for some of the beam-trapping ph.nomena observed in

liquids when they are traversed by a high intensity laser pulse.

We began developing the theory in 1965. Sections II and III pre-

rmt theoretical results obtained by August 19b6. These results were origi-

nally presented in the proposal for this contract. Since that time, some im-

provements have been incorporated in Sections II and III. The trapping 2

threshold was reduced by a factor of n /4 from the formula in the proposal.

Also, new information from the manufacturer of ouc glass samples has been ac-

quired. The graph comparing the theoretical and experimental damage thresholds

for three optical glasses is now based on the new, more accurate values of the

glass bulk parameters. The agreement between theory and experiment is greatly

improved.

Section IV covers an analytical solution to the beam-trapping equa-

tions for the case in which Kerr effect trapping is the dominant trapping mach-

anism and where electrostrictively driven sound waves are weak or absent.

Section V presents some of our computer results, depicting the

acoustic beam-trapping phenomenon. An explanation of the salient features of

the glass damage phenomenon is given.

Section VI presents a brief summary of our work and our plans for

the second half of the contract.

Page 7: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

-: : C

I I I I I [

i I

PERKIN-ELMER Report No. 8995

SECTION II

INTRODUCTION

2.1 A NEW POWER LIMIT

As higher and 'Lgher power lasers are developed, basic linltations

in the power-transmitting capability of materials and propaeation media are be-

ing discovered. Examples of limitations are electric breakdown and beam insta-

bility. During our study of gain saturation and other anomalies in stimulated

Raman effect and in our experimental work in laser damage to glass, we have

identified a new kind of beam instability, acoustic beam trapping.

Acoustic beam trapping is caused by the focusing action of electro-

strictively driven acoustic waves. These waves alter the index of refraction

of the medium by the density changes they cause. All optical materials experi-

ence these electrostriction forces and photoelastic effects.

Acoustic beam trapping sets upper limits to beam power which can be

transmitted in important materials such as glass and air. Furthermore, the

scaling laws for acoustic trapping (power threshold versus beam size and pulse

duration time) differ from the scaling laws for other types of beam Instability,

such as Ken effect or anomalous dispersion trapping. In fact, for many trans-

parent materials, there is a beam size and pulse duration for which the acoustic

trapping threshold is lower than the thresholds fjr other known instabilities.

2.2 THRESHOLD PREDICTABLE

At the present time, we can predict the acoustic trapping power

threshold in glassy materials from a knowledge of material properties, wave-

length, pulse duration time, and beam size. There is a critical power level

for each material which can be calculated from tabulated material properties

and the laser wavelength.

Page 8: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

it »*

i

rii

i

i

I

PERKIN-ELMER Report No. 8995

2.3 RECONCIT.ABLE WITH STEADY STATE TRAPPING THEORY

We have also shown that for any given beam size and type of material,

1 li there is a maximum average rate of increase of power which can be transmitted

without causing acoustic beam trapping. So long as pownr is added to a beam at

a slower rate than this, it is possible in principle to reach the steady state

trapping threshold predicted by Chiao, Garmire, and lovnes (Ref. 1).

2.4 SCALING LAWS

Figure 1 illustrates tha scaling laws for three different effects

which limit the ability of an optical material to transmit a laser beam. In

the graph, total beam power Is plotted versus beam radius on log-log scales.

The three effects are electric field breakdown, Kerr effect trapping, and acous-

tic trapping.

For electric field breakdown, as in laser-Induced "air sparks,"

I i Power Thresho for Breakdown

2 Id f Constant for A r Beam ^ n V, Material • vÄadlus/

On the log-log plot, tais poyer thresnold Is a line with a slope of 2. Thus,

electric field breakdown limits the material to transmission of power levels

and beam sizes in the right-hand Portion of the graph.

For Kerr effect trapping, as shown by Chlao, Garmlre, and Townes

J CRef. 1)

Power Threshold for .„ , .. N /'Constant for "N Kerr Effec. Trapping = (Wavelength) ^ Material )

Since the threshold is independent ot beam size. It appears as a horizontal line

on the graph. For transmission without Kerr effect trapping, the beam power

level must be below the Kerr effect threshold.

U Together, Kerr effect trapping and electric field breakdown limit

beam transmission to beam sizes and power levels In the lower right portion of

11 the graph. Both of the effects are virtually Independent of the pulse duration

time.

Page 9: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

I

I I I I I

107p

106p

Kerr Effect

ID-

10

I o Ml 0) n

Acoustic Trapping

103p

2 10p

10p

Critical

103^t Power

Electric Field Breakdovin

lOw

Beam Size

10^w 103w io\ 105w

Figure 1 . Scaling Laws for Kerr Effect Trapping, Breakdown, and Acoustic Trapping.

Page 10: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKHM-ELMER Report No. 8995

The acoustic trapping threshold can be calculated as in Section III

as long as the laser pulse duration time is less than the time required for

sound to cross the beam radius. The scaling luw for acoustic trappiug is

9 Pover Threshold for 2 / Constant for "> ^Beam Radius) ' Acoustic Trapping = (Wavelength) ^ Material • (Pulse'Duration)^"

for the limiting case of a very short pulse or a very large beam. Tf the pulse

duration is held constant, the threshold power can ba plotted versus beam size

as shown in the graph. The acoustic trapping power threshold curve flattens

out at the bcttora, at the critical power level. This critical po^er level,

P , is characteristic only of the material and the wavelength.

For the domain in which the pulse duration is longer than the time

T require,-1, for sound to cross the beam, the acoustic trapping threshold remains

at P . This power level can be exceeded only if ths power does not increase

more rapidly than P /T in any time T.

If the pulse duration is increased by a numerical factor, and the

beam radius is increased by the same factor, the power threshold remains con-

scant. Thus, in the leg-log graph, the curve Is shifted to the right by the

log of the numerical factor by which thfe beam radius is increased. Thus the

graph shows three extra curves shifted by factors of 10,100, and 1000.

t In some materials the Kerr effect trapping threshold is lower than

• the critical power level for acoustic trapping. In these materials Kerr effect

Strapping will occur before acoustic trapping. However, most common optical

materials have a critical power level lower than the Kerr affect trapping thresh-

_ old. For these materials there- will always be a domain of beam size and pulse

p duration time in which acoustic beam trapping sets the maximum transroittable

power for the materialt

^ 2.5 TYPICAL RESULTS FOR GLASS

Acoustic trapping is an important cause of laser damage to optical

glass, as shown In experiments perfomwd by Steinberg, Atwood, Lea, and Ward

n(Ref. 2). The theoretical trapping threshold is compared with the experimental

I

II i «

Page 11: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIIM-ELMER Report No. 8995

damage tareshold in Figure 2. Generally, the trapping threshold is below the

experimental damage threshold points, as expected. In particular, note the

experimental damage threshold curve for dense flint. The left, lower portion

has a slope of 2. In that region, the cause of damage is orobably electric

breakdown. The right, upper portion fits the curve for acoustic trappin to

within th.? experimental repeatability, and it scales the same way. For fused

silica and BK-7 the agreement between the accastic trapping threshold and the

measured damage threshold is even better.

The experimental results are not attributable to Kerr effect trap-

ping because of the dependence on beam size. Also, in glass, Kerr effect is

so weak that the power threshold for Kerr effect trapping is /, negawatts, which

is above the to? of the graph.

Page 12: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Rer^rt No. 8995

2000

I I I I I I I I I I I I

T 1 i I I I I

1000

800 -

')00

400

300

200

100

80

60

^ 40 .0

5 30 •H

u 20 0)

5

3 fi 10

!/

/ —• O •— Fused Silica

-■ D--B— BK-7

—A A A— SFS-5 Dense Fl^it

u o 0)

/

(U so

•3 E u vt <u a o. x o

Cd Z

01 60

e Q

J LI I I I I X J I I I I I I 20 30 40 60 80 100 200 300 400

Baam Radius w in Microns

600 800 1000

Figure 2 • Laser Damage Thresholds for Three Optical Glasses

^^-. 5 i-

I

Page 13: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER

SECTION III

DERIVATION OF THE ACOUSTIC BEAM TRAP¥TNG THRESHOLD

Report No. 8995

I 3.1 THE DRIVING FORCE

Electrostriction is the force exerted by an electric field on a

material uedium, when the force is proporcional to the square of the field.

The net body force is th<?.i proportional to the power intensity gradient, in

a lossless medii.m. The -elevant permittivity ratio at optical frequencies is

the square of the refractive index. Thus, the net body force f per unit

volume due to the light beam is (Ref. 3)

I I I i

f = (:/6n c) (n2 + 2) (n2 - 1) 7 I. 00 o

Thus a cylindrical beam will drive a radially proprgating -.ound wave.

3.2 THE PHOTOELASTIC EFFECT

The refractive index change ^n due co small acoustic compressions

is defined as

An = n - n = (p - po) »n/dp

Assuming constant polarizability per molecule, ^n/ftp may be calculated by

differentiating the Clausius-Mosotti relation. Then we obtain

Ln = (1/6 no)(n2 + 2) (n2 - 1) a

where a Is the normalized compression,

cr = (p - p0) /p0

Page 14: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKINELMER Report No. 8995

1

Other effects, such as the Kerr effect, may add to Ar1. They will

not have the same distribttion as acoustic compressi-in, in general, nor will

they vary the same way with beam size. However, they can be added to , n later.

Let us ignore them for the present.

3.3 THE LASER BEAM

For the purpose of this discussion, the laser pulse is represented

by two impulses of equal energy. Each impulse has a Gaussian radial intensity

distribution. The radius is measured to the point where the intensity drops

to 1/e of the peak intensity. The beam is circularly symmetrical and gently

focused. At the focus the intensity distribution is

I(r,t) = (W/TT W2) exp (-r2/w2) 6 (t ± /it/2)

where W is the energy of one impulse. The two impulses are separated by a

time, At, equal to 2/3 the half-height duration time of a laser pulse. Thus

if the physical Lassr pulse has a Gaussian time distribution, the two impulses

will occur at the centroids of the two halves of the physical pulse.

The temporal intensity distribution is shown in Figure 3, and the

spatial intensity dlatribution is sVetched in Figure 4, The resulting force

distribution is given in Figure 5.

3.4 THE ACOUSTIC WAVE

The electrostrictive force drives -y radially propagating acoustic

wave. Usually when acoustic beam trapping occurs the boundaries of the mate-

rial are so far from the beam center that the trapping event occurs sooner

than sound can be reflected from the boundaries. Hence, the acoustic wave

equation applies with the following conditions:

1) Circular symmetry

2) Infinite homogeneous Isotropie medium

3) Solution at the beam center is well behaved

4) Solution at infinite radius is zero

3) Solution is not a function of z (paraxial beam case)

Page 15: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

i

I

PERKIN-ELMER Report No. 8995

I

[ I

I

Gaussian Temporal Intensity Distribution

Two-Impulse Approximation

fc

Half-height Duration

t Figure 3. Laser Pulse Temporal Intensity Distribution,

and the Two-Impulse Approximation.

10

Page 16: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

p

^

Figure 4.

Gaussian Spatial Intensity Distribution

Figure 5.

Initial Force Distribution

Figure 6.

Initial Rate of Change of Compression Distribution

Figure y,

Hankel Transform of Initial Rate of Change of Compression

Figure 8.

Acoustic Compression at T = 0.596. Proportional to Refractive Index Change

11

-ifnm

Page 17: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

I PERKIN-ELMER Report No. 8995

I I

The acoustic wave equation can be solved for the case of an impulse

driving force. Since the equation is linear, it is convenient to use dimen-

sionless variables and a normalized impulse.

Let the dimensionless radial coordinate be x, the radius measu/ed

in units of the characteristic beam radius w;

x = r/w

There is a characteristic acoustic response time, T, equal to the time re-

quired for sound to traverse the beam radius;

T = w/v

where v is the speed of sound for a two-dimensional compression wave.

v = (Y/2po)1/2 (1 - f)"

1/4 (1)

Here Y is Young's modulus and c Is Poisson's ratio for solid media.

Let the dimensionless time variable be T, time measured in units of T.

T = t/T = (v/w) t

This choice of units makes the velocity of sound equal unity in ■

the wave equation. Also, the acoustic velocity ü is normalized by the velocity

ot sound, and the acoustic displacement u is normalized by w. The acoustic

compression a is the negative divergence of the displacement:

O = - 7 * U

w' -e U - u/w. Differentiating with respect to T, we have

ha/bT = - V • ÖU/ÖT

where dU/dT = u/v

i 12

Page 18: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

1 PERKIN-ELMER

Report No. 8995

i! i

i 1

I 1

| i

i

\ i

The initial rate of change of compression distribution appears in

Figure 6,

When the electrostrictive force acts impulsively on a medium ini-

tially at rest, there is no immediate displacement or compression. The ini-

tial velocity and initial rate of change of compression may be deduced from

conservation of momentum,

I f dt P u vo

The compression a is a function of radius and time only. To sum

up, the problem is to solve the acoustic wave equation:

2 V a

2 /^ 2 n (T/öT

subject to

In 111

ID

where

1)

2)

3)

4)

5)

a = (T(X,T)

a(x,o) = o 2 2

aa(x,0)/^T = A(l - x ) exp (-x )

3 2 2 A = (2W/3im cvpow ) (no + 2) (no - 1)

CT (0,T) ^OO

afoc, T) =0

Tb" equation and all of the conditions except (3) are satisfied

by a (x, T) = (B/A)J (yx) sin (yr) where B and y are arbitrary. A linear

superposition of solutions will satisfy condition (3) as well. The correct

combination of solutions g(y) is given by the Hankel transform (not the

Fourier transform because of circular symmetry). The weight function if x

13

Page 19: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

f {

PERKIN-ELMER Report No. 8995

and the eigenfunctions are J (yx) corresponding to 1 and sin (xy) or cos (yx)

for a Fourier transform. Thus,

OB

8(y) = j £ [A(l - £) exp (-^2) 1 J0(Cy)d^

--. *vy2/8) exp (-y2/4)

This distribution is plotted in Figure 7. Now superposition is

applied to obtain the solution.

Cr(x,T) « A f (y2/8) exp (-y2/4) J (yx) sin (yT) dy "o

This difficult integral has been evaluated numerically as described below.

For the important on-axis case, where x = 0, the integration may be performed

analytically, since J (0) = 1.

O(0,T) = A f (y2/8) exp (-y2/4) sin (yr) dy Jo

A T/2 + (A/2) (1 - 2T2) exp (-T2) f exp (^ ) d^ Jo

= A (T- 4T3/3 + - - -)

This function is plotted in Figure 9.

3.5 THE NUMERICAL INTEGRATION

The integral for o(x, T) was computed in our Scientific Computer

Facility. We programmed the Scientific ^ Systems 9300 Computer with 65

Fortran IV statements. Simpson's method WdS used. The Bessel function was "8

generated by a polynomial approximation accurate to 5 x 10 absolute error

from the Handbook of Mathematical Functions (Ref. 4), One hundred and one

values of y were used each time the integral was evaluated. The integral was

computed at. 101 values of x for each of 124 values of r. Time was saved by

storing parts of the kernel which do not change, and by buffering the printer.

Running time was about 40 minutes, or about 1 millisecond for each evaluation

of the Bessel function.

14

Page 20: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

i F

1 =

l

l

01 B

o

s § m 10 0) u

u

s §

I ON

0)

F>4

v/(.t'o)^

15

__—^__ i

Page 21: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. £995

i.

The computed points were saved on magnetic tape. We displayed the

points with the on-line oscilloscope associattd with our SDS 930 Computer.

The Fortran II display program allowed us to show stationary frames continu-

ously, or to show frames sequentially at any comfortable viewing rat»?. The

display was useful in studying the wave motion and ir checking the computation.

Refer to Figures 10 and 11 for plots oif wave amplitude o/A versus

x at various values of T. The compression at maximum on-axis amplitude is

plotted in Figure 8, for comparison with the other functions in the solution.

3.6 THE TRAPPING CONDITION

An unexpected result of the computed solution of the acoustic wave

is the fact that the compression is greatest on axis, for 0 < T ' 0.85. Thus,

it is, from the start, a focusing distribution. There is no latency period of

zero or ne ative focusing before somo positive focusing begins. The on-axis

solution shown in Figure 9 shows that the convergence (reciprocal focal length)

varies linearly with r, for T « 1.

The two-impulse model is valid for T < 0.85. The first impulse

starts an acouscic wave, and the acoustic wave builds up a distributed lens

cf increasing convergence. Diffraction of the second impulse will be defeated

if, at the time the second impulse occurs,

An/n > (l/2)(\/27rnow)2

This on-axis An will be the sum of An caused by the acoustic wave and Kerr

effect. The acoustic An is a function of beam energy and time duration, while

Kerr effect depends on intensity. Thus the two effects do not vary the same

way with beam size.

3.7 THE TRAPPING THRESHOLD

For optical glass and other materials in which Kerr effect is weak,

we may calculate the approximate energy threshold for trapping by equating the

16

Page 22: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

I Id

n I i1

u

...

<

« 4)

■^•^^«•mmfiww

Radial Distance

T = 0.1T/21/2

T - 0.3T/2 1/2

3.54w

T - 0.2T/2 1/2

T • 0.4T/2 1/2

I 1 i

"

0

11

,.1/2 T - 0.517

Figure 10.

T « 0.6T/2

Sound Wave Compression at Various Times After An Initial Impulse

17

1/2

^

Page 23: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

s

-i-.

H

I r

r

i i

i.

PERKIN-ELMER

T = 0.85 T/2 1/2

Report No. 8995

ITS«"" •-. -••-"

T = 0.8T/2 1/2

■*r

T = 0.9T/2 1/2

Figure II. Sound Wave Compression at Various Times Near the First Peak

L 18

Page 24: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN ELMER Report No. 8995

required axial n with the . ,n due to the acoustic wave. The total beam

energy W ^ 2W. Hence

TRAP 4:i o ^o

2 ..27 2 ,? a/A (no . 2) (no - 1

-

For T -^ !. a/A =- T = tX v/w

Thus, for pulses that are short compared with the beam radius transit time,

W, TRAP 4it

2 no 9c\ T _Q o.

' 2.2 _2 ' ^t (n" • 2) (n^ - 1)

we may define the power for the impulse pair as

p = "TOTAL7-1-

Then for shori; duration pulses (iT^l) we have

TRAP 4 it 9c\ r o ^o

1 2 _2 2 ,.2

., 2 n — 2

(n + 2) (n - 1) " At o o

; 1 =

'

For T - 2'1/Z, w2/At2 = v2/2

Thus, when the pulse duration is matched to the beam transit time, PTBAp =

P,RTT1 . . a power level charactenstic only ox the material (and the wavelength).

The quantity,

P /T CRITICAL'

is the maximum average rate at which power may be increased in a material over

a time, T, without acoustic trap formation.

1 19

Page 25: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

w 3.8 NUMERICAL CALCULATIONS

r The acoustic trapping threshold curve for a material can be ob-

** tained with simple calculations. It is only necessary to obtain the sound

j- velocity, the beam size that matches the pulse duration timp, and the critical

H power level. Then simple graphical methods yield the threshold curve.

| Sound Velocity; The relevant sound velocity is that for a radially

propagating compressional wave. For solids, use equation (1) in Section 3.4 i

f| The Young's mo ulus, density, and Poisson's ratio for glasses are found in the

"* newer Schott glass catalogs.

| Example: "or BK-7, Y = 8310 kp/mm , p = 2.51 g/cm3, f = 0.208.

Multiply p by 10 to obtain the density in kilograms per cubic meter, and Y

j by 9.81 x 10 to obtain newtons per square meter from kiloponds per square li 3 millimeter. Hence v = 4.27 x 10 meters per second.

jj Matched Beam Size: Let the pulse duration time p be measured at

the half-peak points. Then w . . = pv. r* r r matched t

I

I I

Critical Power Level: The critical power level for the material

is

^ritical^81 C>^8«> PVV f(no+ 2> V " l>2 1

Example: For BK-7, at \ = 694.3 nm, n = 1.45. The critical

power level is 44C kilowacts.

Plotting the Threshold: Use log-log graph paper having equal siz;

divisions for the horizontal and vertical scales. Draw a horizontal line

corresponding to the critical power level. Mark the matched point on the

critical power line, corresponding to the matched beam size.

For a quick, conservative overestimate of uhe trapping threshold

I curve, simply draw a line with a slope of 2 upward from the matched point.

(The angle is 63.5' to the horizontal.)

20

Page 26: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

..

a. i i'i i ii

PERKIN ELMER Report No. 8995

i

For an accurate graph, plot the curve in Table I on a second

piece of log-log paper. Lay the first pi^ce over the second so the axes are

parallel and the natched point is just over the bottom point on the second

graph. Trace the curve.

TABLE I

NUMERICAL VALUES USEFUL FOR PLOTTING ACOUSTIC TRAPPING THRESHOLD CURVES

-: i

I • Horizontal Axis Vertical Axis

l ■ h

'I . .

I

j

1 n i if

D

0.91

1.00

1.10

1.25

1.43

1.68

2.00

2.50

3.33

5.00

10.0

0.965

1.00

1.05

1.20

1.43

1.78

2.36

3.48

5.90

12.8

55.0

::

21

Page 27: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

o ID

D L L L

[ L

t L

SECTION IV

AN ANALYTICAL SOLUTION OF THE LASER BEAM RADIUS EQUATION FOR KERR EFFECT TRAPPING

A high intensity laser beam passing through a material medium can

focus itself into a long, thin filament and propagate without normal diffrac-

tion spreading. This self-trapping phenomenon arises when the medium's re-

fractive index n is higher along the beam axis than along the beam edges.

Such refractive index distribution acts like a series of thin positive lenses,

as shown in Figure 12,

The high intensity beam can set up a focusing refractive index

distribution by everal physical mechanisms, such as electrostriction, anom-

alous dispersion, and reorientation of molecular dipole moments. The last

effect is Kerr effect, after its discoverer, John Kerr (1824-1907). In

liquids the molecular reorientation can occur in times on the order oi 10

picoseconds. The effect is thus virtually instantaneous compared with the

duration of nanosecond laser pulses, but it is much too slow to follow 500

t rahertvs light wave fields. The local change in index of refraction is

proportional to the local beam intensity. When the beam power is above a

certain threshold power level, the laser beam focuses itself to a smaller

beam size and higher intensity. The smaller, higher intensity beam causes

still stronger self-focusing until the beam has trapped itself at a small

radius, limited by diffraction.

4.1 THE BEAM RADIUS EQUATION

This note sets forth an analysis of Kerr effect trapping. The

analysis is based on a ray tracing equation published by Tien, Gordon, and

Whinnery (Ref. 6). Equation (9) in that article reduces to

2 2-3 2 2 2 3 a/7\z = a (\/?iTno) + (a/n^ n(0,z,t)/^r (2)

t 22

Page 28: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

uoiinqiaqsTa xapui aAijDHJja^

= I i

:-

l Ü

a» tn C

at u c 0) 3 cr W3

c o

a o M

<U

U o

>4-( c o u

-i X! T J-l u O

•,■■ XI *J a. c O >-(

c 0 QJ

•rH J= ■ 4-1 *-< V) 3 01

XI 00 u •H C 3 ^1 •H on

■U 4J --4 in <0 u-i

•H M Q 60 x;

(U u >: u o u s x>

•u ■H c «4-1

c-i 5^ O XI

V tn t C 9)

•rJ ai oo *J (U TJ U (O O) to u <u u

U-l X) x: u 4J OS >.

<T) 00 • ^-* e c >4-i « O >W

•D ai -H o •H U 05 0 B C

^^ 4) T3 n> o -t C J=

XI in nj -u «0 > u ■H tn to nj 3 U -< At o- >; x

U Q to (0

U 41 C so £ £ O c H *J

•r-l e > 00 o « . C -H X to 0 i-i

4) ^H 10 E T» <0 T3 3 ■H U •-I ? /-s «0 T3 Ci N " 0) •a « S e M

to c !0 Q) S-/ Q)

so u c JS Q •rt ui VJ E

c oo c 0) C 01 o _4 01 -H •^ CO 1-1 ■u > <u (d ■M x: u 00 3 u 4) to cr to j3 tt. ii3 a. u o u a. • S a r-l

« P3 ti

u 3 00

Ui

23

Page 29: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

E I I I I

[

PERKIN-ELMER Report No. 8995

where

a(z,t) - beam radius to the point where the intensity is

1/e of the peak intensity,

n(r,z,t) = local index of refraction

n = undisturbed or nominal index of the medium o \ = vacuum wavelength of the laser beam,

r,i5 = cylindrical coordinates for the beam.

For this equation to be valid Tien. Gordon, and Whinnery require

that

1) "the light beam is...launched v;ith a Laguerre-Gaussian

or Hermite-Gaussian field distribution"

2) "the refractive index of the medium varies slowly in space

(negligibly in an optical wavelength)"

3) the variations in n are small compared with n

The last two restrictions are violated when the beam is trapped

to a filament with approximately one wavelength "adius, and tht scattered

light is spread over a wide wavelength range. However, the equation can be

used to study the collapse of the beam toward the trapped condition, anJ to

study the phenomenon when the beam power is below threshold.

U.L BEAM PROPAGATION WITHOUT SELF-FOCUSING

When the beam intensity is weak, the last term in equation (2)

is negligible. Thy ray paths are then hyperbolas. The beam radius solution

is

9 7, 1/2 a(z) - (a + a z) ^ (\z/2jm a )

1 o o J (3)

where a and a. are, respectively, the initial radius and initial slope of the o i r

rays just after entering the medium at z = 0. For example, let \ = l^m,

n = 1.5, a = O.lmm, a, = -0.001839. Then the beam is focused to a radius of o o 1

0.05mm at z = 4.1cm. This is the case of a very gently focused beam. The

24

Page 30: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKINELMER Repent No. 8995

first frame of Figure 14 shows a plot of a(z) for this case . Note that the

vertical, radius scale is greatly exaggerated.

4.3 THE LASER BEAM INTENSITY EQUATION

The lowest order mode for the laser beam has a nearly Gaussian

intensity profile characterized by the radius, a. to the point where the

intensity drops to 1/e of the peak intensity. If the propagation medium is

not too violently inhemogeneous and if the inhomogenelty is radially sym-

metrical, the Gaussian profile is maintained along the beam, although the

beam radius varies ier^nse of diffraction and focusing by the inhomogene!ties.

I at

The laser power may also vary as a function of time. For Q-switched

pulse lasers generally the pulse energy can be measured and some rough idea of

the time distribution of the pulse can be obtained. For purposes of analysis

we can use a simple pulse shape such as the unit quartic pulse shown in

Figure 13. This shape has a continuous derivative and finite extent.

Thus the laser beam intensity is

2 2 2 I(r,z,t) = (W/ita^) exp (-r /a ) q (t;p)

11

I I i SI

where W is the pulse energy and q(t;p) Is the unit quartic pulse of Lime

constant p:

2 2 , , I lO/ifiPl ir:/pj (2-t/v) l(t;p) = f (15/16P) ^^ (^/P) ' 0 s t s 2p

(_0 o therwise.

'■

Of course, a is a(z,t) given by the solution of equation (2). The intensity

distribution is so normalized that

p 27t OD , »

I de r dr I dt I (r,z,t) = W o l o ,! -•

independent of the value of z.

25

Page 31: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

I I I I I t E

t I I I I 1 I I

PERKJN-ELMER

15/32 p

Report No. 8995

, , f (15/16p) (t/p)2(2-t/p)2 for 0 s t < 2p

q^ 'p; I 0 otherwise

J: q(t;p) dt = 1

Figure 13. Unit Quartic Pulse

26

Page 32: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIISHELMER Report No. 8995

11

4.4 VARIATION OF REFRACTIVE INDEX WITH INTENSITY

> The change in refractive index can be expressed in cgs units as

2 2 n = n + — \JE +

• where J ia the high frequency Kerr constant due to molecular rotation. Thus

— the change in index is proportional to the local intensity for either a plane

i polarized or circularly polarized beam (although the constants of proportion-

ality differ for the two cases). Let the constant of proportionality be ¥.,

depending on the wavelei.gth, the material, and the beam polarization.

n(r,z,t) = n + KI (r,z,t)

In equation (2) the last term becomes

f -3 ; -a J(2 YM/nnJ q (t;p)

f Since both terms on the right of equation (2) are proportional to

' a , the solution at any time will be a hyperbola whose vrlue and slope at

[r= 0 are a and a , respectively. However, the focal point and semiaxes of the

hyperbola may vary with time. The equation becomes

I

1 1

1 1 1

a2a/3z2 = f(t)a"3

where

f(t) = (\/2rmo)2 - 2KWq(t;p)/nno.

The «olution is

r 2 -2 2 T 1/2 a(z,t) = j (ao + a1z) + f(t) ao z j

A negative, zero or complex value of a would be a non-physical

. solution. This places a restriction on the minironn value of f(t) which occurs II 2 when t = p. We have a (z,p) equal to a quadratic expression in z. In order

2 that a (z,p) be greater than zero for all real values of z, the discriminant

2 of the quadratic expression must be negative. Then the only zeros of a will

27

Page 33: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

I E

PERKHSHEUS/IER Report No. 8995

be complex. This condition implies that f(t) must always be positive. The

minimum value of f(t) occurs when t-p, i.e., at the peak of the laser pulse.

For a valid solution we must have

1 2 An/nc < ^ (\/2itnoa)

where ^n, the on-axis increase in refractive index, is

2 £,n = KWq(p;p)/na

This condition is equivalent to the Pierce stability criterion L/f < 4. where

L and f are the spacing and focal length of lenses in a sequence (Ref. 5).

[

The threshold value of KW/p becomes (KW/p)

Let k be the fraction

k= (KW/p)/(KW/p)thre8hold

= 2\*/fln threshold

I

Then f(t) becomes

f(t) = (\/2*no)2 [l-kq(t;p)/q(p;p)]

Refer to Figures 14 through 16 for plots of the beam radius when k = 0.999,

0.95, and 0.85.

I I

28

Page 34: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. B995

\ = lum, n0 = 1.5, an = 100 \im, a1 -0.001839

12.5 cm

Figure 14, Kerr Effect Trapping at 99 Percent Threshold

I " - 29

Page 35: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

I I I r

i

i i i i

i

FERKIN-ELMER Report No. 8995

\ = Uun. a = 1.5, a = lOO^m, a. = -0.001839 Time _ ' o • o '1 0.0p '~

lOOuni

0.2p

0.4p

3

'S

§ 0) PQ

0,6p

O.Bp

1. o

12.5 cm

Figure 15. K.err Effect Trapping at 95 Percent Threshold

L 30

Page 36: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

Time \ - l^iva, n = 1.5, a = lOOmn, a O.Op 0 0 1

lOO^m

= - 0.001839

0.2p

0.4p

l.bp

10 3 •o eg ai

so 0.8p

l.Op

Axial Length_ 12. 5 ens

Figure 16. Kerr Effect Trappinc at 85 Percent Threshold

31

Page 37: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

SECTION V

I I D

L

BEAM TRACING IN THE ACOUSTIC TRAP

The theory, developed from the two-impulse model in Sections II

and Hi, yields a threshold for acoustic beam trapping which is below the ex-

perimental damage thresnolds for the glassy materials tested. It was solved

by completely analytical means, using tabulated properties of special functions.

The computer was only used to obtain values for the radial sound wave off the

axis of the beam-. That computation was not essential to the solution, although

it did give the following important insight: The refractive index dis;ribution

set up by the acoustic wave is initially a focusing distribution, although its

strength viries with time. However, the theory did not allow for the continuous

interaction of light and sound, nor did it show how the focusing action develops

with time.

5.1 OTHER FEATURES OF THE GLASS DAMAGE PHENOMENON

A complete theory should also explain the other observed features

of internal, filamentary glass damage, such as the following:

0 L

0

1) The Spectrum of the Side-Scattered Light

A white light flash is seen when the damage event occurs. o o

Ruby leser light may be shifted from 69^JA down to 4000A.

Such large shifts can only be explained by a very .=trong

dynamic refractive index change.

2) Starting Location of the Pamage Track

Even if the laser Is focused at the entrance face of the

glass samole, the damage track never begins at the entrance

face. There is always a short interval between the entrance

face of the sample and the start of the track.

32

Page 38: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

I I I i I I I I

11

PERKIN-ELMER Report No. 8995

3) Exit Face PlttinR

Usually the track ends in a pit on the exit face of the

glass sample. As a matter of fact, the threshold for the

exit face pitting appears to be a little lower than the

threshold for track formation. Nevertheless, the phenomena

is different from the type of surface damage reported by

others.

4) Location of Damaae Stars

The damage track is not always continuous; it may start and

stop several times in the sample. Often there is a damage

star on the upstream end of the track. These damage stars

are localized regions of gross fracture. Occasionally

they show discoloration, indicating possible chemical de-

composition. If the damage track extends upstream of the

damage star it usually only does so for a very shoit distance

compared to the extension downstream from the damage star.

The location of the damage star is intrigu.'ng. If the

damage starts in the damage star and then propagates down-

stream to form the track, why doesn't the damage star cast

a downstream shadow and prevent the beam from concentrating

in the thin filament?

5) Track Propagation Speed

Since the damige event occurs in nanoseconds, its dynamics

are difficult to follow. There is experimental evidence

that the event which forms the track propagates at u^out

10 times the speed of sound in the glass. The experimental

evidence does not show unequivocally that the propagation

is either upstream or downstream.

Most of these features are explained qualitatively in the theory of

electrostricttvely driven acoustic beam trapping.

33

« I i I ii

Page 39: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER eport No. 8995

5.2 DEFINITION OF THE PROBLEM

Our purpose is to understand how the acoustic beam trap forms, not

how the damage occurs. Therefore we only need consider small acoustic pressure

changes in the material. We will see how these lead to a beam instcbility when

the oower exceeds a certain threshold. The nature of the instability is such

that the beam rapidly focuses itself into a thin filament. The damage occurs

when the beam is in this trapped condition. The instability itself can be

studied by means of a geometric rr.' tracing equation which includes first order

diffraction effects. This beam tracing equation is equation (2) of Section IV.

The sound wave satisfies essentially the same conditions as those

given in Section III, except that it is driven continuously by the light wave,

rather than running inertially after an impulse. For a gently focused beam,

axial components of the sound wave field are negligible in comparison with

the radial components.

The desired solution will show a graph of the beam radius plotted

versus axial lengtli for each of many small time intervals during the laser

pulse. The time distribution of the laser pulse can be assumed to be the

quartic pulse distribution described in Section IV, since the actual details

of the shape of the laser pulse do not have too much bearing on the threshold

for damage and on the tire development of the trap.

5.3 THE COUPLED EQUATIONS

The sound wave equation is

2,2 -1 . 2.2 -22 .2

- (l/6vnj) )(n + 2)(n - 1) (WAOqU;?) 0*0 o o

-4,, 2.2. . 2,2. . 4 A (l-r /a ) exp (-r /a )

I I

The beam tracing equation is

a2a/3z2 = a'V^nn^2 + (a/6n2)(n2 + 2) (n* - 1) ft2a(o, z, t)/?r

2

34

Page 40: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

I I I

[ n i •

-.

PERKIN-ELMER Report No. 8995

The general solution of these coupled, nonlinear, partial dif-

ferential equations has been left as an exercise for the computer. One

simplification can be made immediately. The beam tracing equation involves

only the second derivative with respect to r of the on-axis sound wave solution.

Since ehe sound wave equation is linear we T»ay take its Hankel transform ana-

lytically. When z is held fixed, the Hankel transform of Q satisfies an

ordinary differential equation in t. This ordinary equation may be solved using

initial values of the Hankel transform, and of its derivative with respect to

time, and assuming a step-plus-ramp driving function. The solution is then

found by taking the inverse Hankel transform. Since the beam tracing equation

involves only the second derivative of the sound wave amplitude at r = 0, and

the inverse Hankel transform involves r only as part of the argument of the

zeroth-order Bessel function, we may perform the differentiation analytically

under the integral sign. Thus for each time step and each value of z, the

sound wave equation can be solved with only one numerical integration. The

use of a step-plus-ramp driving function allows us to take a much coarser

time step than we would require if we used an impulsively driven equation.

5.4 THE STEP-EY-STEP COMPUTER SOLUTION

The problem is divided up into sections and si icen as shown In

Figure 17. In its present state of development the computer program can hano^e

101 sections. The zeroth section is the entrance face of the solid oi liquid

medium.. The initial beam radius and entrance angle are fixed on the zeroth

section. Since the solution is radially symmetrical, it will be the same in

each radial slice.

The computer program must store values of a for each of the 101

sections. Also for each section, it must store 40 values of the Hanhel trans-

form of the sound wave, and 40 values of the derivative of the Hankel transform

of the sound wave with respect to time. An array of 101 points is also set

aside for the second derivative of the sound wave with respect to r. No

Bessel function storage or Bessel function subprogram is required. Initially

the Hankel transform of the sound wave and its derivative with respect to time

35

■^r:-.:.-.- CT^ --■ ■-

Page 41: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

I m >, t-i n d < 60 C

a, a <a U H E CD 01 M

(U J3

u

4) a

3 00

« « K •4 m •<

36

Page 42: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

are set to zero along with the second derivative of the sound wave with re-

spect to r. The solution then proceeds in steps:

i .

I : -

1) Solve the beam tracing equation using the Runge-Kutta method.

Tie axial step size In the solution is 1/5 of the Interval

between sections. Values of the second derivative of the

sound wave with respect to r are interpolated by Lagrange

cubic fitting.

2) Store 101 values of the beam radius a. Also save these

values on tape for later plotting.

i.

i 11

i U

f?

I.

3) If the time has not yet reached 2.Op, the end of the

laser pulse, advance the time by a snail time step and

obtain the new value of the beam power at that time.

4) For each section of the beam compute the second

derivative of the sound wave with respect to r by

using Simpson intergration. At the same time revise

the value of the Hankcl transform of the sound wave

and its time derivative for the next time step.

5) Go to 1.

The program parameters are the wavelength, the nominal refractive

index, the initial beam radius and the initial sW ^r "' . beam radius at

the entrance face, and the fraction of pulse power divided by theoretical

threshold.

5.5 TYPICAL RESULTS FOR ACOUSTIC BEAM TRAPHNG

Figure 18 depicts acoustic beam trapping when the pulse power is

3007, of the threshold, There are 11 frames showing the beam radius plotted

versus axial length tor 11 different times during the pulse. The first frame

shows the path of the beam when the illumination has just begun. The sixth

frame gives the trace when the pulse power is highest. The 11th frame shows

the path taken by the light in the trailing end of the pulse.

37

Page 43: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

r PERKIN-ELMER Report No. 8995

1 n

L

L L

Time P = 0-7 a/v> ^ = l^> n0 = 1-5^ a

0 = 100^; ^l

O.Op

0.2p

•0.001839

0.4p

0.6p

O.Bp

a i-op at as

n I.2p

1.4p

1.6p

1.8p

2.Op

lOOum

0 Axial Length 12.5 cm

Figure 18• Acoustic Beam Trapping at 300 Percent Threshold

38

Page 44: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Report No. 8995

| The action may be interpreted by considering that each section of

the beam forms an acoustic lens. Initially the lenses have zero power, but

i some of them build up focusing power faster than others. The strongest lens

is formed at the focus of the beam. As its focusing power builds up in time

it causes the beam to be refocused at some distance downstream. As the

strength of the acoustic lens continues to build up, its focal length is

shortened so the second focus moves upstream. The second focus likewise forms

** a third focus and so on. A little after the laser pulse is half over, three

« definite foci have appeared in the 12-1/2cm region plotted. At the same time

li that these acoustic lenses are becoming stronger because of the inertial

properties of the material, the laser pulse power is decreasing. The net

result is regions of sharper and sharper focus interspersed by regions in which

the beam radius is increasing. Eventually the beam will be able to escape

again to large radii. E .

I

I j

I s

F n i.

II:

Note in the last frame that the sharpest focus occurp the farthest

downstream. This is to be expected since the focus there is the result of

two strong lenses upstream. Another important feature is the motion of the

first focus. Notice how it moves upstream in the beam. This is because the

acoustic lenses in the first few sections of the sample are also developing

in strength with time.

Several features of the damage phenomenon are immediately ex-

plained. The fact that white light is scattered largäly from the region of

the damage stars is easily seen to follow from the fact that the highest

focusing and therefore the highest index of refraction occur at places which

become damage stars. Also the track propagates backward, becoming more and

more sharply defined until it ends in a damaged star. If there are more than

one, the damage stars occur in backward sequence. The farthest downstream

occurr first, then the next one upstream, then finally the damage star closest

to the upstream end of the beam. Thus, the damage event does not cast a shadow

that gets in its own way.

39

Page 45: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

I t [ I

PERKIN-ELMER Report No. 8995

The reason why it is difficult to measure the spee;' of propagation

of the damage event is now apparent. The intensity varies widely and more than

one focus region may move past the velocity-sensing optics.

The pit on the exit face of the glass sample is formed by axial

sound wave components which do become large when the beam is sharply focused.

There can never be a damage pit on the entrance face. The entrance conditions

of the beam are fixed; therefore, the bean cannot focus itself sharply at the

entrance face of the glass.

All of the qualitative features of the phenomenon are a sharp con-

trast to the phenomenon of Kerr effect trapping, in which there is one focus

region which moves downstream.

40

*■■*•■*- —^rl —: -

Page 46: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER i 1?

l n i

10

In

D

Report No. 8995

following:

SECTION VI

SUMMARY OF RESULTS AND FUTURE PLANS

During the first half of the contract, we have accomplished the

1) Kerr Effect Trapping Analyzed; An analytical solution was

obtained for Kerr effect trapping. Computer plots were

made showing the time development of this trapping phenomenon.

A threshold for the effect was obtained.

2) Two-Impulse Model Improved; Some corrections and improve-

ments were made in the earlier two-impulse model for the

acoustic beam trapping phenomenon.

3) Acoustic Trapping Equations Derived; Equations were derived

for acoi;stic beam trapping when the light pulse is treated

as having a finite time distiribution.

4) Driven Sound Wave Progranmed: Programs were written and

tested to show the development of the sound wave with a

time distributed forcing function.

5) Beam-Tracing Program Completed; Programs were written to

compute and display the acoustic beam trapping phenomenon

in which the sound-Light interaction is carried out.

6) Combination Program Started: A computer program to study

trapping phenomena when both Kerr effect and acoustic trap-

ping are operative has been written and partly tested.

41

Page 47: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

PERKIN-ELMER Rpport No. 3995

,-s

[

During the second half of the contract, we expect to complete

the development of the combination Kerr effect and acoustic trapping program.

Ve will prepare an article for publication to set forth the theory and dis-

play seve-rai cases of interest. If uine permits, we will also prepare a

compi".i.on article which sets forth the experimental results from our old glass

damage study contract in more detail. Another possibility will be to study

the reduced irapping thresholds expected if the material already has a focusing

index distribution in it due to thermal strains.

A second quarterly letter report and a final report will be

s-hmitted.

i; r L

f

I I L

42

Page 48: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

I 1 I I 1 1 1 I

PERKIN ELMER Report No. 8995

REFERENCES

R.Y. Chiao, E. Garmire, and C.l. Townes, "Self Trapping of Optical Beams," Phys. Rev. Letters 13. 15, (October 12, 1964), pp 479-482.

G.N. Steinberg, J.G, Atwood, P.H. Ld and S.A. Ward, Research into the Causes af Laser Damage to Optical Components; Final Report, the P rkin-Elmer Corporation TR-79A5, February 28, 1965. Prepared for US Army Electronics Command AMSEL-HL-L, Fort Monmouth. Contract #LA28043 AMC 00009.

3 Js'Miu. Adams Stratton, Electromagnetic Theory, McGraw-Hill, 1941, p.140.

4 Milton Abraraowitz an-' Irene A. Stegun, Handbook of Mathematical Functions with Formulas. Graphs, and Mathematical Tables, National Bureau of Standards, Appl. Mathematics Series No. 55, U.S. Government Printing Office, Washington, D.C.. June 1964, Section 9.4, pp 369f.

T

J.R. Pierce, Theory and Design of Electron Beams. D. Van Nostrand, 2nd ed., (1S54,, p. 195.

'P.K. Tien, J.P. Gordon, and J.R. Whinnery, "Focusing of a Light Beam of Gaussian Field Distribution in Continuous and Periodic Lens-Like Media," Proceedings of the IEEE, 53, No. 2, Febi ary 1965, pp 129-136.

Ü

43

Page 49: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

E

I I I I i I i

e

L

ö i.

Unclassified Security Claisificaiion

DOCUMENT CONTROL DATA • R&O (Seturny classificatwn of Hilf, body of abstract and mdrxinf annotation must be entered when the overall report is classified)

I. ORtCINATINC ACTIVITY (Corporate author)

The Perkln-Elmer Corporation Nor-'alk, Connecticut

l «(PORT TITLt

Acoustic Self-Trapping of Laser Beams

J<5. REPOUT SECURITY CLASSIFICATION

Unclassified TS anou>-

*. DESCRIPTIVE NOTEb (Type of •■upon and inclusive ilatei)

First Swülmnaal T«chnlcal Summary. Report 5. AUTHOWSJ (Last name, first name, inU ail

Kerr, Edwin L.

«- REPORT DATE

31 Ortohpr i%7 to. CONTRACT OR ORANT NO.

N00014-67-C-0468 b. OROJECT NO.

c. TASK

ARPA Order 306

7a TOTAL NO. OF PASES

45 TK NO.or REFS

6 Sa ORIGINATOR'S REPORT NUMBERTV

8995

»6. OTHER REPORT HOlSltAny other numbers thai may he assigned this report)

10. AVAILABILITY'^IMITATION NOTICES

"Qualified requestors may obtain copies of this report from Dnc"

H. SUPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY

Advanced Research Projects Agency, The Office of Naval Research, Department of Defense

IS. ABSTRACT

Internal filamentary glass damage jrused by high power Q-switched pulse lasers has been reported in the literature. Also, filamentary trap- ping in liquids has been observed. This report presents a theory of acoustic beam trapping. In the theory, electrostriction is the sound wave driving force and sound wave compressions cause focusing of the light wave fields. The theoretical trapping thresholds are calculated from the laser

wavelength and the density, refractive index. Young's modulus, and Poisson s ratio of e solid material, or the density, refractive index, and speed of sound o.* a liquid medium. The acoustic trapping thresholds agree with ex- periment*' I jlass damage thresholds to within experimental accuracy, and the) varv the same way with initial beam size. Trapping imv/ies showing the formation of strongly focused regions are presented. An explanation is given for most of the salient features observed in the damage phenomenon.

An analysis of Kerr effect trapping is also presented for purposes

of comparison. -^^

DO 'onm 1473 '"' I JAN M *

Uncla^alfled Security Clataification

44

Page 50: PERKIN-ELMER · glass bulk parameters. The agreement between theory and experiment is greatly improved. Section IV covers an analytical solution to the beam-trapping equa- tions for

Unclassified Security Classification

KEY WORDS

Acoustic Trapping Electrostrlctlon Kerr Effect Trapping Glass Damage Beam Trapping Sound-Light Interaction

BOLE

LlMK 6

»T

t

J I INSTRDCTIO;^

I. ORIGINATING ACTIVITY; Enter ihr name and address of the contractor, subcontractor, granlr Department of Defense activity ur other organization (corporate author) issuing the report,

2u. HKPORT SECURITY CLASSIFICATION. Enter the over- all security classification uf the report. Indicate whether "Restricted Data" is included. Marking is to be in accord- ance ^ith appropriate security regulations.

2b. CROUP: Automatic doMMradmg is specified in DoD Directive ^200.10 and Armed ForcesTn-Justrial Manual. Enter the group number. Also, when applicable, show that optional markings have been used for Croup 3 and Croup 4 as authorized. 3. REPORT TITLE; Enter the complete report title in all > apital letters. Titles in alt cases should be unclassified, it a meaningful title cannot lie selected without < lassifica- Hon, shov. title classification .n all capitals in parenthesis immediately billowing the title. 4. DESCRIPTIVE NOTES: II appropriate, enter the type of report, e.g., interim, progress, summr.ry, annual, or final, (>ive the inclusive dates when a sprcilic repotting period is covered. 5. AUTHORtS): Enter the namely) of aulhoris) as shown on or in the report. Enter last name, first name, middle initial, II military, show rank and branch of service. The name of the principal author is an absolute minimum requirement.

6. REPORT DATE; Enter the date of the report as day, month, year, or month, year. If more than one date appears on the report, use date of publication.

7a. TOTAL NUMBER OF PAGES: The total page count should follow normal pagination procedures, i.e., enter the number of pages containing information.

7ft. NUMBER OF REFERENCES: Enter the total number of references cited in the report, 8o. CONTRACT OR GRANT NUMBER: II appropriate, enter the applicable number of the contract or ««rant under which the report was written.

Hh, Be, ft »d. PROJECT NUMBER: Enter the appropriute military department identification, such as project number, subproject numbe,, system numbers, task number, etc.

9a. ORIGINATOR'S REPORT NUMBERfS): Enter the offi- cial report number by which the document will be identified and controlled by the originating activity. This number must be unique to this report, 96. OTHER REPORT NUMBER(S): If the report has been assigned any other report numbers (either by the originator or by the sponsor/, also enter this numberfsi.

10. AVAILABILITY LIMITATION NOTICES: Enter any limi- tations on further dissemination of the report, other than those imposed by security classification, using standard statements such as:

(1) l*0ualilied requesters may obtain copies of this report from DDC,"*

12) "Foreign announcement and dissemination i>f this report by DDC is not .luthorized,"

(3) "U* S* Government agencies mav obtain copies of this report directly from DDC. Other qualified DDC users shall request through

(4) "U. S. military agencies may obtain copies of this report directly from DDC. Other qualified users shall request through

(5) "All distribution of this report is controlled. Quali- fied DDC users •''all request throogh

II the report has been furnished to the Office ol Technical Services, Department of Commerce, tor s.ile to the public, indi- cate this tart and enter the price, if known.

11. SUPPLEMENTARY MOTES: Use for additional explana- tory notes.

12. SPONSORING MILITARY ACTIVITY: Enter the name of the departmental project office or laboratory sponsoring (pay- mg for) he research and development. Include address.

13. AB >TR AC'I. c'.nter an abstract giving a brief and factual summary of the document indicative of the report, even though it may also appear elsewhere in the body of the tech- nical report. If additional apace is required, a continuation sheet shall be attached.

It is highly desirable that th. abstract ol classified re- ports be unclassified. Each «aragraph of the abstract shai' end with an indication of the ;ii.;iti"y security claasilication of the information in the paragraph, represented aa (TS), (Sj, (C), or (V).

There is no limitation on the length of the abs*ract. How- ever, the suggested iragth is from 150 to L'2S words.

14. KEY WORDS: Key words are technically meaningful terms or short phrases that characterize a report and may be used as indei entries for cataloging the report. Key words must be selected so that no security classification is required. Identi- fiers, such ss equipment model designation, trade name, mili- tary project rode name, geographic tocation, may be used aa key words bjt will be lollawed by an indication of technical context. The assignment of links, rules, and wvigk.a is optional,

f f? ünclas sIfled Security CUssification

A5