Dirac delta function and Heaviside step function

pedroobv
Messages
9
Reaction score
0
[SOLVED] Dirac delta function and Heaviside step function

In Levine's Quantum Chemistry textbook the Heaviside step function is defined as:

H(x-a)=1,x>a
H(x-a)=0,x<a
H(x-a)=\frac{1}{2},x=a​

Dirac delta function is:

\delta (x-a)=dH(x-a) / dx​

Now, the integral:

\int ^{\infty}_{-\infty}f(x)\delta (x-a)dx​

Is evaluated using integration by parts considering

u=f(x), du=f'(x)
dv=\delta (x-a)dx, v=H(x-a)​

We have then:
\int ^{\infty}_{-\infty}f(x)\delta (x-a)dx=f(x)H(x-a)|^{\infty}_{-\infty}-\int ^{\infty}_{-\infty}H(x-a)f'(x)dx

\int ^{\infty}_{-\infty}f(x)\delta (x-a)dx=f(\infty)-\int ^{\infty}_{-\infty}H(x-a)f'(x)dx​

Since H(x-a) vanishes for x<a, the integral becomes:

\int ^{\infty}_{-\infty}f(x)\delta (x-a)dx=f(\infty)-\int ^{\infty}_{a}H(x-a)f'(x)dx=f(\infty)-\int ^{\infty}_{a}f'(x)dx​

This is the point where my question arrives. H(x-a) is considered to have a value of unity for all the integral and that's why it is pulled out of the integral as a constant, however the lower bound of the integral is a and in this point H(x-a)=1/2. Could you please tell me if the following explanation is correct?

I think that because in all the integral, except in a, H(x-a)=1 and since the upper bound is infinity the value of the integral at the point a can be ignored.

If I'm wrong, any suggestion for correcting my explanation will be appreciated.
 
Physics news on Phys.org
You can always ignore the value of a function at a single point when you are integrating. If you remember how Riemann sums work, the contribution of the single point can be put into a rectangle of height 1/2 and width 0. So it makes no contribution.
 
Ok, thank you very much. That answers my question.
 
##|\Psi|^2=\frac{1}{\sqrt{\pi b^2}}\exp(\frac{-(x-x_0)^2}{b^2}).## ##\braket{x}=\frac{1}{\sqrt{\pi b^2}}\int_{-\infty}^{\infty}dx\,x\exp(-\frac{(x-x_0)^2}{b^2}).## ##y=x-x_0 \quad x=y+x_0 \quad dy=dx.## The boundaries remain infinite, I believe. ##\frac{1}{\sqrt{\pi b^2}}\int_{-\infty}^{\infty}dy(y+x_0)\exp(\frac{-y^2}{b^2}).## ##\frac{2}{\sqrt{\pi b^2}}\int_0^{\infty}dy\,y\exp(\frac{-y^2}{b^2})+\frac{2x_0}{\sqrt{\pi b^2}}\int_0^{\infty}dy\,\exp(-\frac{y^2}{b^2}).## I then resolved the two...
Hello everyone, I’m considering a point charge q that oscillates harmonically about the origin along the z-axis, e.g. $$z_{q}(t)= A\sin(wt)$$ In a strongly simplified / quasi-instantaneous approximation I ignore retardation and take the electric field at the position ##r=(x,y,z)## simply to be the “Coulomb field at the charge’s instantaneous position”: $$E(r,t)=\frac{q}{4\pi\varepsilon_{0}}\frac{r-r_{q}(t)}{||r-r_{q}(t)||^{3}}$$ with $$r_{q}(t)=(0,0,z_{q}(t))$$ (I’m aware this isn’t...
It's given a gas of particles all identical which has T fixed and spin S. Let's ##g(\epsilon)## the density of orbital states and ##g(\epsilon) = g_0## for ##\forall \epsilon \in [\epsilon_0, \epsilon_1]##, zero otherwise. How to compute the number of accessible quantum states of one particle? This is my attempt, and I suspect that is not good. Let S=0 and then bosons in a system. Simply, if we have the density of orbitals we have to integrate ##g(\epsilon)## and we have...
Back
Top