Bra-ket notation, utvecklad av Dirac, utgör en grundläggande ram för att beskriva kvanttillstånd och operatorer med hjälp av linjär algebra. Den skiljer sig från den klassiska vågfunktionsbeskrivningen genom att använda vektorer (kets |·〉) och duala vektorer (bras 〈·|) i Hilbertrummet. En ket |a〉 representerar en kvanttillståndsvektor i kolonnform medan motsvarande bra 〈a| är dess komplexkonjugerade transponat, alltså radvektorn. Denna dualitet är avgörande för beräkningar och uttryck för inre produkter 〈a|b〉, vilka är komplexa tal men ger alltid reella förväntningsvärden när operatorer representerar mätbara storheter.

Operatorer, som betecknas med hatt (ˆ), är linjära transformationer på dessa tillståndsvektorer och motsvarar observerbara fysiska storheter, såsom energi (Hamiltonoperatorn Ĥ), position (x̂) och rörelsemängd (p̂). Dessa operatorer är Hermitiska, vilket betyder att de är lika med sin egen adjungat († = Â), och garanterar att deras egenvärden, vilka representerar möjliga mätvärden, är reella tal. Kommutatorn [Â, B̂] = ÂB̂ − B̂ är en central konstruktion som anger i vilken mån två operatorer är kompatibla; speciellt är position och rörelsemängd icke-kommutativa med kommutatorn [x̂, p̂] = iħ, vilket är grundläggande för osäkerhetsprincipen.

Inom en given bas av ortonormerade eigenkets |a_i〉 gäller slutehetsrelationen ∑_i |a_i〉〈a_i| = 1, vilket innebär att varje tillstånd kan uttryckas som en linjärkombination av dessa basvektorer. Detta förenklar representation och beräkning av operatorers verkan.

Den harmoniska oscillatorn, en fundamental modell inom kvantmekanik, beskrivs av Hamiltonoperatorn Ĥ = p̂²/(2m) + (1/2)mω²x̂², där m är massan och ω den vinkelfrekvens som karakteriserar systemets bindning. Lösningen till Schrödingerekvationen för denna oscillator erhålls elegant genom att införa annihilations- (â) och skapandeoperatorer (â†), vilka kombinerar position och rörelsemängdsoperatorer på ett sådant sätt att de höjer eller sänker oscillatorns vibrationskvanttal v.

Dessa operatorer uppfyller kommutationsrelationen [â, â†] = 1, och Hamiltonian kan skrivas som Ĥ = ħω(â†â + 1/2), vilket tydligt visar att energin är kvantiserad i diskreta nivåer Ev = ħω(v + 1/2). Genom användning av â och ↠kan man konstruera energinivåernas egenvektorer |v〉, där â sänker nivån och ↠höjer den, vilket underlättar förståelsen av tillståndens struktur och övergångar utan att lösa differentialekvationer direkt.

Det är avgörande att förstå att mätresultat i kvantmekanik alltid är reella tal som framkommer från egenvärden till Hermitiska operatorer, och inte själva operatorerna eller tillstånden. Bra-ket notationens komplexa tal och operatorers icke-kommutativa egenskaper ger upphov till fenomen som osäkerhetsprincipen, som i sin tur formar det fundamentala skiljelinjen mellan klassisk och kvantmekanik.

Det är också viktigt att inse att det praktiska arbetet ofta sker inom ett begränsat delrum av hela Hilbertrummet, där man väljer en ändlig uppsättning basvektorer för att approximera tillstånd. Denna begränsning gör det möjligt att hantera komplexa system på ett hanterbart sätt samtidigt som slutehetsrelationen säkerställer att alla tillstånd kan uttryckas inom detta delrum.

För en djupare förståelse krävs insikten att operatorernas algebra och deras kommutationsrelationer inte bara är matematiska konstruktioner utan speglar de fysiska begränsningarna och möjligheterna i kvantvärlden. Till exempel är det genom förståelsen av annihilations- och skapandeoperatorernas funktion som moderna kvantfältteorier och avancerade metoder inom kvantkemin utvecklats, vilket visar att denna abstrakta notation har direkt koppling till experiment och teknologi.

Hur representeras det vibrationella Hamiltonian för vatten i kondenserade faser?

