26 svar
388 visningar
H-lena 24 – Fd. Medlem
Postad: 2 mar 2017 18:45

Integrera med deltafunktion

Hej!

Jag vill integrera följande funktion 

 

0δ(r)rsin(qr)dr

där q är en konstant och δ(r) är deltafunktionen. Jag vet att δ(r)dr=1 

men nu har jag ju undre gränsen som 0, hur blir det då?

 

Tacksam för svar!

Freewheeling 220 – Fd. Medlem
Postad: 2 mar 2017 19:41

I allmänhet gäller det att:

0f(x)δ(x)dx = - f(x)δ(x)dx

H-lena 24 – Fd. Medlem
Postad: 2 mar 2017 19:56

jaha okej tack. Så då kan jag räkna integralen med gränserna -rδ(r)sin(qr)dr ?

Freewheeling 220 – Fd. Medlem
Postad: 2 mar 2017 21:43 Redigerad: 2 mar 2017 21:45

Hej,

När jag tänker efter så har jag även vid något tillfälle stött på definitionen:

0f(x)δ(x)dx = f(0)2

Det har att göra med att man kan betrakta Dirac-pulsen som gränsvärdet av en jämn funktion vn(x), vilket skulle ge:

f(0) =-f(x)δ(x)dx = limn-f(x)vn(x)dx = limn20f(x)vn(x)dx = 20f(x)δ(x)dx

Det samband jag gav ovan i inlägget bygger på någon slags generalisering tror jag, där man definierar

Sf(x)δ(x)dx = f(0) om 0S0 annars

Min poäng är att Dirac delta-funktionen verkar vara någonting knepigt i grund och botten och det kan vara värt att kolla upp exakt hur ni definierar saker och ting i kursen.

Ett ytterligare klargörande är att jag i båda fallen antar jag att vi integrerar över intervallet [0,∞) och inte över (0,∞) vilket skulle innebära att integralen är noll. Så om detta på något vis framgår av uppgiften så är det värt att ta hänsyn till.

Edit: Usch, vad jag sabbade till det med tecknen innanför integralerna. Hoppas det går att läsa...

emmynoether 663 – Fd. Medlem
Postad: 2 mar 2017 22:31 Redigerad: 2 mar 2017 22:37

Du måste göra ett koordinatbyte av något slag, du får inte använda dig av diracs deltafunktion i en ändpunkt precis som du säger eftersom den måste integreras över ett område. Deltafunktionen är egentligen ingen funktion, det är en distribution. Testa gör ett byte till kartesiska koordinater och utför integralen där.

emmynoether 663 – Fd. Medlem
Postad: 2 mar 2017 22:34 Redigerad: 2 mar 2017 23:25
H-lena skrev :

jaha okej tack. Så då kan jag räkna integralen med gränserna -rδ(r)sin(qr)dr ?

 Nej så kan du inte göra.

 

EDIT: Jag tänkte på det lite till, det finns tillfällen då du kan utvidga intervallet och dela med två och sedan göra på det där viset men det fungerar inte allmänt. Det beror helt på vilket problem du försöker lösa.

H-lena 24 – Fd. Medlem
Postad: 3 mar 2017 09:23

Tack för svar. Jag känner mig lite förvirrad över den här integralen just nu måste jag säga. I kartesiska koordinater ser den ut såhär:

 

-δ(x-r)eiqxd3x med lite konstanter framför. Det är en Born approximation av första ordningen om det säger nånting. Det ges inga gränser i uppgiften så därför antar jag att det är från - oändligheten till +oändligheten. 

emmynoether 663 – Fd. Medlem
Postad: 3 mar 2017 11:24

Jobbar du alltså i sfäriska koordinater från början? Är det en spridning mot en Dirac potential som ligger i origo i sfäriska koordinater? 

H-lena 24 – Fd. Medlem
Postad: 3 mar 2017 14:12 Redigerad: 3 mar 2017 14:12

Det är en spridning mot en dirac potential antar jag den ser ut på följande vis

