U području molekulske dinamike često se susreću dvije naizgled suprotstavljene tvrdnje koje, paradoksalno, koegzistiraju među ozbiljnim istraživačima. S jedne strane stoji uvjerenje da su simulacije dovoljno precizne za predviđanje ponašanja molekula u gotovo svim kemijskim uvjetima, dok se s druge strane zagovara stav da su takve simulacije uvijek kompromis između točnosti i računalne izvedivosti te stoga nikada ne mogu u potpunosti zamijeniti eksperiment. Tijekom vlastitog iskustva svjedočio sam ispreplitanju ovih pozicija bilo je istraživača koji su bez zadrške vjerovali u rezultate simulacija bez dodatne validacije, što je znalo dovesti do suptilnih pogrešaka u tumačenju reakcijskih mehanizama ili termodinamičkih svojstava, osobito u složenim sustavima poput proteinskih sklopova ili katalitičkih površina.
Na molekulskoj razini, simulacije molekulske dinamike (MD) temelje se na numeričkom rješavanju Newtonovih jednadžbi gibanja za skupinu atoma ili molekula. Svaka čestica tretira se kao klasični objekt kojim upravljaju sile određene potencijalnim funkcijama, koje modeliraju različite vrste interakcija: kovalentne veze, kutove, dihedrale, ali i van der Waalsove i Coulombove sile između ne-kovalentno vezanih atoma. Za razliku od klasične kemijske kinetike ili termodinamike koje opisuju prosječna svojstva sustava u ravnoteži ili ukupnu brzinu reakcije, MD simulacije omogućuju praćenje putanja pojedinačnih atoma s vremenskom rezolucijom reda femtosekundi ($10^{-15}$ s). To otkriva detalje o prijelaznim stanjima i međumolekulskim interakcijama koje bi inače ostale neprimijećene.
Ono što posebno izaziva pažnju jest koliko stvarnost odstupa od idealiziranih modela koje nalazimo u udžbenicima. U jednom projektu gdje smo proučavali adsorpciju malih organskih molekula na metalnu površinu zlatnog katalizatora, standardni force field modeli nisu uspjeli zadovoljavajuće opisati polarizabilnost elektronske gustoće na granici faza. Posljedica toga bilo je precjenjivanje afiniteta adsorpcije za gotovo 20 %, što nije zanemarivo ako se cilja dizajn katalizatora s optimalnom selektivnošću. Tek nakon uključivanja kvantnomehaničkih parametara iz ab initio simulacija kao korekcije u klasični MD model (tzv. QM/MM pristup), dobili smo rezultate koji se slažu s eksperimentalnim podacima o adsorpciji pri temperaturi od $298\,K$ i koncentraciji organske tvari od $0.1\,mol/L$.
Važno je također imati na umu da povezivanje strukture i svojstava unutar MD simulacija zahtijeva razumijevanje kako promjene konformacije utječu na kinetiku i termodinamiku procesa. Primjerice, fleksibilnost lanca polimera može direktno utjecati na njegovu toplinsku stabilnost i mehanička svojstva. To povezujemo putem analize slobodne energije promjene konformacije $$\Delta G = \Delta H - T \Delta S,$$ gdje entalpijski doprinos $\Delta H$ proizlazi iz vodikovih veza i drugih interakcija unutar molekule, dok entropijski dio $\Delta S$ odražava stupanj slobode gibanja lanaca.
Zanimljivu kemijsku anomaliju primijetio sam tijekom simulacija vode u uskim nano-poredima grafenskih membrana: iako tradicionalna termodinamika sugerira da bi voda spontano prolazila kroz ovakve kanale zbog hidrofobnih efekata smanjenja površinske napetosti, naše MD simulacije pokazuju da specifične međumolekulske orijentacije i formiranje lanaca vodenih mostova stvaraju energijsku barijeru kretanju koja značajno smanjuje propusnost za red veličine nižu nego što bismo očekivali prema pojednostavljenim modelima.
Da ilustriram praktičnu primjenu MD simulacija kvantitativnim podatcima: razmotrimo reakciju prijenosa protona između dviju karboksilnih skupina u vodenom okruženju:
$$\mathrm{R-COOH} + \mathrm{R'-COO^-} \rightleftharpoons \mathrm{R-COO^-} + \mathrm{R'-COOH}.$$
Koristeći MD simulaciju parametriranu eksperimentalnim podacima pri pH 7 i temperaturi $310\,K$, moguće je izračunati slobodnu energiju prijenosa $\Delta G^\ddagger$, odnosno aktivacijsku barijeru reakcije. Primjerice, metode slobodne energije poput umbrella samplinga mogle su nam dati vrijednost $\Delta G^\ddagger = 35\,kJ/mol$. Ovaj rezultat ukazuje na umjerenu kinetičku barijeru za protonski prijenos unutar vodene okoline dovoljnu za relativno brzu izmjenu protona, ali ne trenutačnu što odgovara eksperimentalnim mjerenjima sličnih sustava.
S druge strane konstanta ravnoteže $K$ može se izraziti termodinamički kao
$$K = e^{-\frac{\Delta G}{RT}},$$
gdje $\Delta G$ označava razliku slobodnih energija između produkata i reaktanata. Time dobivamo kemijski smisao predviđanja smjera reakcije pod određenim uvjetima koristeći mikroskopsku sliku iz simulacija.
Sve navedeno pokazuje da matematički opisi interakcija sami po sebi nisu dovoljni; nužno je rigorozno procijeniti koliko naši modeli reflektiraju stvarne kemijske uvjete uključujući temperaturu, koncentraciju iona te lokalnu polarizabilnost medija. Samo takva integrirana perspektiva može pružiti pouzdane predikcije.
Pitanje koje leži duboko ispod ove problematike jest ono što povijest znanosti često ističe: čak ni najsofisticiraniji modeli teško hvataju složenost stvarnosti jer ona nije samo kompleksna već otpornija na formalizaciju. Kako onda kontinuirano unaprijediti metode tako da one ne budu tek preslika prošlosti nego pouzdan alat za inovacije u kemijskoj industriji? To zahtijeva neprestano preispitivanje temelja naših pristupa i otvorenost novim paradigmama lekcija koju sam naučio kroz spoj teorijske fizike i eksperimentalne kemije tijekom svoje akademske formacije.
Generira se sažetak…