Hydrogen-disordered ice generator
Project description
GenIce
======
Generate proton-disordered ice structures for GROMACS.
Requirements
------------
- Python 3
- NetworkX
- numpy
Note: WinPython includes all of these requirements. ##Installation
GenIce is registered to PyPI (Python Package Index). Install with pip3.
::
pip3 install genice
Uninstallation
--------------
::
pip3 uninstall genice
Usage
-----
::
usage: genice [-h] [--rep REP REP REP] [--dens DENS] [--seed SEED]
[--format gmeqdXoc] [--water model] [--guest D=empty] [--debug]
[--quiet]
Type
positional arguments:
Type Crystal type (1c,1h,etc.)
optional arguments:
-h, --help show this help message and exit
--rep REP REP REP, -r REP REP REP
Repeat the unit cell in x,y, and z directions. [2,2,2]
--dens DENS, -d DENS Specify the ice density in g/cm3
--seed SEED, -s SEED Random seed [1000]
--format gmeqdXoc, -f gmeqdXoc
Specify file format [g(romacs)|m(dview)|e(uler)|q(uate
rnion)|d(igraph)|o(penScad)|c(entersofmass)]
--water model, -w model
Specify water model. (tip3p, tip4p, etc.)
--guest D=empty, -g D=empty
Specify guest in the cage. (D=empty, T=co2, etc.)
--debug, -D Output debugging info.
--quiet, -q Do not output progress messages.
Example
-------
- To make a 3x3x3 units of a hydrogen-disordered ice IV (4) of TIP4P
water in GROMACS .gro format:
::
genice --water tip4p --rep 3 3 3 4 > ice4.gro
- To make a CS1 clathrate hydrate structure of TIP4P water containing
CO2 in GROMACS .gro format:
::
genice -g 12=co2 -g 14=co2 --water tip4p CS1 > cs1.gro
- To make a 2x2x4 units of CS2 clathrate hydrate structure of TIP4P
water containing THF (united atom with a dummy site) in the large
cage in GROMACS .gro format:
::
genice -g 16=uathf6 --water tip4p --rep 2 2 4 CS2 > cs2-224.gro
Structure generation
--------------------
The program generates various ice lattice with proton disorder and
without defect. Total dipole moment is always set to zero. The minimal
structure (with --rep 1 1 1 option) is not always the unit cell of the
lattice because it is difficult to deal with the hydrogen bond network
topology of tiny lattice under periodic boundary condition. Note that
the generated structure is not optimal according to the potential
energy.
Ice structures
--------------
+----------+---------------+-------------+
| Symbol | Description | References |
+==========+===============+=============+
| 1h, 1c | Most popular | |
| | Ice I | |
| | (hexagonal or | |
| | cubic) | |
+----------+---------------+-------------+
| 2d | Hypothetical | Nakamura, |
| | Proton-disord | Tatsuya et |
| | ered | al. |
| | Ice II. | “Thermodyna |
| | | mic |
| | | Stability |
| | | of Ice II |
| | | and Its |
| | | Hydrogen-Di |
| | | sordered |
| | | Counterpart |
| | | : |
| | | Role of |
| | | Zero-Point |
| | | Energy.” |
| | | The Journal |
| | | of Physical |
| | | Chemistry B |
| | | 120.8 |
| | | (2015): |
| | | 1843–1848. |
| | | Web. |
+----------+---------------+-------------+
| 3, 4, 6, | Conventional | |
| 7, 12 | high-pressure | |
| | ices III, IV, | |
| | VI, VII, and | |
| | XII. | |
+----------+---------------+-------------+
| 5 | Monoclinic | |
| | ice V | |
| | (testing). | |
+----------+---------------+-------------+
| 16 | Negative-pres | Falenty, |
| | sure | A., Hansen, |
| | ice XVI(16). | T. C. & |
| | | Kuhs, W. F. |
| | | Formation |
| | | and |
| | | properties |
| | | of ice XVI |
| | | obtained by |
| | | emptying a |
| | | type sII |
| | | clathrate |
| | | hydrate. |
| | | Nature 516, |
| | | 231-233 |
| | | (2014). |
+----------+---------------+-------------+
| 17 | Negative-pres | del Rosso, |
| | sure | Leonardo, |
| | ice XVII(17). | Milva |
| | | Celli, and |
| | | Lorenzo |
| | | Ulivi. “Ice |
| | | XVII as a |
| | | Novel |
| | | Material |
| | | for |
| | | Hydrogen |
| | | Storage.” |
| | | Challenges |
| | | 8.1 (2017): |
| | | 3. |
+----------+---------------+-------------+
| 0 | Hypothetical | Russo, J., |
| | ice "0". | Romano, F. |
| | | & Tanaka, |
| | | H. New |
| | | metastable |
| | | form of ice |
| | | and its |
| | | role in the |
| | | homogeneous |
| | | crystalliza |
| | | tion |
| | | of water. |
| | | Nat Mater |
| | | 13, 733-739 |
| | | (2014). |
+----------+---------------+-------------+
| i | Hypothetical | Fennell, C. |
| | ice "i". | J. & |
| | | Gezelter, |
| | | J. D. |
| | | Computation |
| | | al |
| | | Free Energy |
| | | Studies of |
| | | a New Ice |
| | | Polymorph |
| | | Which |
| | | Exhibits |
| | | Greater |
| | | Stability |
| | | than Ice I |
| | | h. J. Chem. |
| | | Theory |
| | | Comput. 1, |
| | | 662-667 |
| | | (2005). |
+----------+---------------+-------------+
| C0-II | Filled ice C0 | Smirnov, G. |
| | (Alias of | S. & |
| | 17). | Stegailov, |
| | | V. V. |
| | | Toward |
| | | Determinati |
| | | on |
| | | of the New |
| | | Hydrogen |
| | | Hydrate |
| | | Clathrate |
| | | Structures. |
| | | J Phys Chem |
| | | Lett 4, |
| | | 3560-3564 |
| | | (2013). |
+----------+---------------+-------------+
| C1 | Filled ice C1 |
| | (Alias of |
| | 2d). |
+----------+---------------+-------------+
| C2 | Filled ice C2 |
| | (Alias of |
| | 1c). |
+----------+---------------+-------------+
| sTprime | Filled ice | Smirnov, G. |
| | sT' | S. & |
| | | Stegailov, |
| | | V. V. |
| | | Toward |
| | | Determinati |
| | | on |
| | | of the New |
| | | Hydrogen |
| | | Hydrate |
| | | Clathrate |
| | | Structures. |
| | | J Phys Chem |
| | | Lett 4, |
| | | 3560-3564 |
| | | (2013). |
+----------+---------------+-------------+
| CS1, | Clathrate | Matsumoto, |
| CS2, | hydrates CS1 | M. & |
| TS1, HS1 | (sI), CS2 | Tanaka, H. |
| | (sII), TS1 | On the |
| | (sIII), and | structure |
| | HS1 (sIV). | selectivity |
| | | of |
| | | clathrate |
| | | hydrates. |
| | | J. Phys. |
| | | Chem. B |
| | | 115, |
| | | 8257-8265 |
| | | (2011). |
+----------+---------------+-------------+
| RHO | Hypothetical | Huang, Y et |
| | ice at | al. “A New |
| | negative | Phase |
| | pressure ice | Diagram of |
| | 'sIII'. | Water Under |
| | | Negative |
| | | Pressure: |
| | | the Rise of |
| | | the |
| | | Lowest-Dens |
| | | ity |
| | | Clathrate |
| | | S-III.” |
| | | Science |
| | | Advances |
| | | 2.2 (2016): |
| | | e1501010–e1 |
| | | 501010. |
+----------+---------------+-------------+
Please ask vitroid@gmail.com to add new ice structures. ##Water models
+-------------+-------------------+
| symbol | type |
+=============+===================+
| ``tip3p`` | TIP3P (default) |
+-------------+-------------------+
| ``tip4p`` | TIP4P |
+-------------+-------------------+
| ``tip5p`` | TIP5P |
+-------------+-------------------+
Guest molecules
---------------
+-----------------------------------------+--------------------------+
| symbol | type |
+=========================================+==========================+
| ``co2`` | CO2 |
+-----------------------------------------+--------------------------+
| ``uathf`` | United atom 5-site THF |
+-----------------------------------------+--------------------------+
| ``g12``,\ ``g14``,\ ``g15``,\ ``g16`` | A monatomic dummy site |
+-----------------------------------------+--------------------------+
======
Generate proton-disordered ice structures for GROMACS.
Requirements
------------
- Python 3
- NetworkX
- numpy
Note: WinPython includes all of these requirements. ##Installation
GenIce is registered to PyPI (Python Package Index). Install with pip3.
::
pip3 install genice
Uninstallation
--------------
::
pip3 uninstall genice
Usage
-----
::
usage: genice [-h] [--rep REP REP REP] [--dens DENS] [--seed SEED]
[--format gmeqdXoc] [--water model] [--guest D=empty] [--debug]
[--quiet]
Type
positional arguments:
Type Crystal type (1c,1h,etc.)
optional arguments:
-h, --help show this help message and exit
--rep REP REP REP, -r REP REP REP
Repeat the unit cell in x,y, and z directions. [2,2,2]
--dens DENS, -d DENS Specify the ice density in g/cm3
--seed SEED, -s SEED Random seed [1000]
--format gmeqdXoc, -f gmeqdXoc
Specify file format [g(romacs)|m(dview)|e(uler)|q(uate
rnion)|d(igraph)|o(penScad)|c(entersofmass)]
--water model, -w model
Specify water model. (tip3p, tip4p, etc.)
--guest D=empty, -g D=empty
Specify guest in the cage. (D=empty, T=co2, etc.)
--debug, -D Output debugging info.
--quiet, -q Do not output progress messages.
Example
-------
- To make a 3x3x3 units of a hydrogen-disordered ice IV (4) of TIP4P
water in GROMACS .gro format:
::
genice --water tip4p --rep 3 3 3 4 > ice4.gro
- To make a CS1 clathrate hydrate structure of TIP4P water containing
CO2 in GROMACS .gro format:
::
genice -g 12=co2 -g 14=co2 --water tip4p CS1 > cs1.gro
- To make a 2x2x4 units of CS2 clathrate hydrate structure of TIP4P
water containing THF (united atom with a dummy site) in the large
cage in GROMACS .gro format:
::
genice -g 16=uathf6 --water tip4p --rep 2 2 4 CS2 > cs2-224.gro
Structure generation
--------------------
The program generates various ice lattice with proton disorder and
without defect. Total dipole moment is always set to zero. The minimal
structure (with --rep 1 1 1 option) is not always the unit cell of the
lattice because it is difficult to deal with the hydrogen bond network
topology of tiny lattice under periodic boundary condition. Note that
the generated structure is not optimal according to the potential
energy.
Ice structures
--------------
+----------+---------------+-------------+
| Symbol | Description | References |
+==========+===============+=============+
| 1h, 1c | Most popular | |
| | Ice I | |
| | (hexagonal or | |
| | cubic) | |
+----------+---------------+-------------+
| 2d | Hypothetical | Nakamura, |
| | Proton-disord | Tatsuya et |
| | ered | al. |
| | Ice II. | “Thermodyna |
| | | mic |
| | | Stability |
| | | of Ice II |
| | | and Its |
| | | Hydrogen-Di |
| | | sordered |
| | | Counterpart |
| | | : |
| | | Role of |
| | | Zero-Point |
| | | Energy.” |
| | | The Journal |
| | | of Physical |
| | | Chemistry B |
| | | 120.8 |
| | | (2015): |
| | | 1843–1848. |
| | | Web. |
+----------+---------------+-------------+
| 3, 4, 6, | Conventional | |
| 7, 12 | high-pressure | |
| | ices III, IV, | |
| | VI, VII, and | |
| | XII. | |
+----------+---------------+-------------+
| 5 | Monoclinic | |
| | ice V | |
| | (testing). | |
+----------+---------------+-------------+
| 16 | Negative-pres | Falenty, |
| | sure | A., Hansen, |
| | ice XVI(16). | T. C. & |
| | | Kuhs, W. F. |
| | | Formation |
| | | and |
| | | properties |
| | | of ice XVI |
| | | obtained by |
| | | emptying a |
| | | type sII |
| | | clathrate |
| | | hydrate. |
| | | Nature 516, |
| | | 231-233 |
| | | (2014). |
+----------+---------------+-------------+
| 17 | Negative-pres | del Rosso, |
| | sure | Leonardo, |
| | ice XVII(17). | Milva |
| | | Celli, and |
| | | Lorenzo |
| | | Ulivi. “Ice |
| | | XVII as a |
| | | Novel |
| | | Material |
| | | for |
| | | Hydrogen |
| | | Storage.” |
| | | Challenges |
| | | 8.1 (2017): |
| | | 3. |
+----------+---------------+-------------+
| 0 | Hypothetical | Russo, J., |
| | ice "0". | Romano, F. |
| | | & Tanaka, |
| | | H. New |
| | | metastable |
| | | form of ice |
| | | and its |
| | | role in the |
| | | homogeneous |
| | | crystalliza |
| | | tion |
| | | of water. |
| | | Nat Mater |
| | | 13, 733-739 |
| | | (2014). |
+----------+---------------+-------------+
| i | Hypothetical | Fennell, C. |
| | ice "i". | J. & |
| | | Gezelter, |
| | | J. D. |
| | | Computation |
| | | al |
| | | Free Energy |
| | | Studies of |
| | | a New Ice |
| | | Polymorph |
| | | Which |
| | | Exhibits |
| | | Greater |
| | | Stability |
| | | than Ice I |
| | | h. J. Chem. |
| | | Theory |
| | | Comput. 1, |
| | | 662-667 |
| | | (2005). |
+----------+---------------+-------------+
| C0-II | Filled ice C0 | Smirnov, G. |
| | (Alias of | S. & |
| | 17). | Stegailov, |
| | | V. V. |
| | | Toward |
| | | Determinati |
| | | on |
| | | of the New |
| | | Hydrogen |
| | | Hydrate |
| | | Clathrate |
| | | Structures. |
| | | J Phys Chem |
| | | Lett 4, |
| | | 3560-3564 |
| | | (2013). |
+----------+---------------+-------------+
| C1 | Filled ice C1 |
| | (Alias of |
| | 2d). |
+----------+---------------+-------------+
| C2 | Filled ice C2 |
| | (Alias of |
| | 1c). |
+----------+---------------+-------------+
| sTprime | Filled ice | Smirnov, G. |
| | sT' | S. & |
| | | Stegailov, |
| | | V. V. |
| | | Toward |
| | | Determinati |
| | | on |
| | | of the New |
| | | Hydrogen |
| | | Hydrate |
| | | Clathrate |
| | | Structures. |
| | | J Phys Chem |
| | | Lett 4, |
| | | 3560-3564 |
| | | (2013). |
+----------+---------------+-------------+
| CS1, | Clathrate | Matsumoto, |
| CS2, | hydrates CS1 | M. & |
| TS1, HS1 | (sI), CS2 | Tanaka, H. |
| | (sII), TS1 | On the |
| | (sIII), and | structure |
| | HS1 (sIV). | selectivity |
| | | of |
| | | clathrate |
| | | hydrates. |
| | | J. Phys. |
| | | Chem. B |
| | | 115, |
| | | 8257-8265 |
| | | (2011). |
+----------+---------------+-------------+
| RHO | Hypothetical | Huang, Y et |
| | ice at | al. “A New |
| | negative | Phase |
| | pressure ice | Diagram of |
| | 'sIII'. | Water Under |
| | | Negative |
| | | Pressure: |
| | | the Rise of |
| | | the |
| | | Lowest-Dens |
| | | ity |
| | | Clathrate |
| | | S-III.” |
| | | Science |
| | | Advances |
| | | 2.2 (2016): |
| | | e1501010–e1 |
| | | 501010. |
+----------+---------------+-------------+
Please ask vitroid@gmail.com to add new ice structures. ##Water models
+-------------+-------------------+
| symbol | type |
+=============+===================+
| ``tip3p`` | TIP3P (default) |
+-------------+-------------------+
| ``tip4p`` | TIP4P |
+-------------+-------------------+
| ``tip5p`` | TIP5P |
+-------------+-------------------+
Guest molecules
---------------
+-----------------------------------------+--------------------------+
| symbol | type |
+=========================================+==========================+
| ``co2`` | CO2 |
+-----------------------------------------+--------------------------+
| ``uathf`` | United atom 5-site THF |
+-----------------------------------------+--------------------------+
| ``g12``,\ ``g14``,\ ``g15``,\ ``g16`` | A monatomic dummy site |
+-----------------------------------------+--------------------------+
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
GenIce-0.8.5.tar.gz
(33.6 kB
view details)
Built Distribution
File details
Details for the file GenIce-0.8.5.tar.gz
.
File metadata
- Download URL: GenIce-0.8.5.tar.gz
- Upload date:
- Size: 33.6 kB
- Tags: Source
- Uploaded using Trusted Publishing? No
File hashes
Algorithm | Hash digest | |
---|---|---|
SHA256 | 1be1ee2ebd4119664dc03276edea178506c6a35f7b8f99a426dfc740c30b58a2 |
|
MD5 | 9b71009c5c0d74f86c0f8d32edd7d4b2 |
|
BLAKE2b-256 | 76bfecd49a076855cc9d50e869c7e451f24886872fc655db97cbdb6c1c1fba8f |
File details
Details for the file GenIce-0.8.5-py2.py3-none-any.whl
.
File metadata
- Download URL: GenIce-0.8.5-py2.py3-none-any.whl
- Upload date:
- Size: 64.5 kB
- Tags: Python 2, Python 3
- Uploaded using Trusted Publishing? No
File hashes
Algorithm | Hash digest | |
---|---|---|
SHA256 | f4669f6f0e3493c933cd06d755a8ca561b96ad705e8ba9baba7152409446e6ae |
|
MD5 | 900512b3f40c7589d85cbfd83e4e670b |
|
BLAKE2b-256 | 125bcb841fd23145e8d27d547752c3190133dbfc85a631e8ae7b1720f39cb901 |