V(r)=αδ(r), α>0 det är allt som ges i uppgiften, man ska hitta scattering cross section, så jag använde Born-approx. första ordningen med ekvationen och tänkte använda absolutbeloppet av kvadraten f(k,k') efter att jag räknat ut den. 

 

f(k,k')=-m2πhd3x'V(x')ei(k-k')x' där k-k'=q (Det ska stå h streck i kvadrat inte h)

Jag vet inte hur jag ska lösa den integralen med den lilla information jag har, därför ville jag skriva om den med sfäriska koordinater som man kan göra enligt

f(θ)=-m2πh2q0rV(r)sin(qr)dr (h=h streck)

men det gick ju inte så lätt som jag trodde. 

Freewheeling 220 – Fd. Medlem
Postad: 3 mar 2017 17:26
emmynoether skrev :

Du måste göra ett koordinatbyte av något slag, du får inte använda dig av diracs deltafunktion i en ändpunkt precis som du säger eftersom den måste integreras över ett område. Deltafunktionen är egentligen ingen funktion, det är en distribution. Testa gör ett byte till kartesiska koordinater och utför integralen där.

Man kan inte i allmänhet resonera som så att:

-f(x)δ(x)dx = Af(x)δ(x)dx + Bf(x)δ(x)dx

Där A = (-∞,0) och B = [0,∞) och utnyttja att integralen över A är noll?

emmynoether 663 – Fd. Medlem
Postad: 3 mar 2017 22:45
Freewheeling skrev :
emmynoether skrev :

Du måste göra ett koordinatbyte av något slag, du får inte använda dig av diracs deltafunktion i en ändpunkt precis som du säger eftersom den måste integreras över ett område. Deltafunktionen är egentligen ingen funktion, det är en distribution. Testa gör ett byte till kartesiska koordinater och utför integralen där.

Man kan inte i allmänhet resonera som så att:

-f(x)δ(x)dx = Af(x)δ(x)dx + Bf(x)δ(x)dx

Där A = (-∞,0) och B = [0,∞) och utnyttja att integralen över A är noll?

 Nej så funkar det inte. Problemet är att deltafunktionen är en distribution som agerar på funktioner med kompakt stöd, det funkar alltså inte att bara slänga ut den på en ändpunkt för en ändpunkt har inte kompakt stöd. Det jag själv brukar göra är att skriva om integralen på koordinatoberoende form och sedan gå över till kartesiska koordinater där deltafunktionen har kompat stöd. När jag försöker göra det för det här problemet får jag dock bara noll som svar hela tiden. Jag ska kika lite på det igen sen och återkommer om jag hittar någon lösning.

Freewheeling 220 – Fd. Medlem
Postad: 4 mar 2017 08:05

Tack, emmynoether.

H-lena 24 – Fd. Medlem
Postad: 4 mar 2017 09:53 Redigerad: 4 mar 2017 09:53

Jag fick tipset att integrera över kartesiska koordinater och evaluera i 0 för att då får man ut allt förutom deltafunktionen (antar att den är 1 i punkten 0). Jag kom inte direkt någon vart med den ledtråden ändå dock. 

Albiki 5096 – Fd. Medlem
Postad: 4 mar 2017 18:00

Hej!

Jag har förstått frågan såhär. Du vill bestämma spridningstvärsnittet f(k,k')f(\mathbf{k},\mathbf{k}') där k\mathbf{k} och k' \mathbf{k}' är två vågvektorer i rummet 3.\mathbb{R}^3. Detta är ett tal som bestäms av följande integral.

    f(k,k')=-m2π3V(x)ei(k-k')·xdx\displaystyle f(\mathbf{k},\mathbf{k}') = -\frac{m}{2\pi \hbar}\iiint_{\mathbb{R}^3}V(\mathbf{x})e^{i(\mathbf{k}-\mathbf{k}')\cdot\mathbf{x}}\,d\mathbf{x}.

Albiki

Albiki 5096 – Fd. Medlem
Postad: 4 mar 2017 18:03

Hej!

Potentialen V V ges av en Dirac-deltafunktion i origo.

    V(x)=αδ(x-0) \displaystyle V(\mathbf{x})=\alpha\delta(\mathbf{x}-\mathbf{0}) ,

där α \alpha är ett negativt tal.

Albiki

Albiki 5096 – Fd. Medlem
Postad: 4 mar 2017 18:06

Hej!

Med en sådan potential blir spridningstvärsnittet följande integral.

    f(k,k')=m|α|2π3δ(x-0)ei(k-k')·xdx=m|α|2πei(k-k')·0=m|α|2π . \displaystyle f(\mathbf{k},\mathbf{k}')=\frac{m|\alpha|}{2\pi\hbar}\iiint_{\mathbb{R}^3}\delta(\mathbf{x}-\mathbf{0})e^{i(\mathbf{k}-\mathbf{k}')\cdot \mathbf{x}}\,d\mathbf{x} = \frac{m|\alpha|}{2\pi\hbar}e^{i(\mathbf{k}-\mathbf{k}')\cdot \mathbf{0}} = \frac{m|\alpha|}{2\pi\hbar}\ .

Albiki

Albiki 5096 – Fd. Medlem
Postad: 4 mar 2017 18:08

Hej!

Med en sådan potential blir alltså spridningstvärsnittet konstant för alla vågvektorer.

Vad betyder detta i ett kvantmekaniskt sammanhang?

Albiki

Albiki 5096 – Fd. Medlem
Postad: 4 mar 2017 18:11

Hej!

Förutsatt att jag har förstått frågan rätt, förstås.

Albiki

H-lena 24 – Fd. Medlem
Postad: 4 mar 2017 18:47

Ja det stämmer att det är den integralen men jag vet inte hur jag ska räkna ut den vilket är uppgiften. Jag vet inte hur gränserna ska bli eller hur jag ska göra med deltafunktionen.

Albiki 5096 – Fd. Medlem
Postad: 4 mar 2017 19:08

Hej!

Om du tar en titt på ett av mina inlägg så ser du hur jag beräknande integralen när jag satte in deltafunktionen.

Albiki

H-lena 24 – Fd. Medlem
Postad: 4 mar 2017 19:35

Oj förlåt jag såg inte. Det var sådär jag räknade ut den först men var inte säker på om jag hade gjort rätt!

Guggle 1364
Postad: 4 mar 2017 20:50 Redigerad: 4 mar 2017 21:12

H-lena och Albiki, någonstans på vägen slarvar ni bort en kvadrering av \hbar vilket stör mig.


Första ordningens Bornapprox ges ju av

f(1)(k',k)=-14π2m2d3x'ei(k-k')·x'V(x') f^{(1)}(\mathbf{k}',\mathbf{k})=-\frac{1}{4\pi}\frac{2m}{\hbar^2}\int d^3 x' e^{i(\mathbf {k}-\mathbf{k'})\cdot\mathbf{x}'}V(\mathbf{x}')

