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 15: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	11446	14. Jul 2023 17:24 Mispel
	Zustandsdichte in verschiedenen Dimensionen	0	RogerKlotz	2652	10. Jul 2022 13:04 RogerKlotz
	Pauli-Matrizen als Basis der hermiteschen 2x2 Matrizen.	3	manuel459	2885	15. März 2022 21:51 TomS
	Kombinierte Zustandsdichte Quantenpunkt	0	Camina Druma	1257	17. Dez 2021 00:26 Camina Druma
	Kombinierte Zustandsdichte parabolischer Bänder	0	Marlene L.	997	28. Nov 2021 14:05 Marlene L.

Themen		Antworten	Autor	Aufrufe	Letzter Beitrag
Verwandte Themen - die Größten
	Phononen	13	Phonon	2685	07. Nov 2012 01:21 TomS
	Zustandsdichte berechnen	10	TruEnemy	4677	08. Feb 2013 22:36 TruEnemy
	Herleitung des relativistischen Impuls	10	max_doering	5972	19. Aug 2011 20:52 franz
	Messung eines Quantensystems zur Basis	9	Tetrix	1971	26. Jan 2021 21:24 Tetrix
	Zeeman-Effekt: experimentelle Basis?	8	rüdisüli	1506	20. Jul 2016 11:09 TomS

Themen		Antworten	Autor	Aufrufe	Letzter Beitrag
Verwandte Themen - die Beliebtesten
	Zustandsdichte endliche/unendliche lineare Kette	0	Mispel	11446	14. Jul 2023 17:24 Mispel
	Herleitung des relativistischen Impuls	10	max_doering	5972	19. Aug 2011 20:52 franz
	Zustandsdichte berechnen	10	TruEnemy	4677	08. Feb 2013 22:36 TruEnemy
	Phononen	7	nekros7	4284	19. Aug 2006 18:16 as_string
	Energie/Anzahl Phononen	6	PhysikFricke	3464	12. März 2012 16:26 PhysikFricke