gmx-tcaf - Online in de Cloud

PROGRAMMA:

NAAM

gmx-tcaf - Bereken viscositeiten van vloeistoffen

KORTE INHOUD

gmx tcaf [-f [<.trr/.cpt/...>]] [-s [<.tpr/.gro/...>]] [-n [<.ndx>]]
[-ot [<.xvg>]] [-oa [<.xvg>]] [-o [<.xvg>]] [-van [<.xvg>]]
[-oc [<.xvg>]] [-ov [<.xvg>]] [-b ] [-e ]
[-DT ] [-[nu] [-xvg ] [-[geen]mol] [-[geen]k34]
[-gew ] [-acflen ] [-[nee]normaliseren] [-P ]
[-fitfn ] [-beginnen ] [-eindfit ]

PRODUCTBESCHRIJVING

gmx tcaf berekent transversale huidige autocorrelaties. Deze worden gebruikt om de
schuifviscositeit, eta. Voor details zie: Palmer, Phys. Rev. E 49 (1994) pp 359-366.

Dwarsstromen worden berekend met behulp van de k-vectoren (1,0,0) en (2,0,0), elk ook in
the y- En z-richting, (1,1,0) en (1,-1,0) elk ook in de 2 andere vlakken (deze
vectoren zijn niet onafhankelijk) en (1,1,1) en de 3 andere doosdiagonalen (ook niet
onafhankelijk). Voor elke k-vector worden de sinus en cosinus gebruikt, in combinatie met de
snelheid in 2 loodrechte richtingen. Dit geeft een totaal van 16*2*2=64 transversale
stromen. Voor elke k-vector wordt één autocorrelatie berekend, wat 16 oplevert
TCAF's. Elk van deze TCAF's is aangepast aan f(t) = exp(-v)(cosh(Wv) + 1/W sinh(Wv)), v = -t/(2
tau), W = sqrt(1 - 4 tau eta/rho k^2), wat 16 waarden van tau en eta oplevert. De fit
gewichten vervallen exponentieel met tijdconstante w (gegeven met -gew) als exp(-t/w), en de
TCAF en fit worden berekend tot tijd 5*w. De eta-waarden moeten worden aangepast aan 1 - a
eta(k) k^2, waaruit men de schuifviscositeit bij k=0 kan schatten.

Als de doos kubusvormig is, kan men de optie gebruiken -oc, die de TCAF's over alle
k-vectoren met dezelfde lengte. Dit resulteert in nauwkeurigere TCAF's. Zowel de kubieke
TCAF's en fits zijn geschreven om -oc De kubieke eta-schattingen worden ook geschreven naar -ov.

Met optie -mol, de transversale stroom wordt bepaald door moleculen in plaats van atomen. In
In dit geval zou de indexgroep moeten bestaan ​​uit molecuulnummers in plaats van atoomnummers.

De k-afhankelijke viscositeiten in de -ov bestand moet worden aangepast aan eta(k) = eta_0 (1 - ak^2)
om de viscositeit bij oneindige golflengte te verkrijgen.

Opmerking: zorg ervoor dat je coördinaten en snelheden vaak genoeg schrijft. De beginletter,
niet-exponentieel, een deel van de autocorrelatiefunctie is erg belangrijk voor het verkrijgen van een
past goed.

OPTIES

Opties om invoerbestanden op te geven:

-f [<.trr/.cpt/...>] (traj.trr)
Volledig precisietraject: TRR cpt tng

-s [<.tpr/.gro/...>] (topol.tpr) (Optioneel)
Structuur+massa(db): TPR gro g96 pdb br ent

-n [<.ndx>] (index.ndx) (Optioneel)
Indexbestand

Opties om uitvoerbestanden te specificeren:

-ot [<.xvg>] (transcur.xvg) (Optioneel)
xvgr/xmgr-bestand

-oa [<.xvg>] (tcaf_all.xvg)
xvgr/xmgr-bestand

-o [<.xvg>] (tcaf.xvg)
xvgr/xmgr-bestand

-van [<.xvg>] (tcaf_fit.xvg)
xvgr/xmgr-bestand

-oc [<.xvg>] (tcaf_cub.xvg) (Optioneel)
xvgr/xmgr-bestand

-ov [<.xvg>] (visc_k.xvg)
xvgr/xmgr-bestand

Andere opties:

-b (0)
Eerste frame (ps) om van traject te lezen

-e (0)
Laatste frame (ps) om van traject te lezen

-DT (0)
Gebruik frame alleen als t MOD dt = eerste keer (ps)

-[nu (Nee)
Uitvoer bekijken .xvg, .xpm, .eps en .pdb bestanden

-xvg
xvg-plotopmaak: xmgrace, xmgr, geen

-[geen]mol (Nee)
Bereken TCAF van moleculen

-[geen]k34 (Nee)
Gebruik ook k=(3,0,0) en k=(4,0,0)

-gew (5)
Exponentiële vervaltijd voor de TCAF-fitgewichten

-acflen (-1)
Lengte van de ACF, standaard is de helft van het aantal frames

-[nee]normaliseren (Ja)
ACF normaliseren

-P (0)
Order of Legendre polynoom voor ACF (0 geeft niets aan): 0, 1, 2, 3

-fitfn (Geen)
Fit-functie: geen, exp, aexp, exp_exp, exp5, exp7, exp9

-beginnen (0)
Tijd waar de exponentiële aanpassing van de correlatiefunctie moet beginnen

-eindfit (-1)
Tijd waar de exponentiële fit van de correlatiefunctie moet eindigen, -1 is tot de
einde

Gebruik gmx-tcaf online met behulp van onworks.net-services