Différences
Ci-dessous, les différences entre deux révisions de la page.
| Les deux révisions précédentes Révision précédente Prochaine révision | Révision précédente | ||
|
developpement:productions:logiciels:aquasol:implementation [2013/06/17 14:09] cicaluga [Implémentation de nouvelles fonctionnalites dans le code AquaSol] |
developpement:productions:logiciels:aquasol:implementation [2015/01/07 10:04] (Version actuelle) |
||
|---|---|---|---|
| Ligne 3: | Ligne 3: | ||
| ====== Implémentation de nouvelles fonctionnalites dans le code AquaSol ====== | ====== Implémentation de nouvelles fonctionnalites dans le code AquaSol ====== | ||
| - | * **Contact** : Ralf Everaers, Cerasela Calugaru, Sam Meyer\\ | + | {{:developpement:productions:figureaquasol.png?200 |}} |
| + | * **Contact** : Cerasela Calugaru, Ralf Everaers, Sam Meyer\\ | ||
| * **Objectif** : Optimisation (nouveaux algorithmes, parallélisation), implémentation des nouvelles fonctionnalités, analyse du comportement numérique et informatique du code AquaSol | * **Objectif** : Optimisation (nouveaux algorithmes, parallélisation), implémentation des nouvelles fonctionnalités, analyse du comportement numérique et informatique du code AquaSol | ||
| Ligne 14: | Ligne 15: | ||
| Ainsi, nous avons implémenté plusieurs possibilités et notamment une gaussienne tridimensionnelle globale (qui est définie sur tout le domaine) et une gaussienne localisée (en tronquant au seul domaine sphérique de l’atome, suivie d’une rénormalisation). L’écart type sigma peut être choisi fixé pour tous les atomes (des simulations avec la valeur sigma=0.35 ont été effectuées, cette valeur assurant une distribution de la charge d’environ 99% à l’intérieur de la sphère atomique de rayon unitaire) ou calculé en fonction du rayon atomique. | Ainsi, nous avons implémenté plusieurs possibilités et notamment une gaussienne tridimensionnelle globale (qui est définie sur tout le domaine) et une gaussienne localisée (en tronquant au seul domaine sphérique de l’atome, suivie d’une rénormalisation). L’écart type sigma peut être choisi fixé pour tous les atomes (des simulations avec la valeur sigma=0.35 ont été effectuées, cette valeur assurant une distribution de la charge d’environ 99% à l’intérieur de la sphère atomique de rayon unitaire) ou calculé en fonction du rayon atomique. | ||
| Le choix de la méthode (trilinéaire, sphérique ou gaussiennne) est étroitement lié à la manière de modéliser les charges des atomes fixes à l’échelle d’étude : par une distribution de Dirac au centre atomique ou par une distribution constante dans la sphère atomique ou par une distribution gaussienne tronquée à la sphère atomique. Ainsi, des réflexions ont été menées concernant la capacité des trois modèles à représenter les charges des atomes fixes dans le cas de nos problèmes. En effet, à l’échelle de notre étude (de l’ordre de l’Angström), la modélisation des charges fixes par une distribution gaussienne semble bien adaptée dès lors que la discrétisation spatiale est suffisamment fine (ce qui implique un nombre suffisant de points dans la sphère atomique). | Le choix de la méthode (trilinéaire, sphérique ou gaussiennne) est étroitement lié à la manière de modéliser les charges des atomes fixes à l’échelle d’étude : par une distribution de Dirac au centre atomique ou par une distribution constante dans la sphère atomique ou par une distribution gaussienne tronquée à la sphère atomique. Ainsi, des réflexions ont été menées concernant la capacité des trois modèles à représenter les charges des atomes fixes dans le cas de nos problèmes. En effet, à l’échelle de notre étude (de l’ordre de l’Angström), la modélisation des charges fixes par une distribution gaussienne semble bien adaptée dès lors que la discrétisation spatiale est suffisamment fine (ce qui implique un nombre suffisant de points dans la sphère atomique). | ||
| + | |||
| + | ➢ **Charges fixes, charges des ions libres, charges correspondant a la divergence de la polarisation locale** : calcul de toutes les charges intervenant dans l’équation et obtention des fichiers de sortie avec des cartes des toutes ces charges | ||
| + | |||
| + | ➢ **Energie électrostatique** : calcul de l’énergie électrostatique avec la formule basée sur la connaissance du gradient du potentiel électrostatique (en plus de la méthode «self-energy» déjà implémentée) | ||