EM Lasing Medium Green's Function

Consider electromagnetic waves polarized in the direction propagating in one dimension in an active laser medium with current density J ( x ) J(x) . Maxwell's equation for the z z -component of the electric field of these waves then reads

( d 2 d x 2 κ 2 ) E z ( x ) = J ( x ) \left(\dfrac{d^2}{dx^2} - \kappa^2 \right) E_z (x) = J(x)

for κ \kappa some constant. Suppose there is a source at x = 0 x= 0 and that there is a conducting mirror far away so that effectively E z ( 0 ) = 1 E_z (0) = 1 and E z ( ) = 0 E_z (\infty) = 0 are the boundary conditions. Find the Green's function G ( x , y ) G(x,y) for x < y x<y that will allow for determination of the electric field anywhere in space.

e κ x e κ y κ sinh ( κ x ) e^{-\kappa x} - \frac{e^{-\kappa y}}{\kappa} \sinh (\kappa x) e κ x e κ y 2 κ cosh ( κ x ) e^{-\kappa x} - \frac{e^{\kappa y}}{2\kappa} \cosh (\kappa x) sinh ( κ y ) sinh ( κ x ) κ sinh ( κ y ) \frac{\sinh (\kappa y) \sinh (\kappa x)}{\kappa \sinh (\kappa y)} sinh ( κ x ) sinh ( κ y ) κ sinh ( κ x ) \frac{\sinh (\kappa x) \sinh (\kappa y)}{\kappa \sinh (\kappa x)}

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

1 solution

Matt DeCross
Apr 29, 2016

The homogeneous equation is solved by the general expression:

G ( x , y ) = { c 1 e κ x + c 2 e κ x , x < y d 1 e κ x + d 2 e κ x , x > y . G(x,y) = \begin{cases} c_1 e^{\kappa x} + c_2 e^{-\kappa x}, \qquad &x<y \\ d_1 e^{\kappa x} + d_2 e^{-\kappa x}, \qquad &x>y \end{cases}.

The boundary conditions G ( 0 , y ) = 1 , G ( , y ) = 0 G(0,y) = 1, G(\infty,y) = 0 impose the constraints on the coefficients:

c 1 + c 2 = 1 d 1 = 0 \begin{aligned} c_1 + c_2 &= 1 \\ d_1 = 0 \end{aligned}

Solving for c 2 c_2 and d 2 d_2 and plugging into the expression for G ( x , y ) G(x,y) yields:

G ( x , y ) = { 2 c 1 sinh ( κ x ) + e κ x , x < y d 2 e κ x , x > y . G(x,y) = \begin{cases} 2c_1 \sinh (\kappa x) + e^{-\kappa x}, \qquad &x<y \\ d_2 e^{-\kappa x} , \qquad &x>y \end{cases}.

Enforcing continuity at x = y x=y gives:

2 c 1 sinh ( κ y ) + e κ y = d 2 e κ y , 2c_1 \sinh (\kappa y) + e^{-\kappa y}= d_2 e^{-\kappa y} ,

and requiring the derivative to jump by unity yields:

κ d 2 e κ y 2 c 1 κ cosh ( κ y ) + κ e κ y = 1. -\kappa d_2 e^{-\kappa y } - 2c_1 \kappa \cosh (\kappa y) + \kappa e^{-\kappa y} = 1.

Solving the two conditions for c 1 , d 1 c_1, d_1 yields:

c 1 = e κ y 2 κ d 2 = 1 sinh ( κ y ) κ \begin{aligned} c_1 &= -\frac{e^{-\kappa y}}{2 \kappa } \\ d_2 &= 1-\frac{\sinh (\kappa y)}{\kappa } \end{aligned}

Plugging c 1 c_1 into G ( x , y ) G(x,y) for x < y x<y , find the result:

G ( x , y ) = e κ x e κ y κ sinh ( κ x ) , x < y , G(x,y) =e^{-\kappa x}-\frac{e^{-\kappa y}}{\kappa } \sinh (\kappa x) , \qquad x<y,

as claimed.

This is wrong. The green function G(0,y)=0 is the correct selection.

Hamidreza Ramezani - 2 years, 7 months ago

Log in to reply

No, I believe the original solution is correct. We do not have Dirichlet boundary conditions in this problem, unlike the example in the wiki this problem is used in. The boundary condition G(0,y) = 0 is not the correct boundary condition to use. In this problem we have added a source at x=0 instead of a conducting mirror, which changes the boundary conditions as given in the problem statement.

Matt DeCross - 2 years, 7 months ago

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...