I studier av bulk vatten är det viktigt att förstå hur olika typer av vibrational kopplingar påverkar de experimentella spektra som observeras i IR och Raman-spektroskopi. För att förklara detta måste vi först se på hur intramolekylära och intermolekylära kopplingar interagerar med de vibrationsspektra som genereras av vattenmolekyler.

När man tittar på OH-sträckfrekevenser i vätskefas vatten, märker man att starkare vätebindningar leder till lägre frekvenser för OH-stretching. Detta innebär att den starkare vätebindningen får den relativa OH-sträckningen att bli mindre energirik, vilket avspeglas i en sänkning av sträckfrekvensen. I detta sammanhang är det viktigt att notera att icke-Condon-effekter inte verkar påverka övergångspolariserbarheten i det aktuella systemet på ett signifikant sätt, vilket innebär att intramolekylär vibrationskoppling är den främsta faktorn som påverkar spektrala förändringar.

Spektren som ses i figurerna 4.4a och 4.4b i den numeriska simuleringen visar tydligt denna effekt. När endast intramolekylära kopplingar beaktas (fig. 4.4b), blir spektrumet blåförskjutet jämfört med det utan kopplingar (fig. 4.4a), vilket beror på att den antisymmetriska sträckningen har högre frekvens och större övergångsdipol än den symmetriska. Den antisymmetriska sträckningen leder alltså till en högre frekvens eftersom kopplingskonstanten är negativ. Omvänt är den symmetriska sträckningen mer framträdande i Raman-spektroskopi, vilket gör att vi ser en liknande topp i spektrumet i fig. 4.4b och 4.4f. I Raman-spektroskopi är det alltså viktigt att förstå skillnaden mellan hur intramolekylär symmetrisk och antisymmetrisk sträckning påverkar intensiteten i spektrala toppar.

Det intermolekylära vibrationskopplingen har också en påtaglig effekt, vilket kan ses i fig. 4.4g. Här förskjuts det lågfrekventa området i Raman-VV-spektret betydligt. Detta beror på att in-fas sträckningar, som är delokaliserade över ett stort antal OH-grupper (ungefär 10), får en lägre frekvens och därmed en större tvärsnittsarea i Raman-VV-spektrumet. Detta fenomen beror på negativa kopplingskonstanter som leder till ett större bidrag från dessa sträckningar i Raman-VV-spektrumet.

Det är också viktigt att förstå att begreppen "in-phase" och "anti-phase" i sammanhanget av vätebindningar inte bara handlar om fasen av sträckningarna, utan också om hur de orienterade OH-gruppernas övergångspolariserbarheter bidrar till det totala spektrumet. I Raman-VH och IR-spektrumen förstärks "anti-phase" sträckningar, eftersom deras övergångspolariserbarheter inte är lika känsliga för orienteringen av OH-grupperna, vilket inte är fallet för in-phase sträckningar.

Vidare kan Fermi-resonans kopplingar påverka spektra genom att skapa specifika dippar och toppar. I Raman-VV-spektrumet kan man till exempel observera en dip vid 3300 cm−1, som kallas Evans-fönstret. Denna dip uppstår till följd av en resonans mellan OH-sträckningens fundamentala frekvens och HOHbend-övertonen, vilket leder till en splittring av det intermolekylärt breddade OH-sträckbands spektrum och därmed en bimodal funktion i spektrumet. Fermi-resonansen skapar även en svag axel i IR-spektret runt 3200 cm−1.

Det är också väsentligt att förstå att Fermi-resonans inte märkbart påverkar Raman-VH-spektret, vilket gör att detta spektrum bibehåller sin relativt enkla form jämfört med Raman-VV.

Vid en noggrannare undersökning av bulk vatten och dess spektra måste vi erkänna att det är svårt att enkelt koppla OH-sträckningarna till specifika kemiska arter, såsom starkt och svagt vätebindade molekyler. Detta innebär att den klassiska uppdelningen av vatten i starkt och svagt vätebundna grupper inte alltid ger en tillräcklig beskrivning av komplexiteten i intermolekylära interaktioner i vätskor.

För att förstå relationen mellan strukturen hos vätskan och de spektrala signaturerna som observeras i experimenten, används olika metoder, bland annat analys av MD-trajectories. I detta avseende introducerades M-indexet, ett kraftfullt verktyg för att analysera vätebindande molekylsystem på basis av molekylens dimerstatus och vätebindningens karaktär. Genom att definiera MOH-indexet (som ursprungligen föreslogs av Ohno och hans medarbetare) kan man använda en enkel formel för att karakterisera styrkan av vätebindningen för varje OH-grupp i vattensystemet.

