Particle in a box: Difference between revisions
imported>Michael Underwood m (Added See Also section linking to spherical well) |
imported>Johan Förberg |
||
(30 intermediate revisions by 8 users not shown) | |||
Line 1: | Line 1: | ||
{{subpages}} | {{subpages}} | ||
{{Image|Classical particle in a box with velocity.jpg|right|207px|Figure 1: The classical view of a particle (with velocity <math>v</math>) in a 1D box. With constant kinetic energy it bounces back and forth between the walls at <math>x=0</math> and <math>x=L</math>. Classically the particle is allowed to have any non-negative value for its kinetic energy, including zero.}} | |||
The '''particle in a box''' or '''infinite square well''' problem is one of the simplest non-trivial solutions to [[Schrödinger equation|Schrödinger's wave equation]]. As such it is often encountered in introductory [[quantum mechanics]] material as a demonstration of the [[quantization (physics)|quantization]] of energy. | The '''particle in a box''' or '''infinite square well''' problem is one of the simplest non-trivial solutions to [[Schrödinger equation|Schrödinger's wave equation]]. As such it is often encountered in introductory [[quantum mechanics]] material as a demonstration of the [[quantization (physics)|quantization]] of energy. | ||
In its simplest form the problem is one-dimensional, and involves a single [[particle]] living in an infinite [[potential well]]. | In its simplest form the problem is one-dimensional (1D), and involves a single [[particle]] living in an infinite [[potential well]]. | ||
What this means in the 1D case, which we examine here, is that there is a particle that is free to move back and forth along | What this means in the 1D case, which we examine here, is that there is a particle that is free to move back and forth along | ||
a line in between two points (e.g. 0 and <math>L</math>, or <math>-a/2</math> and <math>a/2</math>). The particle is | a line in between two points (e.g. 0 and <math>L</math>, or <math>-a/2</math> and <math>a/2</math>). The particle is | ||
Line 11: | Line 12: | ||
==The 1D square well== | ==The 1D square well== | ||
===Setting up the problem=== | ===Setting up the problem=== | ||
As depicted in Figure 1, the [[potential energy]] for this system can be given by | |||
As depicted in Figure 1, the potential can be given by | |||
:<math> | :<math> | ||
V(x)=\left\{\begin{align} | V(x)=\left\{\begin{align} | ||
Line 21: | Line 21: | ||
That is, in a region of length <math>L</math> (the box or well) the potential is zero; everywhere else it is infinite. | That is, in a region of length <math>L</math> (the box or well) the potential is zero; everywhere else it is infinite. | ||
An immediate consequence of this is that the particle must be located somewhere between <math>x=0</math> and <math>x=L</math> | Since the total [[energy]] of the particle is given by the sum of its potential and kinetic energies, if it is in a region of infinite potential then the particle itself must have infinite energy, which of course cannot happen. | ||
An immediate consequence of this potential is therefore that the particle must be located somewhere between <math>x=0</math> and <math>x=L</math>. | |||
This means that we only have to solve Schrödinger's equation inside the box and since the potential does not depend on time we can use the time-independent version of the Schrödinger equation. With <math>V(x)=0</math> in the box the [[wavefunction]], <math>\psi(x)</math>, that describes the state of the particle must satisfy the [[differential equation]] (DE) | |||
:<math> -\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2}\psi(x)=E\psi(x)\ , </math> | :<math> -\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2}\psi(x)=E\psi(x)\ , </math> | ||
where <math>\hbar</math> is [[Planck's constant]] divided by 2[[pi|π]], and <math>m</math> and <math>E</math> are the [[mass]] and energy of the particle, respectively. | |||
This is nothing but the [[eigenvalue problem]], and our task is to determine the energy eigenvalue(s) <math>E</math> and [[eigenstate]](s) <math>\psi</math> that solve it. | where <math>\scriptstyle\hbar</math> is [[Planck's constant]] divided by 2[[pi|π]], and <math>m</math> and <math>E</math> are the [[mass]] and mechanical energy of the particle, respectively. | ||
:<math> | This is nothing but the [[eigenvalue problem]], and our task is to determine the energy eigenvalue(s) <math>E</math> and [[eigenstate]](s) <math>\psi</math> that solve it. To simplify the expressions, let | ||
where <math>a</math> and <math>b</math> are | |||
The first condition we can use is that the wavefunction must be [[continuous]]. | :<math> k = \sqrt{\frac{2mE}{\hbar^2}}.</math> | ||
which tells us that the wavefunction must be | |||
We shall later see that this choice of letter is justified, since ''k'' corresponds to the wave number of the solutions. We are now faced with a linear, second-order differential equation. The real solutions are most readily written as linear combinations of | |||
:<math> \psi(x) = a\cos (kx) + b \sin (kx)\ ,</math> | |||
where <math>a</math> and <math>b</math> are constants to be determined from conditions on the system. | |||
(The interested reader is invited to find the second derivative of this function and substitute it back into the differential equation above to verify that it is indeed satisfied.) | |||
The first condition we can use is that the wavefunction must be [[continuous]]. | |||
We therefore know that the particle cannot exist at any position <math>x<0</math>, | |||
which tells us that the wavefunction must be zero at negative positions (in order for the probability of finding the particle there to also be zero). | |||
Continuity then implies that ''at'' <math>x=0</math> the wavefunction is also zero, so we have | Continuity then implies that ''at'' <math>x=0</math> the wavefunction is also zero, so we have | ||
:<math> 0=\psi(0)=a\cos(0)+b\sin(0)=a\ | :<math> 0=\psi(0)=a\cos(0)+b\sin(0)=a\ .</math> | ||
With <math>a=0</math> we now know that the wavefunction is | |||
:<math> \psi(x)=b\sin | :<math> \psi(x)=b\sin (kx)\ .</math> | ||
===Energy quantization=== | ===Energy quantization=== | ||
Line 68: | Line 78: | ||
===The solutions=== | ===The solutions=== | ||
[[Image:particle-in-a-box.png|thumb|320px|Figure 2: A plot of <math>\psi_n(x)</math> for the first 7 wavefunctions, <math>n=1</math> (darkest) to <math>n=7</math> (lightest).]] | [[Image:particle-in-a-box.png|thumb|320px|{{#ifexist:Template:Particle-in-a-box.png/credit| | ||
{{Particle-in-a-box.png/credit}}<br/>|}}Figure 2: A plot of <math>\psi_n(x)</math> for the first 7 wavefunctions, <math>n=1</math> (darkest) to <math>n=7</math> (lightest).]] | |||
We have now solved the particle-in-a-box problem. According to the Schrödinger wave equation the particle confined to the (infinite) box can only take on certain energy [[eigenvalues]], namely | We have now solved the particle-in-a-box problem. According to the Schrödinger wave equation the particle confined to the (infinite) box can only take on certain energy [[eigenvalues]], namely | ||
:<math> E_n=\frac{n^2\pi^2\hbar^2}{2mL^2}\ ,</math> | :<math> E_n=\frac{n^2\pi^2\hbar^2}{2mL^2}\ ,</math> | ||
Line 107: | Line 119: | ||
===The free space limit=== | ===The free space limit=== | ||
We can think of 'free space' as an infinitely wide box. Therefore if we take <math>L</math>, the width of the box, | We can think of [[Free space (electromagnetism)|'free space']] as an infinitely wide box. Therefore if we take <math>L</math>, the width of the box, | ||
to infinity we would expect to find that in the limit we recover the result of a [[free particle]]. In fact, | to infinity we would expect to find that in the limit we recover the result of a [[free particle]]. In fact, | ||
this is exactly what happens. | this is exactly what happens. | ||
Line 136: | Line 148: | ||
which is independent of <math>L</math>. This means that we would still find the particle if we looked in | which is independent of <math>L</math>. This means that we would still find the particle if we looked in | ||
''all of space'' after taking the limit <math>L\to\infty.</math> | ''all of space'' after taking the limit <math>L\to\infty.</math> | ||
==Non-stationary states== | |||
Eigenstates of the time-independent Schrödinger equation have a probability distribution that does not change with time, as <math>|\psi(x,t)|^2 = |\psi(x,0)|^2</math>. However, the probability distributions of wavefunctions that are sums of eigenfunctions with different energies do depend on time. In this section we will look at an example of this, and see how the energy eigenstates we found above allow us to solve for the evolution of the particle when it starts in an arbitrary state. | |||
By solving the time-independent Schrödinger equation above we found the energy eigenstates of the particle in a 1D box. As the eigenstates of the total Hamiltonian these states are [[stationary state|stationary]], meaning that they do not change over time. However we also saw that these states form a basis for the space of all possible states of the particle, so what happens when it is in a [[superposition]] of multiple energy states? | |||
Since such a state is no longer an eigenstate of the Hamiltonian we know that it won't be stationary, so we will have to use the time-''de''pendent Schrödinger equation. | |||
As a specific example suppose that at time <math>t=0</math> we know that the particle is in the superposition of the first two energy levels given by | |||
:<math>\Psi(x,t=0) = \frac{1}{\sqrt{2}}\Big(\psi_1(x)+\psi_2(x)\Big)</math> | |||
where the factor of <math>\scriptstyle 1/\sqrt{2}</math> is there to make the state normalized. | |||
Our goal now is to determine how the particle's state changes with time, i.e. to find <math>\scriptstyle\Psi(x,t)</math> for <math>t>0</math>. | |||
Since the Hamiltonian of the system is independent of time and the state is expressed in terms of energy eigenstates we can simply write down the wavefunction of the particle at time <math>t</math> as | |||
:<math> | |||
\Psi(x,t)=\frac{1}{\sqrt{2}}\left(e^{iE_1t/\hbar}\psi_1(x)+e^{iE_2t/\hbar}\psi_2(x)\right)\ . | |||
</math> | |||
Since we have already solved for the <math>E_n</math> and <math>\psi_n</math>, we can plug in the values from above to find the actual form of the solution. But first let's simplify our notation slightly. [[Dimensional analysis]] tells us that <math>\scriptstyle E_n/\hbar</math> must have units of [[frequency]], so let us define | |||
:<math> | |||
\omega_n = \frac{E_n}{\hbar} = \frac{n^2\pi^2\hbar}{2mL^2} | |||
</math> | |||
{{Image|Particle in a box (time evolution).gif|right|320px|Figure 3: The time-evolving probability density for a coherent superposition of the first two energy eigenfunctions (<math>n=1</math> and <math>n=2</math>).}} | |||
from which we can see that <math>\scriptstyle \omega_n=n^2\omega_1</math>. | |||
Now when we substitute the energy values and states into our time-evolving wavefunction we find | |||
:<math> | |||
\begin{align} | |||
\Psi(x,t) | |||
&=\frac{1}{\sqrt{2}}\left( | |||
e^{i\omega_1 t}\sqrt{\frac{2}{L}}\sin\left(\frac{\pi x}{L}\right) | |||
+e^{i\omega_2 t}\sqrt{\frac{2}{L}}\sin\left(\frac{2\pi x}{L}\right) | |||
\right) \\ | |||
&= | |||
\frac{e^{i\omega_1 t}}{\sqrt{L}} | |||
\left( | |||
\sin\left(\frac{\pi x}{L}\right) | |||
+ | |||
e^{3i\omega_1 t} | |||
\sin\left(\frac{2\pi x}{L}\right) | |||
\right) | |||
\end{align}\ . | |||
</math> | |||
To help us visualize what this wavefunction represents we will look at the probability density for the particle to be at a given location in the box as a function of time. This gives us | |||
:<math> | |||
\begin{align} | |||
p(x,t)&=\Psi^*(x,t)\Psi(x,t) \\ | |||
&= | |||
\frac{1}{L}\left( | |||
\sin\left(\frac{\pi x}{L}\right) | |||
+e^{-3i\omega_1 t}\sin\left(\frac{2\pi x}{L}\right) | |||
\right) | |||
\left( | |||
\sin\left(\frac{\pi x}{L}\right) | |||
+e^{3i\omega_1 t}\sin\left(\frac{2\pi x}{L}\right) | |||
\right) \\ | |||
&= | |||
\frac{1}{L} | |||
\left( | |||
\sin^2\left(\frac{\pi x}{L}\right) | |||
+\sin^2\left(\frac{2\pi x}{L}\right) | |||
+2\cos(3\omega_1 t)\sin\left(\frac{\pi x}{L}\right)\sin\left(\frac{2\pi x}{L}\right) | |||
\right)\ . | |||
\end{align} | |||
</math> | |||
This probability density is shown in Figure 3, from which we can see that the most likely place to find the particle oscillates back and forth across the box. This oscillation occurs at frequency <math>\omega_2-\omega_1=3\omega_1</math>. If the particle is charged, the oscillating charge density will produce electromagnetic radiation (light) with photon energy equal to the energy difference between the levels ''n'' = 1 and ''n'' = 2. | |||
The method used here can be generalized to solve for the evolution of ''any'' initial state, due to the completeness of the eigenstates of the Hamiltonian. An arbitrary initial state can be written as | |||
:<math> | |||
\Psi(x,0) | |||
= | |||
\sum_{n=1}^{\infty} | |||
c_n\psi_n(x) | |||
</math> | |||
where the [[expansion coefficients]] <math>c_n</math> are in general complex and satisfy the normalization condition | |||
:<math> | |||
\sum_{n=1}^\infty |c_n|^2=1\ . | |||
</math> | |||
The state of the particle at a later time can then be written simply as | |||
:<math> | |||
\Psi(x,t)=\sum_{n=1}^\infty c_n e^{iE_nt/\hbar}\psi_n(x)\ , | |||
</math> | |||
as was done in the example above. | |||
==Generalization to 3D (to be done)== | ==Generalization to 3D (to be done)== | ||
==References== | ==References== | ||
{{reflist}} | |||
Latest revision as of 11:31, 5 April 2011
![](http://s9.addthis.com/button1-share.gif)
The particle in a box or infinite square well problem is one of the simplest non-trivial solutions to Schrödinger's wave equation. As such it is often encountered in introductory quantum mechanics material as a demonstration of the quantization of energy. In its simplest form the problem is one-dimensional (1D), and involves a single particle living in an infinite potential well. What this means in the 1D case, which we examine here, is that there is a particle that is free to move back and forth along a line in between two points (e.g. 0 and , or and ). The particle is however constrained to this region, and cannot under any circumstances be found outside of it.
In this article we will rigorously define this problem and find the solutions to the corresponding Schrödinger wave equation, and then examine some of the properties and implications of these solutions.
The 1D square well
Setting up the problem
As depicted in Figure 1, the potential energy for this system can be given by
That is, in a region of length (the box or well) the potential is zero; everywhere else it is infinite. Since the total energy of the particle is given by the sum of its potential and kinetic energies, if it is in a region of infinite potential then the particle itself must have infinite energy, which of course cannot happen. An immediate consequence of this potential is therefore that the particle must be located somewhere between and .
This means that we only have to solve Schrödinger's equation inside the box and since the potential does not depend on time we can use the time-independent version of the Schrödinger equation. With in the box the wavefunction, , that describes the state of the particle must satisfy the differential equation (DE)
where is Planck's constant divided by 2π, and and are the mass and mechanical energy of the particle, respectively.
This is nothing but the eigenvalue problem, and our task is to determine the energy eigenvalue(s) and eigenstate(s) that solve it. To simplify the expressions, let
We shall later see that this choice of letter is justified, since k corresponds to the wave number of the solutions. We are now faced with a linear, second-order differential equation. The real solutions are most readily written as linear combinations of
where and are constants to be determined from conditions on the system. (The interested reader is invited to find the second derivative of this function and substitute it back into the differential equation above to verify that it is indeed satisfied.) The first condition we can use is that the wavefunction must be continuous. We therefore know that the particle cannot exist at any position , which tells us that the wavefunction must be zero at negative positions (in order for the probability of finding the particle there to also be zero). Continuity then implies that at the wavefunction is also zero, so we have
With we now know that the wavefunction is
Energy quantization
Using the same continuity argument at tells us that
One solution to this is of course . However, this would mean that the wavefunction vanishes everywhere -- implying that there is no particle! The other way to satisfy this equality is to have the sine term vanish, which will happen if
where mathematically can be any integer. This problem is not purely mathematical though, and we know for physical reasons that , , , and are all greater than zero. This means that must also be greater than zero, so is restricted to being a natural number, . We can now solve for the particle's energy,
where we have labelled the energy by the integer . We have just derived energy quantization! Without the potential well, i.e. if the particle was free, its energy would be allowed to take on any real number. Once inside the box though, only a specific discrete set of energy eigenvalues is permitted.
Substituting back into , we now have an infinite number of possible wavefunctions for the particle,
The final step is to determine the value of . This requires another condition that we can impose upon the problem.
Normalization
The probability of finding the particle somewhere must be . Since the wavefunction vanishes everywhere outside , the normalization condition is
Solving for tells us , where is some overall quantum phase (a real number). Because overall phases do not affect measurable results we are free to choose any value we want for without affecting the physics. For simplicity then, we set .
The solutions
We have now solved the particle-in-a-box problem. According to the Schrödinger wave equation the particle confined to the (infinite) box can only take on certain energy eigenvalues, namely
and the wavefunction for the particle when it is in the th energy eigenstate is
inside the box, and outside it.
Figure 2 shows a plot of the first seven solutions for the wavefunction of a particle in a box. Note how as increases the number of waves that fit in the box increases as well, meaning that the wavelength has decreased and therefore frequency has increased. As with electromagnetic waves a higher frequency corresponds to a higher energy, in agreement with the expression for .
Alternate description
Another common (equivalent) way to describe this problem is to shift the location of the box, so that it runs from to instead of to . The derivation of the allowed wavefunctions is very similar to what we did here, or the results can be obtained via coordinate transformation (specifically, a translation). With this different description of the square well the energy eigenvalues remain unchanged (as is to be expected[1]) while the resulting wavefunctions are given by[2]
if is odd, and
if is even.
Properties of the particle in a box
Completeness
The are complete, so form a basis for the possible wavefunctions, or states, of the particle. This means that an arbitrary state, , can be expressed as
where is the probability of measuring the particle's energy to be , and therefore finding the particle in state .
Non-vanishing ground state energy
The particle is not allowed to have zero energy, because if then , which makes the wavefunction
But if for all then there is zero probability of finding the particle no matter where you look -- which means that there is no particle! What this implies then is that the lowest allowed energy or ground state energy for the particle in the box is
The free space limit
We can think of 'free space' as an infinitely wide box. Therefore if we take , the width of the box, to infinity we would expect to find that in the limit we recover the result of a free particle. In fact, this is exactly what happens.
Mathematically we would find the same result for either of the solutions above, but physically free space doesn't have an infinite potential wall at . It is therefore more intuitive to look at the second solution, in which the potential well extends from to .
Notice that the energy spacing between two consecutive energy levels is
It is easy to see that for any finite value of the value of goes to zero as goes to infinity. This tells us that a free particle can have any (real, positive) energy.
Perhaps surprisingly the wavefunction also appears to go to zero everywhere in the limit, which at first may appear to be a problem. As discussed above if the wavefunction vanishes everywhere then there shouldn't be a particle. The difference in this case is subtle, and has to do with the fact that we are dealing with infinities. The important point to note is that if the particle is in state then the probability to find it if we look in the entire box is given by
which is independent of . This means that we would still find the particle if we looked in all of space after taking the limit
Non-stationary states
Eigenstates of the time-independent Schrödinger equation have a probability distribution that does not change with time, as . However, the probability distributions of wavefunctions that are sums of eigenfunctions with different energies do depend on time. In this section we will look at an example of this, and see how the energy eigenstates we found above allow us to solve for the evolution of the particle when it starts in an arbitrary state.
By solving the time-independent Schrödinger equation above we found the energy eigenstates of the particle in a 1D box. As the eigenstates of the total Hamiltonian these states are stationary, meaning that they do not change over time. However we also saw that these states form a basis for the space of all possible states of the particle, so what happens when it is in a superposition of multiple energy states? Since such a state is no longer an eigenstate of the Hamiltonian we know that it won't be stationary, so we will have to use the time-dependent Schrödinger equation.
As a specific example suppose that at time we know that the particle is in the superposition of the first two energy levels given by
where the factor of is there to make the state normalized. Our goal now is to determine how the particle's state changes with time, i.e. to find for . Since the Hamiltonian of the system is independent of time and the state is expressed in terms of energy eigenstates we can simply write down the wavefunction of the particle at time as
Since we have already solved for the and , we can plug in the values from above to find the actual form of the solution. But first let's simplify our notation slightly. Dimensional analysis tells us that must have units of frequency, so let us define
from which we can see that . Now when we substitute the energy values and states into our time-evolving wavefunction we find
To help us visualize what this wavefunction represents we will look at the probability density for the particle to be at a given location in the box as a function of time. This gives us
This probability density is shown in Figure 3, from which we can see that the most likely place to find the particle oscillates back and forth across the box. This oscillation occurs at frequency . If the particle is charged, the oscillating charge density will produce electromagnetic radiation (light) with photon energy equal to the energy difference between the levels n = 1 and n = 2.
The method used here can be generalized to solve for the evolution of any initial state, due to the completeness of the eigenstates of the Hamiltonian. An arbitrary initial state can be written as
where the expansion coefficients are in general complex and satisfy the normalization condition
The state of the particle at a later time can then be written simply as
as was done in the example above.
Generalization to 3D (to be done)
References
- ↑ The energy of a particle should not depend in any way on our choice of coordinate system, i.e. where we choose to place the origin .
- ↑ A. Messiah (1999) Quantum Mechanics: Two volumes bound as one ISBN 0486409244