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

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

Themen		Antworten	Autor	Aufrufe	Letzter Beitrag
Verwandte Themen - die Beliebtesten
	Herleitung des relativistischen Impuls	10	max_doering	5219	19. Aug 2011 21:52 franz
	Zustandsdichte berechnen	10	TruEnemy	4348	08. Feb 2013 22:36 TruEnemy
	Phononen	7	nekros7	3882	19. Aug 2006 19:16 as_string
	Energie/Anzahl Phononen	6	PhysikFricke	2866	12. März 2012 16:26 PhysikFricke
	Diamant, Phononen	4	neuer	2822	13. Dez 2011 02:05 Braino