M-indexet tillåter oss att ranka vätebindningarnas styrka i fyra nivåer, och det har visat sig vara ett effektivt sätt att koppla spektrala data till de specifika strukturella egenskaperna hos molekyler i vätska. Genom att räkna antalet vätebindningar som en vattenmolekyl är engagerad i, får vi ett enkelt men kraftfullt verktyg för att förstå de komplexa interaktionerna i vätskefas vatten.

Hur kan elektromagnetiska vågor beskrivas och lösas i frekvensdomänen med hjälp av Maxwells ekvationer?

Elektromagnetiska fält, såsom det elektriska fältet, kan analyseras genom att använda Fourier-transformen och byta ut tidsvariabeln mot vinkelfrekvensen ω. Denna transform ger en funktion i frekvensdomänen, vilket förenklar beskrivningen av fältens variation i tid och rum. Genom att Fourier-transformera Maxwells ekvationer erhålls en formel för det komplexa elektriska fältet i frekvensdomänen, vilket möjliggör en mer hanterbar beskrivning av elektromagnetiska vågor i olika material, inklusive dielektrika, metaller och halvledare.

I frekvensdomänen framträder en generell vågekvation som beskriver elektromagnetiska vågor med hänsyn till polarisation och strödtäthet i materialet. Dessa termer kan kombineras i en generaliserad polarisationsvektor och en generaliserad elektrisk flödestäthet, vilket förenklar Maxwells ekvationer genom att eliminera behovet av att explicit hantera verklig laddnings- och strödtäthet. Denna generalisering möjliggör en enhetlig beskrivning av elektromagnetiska vågor i material som annars uppvisar olika elektriska egenskaper, till exempel metaller där den komplexa brytningsindexen kan tolkas inom samma ramverk som för dielektriska material.

Det är viktigt att notera att distinktionen mellan polarisation och elektrisk ström är i många fall godtycklig, och i spektroskopiska experiment är det ofta omöjligt att avgöra vilken komponent som bidrar till det observerade svaret utan förhandsinformation om materialet. Detta understryker att det är fullt möjligt att formulera elektromagnetismens grundläggande ekvationer utan att särskilja mellan dessa bidrag, något som också speglar den mikroskopiska formuleringen av Maxwells ekvationer där polarisationsbegreppet inte alltid är explicit.

Vid analys av gränssnitt mellan olika material, till exempel i tvåskiktsmodellen med komplexa brytningsindex, härleds Fresnelkoefficienterna för reflektion och transmission från vågekvationen och tillhörande randvillkor. Antagandet om plana vågor med en given vinkelfrekvens leder till uttryck för de elektriska och magnetiska fältens komponenter i både infallande, reflekterande och transmitterande vågor. Dessa komponenter kan delas upp i P- och S-polarisation, vilket är centralt för att korrekt beskriva vågornas beteende vid gränsskiktet.

Vid gränssnittet mellan två media tillämpas randvillkor som säkerställer kontinuiteten av de relevanta fältkomponenterna. Genom att analysera ett infinitesimalt tunt skikt kring gränsen kan sambandet mellan amplituderna för infallande, reflekterande och transmitterande vågor fastställas, vilket leder till de välkända Fresnelkoefficienterna. Dessa koefficienter beskriver hur mycket av den infallande vågen som reflekteras respektive transmitteras och är frekvens- och vinkeldependenta, beroende på mediernas optiska egenskaper.

Det är centralt att förstå att de komplexa brytningsindexen, särskilt i metaller, uppstår ur den generaliserade polarisations- och strömbeskrivningen, vilket visar att elektromagnetiska vågor i dessa material inte skiljer sig principiellt från vågor i dielektrika, utan att skillnaden ligger i materialresponsen. Denna insikt möjliggör en sammanhängande teori för elektromagnetiska vågor över olika materialklasser.

Vidare är det av vikt att förstå att frekvensdomänsbeskrivningen inte bara förenklar lösningen av Maxwells ekvationer utan också möjliggör direkt koppling till mätbara storheter i experimentella sammanhang, såsom spektroskopi. Detta kräver dock en noggrann tolkning av materialparametrar och deras frekvensberoende, liksom en förståelse för hur materialets mikroskopiska egenskaper påverkar den makroskopiska responsen.

