I tilfelle av HOD-bøyning er den transisjonelle dipolmomentet ikke parallelt med bisektoren på grunn av asymmetrien. Dette innebærer at formelen i Eq. (3.18) ikke er anvendelig for HOD. I stedet må uttrykket for 〈0|μ̂p|n, τ〉 defineres på en konsistent måte med den transisjonelle polarisabilitetstensoren. Dette betyr at dipolmomentet for HOD-bøyning, så vel som HOH og DOD, vil bli forklart i Seksjon 3.1.7, der den transisjonelle polarisabilitetstensoren diskuteres.

Når man beregner IR-spekteret numerisk ved hjelp av WFP-metoden (Wave Function Propagation), følger prosedyren disse trinnene: først må man oppnå MD-banene. Deretter konstrueres Hamilton-matrisen H̃ ved hvert tidssteg i henhold til kapitlet som omhandler det. Deretter beregnes det transisjonelle dipolmomentet 〈0|μ̂p|n, τ〉 for hvert basisket i hvert tidssteg i samsvar med de relevante ligningene. H̃ diagonaliseres for å få frekvensegenverdien ωξ,(τ−1) og egenvektorkoeffisienten cn,ξ,(τ−1) på hvert tidssteg. Deretter beregnes IR-bølgefunksjonskoeffisienten an,p,τ i hvert tidssteg. For å til slutt oppnå IR-spekteret, utføres en Fourier-transformasjon av tidskorrelasjonsfunksjonen fIR(τ). For at IR-spekteret skal ha tilstrekkelig kvalitet, kreves vanligvis omfattende akkumulering av fIR(τ) langs MD-banen. Dette betyr at jo større antallet MD-tidssteg, varigheten på tidskorrelasjonsfunksjonen, og antallet TIP4P-vannmolekyler er, desto mer tid må brukes på akkumulering for å få et "jevn" IR-spektrum i OH-strekkregionen.

Raman-spekteret, derimot, krever beregning av tidskorrelasjonsfunksjonen for polarisabilitetstensoren. Den totale polarisabilitetstensoren α̂(t) for systemet uttrykkes som summen av polarisabilitetstensorene for hver enkelt molekyl α̂i(t). I tilfelle Raman-spektroskopi, der polarisasjonsavhengige målinger utføres selv for isotrope prøver, introduseres en ny tidskorrelasjonsfunksjon for polarisabilitetstensoren. Denne kan skrives som f_pq(t) = 〈0|α̂pqÛ(t)α̂pq|0〉, hvor p og q representerer de XYZ-laboratorieaksene, og α̂pq er komponenten av polarisabilitetstensoren ved t=0. Tidsutviklingen til denne funksjonen kan behandles på samme måte som for dipolmomentet, men med hensyn til den polarisabilitetsrelaterte bølgefunksjonen.

For Raman-spektroskopi, når man arbeider med et isotropisk system som bulk væske, er de mest populære polarisasjonskombinasjonene VV og VH, der eksitasjons- og spredningslysets polariseringsretninger henholdsvis er parallelle og perpendikulære. Raman-spekteret i disse polarisasjonskombinasjonene kan uttrykkes som integraler som involverer tidskorrelasjonsfunksjonene fXX(t), fYY(t), fZZ(t), og de tverrpolarisasjonene fXY(t), fXZ(t), fYZ(t), avhengig av hvilken type måling som utføres.

Overgangen til Raman-spektroskopi krever at vi har tilgang til den transisjonelle polarisabilitetstensoren, som kan uttrykkes som en funksjon av den initiale polarisabilitetstensoren α̂i(0) og en rotasjonstransformasjon. Den transisjonelle polarisabilitetstensoren i praksis er knyttet til OH-strekkfunksjonene, og dens verdi kan bestemmes gjennom en serie av approksimasjoner som antar at molekylets polarisabilitet er aksialt symmetrisk i forhold til OH-båndet. Når dette er beregnet, kan man relatere de nødvendige komponentene til tensoren for de forskjellige vibrasjonene (for eksempel strekkfundamentals) og beregne overgangen mellom de vibrerende tilstandene.

For OH-strekkfundamentals er det viktig å merke seg at polarisabilitetstensorens avledede komponenter α′n er relevante for de første ordens bidragene til vibrasjonelle overganger. Den transisjonelle polarisabilitetstensoren for strekkene kan dermed uttrykkes ved en roterende komponent som involverer de relevante aksene til molekylet, som gjør det mulig å forstå hvordan polarisabiliteten endres med vibrasjonen i OH-båndet.

