Retour au menu principal

## Mode : fonction de distribution radiale

#### Syntaxe

`atomsk --rdf <fichier> <dR> [options]`

#### Description

Dans ce mode Atomsk calcule la ou les fonction(s) de distribution radiale (radial distribution function, RDF) des systèmes contenus dans le <fichierliste>.

L'analyse est faite pour toutes les "peaux" [R;R+dR] où R varie de 0 jusqu'à <Rmax>, et est moyennée sur tous les atomes (moyenne spatiale). Le rayon maximal <Rmax> et l'épaisseur <dR> de la "peau" doivent être fournis (en Å). Le rayon maximal <Rmax> peut être supérieur à la taille du système.

Si le <fichier> contient un système atomique (dans l'un des formats supportés par Atomsk), alors la RDF est calculée pour ce système.

Si le <fichier> contient une liste de fichiers à traiter, alors la RDF est calculée pour chacun de ces fichiers, et moyennée sur toutes ces configurations, autrement dit il s'agit d'une moyenne spatiale et temporelle. Les fichiers devraient tous correspondre au même système (par exemple, des instantanés d'une simulation en dynamique moléculaire) pour que le calcul ait du sens.

Si les systèmes contiennent plusieurs espèces chimiques alors les fonctions de distribution radiale partielles sont calculées. Par exemple si des atomes A et B existent, alors Atomsk calculera les fonctions de distribution des atomes A autour des atomes A, de B autour de A, et B autour de B.

Après le calcul les fichiers suivants sont produits (<A> et <B> sont remplacés par les symboles atomiques) :

- *_rdf_<A><B>.dat: fichiers contenant les fonctions de ditribution radiale des atomes B autour des atomes A en fonction de la distance.

Si ce mode est appelé avec une ou plusieurs options alors elles seront appliquées au système avant le calcul des fonctions de distribution.

#### Performance

Ce mode utilise un algorithme de recherche de voisins, qui n'est malheureusement pas très optimisé dans Atomsk. Le temps nécessaire pour la recherche de voisin augmente linéairement avec le volume (donc en R3). Il est donc conseillé d'utiliser une valeur de Rmax aussi petite que possible, typiquement entre 10 et 20 Å. Pour gagner du temps, la liste de voisins n'est reconstruite que si le déplacement maximal d'un atome dépasse la valeur de <dR>. Après cela, le temps nécessaire au calcul de la RDF en elle-même augmente linéairement avec le pas dR, et avec le nombre de particules.

#### Exemples

- `atomsk --rdf systeme.cfg 10 0.1`

  Ceci calculera la fonction de distribution radiale pour le système contenu dans `systeme.cfg`, de 0 à 10 Å, en utilisant un pas dR=0.1 Å.

- `atomsk --rdf liste.lst 10 0.1`

  Ceci calculera les fonctions de distribution radiale pour tous les systèmes contenus dans `liste.lst`, de 0 à 10 Å, en utilisant un pas dR=0.1 Å.

Retour au menu principal