Slutligen är det viktigt att inse att formuleringen i frekvensdomänen och användandet av generaliserade fält och polarisationsbegrepp utgör ett kraftfullt verktyg för att analysera och designa optiska och elektromagnetiska system, där komplexa material och gränssnitt spelar en avgörande roll. Denna teoretiska grund är oumbärlig för att förstå och utnyttja ljusets interaktion med materia på djupet.

Hur härleds gränsvillkoren för elektromagnetiska fält vid gränsytan mellan två medier?

När vi betraktar ett område V i rummet, definierat av koordinaterna X,Y,ZX, Y, Z med villkoren X0<X<X0+lXX_0 < X < X_0 + l_X, Y0<Y<Y0+lYY_0 < Y < Y_0 + l_Y och lZ2<Z<lZ2-\frac{l_Z}{2} < Z < \frac{l_Z}{2}, kan vi analysera de elektromagnetiska fälten och deras gränsvillkor vid gränsytan mellan två medier. Detta område omsluts av en yta SS som fungerar som gränsen mellan medierna. Genom att använda divergensteoremet (Gauss sats) kan volymintegralen av divergensen av en vektor, såsom elektrisk flödestäthet D\mathbf{D}, omvandlas till en ytintegral över SS.

Utifrån D=0\nabla \cdot \mathbf{D} = 0 gäller att

VDdV=0,\int_V \nabla \cdot \mathbf{D} \, dV = 0,

vilket via Gauss teorem motsvaras av

SDdS=0.\int_S \mathbf{D} \cdot d\mathbf{S} = 0.

Ytintegralen kan delas upp i tre delar över ytorna vinkelräta mot XX, YY, och ZZ-axlarna. För fall där vågorna propageterar i XZXZ-planet och fälten är oberoende av YY, försvinner bidraget från ytorna i YY-riktningen. När tjockleken i ZZ-riktningen, lZl_Z, närmar sig noll, kvarstår endast ytintegralen över gränsytan Z=0Z=0, vilket leder till villkoret:

DZZ=+0=DZZ=0.D_Z|_{Z=+0} = D_Z|_{Z=-0}.

Detta uttrycker kontinuiteten hos den normala komponenten av den elektriska flödestätheten över gränsytan. Med andra ord, den elektriska flödestätheten är kontinuerlig i riktningen vinkelrät mot gränsytan.

I det övre mediet finns både infallande och reflekterade vågor, medan i det nedre endast transmissionsvågor existerar. Genom att analysera dessa vågor och sätta lika villkoren för den normala komponenten av D\mathbf{D} i båda medierna, erhålls reflektions- och transmissionslagen. Dessa leder till lagarna för reflektion och brytning, där infallsvinkeln är lika med reflektionsvinkeln och brytningsvinkeln styrs av Snells lag:

θi=θr,nasinθi=nbsinθt.\theta_i = \theta_r, \quad n_a \sin \theta_i = n_b \sin \theta_t.

Vidare kan relationer mellan amplituderna för infallande, reflekterade och transmitterade vågor uttryckas, särskilt för de olika polarisationskomponenterna, PP- och SS-polarisation.

Liknande resonemang kan göras för magnetisk flödestäthet B\mathbf{B}, vilket leder till att dess normala komponent också är kontinuerlig över gränsytan:

BZZ=+0=BZZ=0.B_Z|_{Z=+0} = B_Z|_{Z=-0}.

När man betraktar virvlar (curl) av elektriska och magnetiska fält i infinitesimalt tunna områden kring gränsytan och applicerar Stokes sats, härleds ytterligare gränsvillkor för de tangentiella komponenterna av elektriska och magnetiska fält. Till exempel, kontinuiteten av den tangentiella komponenten av det elektriska fältet i XX-riktningen:

EXZ=+0=EXZ=0,E_X|_{Z=+0} = E_X|_{Z=-0},

och motsvarande för magnetfältet. Dessa villkor kompletterar de normala komponenternas kontinuitet och utgör tillsammans de fyra grundläggande gränsvillkoren för elektromagnetiska fält vid en gränsyta.

Utifrån dessa villkor kan Fresnelkoefficienterna för reflektion och transmission formuleras. Dessa koefficienter, som beror på materialens brytningsindex och infallsvinklar, beskriver amplitudförhållandena mellan infallande, reflekterade och transmitterade vågor för både PP- och SS-polarisation. Särskilt är Brewsters vinkel den infallsvinkel där reflektionskoefficienten för PP-polariserad ljus går till noll, vilket innebär fullständig transmission utan reflektion.