Hva er viktig å forstå i tillegg til det som allerede er skrevet? Når man beregner spektra som disse, er det viktig å merke seg at nøyaktigheten og oppløsningen i beregningene er avgjørende for å få et realistisk og pålitelig spektrum. For IR-spektroskopi, for eksempel, kreves det at tidseriene av korrelasjonsfunksjoner blir tilstrekkelig langt for å oppnå en jevn bølgeform som representerer de fysiske spektrene korrekt. For Raman-spektroskopi, er den spesifikke polarisasjonsavhengigheten kritisk for å få en korrekt beskrivelse av spektrene, spesielt når man arbeider med isotropiske systemer der orienteringen av molekylene ikke er definert på forhånd. I tillegg er det viktig å vurdere hvordan molekylenes strukturelle egenskaper kan påvirke de eksperimentelle resultatene og hvordan vi kan bruke den teoretiske modellen for å sammenligne med eksperimentelle data for å validere de antatte teoriene.

Er strukturen til THF-klatrathydrat ved høyt trykk OS-III eller CS-II?

Studier av tetrahydrofuran (THF) klatrathydrat har avslørt at dens struktur under normale forhold ved atmosfærisk trykk og temperatur over 0 °C følger CS-II-strukturen. Denne strukturen er en av de mest kjente for gasshydrater og har blitt brukt som en modell for bedre å forstå hvordan andre klatrathydrater dannes. Imidlertid ble det i 2004 fremmet en ny hypotese av Kurnosov og hans team, som gjennom bruk av nøytrondiffraksjon i høyt trykk, nærmere bestemt rundt 0.8 GPa ved romtemperatur, foreslo at strukturen til THF-klatrathydrat i dette miljøet var en annen, nærmere bestemt OS-III (ortorombisk system). Denne nye strukturen har siden den gang ikke blitt videre undersøkt eller bekreftet på en bredere skala.

Eksperimentelle og beregningsmessige metoder, som Raman-spektroskopi, kan være avgjørende for å bekrefte OS-III-strukturen. I mitt laboratorium ble THF-klatrathydrat under trykk i en diamantanvilselle (DAC) forberedt, og en enkeltkrystall ble dyrket. Dette eksperimentet viste anisotrop vekst, der krystallene vokste i forskjellige retninger. Dette peker på at strukturen ikke er kubisk, og dette er et viktig tegn på at THF under høyt trykk kan danne en ortorombisk struktur i stedet for den tradisjonelt betraktede kubiske CS-II-strukturen. Veksten i krystallen ble overvåket fra et lavt trykk på omtrent 102 kPa, økt til 1 GPa, og deretter redusert til 0.8 GPa, der krystallen ble stabilisert og analysene ble utført.

Raman-spektroskopiske målinger på enkeltkrystaller viste en rekke karakteristiske band i OH- og CH-strekkbåndene, spesielt rundt 3200-3400 cm−1 og 2900-3000 cm−1, som gir verdifull informasjon om krystallens struktur. Ved å sammenligne eksperimentelle data med beregnede spektra for både CS-II og OS-III, fant vi at de eksperimentelle spektrene bedre stemte overens med den beregnede spektren for OS-III-strukturen. Beregnede Raman-spektre for CS-II viste stor avvik fra de eksperimentelle observasjonene, noe som ytterligere støtter hypotesen om at strukturen i dette høytrykksområdet ikke er CS-II, men faktisk OS-III.

Det er viktig å merke seg at det fortsatt ikke er klart hvordan vekstretningen i krystallen, observert under eksperimentene, forholder seg til de krystallakse. Denne informasjonen kan være avgjørende for å fullt ut forstå strukturen til hydratet under høyt trykk. Imidlertid er det eksperimentelle og beregningsmessige resultatene en tydelig indikasjon på at OS-III er den riktige strukturen for THF-klatrathydrat ved 0.8 GPa, og gir en bekreftelse på de opprinnelige funnene til Kurnosov og hans team.

I tillegg til å bekrefte den nye strukturen, viser disse resultatene også kraften i kombinasjonen av eksperimentelle metoder og beregningsbaserte tilnærminger. Det gir oss en bedre forståelse av hvordan molekylære modeller og eksperimentelle teknikker kan samarbeide for å avdekke detaljer om komplekse systemer som gasshydrater.

Et viktig aspekt å forstå er at disse nye resultatene kan ha bredere implikasjoner for studier av andre gasshydrater og deres strukturer under ekstreme forhold. De kan gi oss dypere innsikt i hvordan andre systemer med lignende molekylære strukturer reagerer på trykk og temperatur. Dette kan igjen være relevant for fremtidige anvendelser, både i industriell sammenheng og i studier relatert til planetarisk vitenskap, hvor hydrater kan finnes under ekstreme forhold på fjerne planeter og måner.

