Substitution I don't get

In summary: \frac{-\alpha q'^2}{2}+iq'x - \frac{x^2}{2\alpha} + \frac{x^2}{2\alpha} =... \frac{-\alpha q'^2}{2}-\frac{\alpha x^2}{2\alpha} + iq'(x + \frac{x}{\alpha}) + \frac{x^2}{2\alpha} =...\frac{-\alpha q'^2}{2}-\frac{\alpha x^2}{2\alpha} + iq'(x + \frac{x}{\alpha} + \frac{1}{2\alpha} x^2 - \frac{1}{2\alpha}x^
  • #1
quasar987
Science Advisor
Homework Helper
Gold Member
4,807
32
Yes, it's me and the wave packets... again!

This is taken from the text of Gasiorowicz's Quantum Physics 3rd ed. pp.26.

We have a gaussian wave packet at t=0 that is is described by

[tex]\psi(x,0)=\int_{-\infty}^{\infty}dke^{-\alpha (k-k_0)^2/2}e^{ikx}[/tex]

and we apply the change of variable [itex]q' = k-k_0[/itex]. Suposedly the wave packet becomes

[tex]\psi(x,0)=e^{ik_0x}e^{-x^2/2\alpha}\int_{-\infty}^{\infty}dq'e^{-\alpha q'^2/2}[/tex]

How does one gets to that? When I make the substitution k = q' + k_0, I get

[tex]e^{ik_0x}\int_{-\infty}^{\infty}dq' e^{-\alpha q'^2/2}e^{iq'x}[/tex]
 
Physics news on Phys.org
  • #2
Combine the exponentials under the integral and complete the square (for q') in the resulting exponent.

Regards,
George
 
  • #3
I really don't like that Gasiorowicz book. It's the one that's purple and yellow on the front right? When substituting I get the same thing you do. It seems that he's "pulling out" the [itex]e^{iq'x}[/itex] from the integral and it "becomes" [itex]e^{-x^2/2\alpha }[/itex]. Let's simplify some thing. Let b = -a/2, and let y be used where x was being used, and let x be used where q' was being used. The integrand is then:

dx ebx2 + ixy

We want to see why:

[tex]\int _{-\infty } ^{\infty } dx\, \exp (bx^2 + ixy) = \exp \left (\frac{y^2}{4b}\right ) \int _{-\infty } ^{\infty }dx\, \exp (bx^2)[/tex]

Using Wolfram's "The Integrator" we get that the left side of the equation in indefinite form is the top thing in the attached image, and the right side is the bottom thing.

This is exactly what we want if

[tex]\lim _{x \to \infty} \left [\mbox{Erfi}\left (\frac{2bx + iy}{2\sqrt{b}}\right ) - \mbox{Erfi}\left (\frac{2b(-x) + iy}{2\sqrt{b}}\right )\right ] = \lim _{x \to \infty} \left [\mbox{Erfi}(\sqrt{b}x ) - \mbox{Erfi} (\sqrt{b}(-x))\right ][/tex]

[tex]\lim _{x \to \infty} \left [\mbox{Erfi}\left (\sqrt{b}x + \frac{iy}{2\sqrt{b}}\right ) - \mbox{Erfi}\left (\sqrt{b}(-x) + \frac{iy}{2\sqrt{b}}\right )\right ] = \lim _{x \to \infty} \left [\mbox{Erfi}(\sqrt{b}x ) - \mbox{Erfi} (\sqrt{b}(-x))\right ][/tex]

This would seem plausible because iy/2√b would just be a constant, but it's a complex number, and from what I could tell, Erfi takes real arguments only. I have no idea why the Integrator would put the "i" in the argument, I don't always trust that application.
 

Attachments

  • integral.JPG
    integral.JPG
    4.6 KB · Views: 376
Last edited:
  • #4
Ouch, yes, completing the square is much less involved. However, I still run into a problem, namely that I end up with:

[tex]\psi (x, 0) = e^{ik_0x}e^{-x^2/2\alpha }\int _{-\infty } ^{\infty} dq' \, e^{-\alpha \mathbf{(q' - xi/\alpha )}^2/2}[/tex]

instead of:

[tex]\psi (x, 0) = e^{ik_0x}e^{-x^2/2\alpha }\int _{-\infty } ^{\infty} dq' \, e^{-\alpha \mathbf{q}'^2/2}[/tex]

If [itex]i[/itex] weren't imaginary, then a simple change of variables from q' to q' - xi/a (or something like that) would make everything okay. I haven't done any complex analysis - is such a substitution valid? If q' goes from [itex]-\infty [/itex] to [itex]\infty [/itex], would it make sense to say that q' - xi/a ranges over the line in the complex plane {c + di : d = -x/a}?
 
  • #6
AKG said:
Ouch, yes, completing the square is much less involved. However, I still run into a problem, namely that I end up with:

[tex]\psi (x, 0) = e^{ik_0x}e^{-x^2/2\alpha }\int _{-\infty } ^{\infty} dq' \, e^{-\alpha \mathbf{(q' - xi/\alpha )}^2/2}[/tex]