Det är viktigt att inse att dessa gränsvillkor bygger på antagandet om idealiserade, skarpa gränsytor och linjära, isotropa material. Vid verkliga gränsytor kan ytråhet, materialanisotropi, och icke-linjära effekter påverka fältens beteende och leda till avvikelser från den ideala beskrivningen. Dessutom förutsätts att tidsberoendet i fälten är harmoniskt och stationärt, vilket gör att tidsderivator kan hanteras på ett förenklat sätt. Vid dynamiska eller komplexa materialegenskaper krävs mer avancerade modeller.

För att fullt ut förstå vågens beteende vid gränsytan är det också avgörande att betrakta energiöverföring och effektflöden, vilka kan analyseras med hjälp av Poyntingvektorn. Denna analys visar att medan elektriska och magnetiska fältens komponenter är kontinuerliga enligt de härledda villkoren, kan energi reflekteras och transmitteras i olika proportioner beroende på materialens egenskaper och vågens polarisation.

Hur man beräknar den effektiva andra ordningens icke-linjära optiska susceptibiliteten

För att bestämma den effektiva andra ordningens icke-linjära optiska susceptibiliteten i ett heterodynt detekterat SFG-mätning (Sum Frequency Generation) med olika polarisationskombinationer, är det nödvändigt att använda specifika tensorer och deras relationer för att beskriva och analysera materialets respons vid olika frekvenser. Genom att använda de relevanta tensorer som χXXZ, χYYZ, och χZZZ kan vi härleda uttrycken för de olika effektiva susceptibiliteterna.

För det första, genom att använda de icke-noll tensor elementen i uttryck för polarisering PP enligt ekvationerna (8.104) och (8.105), kan man erhålla de effektiva andra ordningens icke-linjära optiska susceptibiliteter. Dessa uttryck, beroende på polarisationen, kan skrivas som:

  • χeff,SSP(2)=L3YL1YL2Zsin(2θ2i)χYYZ\chi_{\text{eff,SSP}}^{(2)} = L_3Y L_1Y L_2Z \sin(2 \theta_{2i}) \chi_{YYZ}

  • χeff,SPS(2)=L3YL1ZL2Ysin(θ1i)χYZY\chi_{\text{eff,SPS}}^{(2)} = L_3Y L_1Z L_2Y \sin(\theta_{1i}) \chi_{YZY}

  • χeff,PSS(2)=L3ZL1YL2Ysin(θ3r)χZYY\chi_{\text{eff,PSS}}^{(2)} = L_3Z L_1Y L_2Y \sin(\theta_{3r}) \chi_{ZYY}

  • χeff,PPP(2)=L3XL1XL2Zcos(θ3r)cos(θ1i)sin(2θ2i)χXXZ\chi_{\text{eff,PPP}}^{(2)} = -L_3X L_1X L_2Z \cos(\theta_{3r}) \cos(\theta_{1i}) \sin(2 \theta_{2i}) \chi_{XXZ}

Vid vidare analys kan man få ytterligare termer som involverar både refraktionsindex och Fresnel-koefficienter för att bestämma den exakta responsen i olika material, till exempel kvarts, där specifika parametrar som χq\chi_q och npn_p spelar en avgörande roll. Dessa relationer är nödvändiga för att exakt kunna beskriva och förstå hur SFG-fältet E3rSE_{3rS} genereras vid provets yta i förhållande till referensytan av kvarts.

Vid mätningar av den andra ordningens optiska respons i heterodyna SFG är det också viktigt att förstå hur brytningsvinkeln, eller Brewster-vinkeln, påverkar reflektioner och transmission vid olika gränssnitt. Genom att använda de generella formlerna för reflektion vid olika vinklar kan man härleda relationer för materialets optiska egenskaper, exempelvis genom att sätta upp förhållandet mellan vågvektorer för olika lager i en etalon.

Vid beräkning av denna respons måste vi ta hänsyn till de komplexa brytningsindexen för materialet, inklusive både real- och imaginär del, som beskriver hur ljuset dämpas och förändras när det passerar genom flera lager, såsom i fallet med etalonstrukturer bestående av guld och vatten. Den resulterande transmittanskurvan visar skarpa band som är kopplade till resonansfenomen i systemet.

För att kunna utföra dessa beräkningar är det också viktigt att förstå hur man löser uppsättningar av linjära ekvationer som beskriver gränsvillkoren mellan de olika materialen. Detta görs genom att noggrant beakta de elektriska och magnetiska fälten i varje lager och sedan tillämpa kontinuitetskraven vid varje materialgräns.

