2005
DOI: 10.1021/jp050690z
|View full text |Cite
|
Sign up to set email alerts
|

Topological Defects and Bulk Melting of Hexagonal Ice

Abstract: We use classical molecular dynamics combined with the recently developed metadynamics method [Laio, A.; Parrinello, M. Proc. Natl. Acad. Sci. U.S.A. 2002, 99, 20] to study the process of bulk melting in hexagonal ice. Our simulations show that bulk melting is mediated by the formation of topological defects which preserve the coordination of the tetrahedral network. Such defects cluster to form a defective region involving about 50 molecules with a surprisingly long lifetime. The subsequent formation of coordi… Show more

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
1
1
1
1

Citation Types

4
57
0

Year Published

2006
2006
2015
2015

Publication Types

Select...
7
3

Relationship

2
8

Authors

Journals

citations
Cited by 69 publications
(61 citation statements)
references
References 27 publications
4
57
0
Order By: Relevance
“…For simplicity our implementation restricts choices of r 2 such that this cannot occur. The use of the potential energy as a biasing order parameter was first demonstrated by Donadio et al [20] and has since been employed in a study of Lennard-Jones nucleation [6], and by the present authors in the freezing of water [7] and crystallisation of calcium carbonate nanoparticles [8].…”
Section: A Steinhardt Order Parametersmentioning
confidence: 76%
“…For simplicity our implementation restricts choices of r 2 such that this cannot occur. The use of the potential energy as a biasing order parameter was first demonstrated by Donadio et al [20] and has since been employed in a study of Lennard-Jones nucleation [6], and by the present authors in the freezing of water [7] and crystallisation of calcium carbonate nanoparticles [8].…”
Section: A Steinhardt Order Parametersmentioning
confidence: 76%
“…The rings were constructed using the shortest path ring definition (43), which considers a water molecule (using the oxygen atoms as centers) and two of its nearest neighbors with a distance cutoff criterion of 3.3 Å. This cutoff corresponds to the end of the first hydration shell given by the g(r) in BLYP water.…”
Section: Methodsmentioning
confidence: 99%
“…The H-bond estimation is based on the definition described in Ref. 42 . For the presence of a H-bond between two water molecules we calculated the h i j function as a product of two terms:…”
Section: Characterization Of the Watermentioning
confidence: 99%