Lab 3: measurement of Laser Gaussian Beam Profile Lab 3: basic experience working with laser (1) To create a beam expander for the Argon laser () To measure the spot size and profile of the Argon laser Measure before and after the beam expansion Do this by moving a knife edge through the beam Have a computer controlled knife that moves through beam
Knife Edge measurement of Gaussian Beam Consider a Gaussian shaped beam I( r ) r exp w P w r exp w I0 Where P = total power in the beam w = 1/e beam radius at point w(z) This is in cylindrical coordinates r is the radius of the central area
Knife Edge and Gaussian Straight knife edge cutting into a Gaussian shaped beam Measure the total power seen when knife move in x direction Must convert to Cartesian coordinates & integrate Assume - is when the knife fully below the beam x P x y I( x) dx dy w exp w exp w Where P is the total power of the beam I(x) is the intensity measured at position x In x direction the beam is cut: Integrate from x to - In y direction get full beam: integrate from - to + To solve this use the error function or integral of the normal erf ( x ) x e s ds Two ways of fitting this: Fit the power measured: that is the integral Fit the derivative 0
Fit the power measured That is the integral x 0 = centre of beam Then the power measured is given by for x >x 0 P x x 0 I( x) erf w For x<x 0 P x x0 I( x) 1 erf w Must also assume some background light level B In Excel use the Normdist (Normal distribution function) This is slightly different from erf function Fit with the excel function of the following formula I(x)=P*normdist(-x,-x 0,w/1.414,1)+B Where x is the position (starting with x below x 0 ) x 0 is the fitted centre point of the beam w = 1/e point you fit 1 is to make it the integration of the normal distribution B is the background or offset level Set up a spreadsheet with initial estimates of each parameter Have columns with z, I(x n ), fitted I(x n ), (fit-i(x)), error Set a column to sum the error (sum of the squares) Plot the difference (gives an estimate of w and x 0 ) P estimate is max power, B is background level Then use solver under tools to fit with minimizing sum of squares Use sum of squares as fit, others as variables Useful to plot the errors against position (called residuals) Should be on both sides of plots See sample excel layout in appendix B
Fit the difference of Power Measured That is the derivative z 0 = centre of beam Take a derivative of the measurements Best if take a simple derivative di0 dx I x 1 I x0 x x 1 0 Then the plot is Gaussian with the formula: x P x x di 0 exp dx w w Note need to be careful with the derivatives units you use Suggest you plot the difference but not fit it Plotting shows if the curve you are gettings Does it really look like a Gaussian check that See the plots in Appendix B
Appendix 1: Solver in excel Adding solver to excel For the lab 3 you will probably want to use the excel solver add on. Students who are using MS office 007 may find that the excel solver was not loaded into their excel and does not appear in the data menu. In earlier versions a pull down tab showed it still needed to be installed but 007 (and 010) does not. To get the solver. The instructions are hard to find in help also. Attached is a screen shot of the help instructions and the add-ins list showing where the solver is. Using Solver: Here are links to several good sites that give good examples of how to use the solver for a problem where you are adjusting several variables. http://chandoo.org/wp/011/05/11/using-solver-to-assign-item/ http://archives.math.utk.edu/ictcm/vol3/c006/paper.pdf http://web.chem.ucsb.edu/~laverman/chem116/pdf116cl/solver.p df
Appendix B Printout of excel sheet for laser profile fitting.
Nd-YAG Glenn first round: (this is used to roughly profile the laser beam and decide on the resolution for the z-axis needed) 150 91.7 5 84.6 300 70. 375 47.5 450 3.1 55 8.01 600.08 675 0.473 750 0.35 85 0.15 900 0.075 975 0.039 3050 0.015 without laser: 0.03 mw Z0 W/1.414 B P W second round: fitted 370.97 111.4439 0.067357 94.3335 157.5817 step I(z) diff I(z)fitted 150 93.4.4 9.16568 1.53558 9 175 91 3 90.69044 0.09589 8 00 88 3.4 88.50537 0.55396 5 84.6 4.1 85.4711 0.684115 7 50 80.5 5.1 81.3053 0.644056 6 75 75.4 5.4 76.0461 0.41744 300 70 7.5 69.67464 0.105861 5 35 6.5 7.7 6.3911 0.0903 4 350 54.8 8.3 54.745 0.7615 375 46.5 8.4 45.8741 0.39174 3 400 38.1 8.4 37.5417 0.31174 45 9.7 7.5 9.67945 0.0004 450. 6.6458 0.18067 1 475 16. 4.9 16.6036 0.16615 0 500 11.3 3.63 11.7156 0.17441 55 7.67.69 7.9415 0.073563 0 500 1000 1500 000 500 3000 3500 550 4.98 1.8 5.169711 0.03599 575 3.16 1.18 3.33894 0.00546 600 1.98 0.84 1.947875 0.00103 100 65 1.14 0.443 1.13596.1E-05 90 650 0.697 0.37 0.646959 0.00504 675 0.46 0.15 0.36787 0.008496 80 700 0.335 0.065 0.16076 0.014143 70 75 0.7 0.038 0.137608 0.01758 60 750 0.3 0.09 0.099017 0.017684 775 0.03 0.06 0.080966 0.01489 50 800 0.177 0.07 0.07935 0.01083 40 85 0.15 0.08 0.069536 0.006474 30 850 0.1 0.057 0.068169 0.00898 875 0.0963 0.03 0.067645 0.00081 0 900 0.074 0.019 0.067454 4.8E-05 10 95 0.055 0.0145 0.067388 0.000153 0 950 0.0405 0.0103 0.067367 0.0007 975 0.030 0.0073 0.06736 0.001381 0 500 1000 1500 000 500 3000 3500 3000 0.09 0.0047 0.067358 0.001977 305 0.018 0.018 0.067357 0.00416 sum sq 5.47055 Series1 data fitted