Vilket formellt sett bara är en Fouriertransformation av potentialen V m a p q= k-k' \mathbf{k}-\mathbf{k}' (med konstanten -14π2m2 -\frac{1}{4\pi}\frac{2m}{\hbar^2} framför)

emmynoether 663 – Fd. Medlem
Postad: 4 mar 2017 23:00
Albiki skrev :

Hej!

Med en sådan potential blir spridningstvärsnittet följande integral.

    f(k,k')=m|α|2π3δ(x-0)ei(k-k')·xdx=m|α|2πei(k-k')·0=m|α|2π . \displaystyle f(\mathbf{k},\mathbf{k}')=\frac{m|\alpha|}{2\pi\hbar}\iiint_{\mathbb{R}^3}\delta(\mathbf{x}-\mathbf{0})e^{i(\mathbf{k}-\mathbf{k}')\cdot \mathbf{x}}\,d\mathbf{x} = \frac{m|\alpha|}{2\pi\hbar}e^{i(\mathbf{k}-\mathbf{k}')\cdot \mathbf{0}} = \frac{m|\alpha|}{2\pi\hbar}\ .

Albiki

 Jag har rådfrågat runt lite bland vänenr och sen satt jag mig ner och räknade och jag fick precis samma idag när jag tog det från början, dock glömmer du ett minus och \hbar i nämnaren tror jag? I alla fall enlight Griffith's bok. Det H-lena måste missat var att hon inte tog med skalfaktorerna  för deltafunktionen när hon gick över till sfäriska koordinater. 

