Elasticitet av salter

Coulomb-kraft och repulsion

Jonära kristaller hålls ihop av Coulomkraften, men det behövs en repulsiv kraft för att förhindra att avstånden blir mindre och mindre och att kristallen kollapsar till en punkt.
Man kan ta potentialen av en hård sfär, men då finns det inget att räkna vidare med. Realistiska potentialer är lite mjukare. Man kan ta en exponentiell funktion U = B [Graphics:Images/ElastSalt_gr_1.gif], där ρ ≈ 0,3 Å. Eller man kan ta U = B [Graphics:Images/ElastSalt_gr_2.gif], där n ≈ 10. Om man vet [Graphics:Images/ElastSalt_gr_3.gif] och bindningsenergi, kan man bestäma B, n och ρ, som jag visar i en notebook om Born-Mayer teori.

Istället ska jag följa Turton och välja potenslagen med n=10. Vi skriver

[Graphics:Images/ElastSalt_gr_4.gif]

och

[Graphics:Images/ElastSalt_gr_5.gif]
[Graphics:Images/ElastSalt_gr_6.gif]

Eftersom den totala kraften vid jämnvikt är lika med 0, har vi FCoul[[Graphics:Images/ElastSalt_gr_7.gif]] + FRepul[[Graphics:Images/ElastSalt_gr_8.gif]] = 0.
Vi ser att  B = A.

I NaCl har vi laddningar +e och -e, och ett jämnviktsavstånd

[Graphics:Images/ElastSalt_gr_9.gif]

så att

[Graphics:Images/ElastSalt_gr_10.gif]
[Graphics:Images/ElastSalt_gr_11.gif]

Den elektrostatiska kraften är alltså 2,90 nN och den repulsiva kraften är lika stor åt andra hållet.

[Graphics:Images/ElastSalt_gr_12.gif]

[Graphics:Images/ElastSalt_gr_13.gif]

Det var en plot av kraften i newton.
För en plot av den potentiella energi i eV integrerar vi och delar genom e:

[Graphics:Images/ElastSalt_gr_14.gif]

[Graphics:Images/ElastSalt_gr_15.gif]

Vi ser att potentialens minimum ligger vid jämnsviktsavståndet [Graphics:Images/ElastSalt_gr_16.gif].

Kraftkonstanten

Nära kring jämnviktsavståndet kan potentialen approximeras som harmonisk:

[Graphics:Images/ElastSalt_gr_17.gif]

Den återdrivande kraften nära jämnviktsavståndet kan appriximeras linjärt, så att vi får Hookes lag:

[Graphics:Images/ElastSalt_gr_18.gif]

där k är kraftkonstanten som är alltså derivatan av kraften vid jämnviktsavståndet.

Om F[r] = -A [Graphics:Images/ElastSalt_gr_19.gif]/[Graphics:Images/ElastSalt_gr_20.gif] + A[Graphics:Images/ElastSalt_gr_21.gif]/[Graphics:Images/ElastSalt_gr_22.gif], ser vi att  dF/dr = 2A [Graphics:Images/ElastSalt_gr_23.gif]/[Graphics:Images/ElastSalt_gr_24.gif] - 10 A [Graphics:Images/ElastSalt_gr_25.gif]/[Graphics:Images/ElastSalt_gr_26.gif]. Vid [Graphics:Images/ElastSalt_gr_27.gif]har vi alltså k = dF/dr = 8 A/[Graphics:Images/ElastSalt_gr_28.gif]. För NaCl blir det en kraftkonstant

[Graphics:Images/ElastSalt_gr_29.gif]
[Graphics:Images/ElastSalt_gr_30.gif]

Kraftkonstanten mellan en Na- och en Cl-jon är alltså 82 N/m.
I något lämpligare enheter blir det 82 pN/pm.

Elasticitetsmodulen

Makroskopisk definierar man Young's modulus som förhållandet mellan spänning (tryck) och töjning (kompression) av en stav. Spänning är kraft per yta (σ=[Graphics:Images/ElastSalt_gr_31.gif]) och töjning är den relativa längsändringen (ε=[Graphics:Images/ElastSalt_gr_32.gif]). Så vi har

[Graphics:Images/ElastSalt_gr_33.gif]

Om vi nu betraktar vår lilla NaCl-dimer som en stav, ser vi att längden [Graphics:Images/ElastSalt_gr_34.gif] och att tvärsnittsarean [Graphics:Images/ElastSalt_gr_35.gif]. Enligt Hookes lag är [Graphics:Images/ElastSalt_gr_36.gif], så att

[Graphics:Images/ElastSalt_gr_37.gif]
[Graphics:Images/ElastSalt_gr_38.gif]

En enkel teoretisk uppskattning av elasticitetsmodulen ger alltså att elasticitetsmodulen är ungefär 300 GPa. Det experimentella värdet av bulk modulus är 24 GPa, för Young's modulus 40 GPa. Uppskattningen är alltså för hög med en tiopotens.

Experimentell kan man bestämma elasticitetsmodulen ur ljudhastigheten.I en stav gäller:

[Graphics:Images/ElastSalt_gr_39.gif]