Scattering Within Rectangular Waveguide With A Step Obstacle

[skujesco2014]

Dear PF:

I'm currently working in a problem that has had me stranded for several weeks now. The problem reads as follows:

(See attachment)

Consider a beam of quantum particles (that is, the particles are small enough to exhibit non-negligible quantum effects) that propagates through a two-dimensional waveguide of width $H$ from $x=-\infty$ to $x=+\infty$. At $x=0$ the particles encounter a step of height $0<H_0<H$. All walls are impenetrable. Calculate the reflection and transmission coefficients.

Approach

The potential within the waveguide can be described as:

$$V(x,y)=<br /> \begin{cases}<br /> \infty & \text{at} \ AO, OH_0, H_0D, CE, OF\\<br /> 0 & \text{elsewhere}<br /> \end{cases}$$

A particular solution I worked out was:

$$\psi(x,y)=<br /> \begin{cases}<br /> A_1^+\sin k_1x\sin\frac{\pi y}{H}, & x<0, 0<y<H_0 \\<br /> A_1^+e^{ik_1x}\sin \frac{\pi y}{H}+ \sum\limits_mA_m^-e^{-ik_mx}\sin\frac{m\pi y}{H},& x<0, H_0<y<H\\<br /> \sum\limits_mB_m^+e^{ik'_mx}\sin\frac{m\pi(y-H_0)}{H-H_0}, & x>0, H_0<y<H<br /> \end{cases}$$

where $k_1=\sqrt{K^2-(\pi/H)^2}, k_m = \sqrt{K^2-(m\pi/H)^2}, k'_m=\sqrt{K^2-(m\pi/(H-H_0))^2}$, $E=\hbar^2K^2/2m$.

The wavefunction above must satisfy the following boundary conditions:

$$\psi(AO, OH_0, H_0D, CE, OF)=0 ,$$
$$\psi(HH_0^-)=\psi(HH_0^+),$$
$$\psi'(HH_0^-)=\psi'(HH_0^+)$$

From the above, I can calculate the constants $A$ and $B$, but all I get is nonsense. In particular, the above solution is not continuous along the boundary $BH_0$, but as hard as I try, I cannot make a satisfactory modification such that this discontinuity is healed.

What am I doing wrong? Could someone direct me to a similar problem? I'm sure there has to be a treatise for presence of steps in rectangular waveguides, but I can't seem to find any.

Thanks in advance.

 

[Jilang]

Hi skujesco2014!
It looks like for the solutions in regions I and II to be continuous along BHo you need the same solution for both regions with the reflected wave perhaps modified by a exp -(y-Ho) type term. It's just a hunch as I haven't looked at a waveguide problem before.

I'm not sure why you need to introduce the summations where you have though. I thought it was implicit anyway In a general solution.

 

[skujesco2014]

Hi, Jilang. Thanks for your input.

If you notice, the incident waves are the form $A_1^+e^{ik_1x}\sin(\pi y/H)$ but waves in region I are fully reflected whereas in region II they're partially reflected. The sums are introduced because I'm expanding in the basis of the complete sets $\{\sin m\pi y/H\}$ and $\{\sin m\pi (y-H_0)/(H-H_0)\}$ for $x<0$ and $x>0$ respectively. I could do what you say, but the functions would still not match essentially because of the presence of the $e^{ik_1x}$ term.

 

[Jilang]

Yes I see what you are saying, I think both solutions need to include the reflected term. The waves will bounce off the the bottom wall and head upwards, no?

 

[Jilang]

I just reread your last post. Are your first and second expressions the wrong way round then for the segments?

 

[skujesco2014]

The reflected wave from the $0<y<H_0$ region is included: the incoming wave is $A_1^+e^{ik_1x}\sin(\pi y/H)$ and when hitting the step it is 100% reflected; thus, the reflected wave consists of $\sum_m A_m^-e^{-ik_mx}\sin(m\pi y/H)$. Since the wavefunction has to be zero on $OH_0$, we conclude $A_m^-=-A_1^+\delta_{m1}$. Combining exponentials we get $\sin k_1x$.

The wavefunction as it was given satisfies being zero along all the boundaries (i.e., all the walls). If the BCs (boundary conditions) are implemented on the interface $H_0H$ (see attachment) one can determine expressions for the coefficients $A_m^-$ and $B_m^+$ in terms of $A_1^+$ (I can show you this math if you want).

Given a number of particles coming in (so far, there is only one particle incoming in quantum state 1, but this is arbitrary), I'm interested in computing an expression that gives me the number of particles transmitted through the interface $H_0H$ and reflected. I have, in fact, computed these expressions, but I obtain a very strange answer. I discovered that $\psi$ is not continuous along $BH_0$, thus my expressions are wrong. I can't, however, come up with a satisfactory wavefunction for the problem. That's where I'm stuck.

 

[Jilang]

Have you tried including a third term in region II for a backwards reflection at the second frequency?

 

[skujesco2014]

That term is $\sum_m A^−_m e^{−ik_mx}\sin(m\pi y/H)$.

 

[Jilang]

I meant a sin mπ(y-Ho)/(H-Ho) type term.

 

[skujesco2014]

If we do so, then the wavefunction wouldn't be zero on the $AO$ wall.

 

[Jilang]

It would be problematic in region I, but would it work just in region II?

 

[Jilang]

Looking at it holistically, I think there is going to be total reflection unless n= MH/(H-Ho). Only waves of certain frequencies are going to make it through. Is this a filter?

 

[skujesco2014]

It would be problematic in region I, but would it work just in region II?

I don't understand what you meant there.

Looking at it holistically, I think there is going to be total reflection unless n= MH/(H-Ho).

How do you arrive to this conclusion?

 

[Jilang]

What I am thinking is for continuity at the line H Ho you need to fit an integer number of half wavelengths between H and Ho and to also be able to fit an integer number of half wavelengths between H and O. This will only be possible for certain special wavelengths.

 

[skujesco2014]

I see what you mean now.

So, because we have to have an integer number of half-waves in the right-half of the waveguide, only those will that have a node along $BH_0$ will make it through. That would imply then that the wavefunction is forced to have a node along the $BH_0$ line, would you agree?

 

[Jilang]

Yes, that's the way I see it, anything else would be reflected. I can't see how else you avoid the discontinuities.