Zustandsdichte Phononen mit 2 atomiger Basis

Registrieren

FAQ

Suchen

Neue Frage »

Antworten »

Foren-Übersicht -> Mechanik

Autor

Nachricht

DaPhil

Anmeldungsdatum: 24.07.2008
Beiträge: 23

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.

Neue Frage »

Antworten »

Foren-Übersicht -> Mechanik

Themen		Antworten	Autor	Aufrufe	Letzter Beitrag
Verwandte Themen - die Neuesten
	Zustandsdichte endliche/unendliche lineare Kette	0	Mispel	17058	14. Jul 2023 18:24 Mispel
	Zustandsdichte in verschiedenen Dimensionen	0	RogerKlotz	3238	10. Jul 2022 14:04 RogerKlotz
	Pauli-Matrizen als Basis der hermiteschen 2x2 Matrizen.	3	manuel459	3584	15. März 2022 21:51 TomS
	Kombinierte Zustandsdichte Quantenpunkt	0	Camina Druma	1377	17. Dez 2021 00:26 Camina Druma
	Kombinierte Zustandsdichte parabolischer Bänder	0	Marlene L.	1113	28. Nov 2021 14:05 Marlene L.

Themen		Antworten	Autor	Aufrufe	Letzter Beitrag
Verwandte Themen - die Größten
	Phononen	13	Phonon	2783	07. Nov 2012 01:21 TomS
	Zustandsdichte berechnen	10	TruEnemy	4787	08. Feb 2013 22:36 TruEnemy
	Herleitung des relativistischen Impuls	10	max_doering	6092	19. Aug 2011 21:52 franz
	Messung eines Quantensystems zur Basis	9	Tetrix	2119	26. Jan 2021 21:24 Tetrix
	Phononen	8	YO	1883	26. Jun 2021 00:02 Myon

Themen		Antworten	Autor	Aufrufe	Letzter Beitrag
Verwandte Themen - die Beliebtesten
	Zustandsdichte endliche/unendliche lineare Kette	0	Mispel	17058	14. Jul 2023 18:24 Mispel
	Herleitung des relativistischen Impuls	10	max_doering	6092	19. Aug 2011 21:52 franz
	Zustandsdichte berechnen	10	TruEnemy	4787	08. Feb 2013 22:36 TruEnemy
	Phononen	7	nekros7	4345	19. Aug 2006 19:16 as_string
	Pauli-Matrizen als Basis der hermiteschen 2x2 Matrizen.	3	manuel459	3584	15. März 2022 21:51 TomS