Articulo de referencia

haz gaussiano

Valor absoluto instantáneo de la parte real de la amplitud del campo eléctrico de un haz gaussiano TEM 00 , región focal. Mostrando | R mi ( mi ( t 1 ) ) | {\displaystyle |{\mat...

Valor absoluto instantáneo de la parte real de la amplitud del campo eléctrico de un haz gaussiano TEM 00 , región focal. Mostrando|Rmi(mi(t1))|{\displaystyle |{\mathcal {Re}}(E(t_{1}))|}de esta forma, con dos picos para cada frente de onda positivo .
Arriba: perfil de intensidad transversal de un haz gaussiano que se propaga fuera de la página. Curva azul: amplitud del campo eléctrico (o magnético) en función de la posición radial respecto al eje del haz. La curva negra representa la intensidad correspondiente.
Un haz de puntero láser verde de 5 mW, que muestra el perfil TEM 00.

En óptica , un haz gaussiano es un haz idealizado de radiación electromagnética cuya envolvente de amplitud en el plano transversal viene dada por una función gaussiana ; esto también implica un perfil de intensidad (irradiancia) gaussiano. Este modo gaussiano transversal fundamental (o TEM 00 ) describe la salida prevista de muchos láseres , ya que dicho haz diverge menos y puede enfocarse mejor que cualquier otro. Cuando un haz gaussiano se vuelve a enfocar mediante una lente ideal , se produce un nuevo haz gaussiano. Los perfiles de amplitud del campo eléctrico y magnético a lo largo de un haz gaussiano circular de una longitud de onda y polarización dadas están determinados por dos parámetros: la cintura w 0 , que es una medida del ancho del haz en su punto más estrecho, y la posición z relativa a la cintura. [ 1 ]

Dado que la función gaussiana tiene una extensión infinita, los haces gaussianos perfectos no existen en la naturaleza, y los bordes de cualquier haz de este tipo se verían recortados por cualquier lente o espejo finito. Sin embargo, la función gaussiana es una aproximación útil a un haz real en casos donde las lentes o los espejos en el haz son significativamente más grandes que el tamaño del punto w ( z ) del haz.

Fundamentalmente, la función gaussiana es una solución de la ecuación de Helmholtz paraxial , la ecuación de onda para un campo electromagnético. Si bien existen otras soluciones, las familias de soluciones gaussianas son útiles para problemas que involucran haces compactos.

Forma matemática

Las ecuaciones que se muestran a continuación suponen un haz con una sección transversal circular en todos los valores de z ; esto se puede observar al notar que aparece una única dimensión transversal, r . Los haces con secciones transversales elípticas , o con cinturas en diferentes posiciones en z para las dos dimensiones transversales ( haces astigmáticos ) también se pueden describir como haces gaussianos, pero con valores distintos de w 0 y de la ubicación z = 0 para las dos dimensiones transversales x e y .

Perfil de intensidad del haz gaussiano con w 0 = 2 λ .

El haz gaussiano es un modo electromagnético transversal (TEM) . [ 2 ] La expresión matemática para la amplitud del campo eléctrico es una solución de la ecuación de Helmholtz paraxial . [ 1 ] Suponiendo polarización en la dirección x y propagación en la dirección + z , el campo eléctrico en notación fasorial (compleja) viene dado por:

mi(r,z)=mi0incógnita^w0w(z)exp(r2w(z)2)exp(i(kz+kr22R(z)ψ(z))){\displaystyle {\mathbf {E} (r,z)}=E_{0}\,{\hat {\mathbf {x} }}\,{\frac {w_{0}}{w(z)}}\exp \left({\frac {-r^{2}}{w(z)^{2}}}\right)\exp \left(\!-i\left(kz+k{\frac {r^{2}}{2R(z)}}-\psi (z)\right)\!\right)}

donde [ 1 ] [ 3 ]

  • r es la distancia radial desde el eje central de la viga,
  • z es la distancia axial desde el foco (o "cintura") del haz,
  • i es la unidad imaginaria ,
  • k = 2 πn / λ es el número de onda para una longitud de onda en el espacio libre (vacío) λ , y n es el índice de refracción del medio en el que se propaga el haz,
  • E 0 = E (0, 0) , la amplitud del campo eléctrico en el origen ( r = 0 , z = 0 ),
  • w ( z ) es el radio en el que las amplitudes del campo caen a 1/ e de sus valores axiales (es decir, donde los valores de intensidad caen a 1/ e 2 de sus valores axiales), en el plano z a lo largo del haz,
  • w 0 = w (0) es el radio de la cintura ,
  • R ( z ) es el radio de curvatura de los frentes de onda del haz en z , y
  • ψ ( z ) = arctan( z / z R ) es la fase de Gouy en z , un término de fase adicional más allá del atribuible a la velocidad de fase de la luz.

El campo eléctrico físico se obtiene a partir de la amplitud del campo fasorial dada anteriormente, tomando la parte real de la amplitud multiplicada por un factor de tiempo: mifísico(r,z,t)=Re(mi(r,z)miiωt),{\displaystyle \mathbf {E} _{\text{phys}}(r,z,t)=\operatorname {Re} (\mathbf {E} (r,z)\cdot e^{i\omega t}),} dóndeω{\textstyle \omega }es la frecuencia angular de la luz y t es el tiempo. El factor tiempo implica una convención de signos arbitraria , como se analiza en Descripciones matemáticas de la opacidad §  Ambigüedad conjugada compleja .

Dado que esta solución se basa en la aproximación paraxial, no es precisa para haces muy divergentes. La forma anterior es válida en la mayoría de los casos prácticos, donde w 0λ / n .

La distribución de intensidad (o irradiancia ) correspondiente viene dada por

I(r,z)=|mi(r,z)|22η=I0(w0w(z))2exp(2r2w(z)2),{\displaystyle I(r,z)={|E(r,z)|^{2} \over 2\eta }=I_{0}\left({\frac {w_{0}}{w(z)}}\right)^{2}\exp \left({\frac {-2r^{2}}{w(z)^{2}}}\right),}

donde la constante η es la impedancia de onda del medio en el que se propaga el haz. Para el espacio libre, η = η 0 ≈ 377 Ω. I 0 = | E 0 | 2 /2 η es la intensidad en el centro del haz en su cintura.

Si P 0 es la potencia total del haz, I0=2PAG0πw02.{\displaystyle I_{0}={2P_{0} \over \pi w_{0}^{2}}.}

Ancho de haz en evolución

La función gaussiana tiene un diámetro de 1/ ( 2w como se usa en el texto) aproximadamente 1,7 veces el FWHM .

En una posición z a lo largo del haz (medida desde el foco), el parámetro de tamaño del punto w viene dado por una relación hiperbólica : [ 1 ]w(z)=w01+(zzR)2,{\displaystyle w(z)=w_{0}\,{\sqrt {1+{\left({\frac {z}{z_{\mathrm {R} }}}\right)}^{2}}},} donde [ 1 ]zR=πw02norteλ{\displaystyle z_{\mathrm {R} }={\frac {\pi w_{0}^{2}n}{\lambda }}} se denomina rango de Rayleigh, como se explica más adelante, ynorte{\displaystyle n}es el índice de refracción del medio.

El radio del haz w ( z ) , en cualquier posición z a lo largo del haz, está relacionado con el ancho completo a media altura (FWHM) de la distribución de intensidad en esa posición según: [ 4 ]w(z)=FWHM(z)2ln2.{\displaystyle w(z)={\frac {{\text{FWHM}}(z)}{\sqrt {2\ln 2}}}.}

Curvatura del frente de onda

Los frentes de onda tienen curvatura cero (radio de curvatura = ∞) en la cintura. La magnitud de la curvatura del frente de onda aumenta a medida que nos alejamos de la cintura, alcanzando un extremo en la distancia de Rayleigh, z = ± z R (máximo para z = + z R , mínimo para z = - z R ). Más allá de la distancia de Rayleigh, | z | > z R , la magnitud de la curvatura vuelve a disminuir, aproximándose a cero cuando z → ±∞ . La curvatura se suele expresar en términos de su recíproco, R , el radio de curvatura ; para un haz gaussiano fundamental, la curvatura en la posición z viene dada por:

1R(z)=zz2+zR2,{\displaystyle {\frac {1}{R(z)}}={\frac {z}{z^{2}+z_{\mathrm {R} }^{2}}},}

por lo tanto el radio de curvatura R ( z ) es [ 1 ]R(z)=z[1+(zRz)2].{\displaystyle R(z)=z\left[{1+{\left({\frac {z_{\mathrm {R} }}{z}}\right)}^{2}}\right].} Al ser el recíproco de la curvatura, el radio de curvatura invierte su signo y es infinito en el punto más bajo del haz, donde la curvatura pasa por cero.

Haces elípticos y astigmáticos

Muchos haces láser tienen una sección transversal elíptica. También son comunes los haces con puntos de cintura distintos para las dos dimensiones transversales, denominados haces astigmáticos. Estos haces pueden tratarse utilizando las dos ecuaciones de evolución anteriores, pero con valores distintos para cada parámetro en x e y , y definiciones distintas del punto z = 0. La fase de Gouy es un valor único que se calcula correctamente sumando la contribución de cada dimensión, con una fase de Gouy dentro del rango de ± π /4 aportada por cada dimensión.

Un haz elíptico invierte su relación de elipticidad a medida que se propaga desde el campo lejano hasta el centro. La dimensión que era mayor lejos del centro se vuelve menor cerca del centro.

Gaussiana como una descomposición en modos

Las soluciones arbitrarias de la ecuación de Helmholtz paraxial pueden descomponerse como la suma de modos de Hermite-Gauss (cuyos perfiles de amplitud son separables en x e y usando coordenadas cartesianas ), modos de Laguerre-Gauss (cuyos perfiles de amplitud son separables en r y θ usando coordenadas cilíndricas ) o de manera similar como combinaciones de modos de Ince-Gauss (cuyos perfiles de amplitud son separables en ξ y η usando coordenadas elípticas ). [ 5 ] [ 6 ] [ 7 ] En cualquier punto a lo largo del haz z, estos modos incluyen el mismo factor gaussiano que el modo gaussiano fundamental multiplicando los factores geométricos adicionales para el modo especificado. Sin embargo, los diferentes modos se propagan con una fase de Gouy diferente , razón por la cual el perfil transversal neto debido a una superposición de modos evoluciona en z , mientras que la propagación de cualquier modo Hermite-Gaussiano (o Laguerre-Gaussiano) individual conserva la misma forma a lo largo de un haz.

Existen otras descomposiciones modales , pero las gaussianas son útiles para problemas que involucran haces compactos, es decir, donde la potencia óptica está confinada a lo largo de un eje. Incluso cuando un láser no opera en el modo gaussiano fundamental, su potencia generalmente se encuentra entre los modos de orden más bajo utilizando estas descomposiciones, ya que la extensión espacial de los modos de orden superior tiende a exceder los límites del resonador (cavidad) del láser . "Haz gaussiano" normalmente implica radiación confinada al modo gaussiano fundamental (TEM 00 ).

Parámetros del haz

La dependencia geométrica de los campos de un haz gaussiano está regida por la longitud de onda de la luz λ ( en el medio dieléctrico, si no en el espacio libre) y los siguientes parámetros del haz , todos los cuales están conectados como se detalla en las siguientes secciones.

cintura del haz

Ancho del haz gaussiano w ( z ) en función de la distancia z a lo largo del haz, que forma una hipérbola . w₀ : cintura del haz; b : profundidad de foco; zR : rango de Rayleigh ; Θ : dispersión angular total

La forma de un haz gaussiano de una longitud de onda λ determinada está regida únicamente por un parámetro: el radio del haz w₀ . Este valor mide el tamaño del haz en su punto focal ( z = 0 en las ecuaciones anteriores), donde el ancho del haz w ( z ) (definido anteriormente) es mínimo (y, de igual modo, donde la intensidad en el eje ( r = 0 ) es máxima). Como se indicó anteriormente, w ( z ) y w₀ miden el radio del haz , no su diámetro, como se muestra en el diagrama.

A partir de este parámetro se determinan los demás parámetros que describen la geometría del haz. Esto incluye el rango de Rayleigh z R y la divergencia asintótica del haz θ , como se detalla a continuación.

Rango de Rayleigh y parámetro confocal

La distancia de Rayleigh o el alcance de Rayleigh z R se determina en función del tamaño de la cintura de un haz gaussiano:

zR=πw02norteλ.{\displaystyle z_{\mathrm {R} }={\frac {\pi w_{0}^{2}n}{\lambda }}.}

Aquí λ es la longitud de onda de la luz y n es el índice de refracción. A una distancia de la cintura igual al rango de Rayleigh z R , el ancho w del haz es 2 mayor que en el foco, donde w = w 0 , la cintura del haz. Esto también implica que la intensidad en el eje ( r = 0 ) allí es la mitad de la intensidad máxima (en z = 0 ). Ese punto a lo largo del haz también resulta ser donde la curvatura del frente de onda ( 1/ R ) es máxima. [ 1 ]

La distancia entre los dos puntos z = ± z R se denomina parámetro confocal o profundidad de foco del haz. [ 8 ]

divergencia del haz

Las colas de una función gaussiana nunca llegan a cero, pero para los fines de la siguiente discusión, el "borde" de un haz es el radio donde r = w ( z ) . Es ahí donde la intensidad ha disminuido a 1/ de su valor en el eje. Ahora bien, para zzR , el parámetro w ( z ) aumenta linealmente con z . Esto significa que lejos de la cintura, el "borde" del haz (en el sentido anterior) tiene forma de cono. El ángulo entre ese cono (cuyo r = w ( z ) ) y el eje del haz ( r = 0 ) define la divergencia del haz: θ=límitezarctan(w(z)z).{\displaystyle \theta =\lim _{z\to \infty }\arctan \left({\frac {w(z)}{z}}\right).}

En el caso paraxial, como hemos estado considerando, θ (en radianes) es entonces aproximadamente [ 1 ].θ=λπnortew0{\displaystyle \theta ={\frac {\lambda }{\pi nw_{0}}}}

donde n es el índice de refracción del medio a través del cual se propaga el haz, y λ es la longitud de onda en el espacio libre. La dispersión angular total del haz divergente, o ángulo del vértice del cono descrito anteriormente, viene dada por Θ=2θ.{\displaystyle \Theta =2\theta \,.}

Ese cono contiene entonces el 86% de la potencia total del haz gaussiano.

Debido a que la divergencia es inversamente proporcional al tamaño del punto, para una longitud de onda λ dada , un haz gaussiano enfocado en un punto pequeño diverge rápidamente a medida que se propaga lejos del foco. Por el contrario, para minimizar la divergencia de un haz láser en el campo lejano (y aumentar su intensidad máxima a grandes distancias) debe tener una gran sección transversal ( w₀ ) en la cintura (y por lo tanto un gran diámetro donde se lanza, ya que w ( z ) nunca es menor que w₀ ). Esta relación entre el ancho del haz y la divergencia es una característica fundamental de la difracción y de la transformada de Fourier que describe la difracción de Fraunhofer . Un haz con cualquier perfil de amplitud especificado también obedece esta relación inversa, pero el modo gaussiano fundamental es un caso especial donde el producto del tamaño del haz en el foco y la divergencia en el campo lejano es menor que para cualquier otro caso.

Dado que el modelo de haz gaussiano utiliza la aproximación paraxial, falla cuando los frentes de onda están inclinados más de unos 30° con respecto al eje del haz. [ 9 ] A partir de la expresión anterior para la divergencia, esto significa que el modelo de haz gaussiano solo es preciso para haces con cinturas mayores que aproximadamente 2 λ / π .

La calidad del haz láser se cuantifica mediante el producto de parámetros del haz (BPP). Para un haz gaussiano, el BPP es el producto de la divergencia del haz y el tamaño de cintura w₀ . El BPP de un haz real se obtiene midiendo el diámetro mínimo del haz y la divergencia en campo lejano, y calculando su producto. La relación entre el BPP del haz real y el de un haz gaussiano ideal a la misma longitud de onda se conoce como (" "). El para un haz gaussiano es uno. Todos los haces láser reales tienen valores de M² mayores que uno, aunque los haces de muy alta calidad pueden tener valores muy cercanos a uno.

La apertura numérica de un haz gaussiano se define como NA = n sin θ , donde n es el índice de refracción del medio a través del cual se propaga el haz. Esto significa que el rango de Rayleigh está relacionado con la apertura numérica mediante zR=nortew0norteA.{\displaystyle z_{\mathrm {R} }={\frac {nw_{0}}{\mathrm {NA} }}.}

Fase de Gouy

La fase de Gouy es un desplazamiento de fase adquirido gradualmente por un haz alrededor de la región focal. En la posición z, la fase de Gouy de un haz gaussiano fundamental viene dada por [ 1 ].ψ(z)=arctan(zzR).{\displaystyle \psi (z)=\arctan \left({\frac {z}{z_{\mathrm {R} }}}\right).}

Fase de Gouy.

La fase de Gouy produce un aumento en la longitud de onda aparente cerca de la cintura ( z ≈ 0 ). Por lo tanto, la velocidad de fase en esa región supera formalmente la velocidad de la luz . Este comportamiento paradójico debe entenderse como un fenómeno de campo cercano donde la desviación de la velocidad de fase de la luz (como ocurriría exactamente con una onda plana ) es muy pequeña, excepto en el caso de un haz con una gran apertura numérica , en cuyo caso la curvatura de los frentes de onda (véase la sección anterior) cambia sustancialmente en la distancia de una sola longitud de onda. En todos los casos, la ecuación de onda se satisface en cada posición.

El signo de la fase de Gouy depende de la convención de signos elegida para el fasor del campo eléctrico. [ 10 ] Con dependencia e iωt , la fase de Gouy cambia de - π /2 a + π /2 , mientras que con dependencia e - iωt cambia de + π /2 a - π /2 a lo largo del eje.

Para un haz gaussiano fundamental, la fase de Gouy produce una discrepancia de fase neta con respecto a la velocidad de la luz de π radianes (es decir, una inversión de fase) al pasar del campo lejano de un lado de la cintura al campo lejano del otro lado. Esta variación de fase no es observable en la mayoría de los experimentos. Sin embargo, tiene importancia teórica y adquiere un mayor alcance para los modos gaussianos de orden superior . [ 10 ]

Poder a través de una abertura

Si un haz gaussiano está centrado en una apertura circular de radio r a una distancia z de la cintura del haz, la potencia P que pasa a través de la apertura es [ 11 ].PAG(r,z)=PAG0[1mi2r2/w2(z)],{\displaystyle P(r,z)=P_{0}\left[1-e^{-2r^{2}/w^{2}(z)}\right],}

Para un círculo de radio r = w ( z ) , la fracción de potencia transmitida a través del círculo es PAG(w(z),z)PAG0=1mi20,865.{\displaystyle {\frac {P(w(z),z)}{P_{0}}}=1-e^{-2}\approx 0.865.}

De manera similar, aproximadamente el 90% de la potencia del haz fluye a través de un círculo de radio r = 1,07 × w ( z ) , el 95% a través de un círculo de radio r = 1,224 × w ( z ) , y el 99% a través de un círculo de radio r = 1,52 × w ( z ) . [ 11 ]

Parámetro de haz complejo

El tamaño del punto y la curvatura de un haz gaussiano en función de z a lo largo del haz también se pueden codificar en el parámetro complejo del haz q ( z ) [ 12 ] [ 13 ] dado por: q(z)=z+izR.{\displaystyle q(z)=z+iz_{\mathrm {R} }.}

El recíproco de q ( z ) contiene la curvatura del frente de onda y la intensidad relativa en el eje en sus partes real e imaginaria, respectivamente: [ 12 ]

1q(z)=1R(z)iλnorteπw2(z).{\displaystyle {1 \over q(z)}={1 \over R(z)}-i{\lambda \over n\pi w^{2}(z)}.}

El parámetro de haz complejo simplifica el análisis matemático de la propagación de haces gaussianos, y especialmente en el análisis de cavidades resonadoras ópticas utilizando matrices de transferencia de rayos .

Luego, utilizando esta forma, la ecuación anterior para el campo eléctrico (o magnético) se simplifica enormemente. Si llamamos u a la intensidad de campo relativa de un haz gaussiano elíptico (con los ejes elípticos en las direcciones x e y ), entonces se puede separar en x e y de acuerdo con: (incógnita,y,z)=incógnita(incógnita,z)y(y,z),{\displaystyle u(x,y,z)=u_{x}(x,z)\,u_{y}(y,z),}

dónde incógnita(incógnita,z)=1qincógnita(z)exp(ikincógnita22qincógnita(z)),y(y,z)=1qy(z)exp(iky22qy(z)),{\displaystyle {\begin{aligned}u_{x}(x,z)&={\frac {1}{\sqrt {{q}_{x}(z)}}}\exp \left(-ik{\frac {x^{2}}{2{q}_{x}(z)}}\right),\\u_{y}(y,z)&={\frac {1}{\sqrt {{q}_{y}(z)}}}\exp \left(-ik{\frac {y^{2}}{2{q}_{y}(z)}}\right),\end{aligned}}}

donde q x ( z ) y q y ( z ) son los parámetros complejos del haz en las direcciones x e y .

Para el caso común de un perfil de viga circular , q x ( z ) = q y ( z ) = q ( z ) y x 2 + y 2 = r 2 , lo que produce [ 14 ](r,z)=1q(z)exp(ikr22q(z)).{\displaystyle u(r,z)={\frac {1}{q(z)}}\exp \left(-ik{\frac {r^{2}}{2q(z)}}\right).}

Óptica de haces

Diagrama de un haz gaussiano que atraviesa una lente.

Cuando un haz gaussiano se propaga a través de una lente delgada , el haz saliente también es un haz gaussiano (diferente), siempre que el haz viaje a lo largo del eje de simetría cilíndrica de la lente y que la lente sea más grande que el ancho del haz. La distancia focal de la lenteF{\displaystyle f}, el radio de cintura del hazw0{\displaystyle w_{0}}y la posición de la cintura del hazz0{\displaystyle z_{0}}del haz incidente se puede utilizar para determinar el radio de la cintura del hazw0{\displaystyle w_{0}'}y posiciónz0{\displaystyle z_{0}'}del haz saliente.

ecuación de la lente

Como lo derivaron Saleh y Teich, la relación entre los haces entrantes y salientes se puede encontrar considerando la fase que se agrega a cada punto.(incógnita,y){\displaystyle (x,y)}del haz gaussiano a medida que viaja a través de la lente. [ 15 ] Un enfoque alternativo debido a Self es considerar el efecto de una lente delgada en los frentes de onda del haz gaussiano . [ 16 ]

La solución exacta al problema anterior se expresa simplemente en términos de la magnificación.METRO{\displaystyle M}

w0=METROw0(z0F)=METRO2(z0F).{\displaystyle {\begin{aligned}w_{0}'&=Mw_{0}\\[1.2ex](z_{0}'-f)&=M^{2}(z_{0}-f).\end{aligned}}}

La magnificación, que depende dew0{\displaystyle w_{0}}yz0{\displaystyle z_{0}}(o rango de RayleighzR{\displaystyle z_{R}}), viene dado por

METRO=METROr1+r2{\displaystyle M={\frac {M_{r}}{\sqrt {1+r^{2}}}}}

dónde

r=zRz0F,METROr=|Fz0F|.{\displaystyle r={\frac {z_{R}}{z_{0}-f}},\quad M_{r}=\left|{\frac {f}{z_{0}-f}}\right|.}

Una expresión equivalente para la posición del hazz0{\displaystyle z_{0}'}es

1z0+zR2(z0F)+1z0=1F.{\displaystyle {\frac {1}{z_{0}+{\frac {z_{R}^{2}}{(z_{0}-f)}}}}+{\frac {1}{z_{0}'}}={\frac {1}{f}}.}

Esta última expresión deja claro que la ecuación de la lente delgada de óptica geométrica se recupera en el límite que|(zRz0)(zRz0F)|1{\displaystyle \left|\left({\tfrac {z_{R}}{z_{0}}}\right)\left({\tfrac {z_{R}}{z_{0}-f}}\right)\right|\ll 1}También cabe señalar que si|z0+zR2z0F|F{\displaystyle \left|z_{0}+{\frac {z_{R}^{2}}{z_{0}-f}}\right|\gg f}entonces el haz entrante está "bien colimado" de modo quez0F{\displaystyle z_{0}'\approx f}.

Enfoque del haz

En algunas aplicaciones es deseable utilizar una lente convergente para enfocar un haz láser en un punto muy pequeño. Matemáticamente, esto implica minimizar la magnificación.METRO{\displaystyle M}. Si el tamaño del haz está limitado por el tamaño de la óptica disponible, esto generalmente se logra mejor enviando el haz colimado más grande posible a través de una lente de distancia focal pequeña, es decir maximizandozR{\displaystyle z_{R}}y minimizandoF{\displaystyle f}En esta situación, es justificable hacer la aproximación.zR2/(z0F)21{\displaystyle z_{R}^{2}/(z_{0}-f)^{2}\gg 1}, lo que implica queMETROF/zR{\displaystyle M\approx f/z_{R}}y obteniendo el resultadow0Fw0/zR{\displaystyle w_{0}'\approx fw_{0}/z_{R}}Este resultado se presenta a menudo en la forma

2w04πλF#z0F{\displaystyle {\begin{aligned}2w_{0}'&\approx {\frac {4}{\pi }}\lambda F_{\#}\\[1.2ex]z_{0}'&\approx f\end{aligned}}}

dónde

F#=F2w0,{\displaystyle F_{\#}={\frac {f}{2w_{0}}},}

que se encuentra después de suponer que el medio tiene un índice de refracciónnorte1{\displaystyle n\approx 1}y sustituyendozR=πw02/λ{\displaystyle z_{R}=\pi w_{0}^{2}/\lambda }Los factores de 2 se introducen debido a la preferencia común de representar el tamaño de la viga mediante los diámetros de la cintura de la viga.2w0{\displaystyle 2w_{0}'}y2w0{\displaystyle 2w_{0}}, en lugar de los radios de la cinturaw0{\displaystyle w_{0}'}yw0{\displaystyle w_{0}}.

ecuación de onda

Como caso especial de radiación electromagnética , los haces gaussianos (y los modos gaussianos de orden superior que se detallan a continuación) son soluciones a la ecuación de onda para un campo electromagnético en el espacio libre o en un medio dieléctrico homogéneo, [ 17 ] obtenidas al combinar las ecuaciones de Maxwell para el rotacional de E y el rotacional de H , lo que resulta en: 2U=1do22Ut2,{\displaystyle \nabla ^{2}U={\frac {1}{c^{2}}}{\frac {\partial ^{2}U}{\partial t^{2}}},} donde c es la velocidad de la luz en el medio , y U podría referirse al vector del campo eléctrico o magnético, ya que cualquier solución específica para uno determina el otro. La solución del haz gaussiano es válida solo en la aproximación paraxial , es decir, donde la propagación de la onda se limita a direcciones dentro de un pequeño ángulo de un eje. Sin pérdida de generalidad, tomemos esa dirección como la dirección + z, en cuyo caso la solución U generalmente se puede escribir en términos de u, que no tiene dependencia del tiempo y varía de forma relativamente suave en el espacio, con la variación principal espacialmente correspondiente al número de onda k en la dirección z : [ 17 ]U(incógnita,y,z,t)=(incógnita,y,z)mii(kzωt)incógnita^.{\displaystyle U(x,y,z,t)=u(x,y,z)e^{-i(kz-\omega t)}\,{\hat {\mathbf {x} }}\,.}

Utilizando esta forma junto con la aproximación paraxial, 2 u /∂ z 2 puede entonces despreciarse esencialmente. Dado que las soluciones de la ecuación de onda electromagnética solo son válidas para polarizaciones que son ortogonales a la dirección de propagación ( z ), sin pérdida de generalidad, hemos considerado que la polarización está en la dirección x, de modo que ahora resolvemos una ecuación escalar para u ( x , y , z ) .

Sustituyendo esta solución en la ecuación de onda anterior se obtiene la aproximación paraxial a la ecuación de onda escalar: [ 17 ]2incógnita2+2y2=2ikz.{\displaystyle {\frac {\partial ^{2}u}{\partial x^{2}}}+{\frac {\partial ^{2}u}{\partial y^{2}}}=2ik{\frac {\partial u}{\partial z}}.} Escribir las ecuaciones de onda en las coordenadas del cono de luz devuelve esta ecuación sin utilizar ninguna aproximación. [ 18 ] Los haces gaussianos de cualquier cintura de haz w 0 satisfacen la aproximación paraxial a la ecuación de onda escalar; esto se verifica más fácilmente expresando la onda en z en términos del parámetro de haz complejo q ( z ) como se definió anteriormente. Hay muchas otras soluciones. Como soluciones a un sistema lineal , cualquier combinación de soluciones (usando suma o multiplicación por una constante) también es una solución. La gaussiana fundamental resulta ser la que minimiza el producto del tamaño mínimo del punto y la divergencia de campo lejano, como se señaló anteriormente. Al buscar soluciones paraxiales, y en particular aquellas que describirían la radiación láser que no está en el modo gaussiano fundamental, buscaremos familias de soluciones con productos gradualmente crecientes de sus divergencias y tamaños mínimos de punto. Dos importantes descomposiciones ortogonales de este tipo son los modos de Hermite-Gauss y Laguerre-Gauss, que corresponden a la simetría rectangular y circular respectivamente, como se detalla en la siguiente sección. En ambos casos, el haz gaussiano fundamental que hemos estado considerando es el modo de orden más bajo.

Modos de orden superior

modos de Hermite-Gauss

Doce modos de Hermite-Gauss

Es posible descomponer un haz paraxial coherente utilizando el conjunto ortogonal de los llamados modos de Hermite-Gauss , cada uno de los cuales viene dado por el producto de un factor en x y un factor en y . Esta solución es posible debido a la separabilidad en x e y en la ecuación de Helmholtz paraxial escrita en coordenadas cartesianas . [ 19 ] Así, dado un modo de orden ( l , m ) referido a las direcciones x e y , la amplitud del campo eléctrico en x , y , z puede expresarse como: mi(incógnita,y,z)=l(incógnita,z)metro(y,z)exp(ikz),{\displaystyle E(x,y,z)=u_{l}(x,z)\,u_{m}(y,z)\,\exp(-ikz),} donde los factores de dependencia de x e y vienen dados cada uno por: J(incógnita,z)=(2/π2JJ¡w0)1/2(q0q(z))1/2(q(z)q(z))J/2HJ(2incógnitaw(z))exp(ikincógnita22q(z)),{\displaystyle u_{J}(x,z)=\left({\frac {\sqrt {2/\pi }}{2^{J}\,J!\;w_{0}}}\right)^{\!\!1/2}\!\!\left({\frac {{q}_{0}}{{q}(z)}}\right)^{\!\!1/2}\!\!\left(-{\frac {{q}^{\ast }(z)}{{q}(z)}}\right)^{\!\!J/2}\!\!H_{J}\!\left({\frac {{\sqrt {2}}x}{w(z)}}\right)\,\exp \left(\!-i{\frac {kx^{2}}{2{q}(z)}}\right),} donde hemos empleado el parámetro de haz complejo q ( z ) (como se definió anteriormente) para un haz de cintura w 0 en z desde el foco. En esta forma, el primer factor es simplemente una constante de normalización para hacer que el conjunto de u J sea ortonormal . El segundo factor es una normalización adicional que depende de z , que compensa la expansión de la extensión espacial del modo según w ( z )/ w 0 (debido a los dos últimos factores). También contiene parte de la fase de Gouy. El tercer factor es una fase pura que mejora el desplazamiento de fase de Gouy para órdenes J superiores .

Los dos últimos factores explican la variación espacial en x (o y ). El cuarto factor es el polinomio de Hermite de orden J ("forma de los físicos", es decir, H 1 ( x ) = 2 x ), mientras que el quinto explica la disminución de amplitud gaussiana exp(− x 2 / w ( z ) 2 ) , aunque esto no es obvio usando el complejo q en el exponente. La expansión de esa exponencial también produce un factor de fase en x , que explica la curvatura del frente de onda ( 1/ R ( z ) ) en z a lo largo del haz.

Los modos de Hermite-Gauss se designan típicamente como "TEM lm "; el haz gaussiano fundamental puede denominarse TEM 00 (donde TEM es electromagnético transversal ). Multiplicando u l ( x , z ) y u m ( y , z ) para obtener el perfil del modo 2D, y eliminando la normalización de modo que el factor principal se denomine simplemente E 0 , podemos escribir el modo ( l , m ) en la forma más accesible:

mil,metro(incógnita,y,z)=mi0w0w(z)Hl(2incógnitaw(z))Hmetro(2yw(z))×exp(incógnita2+y2w2(z))exp(ik(incógnita2+y2)2R(z))×exp(iψ(z))exp(ikz).{\displaystyle {\begin{aligned}E_{l,m}(x,y,z)={}&E_{0}{\frac {w_{0}}{w(z)}}\,H_{l}\!{\Bigg (}{\frac {{\sqrt {2}}\,x}{w(z)}}{\Bigg )}\,H_{m}\!{\Bigg (}{\frac {{\sqrt {2}}\,y}{w(z)}}{\Bigg )}\times {}\exp \left({-{\frac {x^{2}+y^{2}}{w^{2}(z)}}}\right)\exp \left({-i{\frac {k(x^{2}+y^{2})}{2R(z)}}}\right)\times {}\exp {\big (}i\psi (z){\big )}\exp(-ikz).\end{aligned}}}

En esta forma, el parámetro w 0 , como antes, determina la familia de modos, en particular escalando la extensión espacial de la cintura del modo fundamental y todos los demás patrones de modo en z = 0 . Dado que w 0 , w ( z ) y R ( z ) tienen las mismas definiciones que para el haz gaussiano fundamental descrito anteriormente . Se puede ver que con l = m = 0 obtenemos el haz gaussiano fundamental descrito anteriormente (ya que H 0 = 1 ). La única diferencia específica en los perfiles x e y en cualquier z se debe a los factores del polinomio de Hermite para los números de orden l y m . Sin embargo, hay un cambio en la evolución de la fase de Gouy de los modos sobre z : ψ(z)=(norte+1)arctan(zzR),{\displaystyle \psi (z)=(N+1)\,\arctan \left({\frac {z}{z_{\mathrm {R} }}}\right),}

donde el orden combinado del modo N se define como N = l + m . Mientras que el desplazamiento de fase de Gouy para el modo gaussiano fundamental (0,0) solo cambia en ± π /2 radianes en todo z (y solo en ± π /4 radianes entre ± z R ), este se incrementa por el factor N + 1 para los modos de orden superior. [ 10 ]

Los modos gaussianos de Hermite, con su simetría rectangular, son especialmente adecuados para el análisis modal de la radiación de láseres con cavidades de diseño asimétrico rectangular. Por otro lado, los láseres y sistemas con simetría circular se pueden analizar mejor utilizando el conjunto de modos gaussianos de Laguerre que se presenta en la siguiente sección.

modos de Laguerre-Gauss

Perfiles de intensidad de los primeros 12 modos de Laguerre-Gauss

Los perfiles de haz con simetría circular (o láseres con cavidades con simetría cilíndrica) suelen resolverse mejor mediante la descomposición modal de Laguerre-Gauss. [ 6 ] Estas funciones se escriben en coordenadas cilíndricas utilizando polinomios de Laguerre generalizados . Cada modo transversal se etiqueta nuevamente con dos números enteros, en este caso el índice radial p ≥ 0 y el índice azimutal l , que puede ser positivo o negativo (o cero): [ 20 ] [ 21 ]

Un haz de Laguerre-Gauss con l=1 y p=0. El rojo y el azul indican la intensidad del campo eléctrico con fase positiva y negativa.

(r,ϕ,z)=dolpagLGRAMO1w(z)(r2w(z))|l|exp(r2w2(z))Lpag|l|(2r2w2(z))×exp(ikr22R(z))exp(ilϕ)exp(iψ(z)),{\displaystyle {\begin{aligned}u(r,\phi ,z)={}&C_{lp}^{LG}{\frac {1}{w(z)}}\left({\frac {r{\sqrt {2}}}{w(z)}}\right)^{\!|l|}\exp \!\left(\!-{\frac {r^{2}}{w^{2}(z)}}\right)L_{p}^{|l|}\!\left({\frac {2r^{2}}{w^{2}(z)}}\right)\times {}\\&\exp \!\left(\!-ik{\frac {r^{2}}{2R(z)}}\right)\exp(-il\phi )\,\exp(i\psi (z)),\end{aligned}}}

donde L p l son los polinomios de Laguerre generalizados . C LG lp es una constante de normalización requerida: [ 22 ]dolpagLGRAMO=2pag¡π(pag+|l|)¡02πdϕ0drr|(r,ϕ,z)|2=1,{\displaystyle C_{lp}^{LG}={\sqrt {\frac {2p!}{\pi (p+|l|)!}}}\Rightarrow \int _{0}^{2\pi }d\phi \int _{0}^{\infty }dr\;r\,|u(r,\phi ,z)|^{2}=1,}.

w ( z ) y R ( z ) tienen las mismas definiciones que las anteriores . Al igual que con los modos Hermite-Gaussianos de orden superior, la magnitud del desplazamiento de fase de Gouy de los modos Laguerre-Gaussianos se exagera por el factor N + 1 : ψ(z)=(norte+1)arctan(zzR),{\displaystyle \psi (z)=(N+1)\,\arctan \left({\frac {z}{z_{\mathrm {R} }}}\right),} donde en este caso el número de modo combinado N = | l | + 2 p . Como antes, las variaciones de amplitud transversal están contenidas en los dos últimos factores en la línea superior de la ecuación, que nuevamente incluye la caída gaussiana básica en r pero ahora multiplicada por un polinomio de Laguerre. El efecto del número de modo rotacional l , además de afectar al polinomio de Laguerre, está contenido principalmente en el factor de fase exp(− ilφ ) , en el que el perfil del haz se adelanta (o retrasa) en l fases completas de 2 π en una rotación alrededor del haz (en φ ). Este es un ejemplo de un vórtice óptico de carga topológica l , y puede asociarse con el momento angular orbital de la luz en ese modo.

Modos Ince-Gaussianos

Perfil de amplitud transversal de los modos Ince-Gaussianos pares de orden más bajo

La tercera familia completa de soluciones para la ecuación de onda paraxial son los modos de Ince-Gauss. Describen haces con geometría transversal elíptica caracterizada por la elipticidad.ε{\displaystyle \varepsilon }Los modos Hermite-Gaussiano y Laguerre-Gaussiano son un caso especial de los modos Ince-Gaussiano paraε={\displaystyle \varepsilon =\infty }yε=0{\displaystyle \varepsilon =0}respectivamente. Los modos Ince-Gaussianos se pueden escribir utilizando coordenadas elípticas y polinomios de Ince . Los modos Ince-Gaussianos pares e imparesIGRAMOpag,metromi{\displaystyle \mathrm {IG} _{p,m}^{\mathrm {e} }}yIGRAMOpag,metroo{\displaystyle \mathrm {IG} _{p,m}^{\mathrm {o} }}están dados por: [ 7 ]

IGRAMOpag,metromi(r,ε)=dow0w(z)dopagmetro(iξ,ε)dopagmetro(η,ε)exp(r2w2(z))×expi(kz+kr22R(z)(pag+1)ψ(z)),IGRAMOpag,metroo(r,ε)=Sw0w(z)Spagmetro(iξ,ε)Spagmetro(η,ε)exp(r2w2(z))×expi(kz+kr22R(z)(pag+1)ψ(z)),{\displaystyle {\begin{aligned}\mathrm {IG} _{p,m}^{\mathrm {e} }\left(\mathbf {r} ,\varepsilon \right)&={\frac {Cw_{0}}{w\left(z\right)}}C_{p}^{m}\left(i\xi ,\varepsilon \right)C_{p}^{m}\left(\eta ,\varepsilon \right)\exp \left({\frac {-r^{2}}{w^{2}(z)}}\right)\\&\qquad \qquad \times \exp \mathrm {i} \left(kz+{\frac {kr^{2}}{2R(z)}}-(p+1)\psi (z)\right),\\\mathrm {IG} _{p,m}^{\mathrm {o} }\left(\mathbf {r} ,\varepsilon \right)&={\frac {Sw_{0}}{w\left(z\right)}}S_{p}^{m}\left(i\xi ,\varepsilon \right)S_{p}^{m}\left(\eta ,\varepsilon \right)\exp \left({\frac {-r^{2}}{w^{2}(z)}}\right)\\&\qquad \qquad \times \exp \mathrm {i} \left(kz+{\frac {kr^{2}}{2R(z)}}-(p+1)\psi (z)\right),\end{aligned}}}dóndedo{\displaystyle C},S{\displaystyle S}son constantes de normalización ydopagmetro{\displaystyle C_{p}^{m}},Spagmetro{\displaystyle S_{p}^{m}}los polinomios de Ince pares e impares de ordenpag{\displaystyle p}ymetro{\displaystyle m}respectivamente.ψ(z)=arctan(z/zR){\displaystyle \psi (z)=\arctan(z/z_{\mathrm {R} })}es la fase de Gouy yξ{\displaystyle \xi },η{\displaystyle \eta }son las coordenadas elípticas radiales y angulares definidas por  :incógnita=ε/2w(z)aporrearξporqueη,y=ε/2w(z)sinhξpecadoη.{\displaystyle {\begin{aligned}x&={\sqrt {\varepsilon /2}}\;w(z)\cosh \xi \cos \eta \,,\\y&={\sqrt {\varepsilon /2}}\;w(z)\sinh \xi \sin \eta \,.\end{aligned}}}

modos hipergeométricos-gaussianos

Existe otra clase importante de modos de onda paraxial en coordenadas cilíndricas , en la que la amplitud compleja es proporcional a una función hipergeométrica confluente .

Estos modos presentan un perfil de fase singular y son autofunciones del momento angular orbital del fotón . Sus perfiles de intensidad se caracterizan por un único anillo brillante; al igual que los modos de Laguerre-Gauss, sus intensidades se reducen a cero en el centro (sobre el eje óptico ), excepto para el modo fundamental (0,0). La amplitud compleja de un modo puede expresarse en términos de la coordenada radial normalizada (adimensional).ρ=r/w0{\displaystyle \rho =r/w_{0}}y la coordenada longitudinal normalizadaZ=z/zR{\displaystyle Z=z/z_{\mathrm {R} }}de la siguiente manera: [ 23 ]

pagmetro(ρ,ϕ,Z)=2pag+|metro|+1πΓ(pag+|metro|+1)Γ(pag2+|metro|+1)Γ(|metro|+1)i|metro|+1×Zpag2(Z+i)(pag2+|metro|+1)ρ|metro|×exp(iρ2Z+i)miimetroϕ1F1(pag2,|metro|+1;ρ2Z(Z+i)){\displaystyle {\begin{aligned}u_{{\mathsf {p}}m}(\rho ,\phi ,\mathrm {Z} ){}={}&{\sqrt {\frac {2^{{\mathsf {p}}+|m|+1}}{\pi \Gamma ({\mathsf {p}}+|m|+1)}}}\;{\frac {\Gamma \left({\frac {\mathsf {p}}{2}}+|m|+1\right)}{\Gamma (|m|+1)}}\,i^{|m|+1}\times {}\\&\mathrm {Z} ^{\frac {\mathsf {p}}{2}}\,(\mathrm {Z} +i)^{-\left({\frac {\mathsf {p}}{2}}+|m|+1\right)}\,\rho ^{|m|}\times {}\\&\exp \left(-{\frac {i\rho ^{2}}{\mathrm {Z} +i}}\right)\,e^{im\phi }\,{}_{1}F_{1}\left(-{\frac {\mathsf {p}}{2}},|m|+1;{\frac {\rho ^{2}}{\mathrm {Z} (\mathrm {Z} +i)}}\right)\end{aligned}}}

donde el índice de rotación m es un número entero, ypag|metro|{\displaystyle {\mathsf {p}}\geq -|m|}es de valor real, Γ( x ) es la función gamma y 1 F 1 ( a , b ; x ) es una función hipergeométrica confluente.

Algunas subfamilias de modos hipergeométricos-gaussianos (HyGG) se pueden enumerar como los modos Bessel-Gaussianos modificados, los modos gaussianos exponenciales modificados, [ 23 ] y los modos Laguerre-Gaussianos modificados.

El conjunto de modos hipergeométricos-gaussianos es sobrecompleto y no es un conjunto ortogonal de modos. A pesar de su perfil de campo complicado, los modos HyGG tienen un perfil muy simple en la cintura del haz ( z = 0 ): (ρ,ϕ,0)ρpag+|metro|miρ2+imetroϕ.{\displaystyle u(\rho ,\phi ,0)\propto \rho ^{{\mathsf {p}}+|m|}e^{-\rho ^{2}+im\phi }.}

Véase también

Notas

  1. ^ Svelto , págs . 153-5 .
  2. Svelto, pág. 158.
  3. Yariv, Amnon; Yeh, Albert Pochi (2003). Ondas ópticas en cristales: propagación y control de la radiación láser . J. Wiley & Sons. ISBN 0-471-43081-1OCLC 492184223 
  4. Hill, Dan (4 de abril de 2007). "Cómo convertir medidas FWHM a semianchos de 1/e al cuadrado" . Base de conocimientos de Radiant Zemax . Archivado del original el 4 de marzo de 2016. Recuperado el 7 de junio de 2016 .
  5. Siegman, pág. 642.
  6. 1 2 Probablemente fue considerado por primera vez por Goubau y Schwering (1961).
  7. 1 2 Bandrés y Gutiérrez-Vega (2004)
  8. Brorson, SD (1988). "¿Qué es el parámetro confocal?" . IEEE Journal of Quantum Electronics . 24 (3): 512– 515. Bibcode : 1988IJQE...24..512B . doi : 10.1109/3.155 .
  9. Siegman (1986) pág. 630.
  10. 1 2 3 Paschotta, Rüdiger (12 de diciembre de 2006). "Desplazamiento de fase de Gouy" . Enciclopedia de física y tecnología láser . RP Photonics . Recuperado el 2 de mayo de 2014 .
  11. 1 2 "Melles Griot. Gaussian Beam Optics" (PDF) . Archivado del original (PDF) el 4 de marzo de 2016. Consultado el 7 de abril de 2015 .
  12. 1 2 Siegman, págs. 638–40.
  13. Garg, págs. 165–168.
  14. Véase Siegman (1986), pág. 639, ecuación 29.
  15. ^ Saleh, Bahaa EA; Teich, Malvin Carl (1991). Fundamentos de Fotónica . Nueva York: John Wiley & Sons. ISBN 0-471-83965-5.Capítulo 3, "Óptica de haces"
  16. Self, Sidney (1 de marzo de 1983). "Enfoque de haces gaussianos esféricos". Applied Optics . 22 (5): 658– 661. Bibcode : 1983ApOpt..22..658S . doi : 10.1364/AO.22.000658 . PMID 18195851 . 
  17. 1 2 3 Svelto, págs. 148–9.
  18. Esarey, E.; Sprangle, P.; Pilloff, M.; Krall, J. (1995-09-01). "Teoría y velocidad de grupo de pulsos láser ultracortos y fuertemente enfocados" . JOSA B. 12 ( 9): 1695– 1703. Bibcode : 1995JOSAB..12.1695E . doi : 10.1364/JOSAB.12.001695 . ISSN 1520-8540 . 
  19. ^ Siegman (1986), p645, ecuación. 54
  20. Vallone, G. (8 de abril de 2015). "Sobre las propiedades de los haces circulares: normalización, expansión de Laguerre-Gauss y divergencia en el espacio libre". Optics Letters . 40 (8): 1717– 1720. arXiv : 1501.07062 . Bibcode : 2015OptL...40.1717V . doi : 10.1364/OL.40.001717 . PMID 25872056. S2CID 36312938 .  
  21. Miatto, Filippo M.; Yao, Alison M.; Barnett, Stephen M. (2011-03-15). "Caracterización completa del ancho de banda espiral cuántico de bifotones entrelazados" . Physical Review A. 83 ( 3) 033816. arXiv : 1011.5970 . Bibcode : 2011PhRvA..83c3816M . doi : 10.1103/PhysRevA.83.033816 . ISSN 1050-2947 . 
  22. Nótese que la normalización utilizada aquí (intensidad total para un z fijo igual a la unidad) difiere de la utilizada en la sección #Forma matemática para el modo gaussiano. Para l = p = 0, el modo de Laguerre-Gauss se reduce al modo gaussiano estándar, pero debido a las diferentes condiciones de normalización, las dos fórmulas no coinciden.
  23. 1 2 Karimi et al. (2007)

Referencias

  • Bandrés, Miguel A.; Gutiérrez-Vega, Julio C. (2004). "Ince haces gaussianos" . Letras de Óptica . 29 (2). OSA: 144– 146. Código bibliográfico : 2004OptL...29..144B . doi : 10.1364/OL.29.000144 . PMID 14743992 . 
  • Garg, Anupam (2012). Electromagnetismo clásico en pocas palabras . Princeton, NJ: Princeton University Press. ISBN 978-0-691-13018-7.
  • Goubau, G.; Schwering, F. (1961). "Sobre la propagación guiada de haces de ondas electromagnéticas". IRE Transactions on Antennas and Propagation . 9 (3): 248– 256. Bibcode : 1961ITAP....9..248G . doi : 10.1109/TAP.1961.1144999 . MR 0134166 . 
  • Karimi, E.; Zito, G.; Piccirillo, B.; Marrucci, L.; Santamato, E. (2007). "Vigas hipergeométricas-gaussianas" . Letras de Óptica . 32 (21). OSA: 3053– 3055. arXiv : 0712.0782 . Código Bib : 2007OptL...32.3053K . doi : 10.1364/OL.32.003053 . PMID 17975594 . S2CID 46526713 .  
  • Mandel, Leonard; Wolf, Emil (1995). Coherencia óptica y óptica cuántica . Cambridge: Cambridge University Press. ISBN 0-521-41711-2.Capítulo 5, "Haces ópticos", pág.  267.
  • Pampaloni, F.; Enderlein, J. (2004). "Haces gaussianos, de Hermite-Gauss y de Laguerre-Gauss: una introducción". arXiv : physics/0410021 .
  • Sakpal, S.; Milione, G.; Li, M.; Nouri, M.; Shahoei, H.; LaFave, T.; Ashrafi, S.; MacFarlane, D. (2018). "Estabilidad de haces Ince-Gaussianos en fibras de pocos modos con núcleo elíptico" . Optics Letters . 43 (11): 2656– 2659. Bibcode : 2018OptL...43.2656S . doi : 10.1364/ OL.43.002656 . PMID 29856389. S2CID 46921059 .  
  • Saleh, Bahía EA; Teich, Malvin Carl (1991). Fundamentos de Fotónica . Nueva York: John Wiley & Sons. ISBN 0-471-83965-5.Capítulo 3, "Óptica de haces", págs.  80-107.
  • Siegman, Anthony E. (1986). Láseres . Libros de ciencias universitarias. ISBN 0-935702-11-3.Capítulo 16.
  • Svelto, Orazio (2010). Principios de los láseres (5ª  ed.).
  • Yariv, Amnón (1989). Electrónica cuántica (3ª  ed.). Wiley. ISBN 0-471-60997-8.
  • Tutorial de óptica de haces gaussianos, Newport
  • Simulador científico interactivo de enfoque de haz gaussiano, ephotonics