Ising models of undercooled binary system crystallization; comparison with experimental and pegmatite textures

Simulations of crystal growth in hydrous albite-quartz and albite-orthoclase systems were performed with the Jackson, Gilmer, and Temkin formulation of the Ising model. These simulations demonstrate that at 873 K (approximately 100 K undercooling), comb textures are produced when growth and diffusion probabilities are equivalent. If the growth probability is decreased to 0.1 times that of diffusion, discrete alternating zones of albite and quartz are produced; however, increasing the growth probability of 10 times that of diffusion results in an intergrowth of small domains of albite set in a quartz matrix. At near-liquidus temperatures (973 K), textures similar to those at 100 K undercooling are produced although the crystal-melt interfaces of the comb texture are smoother and the absolute value of the diffusion probability exerts a strong control on the texture, unlike results at 873 K. The textures produced in the simulations with a scale of 10 (super -7) m are similar to experimentally produced textures with scales of 10 (super -6) to 10 (super -4) . The fractal dimension of a comb-textured simulation has been measured and shown to be similar to the fractal dimension of a natural pegmatite with a meter-scale texture. This scaling similarity suggests that these simulations may provide insight into the formation of natural pegmatites despite the 5-7 orders of magnitude difference in scale.

[1]  D. Baker,et al.  NA-K INTERDIFFUSION IN ALKALI FELDSPAR MELTS , 1998 .

[2]  W. Simmons,et al.  The role of Diffusion-controlled Oscillatory Nucleation in the Formation of Line Rock in Pegmatite-plite Dikes , 1997 .

[3]  K. A. Jackson,et al.  Orientation dependence of the distribution coefficient obtained from a spin-1 Ising model , 1997 .

[4]  G. Gilmer,et al.  Microsegregation far from equilibrium , 1996 .

[5]  Temkin,et al.  Monte Carlo simulation of the rapid crystallization of bismuth-doped silicon. , 1995, Physical Review Letters.

[6]  A. Navrotsky Physics and chemistry of earth materials , 1994 .

[7]  David London,et al.  The application of experimental petrology to the genesis and crystallization of granitic pegmatites , 1992 .

[8]  D. Baker Interdiffusion of hydrous dacitic and rhyolitic melts and the efficacy of rhyolite contamination of dacitic enclaves , 1991 .

[9]  C. Ross Ising Models and Geological Applications , 1991 .

[10]  B. Chakoumakos,et al.  Pressure-temperature constraints on the crystallization of the Harding Pegmatite, Taos County, New Mexico , 1990 .

[11]  A. Lasaga,et al.  Crystal-growth kinetics of plagioclase in igneous systems; isothermal H 2 O-saturated experiments and extension of a growth model to complex silicate melts , 1988 .

[12]  A. Navrotsky,et al.  Some Thermodynamic and Experimental Constraints on the Melting of Albite at Atmospheric and High Pressure , 1982, The Journal of Geology.

[13]  D. Eric Crystal growth and nucleation theory and the numerical simulation of igneous crystallization , 1980 .

[14]  G. Lofgren Experimental studies on the dynamic crystallization of silicate melts , 1980 .

[15]  P. M. Fenn The nucleation and growth of alkali feldspars from hydrous melts , 1977 .

[16]  E. G. Ehlers The interpretation of geological phase diagrams , 1972 .

[17]  R. H. Jahns,et al.  Experimental studies of pegmatite genesis; l, A model for the derivation and crystallization of granitic pegmatites , 1969 .

[18]  A.,et al.  Mineral Crystallinity in Igneous Rocks Fractal Method , 2022 .