Delta function potential; Schrodinger Equation

In summary, the TISE for a particle of mass m moving along the x-axis and interacting with an attractive delta function potential at the origin is given by the following parts: Part (a) discusses the difference between a bound state particle and a free particle, where a bound particle has energy E<0 while a continuum eigenfunction has energy E>0. Part (b) shows that the solution for the TISE is given by ##\psi_{(x)} = exp(-|k|x)##. Part (c) explains why the symmetric solution takes the form of ##\psi = exp(-k|x|)## for |x|>a and ##\psi \propto exp(kx) + exp(-kx
  • #1
unscientific
1,734
13

Homework Statement



a0cmrk.png


Consider the TISE for a particle of mass m moving along the x-axis and interacting an attractive delta function potential at origin:
Part(a): What is the difference between a bound state particle and a free particle?

Part(b): Show ##\psi _{(x)} = exp (-|k|x)## is a solution of TISE, given conditions:

Part(c): Explain why symmetric solution is of form:

Part(d): Use boundary conditions to show the following:

Part(e): Investigate this solution graphically

Part(f): Find the significance of this result

Homework Equations


The Attempt at a Solution



Part(a)

A bound particle as energy E<0 while a continuum eigenfunction gives energy E>0.

Part (b)

TISE is given by:

[tex]\hat {H}\psi = E\psi[/tex]

Using first boundary condition provided we get

[tex] 2k = \frac{\lambda}{a}[/tex]

Substituting ##\psi = e^{-k|x|}## into TISE and using what we have just found:

[tex]-\frac{\hbar ^2}{2m}(\frac{\lambda}{2a})^2 - \frac{\hbar ^2 \lambda}{2ma}\delta_{(x)} = E[/tex]

I'm not sure why the delta function goes to zero here?

[tex]E = -\frac{\hbar ^2}{2m}(\frac{\lambda}{2a})^2 [/tex]

Part (c)

For |x|>a, V = 0:
[tex]-\frac{\hbar^2}{2m}\frac{\partial ^2 \psi}{\partial x^2} = E\psi[/tex]

Solving this equation gives an exponential ##\psi = exp (-k|x|)##.

For -a<x<a, the particle is bounded and symmetric about origin so ##\psi \propto exp(kx) + exp(-kx)## which is ##\propto cosh(kx) = A cosh(kx)##.

Part (d)

Boundary Conditions:

At x = a, ##\psi_{|x|<a} = \psi_{|x|>a}##.
Thus we get:

[tex] e^{-ka} = A cosh (ka)[/tex]

Second boundary condition:
At x = a, ##\frac{\partial \psi}{\partial x}_{left} = - \frac{\partial \psi}{\partial x}_{right}##

The delta function gives opposite gradients about a.

[tex]e^{-ka} = A sinh (ka) [/tex]

Using boundary condition earlier,

[tex] tanh (ka) = 1[/tex]
[tex] tanh(ka) = 2 - 1 [/tex]

Can I simply substitute in ##2k = \frac{\lambda}{a}## here?

[tex] tanh(ka) = \frac{2ka}{ka} - 1 [/tex]
[tex] tanh (ka) = \frac{\lambda}{y} - 1 [/tex]

Part (e)

i71miw.png


We can see that ##y_0 < \lambda##.

Since ##tanh(y_0) < 1##, thus ##y_0 > \frac{\lambda}{2}##.

Since for x≠a both delta functions are positive, it implies that the second term corresponding to the two delta functions ##\delta_{(x-a)}## and ##\delta_{(x+a)}## in the schrodinger equation are >0. Thus the total energy subtract that is less than ##-\frac{\hbar ^2\lambda^2}{8ma^2}##.

Part(f)

Two delta function corresponds to electron being attracted to the two protons in nucleus (which is more massive so can treat as stationary). This is a more stable (favourable) configuration, so

[tex]H + H^+ \rightarrow H_2^+[/tex]
 
Physics news on Phys.org
  • #2
unscientific said:

Homework Statement



a0cmrk.png


Consider the TISE for a particle of mass m moving along the x-axis and interacting an attractive delta function potential at origin:
Part(a): What is the difference between a bound state particle and a free particle?

Part(b): Show ##\psi _{(x)} = exp (-|k|x)## is a solution of TISE, given conditions:

Part(c): Explain why symmetric solution is of form:

Part(d): Use boundary conditions to show the following:

Part(e): Investigate this solution graphically

Part(f): Find the significance of this result

Homework Equations


The Attempt at a Solution



Part(a)

A bound particle as energy E<0 while a continuum eigenfunction gives energy E>0.

Part (b)

TISE is given by:

[tex]\hat {H}\psi = E\psi[/tex]

Using first boundary condition provided we get

[tex] 2k = \frac{\lambda}{a}[/tex]

Substituting ##\psi = e^{-k|x|}## into TISE and using what we have just found:

[tex]-\frac{\hbar ^2}{2m}(\frac{\lambda}{2a})^2 - \frac{\hbar ^2 \lambda}{2ma}\delta_{(x)} = E[/tex]

I'm not sure why the delta function goes to zero here?

[tex]E = -\frac{\hbar ^2}{2m}(\frac{\lambda}{2a})^2 [/tex]
You're verifying the solution for x>0 and x<0. At x=0, there's a cusp, so ##\psi''## isn't defined, so the Schrodinger equation doesn't really apply.