instead of:

[tex]\psi (x, 0) = e^{ik_0x}e^{-x^2/2\alpha }\int _{-\infty } ^{\infty} dq' \, e^{-\alpha \mathbf{q}'^2/2}[/tex]

I was going to post that.

Edit: Actually, no. Instead of

[tex]\psi (x, 0) = e^{ik_0x}e^{-x^2/2\alpha }\int _{-\infty } ^{\infty} dq' \, e^{-\alpha \mathbf{(q' - xi/\alpha )}^2/2}[/tex]

I got

[tex]\psi (x, 0) = e^{ik_0x}e^{-x^2/2\alpha }\int _{-\infty } ^{\infty} dq' \, e^{-\alpha \mathbf{(iq' + x/\alpha )}^2/2}[/tex]
 
Last edited:
  • #7
-a(iq + x/a)²/2 = -a(iq + ix/ia)²/2 = -a(i(q - ix/a))²/2 = -ai²(q - ix/a)²/2 = a(q - ix/a)²/2, so your exponent is just the negative of mine. However, I'm pretty sure I completed the square correctly. Also, when you complete the square, say, in the variable q, then you want to end up with something of the form A(q - B)² + C, you shouldn't have iq or anything like that (you should have factored the i out, but there was no i to factor). Actually, you should be able to see that yours is wrong, because when you square it, the q² term will be -a(iq)² = aq², but we know that the original thing we were working with had -aq² as the q² term.
 
  • #8
  • #9
AKG said:
Is that supposed to be a reason for why the "i" is in there, or is that supposed to be reason to distrust Integrator?

About the trust part.
 
  • #10
AKG said:
-a(iq + x/a)²/2 = -a(iq + ix/ia)²/2 = -a(i(q - ix/a))²/2 = -ai²(q - ix/a)²/2 = a(q - ix/a)²/2, so your exponent is just the negative of mine. However, I'm pretty sure I completed the square correctly. Also, when you complete the square, say, in the variable q, then you want to end up with something of the form A(q - B)² + C, you shouldn't have iq or anything like that (you should have factored the i out, but there was no i to factor). Actually, you should be able to see that yours is wrong, because when you square it, the q² term will be -a(iq)² = aq², but we know that the original thing we were working with had -aq² as the q² term.

I wrote it wrong, sorry! Forget the negative sign.

The way I completed the square is