För vi får då

f(θ)=-2m2q0aδ(r)4πr2rsin(qr)dr=-ma2π20δ(r)sin(qr)qrdr=-ma2π2 f(\theta)= \frac{-2 m}{ \hbar^2 q } \int_0^\infty \frac{a \delta(r)}{4 \pi r^2} r \sin(q r) dr = \frac{-m a}{ 2 \pi \hbar^2} \int_0^\infty \delta(r) \frac{\sin(q r)}{q r} dr = \frac{-ma}{2 \pi \hbar^2}

Här ser vi alltså att det funkar att sätta deltafunktion vid gränsen även fast det inte rent matematiskt är väldefinierat. En matematiker skulle skaka på huvudet och säga att det inte är helt korrekt men eftersom vi är ute efter fysiska resultat och allt vi pysslar med är kontinuerligt deriverbart tillräckligt många gånger kommer vi ändå undan i många fall. Nu har vi i alla fall visat det här på två sätt.

emmynoether 663 – Fd. Medlem
Postad: 4 mar 2017 23:05

I sista steget kanske jag ska påpeka att

0δ(r)sin(qr)qrdr=limr0sin(qr)qr=1 \int_0^\infty \delta(r) \frac{\sin(qr)}{qr} dr = \lim_{r \to 0} \frac{\sin(qr)}{qr} =1

Albiki 5096 – Fd. Medlem
Postad: 5 mar 2017 01:56 Redigerad: 5 mar 2017 01:58

Hej!

Jag reflekterade inte över enheterna, eftersom det var hanteringen av deltafunktionen som var i fokus i denna tråd. Följaktligen blev det en miss med konstanten \hbar istället för 2 \hbar^2 , något som tydligen är stötande, vilket jag ber om ursäkt för.

Däremot hade jag svårt att förstå vad som menas med ett negativt spridningstvärsnitt, vilket fick mig att sätta potentialen till negativ ( α \alpha är ett negativt tal); inom fysik är det vanligt att låta potentialer vara negativa. Vad säger Guggle och emmynoether om negativa spridningstvärsnitt? Är det okej, eller är det stötande?

Albiki

H-lena 24 – Fd. Medlem
Postad: 5 mar 2017 09:48

Guggle: Tack, ja jag märkte också att ett h streck i uträkningen. 

 

Tack allihopa för alla svar och all hjälp!

emmynoether 663 – Fd. Medlem
Postad: 5 mar 2017 15:18
Albiki skrev :

Hej!

Jag reflekterade inte över enheterna, eftersom det var hanteringen av deltafunktionen som var i fokus i denna tråd. Följaktligen blev det en miss med konstanten \hbar istället för 2 \hbar^2 , något som tydligen är stötande, vilket jag ber om ursäkt för.

Däremot hade jag svårt att förstå vad som menas med ett negativt spridningstvärsnitt, vilket fick mig att sätta potentialen till negativ ( α \alpha är ett negativt tal); inom fysik är det vanligt att låta potentialer vara negativa. Vad säger Guggle och emmynoether om negativa spridningstvärsnitt? Är det okej, eller är det stötande?

Albiki

 Generellt i spridningsteorin så skriver man, för en inkommande plan våg i z-led, att 

 

ψ(r,θ)=A(eikz+f(θ)eikrr) \psi(r,\theta) = A ( e^{ikz} + f(\theta) \frac{e^{ikr}}{r})

för stora r r . f(θ) f(\theta) är ju spridningsamplituden, jag tror inte att det innebär problem att den är negativa för sannolikheten ges ju alltid av beloppet i kvadrat, och när man normerar så kan man ju välja A så att allting blir riktigt ändå? Jag är dock lite osäker, jag har inte gått igenom detta i detalj någon gång faktiskt så jag låter det vara lite osagt. Svaret bör vara rätt ändå för det är ju vad vi får, och jag utgick ifrån uttrycket i Griffith's bok när jag räknade.

Svara
Close