Part (c)

For |x|>a, V = 0:
[tex]-\frac{\hbar^2}{2m}\frac{\partial ^2 \psi}{\partial x^2} = E\psi[/tex]

Solving this equation gives an exponential ##\psi = exp (-k|x|)##.

For -a<x<a, the particle is bounded and symmetric about origin so ##\psi \propto exp(kx) + exp(-kx)## which is ##\propto cosh(kx) = A cosh(kx)##.

Part (d)

Boundary Conditions:

At x = a, ##\psi_{|x|<a} = \psi_{|x|>a}##.
Thus we get:

[tex] e^{-ka} = A cosh (ka)[/tex]

Second boundary condition:
At x = a, ##\frac{\partial \psi}{\partial x}_{left} = - \frac{\partial \psi}{\partial x}_{right}##

The delta function gives opposite gradients about a.

[tex]e^{-ka} = A sinh (ka) [/tex]
This isn't correct. Integrate the Schrodinger equation between ##x=a-\varepsilon## and ##x=a+\varepsilon## and take the limit as ##\varepsilon \to 0^+##. This will yield the proper boundary condition on the derivatives.

Using boundary condition earlier,

[tex] tanh (ka) = 1[/tex]
[tex] tanh(ka) = 2 - 1 [/tex]

Can I simply substitute in ##2k = \frac{\lambda}{a}## here?

[tex] tanh(ka) = \frac{2ka}{ka} - 1 [/tex]
[tex] tanh (ka) = \frac{\lambda}{y} - 1 [/tex]

Part (e)

i71miw.png


We can see that ##y_0 < \lambda##.

Since ##tanh(y_0) < 1##, thus ##y_0 > \frac{\lambda}{2}##.

Since for x≠a both delta functions are positive, it implies that the second term corresponding to the two delta functions ##\delta_{(x-a)}## and ##\delta_{(x+a)}## in the schrodinger equation are >0. Thus the total energy subtract that is less than ##-\frac{\hbar ^2\lambda^2}{8ma^2}##.

Part(f)

Two delta function corresponds to electron being attracted to the two protons in nucleus (which is more massive so can treat as stationary). This is a more stable (favourable) configuration, so

[tex]H + H^+ \rightarrow H_2^+[/tex]
 
  • #3
vela said:
You're verifying the solution for x>0 and x<0. At x=0, there's a cusp, so ##\psi''## isn't defined, so the Schrodinger equation doesn't really apply.

So the equation##exp(-k|x|)## here only applies when x≠0. We can't find out what the solution is when x = 0? Is it because due to the nature of the delta function it is a sum of infinite series at x = 0?


vela said:
This isn't correct. Integrate the Schrodinger equation between ##x=a-\varepsilon## and ##x=a+\varepsilon## and take the limit as ##\varepsilon \to 0^+##. This will yield the proper boundary condition on the derivatives.

But graphically you can see that it is symmetric about 0. You can also see that at x-ε the gradient is opposite that of at x+ε.
 
  • #4
unscientific said:
So the equation##exp(-k|x|)## here only applies when x≠0. We can't find out what the solution is when x = 0? Is it because due to the nature of the delta function it is a sum of infinite series at x = 0?
At x=0, the lefthand side of the Schrodinger equation isn't defined, so it doesn't really tell you anything about ##\psi##. The continuity requirement, however, allows you to fill in the hole at x=0.

But graphically you can see that it is symmetric about 0. You can also see that at x-ε the gradient is opposite that of at x+ε.
It may look that way in the picture when you have a single delta function, but this time you have two. You no longer have the symmetry you had before. You need to justify the boundary condition mathematically.
 
  • #5


This result shows that the delta function potential is able to represent the interaction between the electron and protons in the nucleus, leading to the formation of a stable molecule. It also shows that the symmetric solution, with its characteristic exponential behavior, is able to accurately describe the behavior of the particle in this potential.
 

FAQ: Delta function potential; Schrodinger Equation

What is a delta function potential?

A delta function potential is a type of potential energy function in quantum mechanics that describes a point-like interaction between a particle and a potential energy barrier. It is represented by the Dirac delta function, which is a mathematical function that is infinite at a single point and zero everywhere else.

How does the delta function potential affect the Schrodinger equation?

The presence of a delta function potential in the Schrodinger equation causes a discontinuity in the wavefunction at the point of interaction. This results in a change in the probability of finding the particle at that point, as well as a change in the overall behavior of the wavefunction.

What are some applications of the delta function potential?

The delta function potential is commonly used in quantum mechanics to model point-like interactions, such as the interaction between an electron and a nucleus in an atom. It is also used in scattering problems, where it can represent a localized potential barrier or a point-like scattering center.

Can the delta function potential be solved analytically?

Yes, in some cases the delta function potential can be solved analytically. This is particularly true for simple potential configurations, such as a delta function potential well or a delta function potential barrier. In more complex systems, numerical methods may be needed to solve the Schrodinger equation.

How does the strength of the delta function potential affect the behavior of the wavefunction?

The strength of the delta function potential directly affects the amplitude of the discontinuity in the wavefunction. A stronger potential will result in a larger discontinuity and a greater change in the behavior of the wavefunction. This can lead to changes in the energy levels and probabilities of finding the particle in different regions of space.

Back
Top