En viktig aspekt som måste beaktas är skillnaden mellan ytbeteende och bulkrespons i material. I vissa fall, som med kvarts, kommer det genererade SFG-fältet att vara beroende av både ytegenskaper och bulkens optiska egenskaper, vilket påverkar den uppmätta signalens intensitet och spektrala fördelning.

För att exakt kunna modellera och tolka resultaten från sådana mätningar, är det också nödvändigt att förstå hur den effektiva andra ordningens icke-linjära susceptibiliteten varierar beroende på materialets kristallstruktur och geometri, och hur dessa faktorer påverkar den optiska responsen vid olika frekvenser.

Hur härleds den molekylära polariserbarhetstensoren från Liouville-ekvationen?

Interaktionen mellan ett ljusfält och en molekyls elektriska dipolmoment kan beskrivas med hjälp av Liouville-ekvationen, vilket möjliggör en analytisk framställning av molekylens polariserbarhet. För enkelhetens skull betraktas här ett tvånivåsystem med Hamiltonoperatorn definierad som summan av det obearbetade Hamiltonianet och en perturbation som representerar interaktionen mellan ljusets elektriska fält och molekylens elektriska dipolmoment.

Hamiltonoperatorn består av termerna för grund- och exciterade tillstånd och en term som innehåller den elektriska dipolmomentoperatorns komponenter multiplicerade med det lokala elektriska fältet. Det antas att fältets tidsberoende kan beskrivas som en exponentiell funktion i form av en komplex fas, vilket är förenligt med Fouriertransformens principer. Under den elektriska dipolapproximationen behandlas molekylen som en punktdipol eftersom molekylens dimensioner är mycket mindre än ljusets våglängd.

Liouville-ekvationen löses perturbativt genom att expandera både Hamiltonoperatorn och densitetsoperatorn i ordningar av perturbationen. För nollte ordningen, då ingen interaktion med ljuset sker, antas att alla molekyler befinner sig i grundtillståndet. Detta innebär att densitetsmatrisens diagonalpost är 1 för grundtillståndet och 0 för det exciterade, samt att inga koherenser förekommer.

Vid första ordningen framkommer en differentialekvation för densitetsmatrisens icke-diagonala element som beskriver koherens mellan grund- och exciterade tillstånd. Utan ljusfält kvarstår koherensen utan avklingning, vilket är orealistiskt då de verkliga systemen påverkas av avfasi. Detta korrigeras genom införandet av en fenemenologisk dämpningsterm, vilken representerar avfasningshastigheten. Denna term gör att koherensen avklingar exponentiellt med tiden, vilket överensstämmer med experimentellt observerad relaxationsdynamik.

När ljusfältet inte är noll löses differentialekvationen genom en konvolution av en responsfunktion med fältets tidsberoende, vilket i frekvensdomänen ger ett uttryck för densitetsmatrisens koherenta komponenter som är proportionella mot det elektriska fältets Fouriertransformerade spektrum.

Genom att beräkna förväntningsvärdet av dipolmomentet med hjälp av densitetsmatrisens första ordningens komponenter erhålls den inducerade dipolmomentets uttryck i termer av elektriska fältets frekvensinnehåll. Den molekylära polariserbarhetstensoren definieras som sambandet mellan det inducerade dipolmomentet och det elektriska fältet. Resultatet är en tensorfunktion av frekvens, vilken i allmänhet måste summeras över alla exciterade tillstånd i systemet.

Denna tensor innehåller information om hur molekylen svarar på det elektriska fältet i olika riktningar och vid olika frekvenser, och är central för förståelsen av molekylers optiska egenskaper. Det är avgörande att inse att uttrycken bygger på approximationer, såsom tvånivåmodellen och dipolapproximationen, men att de kan generaliseras till mer komplexa system med flera exciterade tillstånd.

För att fullt ut förstå de fysikaliska implikationerna av polariserbarhetstensoren krävs insikt i att den beskriver både den resonanta och icke-resonanta responsen hos molekylen, där resonanser uppstår nära övergångsfrekvenser mellan energinivåer. Avfasningsparametern är avgörande för bredden på dessa resonanser och därmed för linjebredden i spektroskopiska mätningar.

Dessutom är det viktigt att förstå att Liouville-ekvationen inte bara hanterar populationsfördelningar utan även kvantkoherenser, vilket möjliggör en kvantmekanisk beskrivning av ljus-materiainteraktioner som inkluderar både absorptions- och emissionsprocesser. Fenomenologiska dämpningsparametrar representerar omgivningens inverkan på systemet, vilket ofta kräver experimentell bestämning eller mer avancerade modeller för att beskriva relaxation och avfasi i realistiska miljöer.