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

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

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