I kvantmekanisk beskrivning av molekylära vibrationer är val av basfunktioner en grundläggande aspekt. För vatten, i kondenserad fas, används en representation där varje basfunktion (ket) specificeras av vibrationskvanttalen för varje lokal mode – stretching och bending. Stretch-fundamentaler betecknas som tillstånd där en enskild OH-stretch är exciterad medan alla andra lokalmoderna är i sitt grundtillstånd. Motsvarande gäller för bend-fundamentaler och bend-overtone. Det totala antalet basfunktioner som behövs för att beskriva systemet beror på hur många vattenmolekyler (N) som inkluderas i simuleringen, där varje vattenmolekyl ger upphov till två stretch- och en bend-mode. Totalt används 4N basfunktioner om både fundamentaler och overtones inkluderas. Om endast OH-stretchregionen är av intresse, reduceras denna bas till 3N genom att utesluta bend-fundamentaler. Om endast HOH-bendregionen är relevant, reduceras basen ytterligare till N.

Genom att införa ett systematiskt numreringsschema för varje basfunktion – till exempel |1〉 till |2N〉 för stretch-fundamentaler, |2N+1〉 till |3N〉 för bend-overtone och |3N+1〉 till |4N〉 för bend-fundamentaler – möjliggörs en kompakt och praktisk matrisrepresentation av det vibrationella Hamiltonian. Schrödingerekvationen löses numeriskt inom detta basrum, vilket medför att spektroskopiska övertoner och kombinationstillstånd utanför denna bas inte kan reproduceras inom denna modell.

Hamiltonmatrisens element delas in i diagonala och icke-diagonala. De diagonala elementen representerar energi för varje enskilt tillstånd, justerade för att exkludera nollpunktsenergin som inte påverkar IR-, Raman- eller SFG-spektra där endast energiövergångar mellan tillstånd är relevanta. För stretch-fundamentaler beror varje diagonalt element på den specifika frekvensen för den exciterade OH-bindningen, ωₛ, medan bending-tillstånden påverkas av motsvarande ω_b. De icke-diagonala elementen inkluderar kopplingar mellan stretch–stretch, bend–bend samt Fermi-resonanskopplingar mellan stretch och bend.

Speciellt intressant i kondenserade faser är hur omgivningen påverkar frekvensen för OH-stretch. Denna frekvens, i v = 0 → 1-övergången, är starkt beroende av vätebindningsstatusen. En OH-grupp som donerar en vätebindning visar en markant rödskiftning. Denna skiftning är inte enbart beroende av styrkan i vätebindningen, utan även av vätebindningsstatusen hos acceptormolekylen samt elektronfördelningen kring syreatomen i donatormolekylen – en effekt som benämns vätebindningskooperativitet.

Att exakt beräkna ωₛ för varje ögonblick i en molekyldynamisk (MD) simulering genom kvantkemiska metoder är praktiskt ogenomförbart på grund av den enorma beräkningskostnaden. För att kringgå detta utvecklade Skinnergruppen vibrationsspektroskopiska kartor. Dessa kartor härleddes ur DFT-beräkningar av vattenkluster i olika mikromiljöer, där frekvensen ωₛ relaterades till det elektriska fältet E längs varje OH-bindning. Det visade sig att E kan användas som deskriptor i en kvadratisk funktion för att approximerat uttrycka ωₛ. Detta elektriska fält kan i sin tur beräknas från de elektrostatiska bidragen från omgivande atomer i varje MD-tidssteg.

Därmed kan man i varje ögonblick i simuleringen, utifrån atomära koordinater och partialladdningar, beräkna E för varje OH-bindning och därigenom bestämma aktuell ωₛ enligt den spektrala kartan. Resultatet blir en dynamisk Hamiltonmatris där varje element förändras över tiden beroende på molekylens miljö, vilket gör det möjligt att rekonstruera spektra med hög tidsupplösning.

Värt att beakta är att denna metod, trots sin approximativa karaktär, lyckas fånga viktiga fysiska fenomen såsom vätebindningskooperativitet och miljöberoende spektralförskjutningar – något som är kritiskt för att tolka spektra från vatten i biologiska system, is, lösningar eller andra kondenserade faser.

Att förstå dessa vibrationella interaktioner i termen av ett Hamiltonianskt ramverk möjliggör inte enbart reproduktion av spektra, utan även en djupare insikt i de intermolekylära krafter som styr vätebindningsnätverkets struktur och dynamik. Det är alltså inte bara en numerisk förenkling, utan en fysisk modellering av hög generalitet.

Fortsatt utveckling kan inkludera högre ordningens tillstånd, såsom stretch–bend-kombinationer och övertoner, genom att utöka basen och anpassa Hamiltonmatrisen därefter. Även inkluderandet av anharmoniska termer och koppling till elektroniska tillstånd kan utveckla modellen mot mer realistiska spektrala simuleringar.