Hvordan elektromagnetisme påvirker vibrasjonsspektroskopi i molekylære systemer

Elektromagnetisme er et omfattende og dyptgående fagområde. Å forstå det fullt ut kan være utfordrende, selv for eksperter, og det er ikke lett å få en fullstendig innsikt i alle aspektene ved det. Denne delen av boken tar for seg en selektiv tilnærming til elektromagnetisme, med fokus på de mest nødvendige ligningene og fysiske størrelsene som kreves for teoretiske beregninger av IR, Raman og sumfrekvensgenereringsspektra for molekylære systemer, i tråd med kapitlene 2 og 3.

I molekylære simuleringer og beregninger av vibrasjonsspektroskopi er det essensielt å forstå hvordan elektromagnetiske enheter fungerer i atomiske enheter, selv om disse enhetene ikke er de mest vanlige i praktisk bruk. For eksempel, når vi beregner den statiske elektriske feltstyrken som påvirker molekylene, brukes Coulombs lov for å bestemme feltstyrken basert på de delvise ladningene i molekylene og avstanden mellom dem. Denne tilnærmingen omformes til atomiske enheter, som gjør det lettere å håndtere de store mengdene av informasjon som kreves for slike beregninger.

Maxwells ligninger utgjør et grunnlag for forståelsen av de elektromagnetiske fenomenene som ligger til grunn for spektroskopiske beregninger. I en kvantemekanisk og klassisk blandet tilnærming er det viktig å kjenne til de makroskopiske variantene av Maxwells ligninger, da de introduserer polarisering, som er en praktisk størrelse i spektroskopiske beregninger. Polariseringen representerer dipolmomentet per volumsenhet i et materiale, og forståelsen av denne størrelsen er avgjørende for å beskrive og analysere signalene som genereres i spektroskopiske eksperimenter.

De makroskopiske Maxwells ligningene omfatter blant annet forholdet mellom elektrisk fluxdensitet (D), elektrisk feltstyrke (E), og polarisering (P). For molekylære systemer som ikke har magnetisering, kan forholdet mellom magnetisk feltstyrke (H) og magnetisk induksjon (B) forenkles til et lineært forhold. Ved å kombinere disse ligningene får man en generell bølgeligning som beskriver hvordan det elektromagnetiske feltet utvikler seg i tid og rom, og hvordan det samhandler med materien.

Gjennom denne tilnærmingen kan man komme fram til uttrykk som beskriver hvordan elektromagnetiske bølger påvirker molekylære systemer, og hvordan materialer reagerer på disse bølgene. Denne reaksjonen, som er et resultat av polariseringen og elektriske strømmer i materialet, utgjør kilden til de spektroskopiske signalene som vi bruker til å utføre vibrasjonsspektroskopi. Det er viktig å merke seg at selv om Maxwells ligninger gir et rammeverk for å forstå disse prosessene, er det i praksis ikke alltid nødvendig å bruke alle detaljene fra disse ligningene når man utfører konkrete beregninger av vibrasjonsspektra.

For beregningene som presenteres i denne boken, er det spesielt nyttig å bruke Coulombs lov i atomiske enheter, og å være klar over de forskjellige forholdene mellom de relevante størrelsene, som dipol–dipol interaksjonsenergi og elektriske felt. Bruken av atomiske enheter gir en enklere tilnærming som kan håndtere de komplekse beregningene knyttet til vibrasjonsspektra, og tillater mer presise og lett forståelige resultater.

Når man forholder seg til det elektromagnetiske feltet og de makroskopiske Maxwells ligningene, er det viktig å ha en solid forståelse av hvordan disse ligningene relaterer seg til molekylære interaksjoner, spesielt i konteksten av vibrasjonsspektroskopi. Dette er nøkkelen til å kunne utføre de nødvendige teoretiske beregningene og få meningsfulle resultater fra eksperimentelle data. Når du utfører beregningene, er det avgjørende å være bevisst på hvilke enheter og konverteringer som er nødvendige for å gjøre beregningene konsistente og korrekte, enten du bruker SI-enheter eller atomiske enheter.

Endelig er det viktig å merke seg at for å oppnå nøyaktige beregninger av vibrasjonsspektra, er det avgjørende å ha god kontroll på forholdet mellom elektriske felt, dipolmoment og polarisering i molekylene, da disse faktorene direkte påvirker hvordan signalene fra molekylene blir målt i eksperimentene. Det er ikke bare spørsmål om matematikk, men også om forståelse av de fysiske prosessene som skjer på mikroskopisk nivå.