Do these solutions look vaguely correct?

  • Thread starter bumclouds
  • Start date
In summary, the conversation discusses the use of MATLAB code and the finite square well in physics problems. It also mentions the use of the WKB approximation and a numerical solver for solving the problem. The x-axis represents normalized energy and the y-axis represents the determinant of a matrix created from boundary and interface conditions. There may be a mathematical artifact at the origin and the number of bound states and their energies depend on the value of V0.
  • #1
bumclouds
25
0

Homework Statement



assignment.jpg


I wrote a whole heap of MATLAB code, and I have no idea whether my answers, or the plots of my answers are anywhere in the vicinity of looking right. I was hoping someone with a lil bit of experience in this type of physics would be able to tell me if I'm on the right track..

Homework Equations



None.

The Attempt at a Solution



for [tex]\alpha[/tex] = 1
alpha1.jpg


for [tex]\alpha[/tex] = 2
alpha2.jpg


for [tex]\alpha[/tex] = 5
alpha3.jpg


for [tex]\alpha[/tex] = 10
alpha10.jpg
 
Physics news on Phys.org
  • #2
This is a variation of the finite square well.

As a result, the wavefuctions and energies should be solvable analytically (analytically meaning you can form closed form expressions for the variables involved), although the relationships for the boundary conditions may yield transcendental equations.

However, if you have to use matlab, you should bear in mind that it is similar to a finite square well, and the energies should be about the same order of magnitude.

I think the WKB approximation may be used here.

You may find this numerical solver useful:

http://qbx6.ltu.edu/s_schneider/physlets/main/finite_sqr_well.shtml [Broken]

What do the x-axis and y-axis represent?
 
Last edited by a moderator:
  • #3
the x-axis is Ehat - the energy, normalised. So: Ehat = E [tex]\frac{L}{2}[/tex], where L is the length of gap in my "well"

and the y-axis is detS - the determinant of the matrix I created from the boundary and interface conditions. So each x-axis crossing should give me an energy eigenvalue shouldn't it?

if so, why am I only getting one zero-crossing in each graph?
 
Last edited:
  • #4
Could I take a look at your S matrix?
 
  • #5
yep.

S =

[tex]e^{-a} \ \ \ \ \ -e^{-jk} \ \ \ \ \ -e^{jk} \ \ \ \ \ 0[/tex]
[tex]ae^{-a} \ \ \ \ \ -jke^{-jk} \ \ \ \ \ jke^{jk} \ \ \ \ \ 0 [/tex]
[tex]0 \ \ \ \ \ e^{jk} \ \ \ \ \ e^{-jk} \ \ \ \ \ -e^{-a}[/tex]
[tex]0 \ \ \ \ \ jke^{jk} \ \ \ \ \ -jke^{-jk} \ \ \ \ \ ae^{-a}[/tex]


This is the matrix "S" when either side of the potential well are the same height. This is copied from the teachers notes that he wrote on the board.



but when region 3 is heigher than region 1, as in the assignment question i got given..

i tried introducting a new vector, B. So that my matrix S changed thus:

[tex]e^{-a} \ \ \ \ \ -e^{-jk} \ \ \ \ \ -e^{jk} \ \ \ \ \ 0[/tex]
[tex]ae^{-a} \ \ \ \ \ -jke^{-jk} \ \ \ \ \ jke^{jk} \ \ \ \ \ 0 [/tex]
[tex]0 \ \ \ \ \ e^{jk} \ \ \ \ \ e^{-jk} \ \ \ \ \ -e^{-B}[/tex]
[tex]0 \ \ \ \ \ jke^{jk} \ \ \ \ \ -jke^{-jk} \ \ \ \ \ Be^{-B}[/tex]


where a = [tex]\frac{pi}{2}\sqrt{Vo-E}[/tex]
and B = [tex]\frac{pi}{2}\sqrt{\alpha Vo-E}[/tex]

and I think this might be where I'm going wrong..
 
  • #6
Your S-Matrix seems fine.

When you say real, is it |Re[function]| or Abs[function] (ie, Re[x+ iy] = x, Abs[x + iy] = [tex]\sqrt{x^2 + y^2}[/tex])?

I worked out a solution in Mathematica, using a slightly different method (basically equating all boundary conditions) and forming the equation

[tex]Tan(k L) = Something Horrible[/tex]
[tex]k = \sqrt{{\frac{\hbar^2 E}{2m}}}[/tex]

From what I observe, the 0 at the origin is a mathematical artifact that comes from dividing by 0.

[tex]Something Horrible = \frac{Something_a}{\frac{Something_b}{0} + Something_c}[/tex]

Additionally, I found it necessary to choose a value for [tex]V_0[/tex] in the form:
[tex]z_0 = \sqrt{\frac{2 m V_0 L^2}{\hbar^2}}[/tex]

Does the question indicate what [tex]V_0[/tex] should be, or am I mistaken?