[tex]\frac{-\alpha q'^2}{2}+iq'x = \left(i\sqrt{\frac{\alpha}{2}}q'+\frac{x}{\sqrt{2\alpha}}\right)^2-\frac{x^2}{2\alpha}[/tex]
 
  • #11
You normally don't see people complete the square as you have. Normally they do this:

Aq² + Bq + C
= A(q² + B/Aq) + C
= A[q² + B/Aq + (B/2A)² - (B/2A)²] + C
= A(q + B/2A)² + (C - B²/4A)

so inside the bracket, the "q" is not multiplied by any factor (the factor appears outside the bracket). However, what you have now is technically correct, but now what do we do with it? Since it's the same thing I have in post 4, how do we get the "-xi/a" to disappear?
 
  • #12
Sorry - from eyeballing the problem I thought completing the square would do the trick. AKG's substitution can be justified by using the residue theorem and an appropriate contour of integration.

Consider

[tex]
\int dz e^{-\frac{\alpha}{2} \left( z - ib \right)^2}
[/tex]

for complex values [itex]z = x + iy[/itex].

Because the integrand in analytic, the residue theorem give that this integral is zero around any closed path (contour) in the complex plane. Consider the closed contour that consists of a horizontal path [itex]z = -x + 0i[/itex] for [itex]x = -a[/itex] to [itex]x = a[/itex], a vertical path [itex]z = a + iy[/itex] for [itex]y = 0[/itex] to [itex]y = b[/itex], a horizontal path [itex]z = x + ib[/itex] for [itex]x = a[/itex] to [itex]x = -a[/itex], and a vertical path [itex]z = a + iy[/itex] for [itex]y = b[/itex] to [itex]y = 0[/itex].

Along the bottom horizontal path, [itex]dz = dx[/itex], and if we let [itex]a \rightarrow[/itex] \infty[/itex], we get the required integral. Along the to horizontal path, [itex]dz = dx[/itex], and if we let [itex]a \rightarrow[/itex] \infty[/itex], we get AKG's substitution. If it can be shown that the integrals along the 2 vertical paths die as [itex]a \rightarrow[/itex] \infty[/itex], then the "substitution" is valid.

Along the right vertical path, [itex]dz = dy[/itex], and the integral is

[tex]
\int_{0}^{b} dy e^{-\frac{\alpha}{2} \left( a + iy - ib \right)^2}.
[/tex]

Square the brackets and factor out the

[tex]
e^{-\frac{\alpha}{2} a^2}.
[/tex]

This dies fast enough to kill the integral. Same thing on the left.

Regards,
George
 

1. What is substitution in science?

Substitution in science refers to the process of replacing one element or compound with another in a chemical reaction or equation. It is a fundamental concept in chemistry and is often used to predict or explain the behavior of different substances.

2. How does substitution work?

In a substitution reaction, the new element or compound takes the place of the original one in the reaction. This can happen through a single or double displacement, depending on the specific reaction. The end result is a new substance or compounds with different properties.

3. What is the difference between substitution and replacement?

Substitution and replacement are often used interchangeably, but they have slight differences in meaning. Substitution specifically refers to the replacement of one element or compound with another in a chemical reaction, while replacement can also refer to the replacement of a part or component in a larger system.

4. What are some real-life examples of substitution in science?

Substitution reactions are common in many everyday processes, such as the use of bleach to remove stains (a substitution reaction between the active ingredient, sodium hypochlorite, and the stain molecules). In biology, enzymes facilitate substitution reactions in the body to break down and convert molecules for energy.

5. How is substitution used in organic chemistry?

In organic chemistry, substitution reactions are used to create new compounds by replacing a functional group or atom with another. This can result in the creation of new compounds with different properties, such as different solubility or reactivity. Substitution reactions are also used in the synthesis of pharmaceuticals and other organic compounds.

Similar threads

  • Introductory Physics Homework Help
Replies
2
Views
1K
  • Introductory Physics Homework Help
Replies
6
Views
207
  • Introductory Physics Homework Help
Replies
3
Views
781
  • Introductory Physics Homework Help
Replies
4
Views
751
  • Introductory Physics Homework Help
Replies
1
Views
1K
  • Quantum Physics
Replies
1
Views
525
  • Introductory Physics Homework Help
Replies
10
Views
1K
  • Introductory Physics Homework Help
Replies
9
Views
3K
  • Introductory Physics Homework Help
Replies
15
Views
1K
  • Calculus and Beyond Homework Help
Replies
3
Views
665
Back
Top