Startseite
Forum
Fragen
Suchen
Formeleditor
Über Uns
Registrieren
Login
FAQ
Suchen
Foren-Übersicht
->
Mechanik
Antwort schreiben
Benutzername
(du bist
nicht
eingeloggt!)
Titel
Nachrichtentext
Smilies
Weitere Smilies ansehen
Schriftfarbe:
Standard
Dunkelrot
Rot
Orange
Braun
Gelb
Grün
Oliv
Cyan
Blau
Dunkelblau
Indigo
Violett
Weiß
Schwarz
Schriftgröße:
Schriftgröße
Winzig
Klein
Normal
Groß
Riesig
Tags schließen
Schreibt eure Formeln hier im Board am besten mit Latex!
So gehts:
Latex-Kurzbeschreibung
|
Formeleditor
[quote="DaPhil"]Hallo, ich habe mir gerade mal die Kette mit einer 2 atomigen Basis (Approxiamtion für Phononen) angeschaut und auch die Dispersionrelation w(k) hergeleitet. Ich möchte mir nun die Zustandsdichte anschauen. Hierzu: [latex] n(\omega)=a/\pi \cdot (d\omega/dk)^{-1} [/latex] Ableiten und alle k's durch omegas ersetzen schaffe ich beim besten willen nicht. Also wollte ich mir das ganze mal per Matlab plotten. Hierzu rechne ich den Abstand zwischen zwei omegas aus, dieser ist ja proportional zur Ableitung, da ich meine k's äquidistant gewählt habe (sind sie ja auch in natura). also sollte die gesuchte Zustandsdichte in etwa so aussehen wie [latex]1/\Delta k[/latex]. Tut sie aber nicht. Ich habe nach stundenlanger Recherche auch nur ein bild davon online gefunden... Hier mein Matlab code: [code] % w(k) ausrechnen a=4; m1=5; m2=2; D=0.5; steps=1e3; k=linspace(0,1*pi/a,steps); f=D*(1/m1+1/m2); fplus=f+D*sqrt((1/m1+1/m2)^2-4/(m1*m2)*(sin(k.*a/2)).^2); fplus=sqrt(fplus); fmin=f-D*sqrt((1/m1+1/m2)^2-4/(m1*m2)*(sin(k.*a/2)).^2); fmin=sqrt(fmin); % figure; % plot(k,fplus,'b',k,fmin,'k',[0,0],[0,1.5],'k',[pi/a,pi/a],[0,1.5],'r',[-pi/a,-pi/a],[0,1.5],'r',[-1.5,1.5],[sqrt(2*D/m2),sqrt(2*D/m2)],'r',[-1.5,1.5],[sqrt(2*D/m1),sqrt(2*D/m1)],'r',[-1.5,1.5],[sqrt(2*D*(1/m1+1/m2)),sqrt(2*D*(1/m1+1/m2))],'r'); % n(w) ausrechnen deltaplus=[]; deltamin=[]; for(j=2:steps) deltaplus(j)=abs(fplus(j)-fplus(j-1)); deltamin(j)=abs(fmin(j)-fmin(j-1)); end deltaplus(1)=deltaplus(2); deltamin(1)=deltamin(2); figure; subplot(2,1,1); plot(1./deltaplus); subplot(2,1,2); plot(1./deltamin); [/code] w(k) ausrechnen funktioniert. n(w) leider nicht... Bitte helft mir bei der Fehlersuche.[/quote]
Optionen
HTML ist
aus
BBCode
ist
an
Smilies sind
an
BBCode in diesem Beitrag deaktivieren
Smilies in diesem Beitrag deaktivieren
Spamschutz
Text aus Bild eingeben
Alle Zeiten sind GMT + 1 Stunde
Gehe zu:
Forum auswählen
Themenbereiche
----------------
Mechanik
Elektrik
Quantenphysik
Astronomie
Wärmelehre
Optik
Sonstiges
FAQ
Sonstiges
----------------
Off-Topic
Ankündigungen
Thema-Überblick
Autor
Nachricht
DaPhil
Verfasst am: 19. Mai 2010 16:25
Titel: Zustandsdichte Phononen mit 2 atomiger Basis
Hallo,
ich habe mir gerade mal die Kette mit einer 2 atomigen Basis (Approxiamtion für Phononen) angeschaut und auch die Dispersionrelation w(k) hergeleitet. Ich möchte mir nun die Zustandsdichte anschauen. Hierzu:
Ableiten und alle k's durch omegas ersetzen schaffe ich beim besten willen nicht. Also wollte ich mir das ganze mal per Matlab plotten. Hierzu rechne ich den Abstand zwischen zwei omegas aus, dieser ist ja proportional zur Ableitung, da ich meine k's äquidistant gewählt habe (sind sie ja auch in natura). also sollte die gesuchte Zustandsdichte in etwa so aussehen wie
. Tut sie aber nicht. Ich habe nach stundenlanger Recherche auch nur ein bild davon online gefunden... Hier mein Matlab code:
Code:
% w(k) ausrechnen
a=4;
m1=5;
m2=2;
D=0.5;
steps=1e3;
k=linspace(0,1*pi/a,steps);
f=D*(1/m1+1/m2);
fplus=f+D*sqrt((1/m1+1/m2)^2-4/(m1*m2)*(sin(k.*a/2)).^2);
fplus=sqrt(fplus);
fmin=f-D*sqrt((1/m1+1/m2)^2-4/(m1*m2)*(sin(k.*a/2)).^2);
fmin=sqrt(fmin);
% figure;
% plot(k,fplus,'b',k,fmin,'k',[0,0],[0,1.5],'k',[pi/a,pi/a],[0,1.5],'r',[-pi/a,-pi/a],[0,1.5],'r',[-1.5,1.5],[sqrt(2*D/m2),sqrt(2*D/m2)],'r',[-1.5,1.5],[sqrt(2*D/m1),sqrt(2*D/m1)],'r',[-1.5,1.5],[sqrt(2*D*(1/m1+1/m2)),sqrt(2*D*(1/m1+1/m2))],'r');
% n(w) ausrechnen
deltaplus=[];
deltamin=[];
for(j=2:steps)
deltaplus(j)=abs(fplus(j)-fplus(j-1));
deltamin(j)=abs(fmin(j)-fmin(j-1));
end
deltaplus(1)=deltaplus(2);
deltamin(1)=deltamin(2);
figure;
subplot(2,1,1);
plot(1./deltaplus);
subplot(2,1,2);
plot(1./deltamin);
w(k) ausrechnen funktioniert. n(w) leider nicht... Bitte helft mir bei der Fehlersuche.