There will be a minimum of one bound state. The maximum number of bound states (zero crossings) depends on how large [tex]V_0[/tex] is.

I think your plots are alright. The only thing to note is that in your first plot [tex]\alpha = 1[/tex], it appears that there are two crossings. If you reduce your scale, I think there's a chance that the line at 4 approaches the x-axis closely,but doesn't actually touch it.

This would explain why that there's only one crossing for the other alphas.

Furthermore, depending on how you implemented the boundary conditions, the S matrix will not be valid beyond a certain energy for the excited states - that energy being [tex]\alpha V_o > E > V_o[/tex], which represent plane waves extending from -infinity and being reflected off the [tex]\alpha V_o[/tex] barrier, while [tex]E > \alpha V_o[/tex] would just give plane waves from +-infinity
 
Last edited:
  • #7
Must I have a value for Vo to find the eigenenergy values?

In all the examples my teacher has ever given, he always seems to write Vo, but not actually specify what it is.
 
  • #8
bumclouds said:
Must I have a value for Vo to find the eigenenergy values?

In all the examples my teacher has ever given, he always seems to write Vo, but not actually specify what it is.

Vo needs to be specified.

For example, in a finite square well, the number of bound states and their energies depends on Vo. As Vo increases, the number of bound states changes as does their energies (slightly).

When Vo -> infinity, we get the infinite square well which has an infinite number of bound states.
 
  • #9
What do you mean by bound states?

I'm really sorry to be such a newbie.!
 
  • #10
A bound state is a state with total energy Kinetic Energy + Potential Energy (appropriately scaled) = E + V < 0.

Alternatively, it is a state whose motion is restricted to a finite region, ie. bounded.

For example, for planets orbiting around the sun, the Kinetic Energy + Gravitational Potential Energy of the system < 0, so they represent bound states (A good thing for Earth!).

One point of note is that the absolute value of potential is arbitrary - only relative values matter. As such, we typically set Potential Energy to be 0 for points that are very far away, and work backwards to find Potential Energy at the origin or example.

So for gravitationally bound systems, the Potential Energy is set to be 0 at infinity, and we find out that it becomes negative for points not at infinity.

Similarly for this modified finite square well, we would set Vo = 0 at infinity, so the Potential Energy of a particle in the well is -Vo for the purposes of determining whether it is bound or not.

The reason we do not set [tex]\alpha V_o[/tex] = 0 at infinite instead is because a particle with energy = Vo would already be able to travel to infinity, let alone a particle with even more energy = [tex]\alpha V_o[/tex]; ie. not bounded by the second definition.

I get the feeling that you are doing Quantum Mechanics without really understanding what is going on - I would encourage reading up on introductory texts such as http://www.amazon.com/dp/0131911759/?tag=pfamazon01-20, or the plethora of online lecture notes.
 
Last edited:
  • #11
Quantum mechanics without really knowing what is going on - that is an accurate description. I'm studying electronic engineering - but we're all forced to take on this subject which is called Electronc Materials. I'm finding it tough with little to no background in Physics.

Anyway, I've bought the book you suggested on eBay. I'm going into uni tomorrow to see if I can recruit some help from somebody. I'll post whatever progress I've made with my assignment here tomorrow night - either as a JPEG or PDF, so you can see. You might find it amusing because I'm so bad.
 

1. What does it mean for a solution to look "vaguely correct"?

For a solution to look "vaguely correct" means that it appears to be somewhat accurate or acceptable, but may not be entirely precise or perfect.

2. How can I tell if a solution looks vaguely correct?

You can tell if a solution looks vaguely correct by comparing it to the expected outcome or solution. If it is close, but not exactly the same, it may be considered vaguely correct.

3. Is a solution that looks vaguely correct still considered valid?

It depends on the context and the level of accuracy required. In some cases, a solution that looks vaguely correct may still be considered valid, while in others it may not be acceptable.

4. What factors can contribute to a solution looking vaguely correct?

Factors that can contribute to a solution looking vaguely correct include human error, lack of precision in measurements or calculations, and incomplete data or assumptions.

5. How can I improve a solution that looks vaguely correct?

To improve a solution that looks vaguely correct, you can double check your calculations and assumptions, gather more accurate data, and seek feedback from others to identify any potential errors or areas for improvement.

Similar threads

  • Advanced Physics Homework Help
Replies
1
Views
2K
  • Advanced Physics Homework Help
Replies
1
Views
2K
Replies
2
Views
961
  • Advanced Physics Homework Help
Replies
11
Views
2K
  • Advanced Physics Homework Help
Replies
6
Views
1K
  • Math Proof Training and Practice
Replies
9
Views
1K
  • Engineering and Comp Sci Homework Help
Replies
2
Views
820
  • Advanced Physics Homework Help
Replies
1
Views
2K
  • Introductory Physics Homework Help
Replies
4
Views
1K
  • Calculus
Replies
8
Views
1K
Back
Top