We use cookies to show content based on your preferences. If you continue to browse you accept their use and installation. More information. >

FMF - Friends of Minerals Forum, discussion and message board
The place to share your mineralogical experiences


Spanish message board






Newest topics and users posts
24 Oct-21:41:08 Re: don lum collection (Don Lum)
24 Oct-21:40:15 Re: the mizunaka collection - quartz (Am Mizunaka)
24 Oct-20:05:46 Silvia's collection - huge siderite (Silvia)
24 Oct-19:57:24 Re: the mizunaka collection - quartz (Silvia)
24 Oct-18:44:45 The mizunaka collection - quartz (Am Mizunaka)
24 Oct-15:15:23 Re: collection of sante celiberti (Sante Celiberti)
24 Oct-14:52:30 Re: collection of michael shaw (Michael Shaw)
24 Oct-13:43:08 Re: collection of michael shaw (Silvia)
24 Oct-11:17:10 Re: collection of michael shaw (Michael Shaw)
24 Oct-09:37:12 Re: collection of volkmar stingl (Volkmar Stingl)
24 Oct-09:28:16 Re: collection of michael shaw (Volkmar Stingl)
24 Oct-09:04:01 Re: collection of michael shaw (Michael Shaw)
24 Oct-03:54:43 Re: collection of tobi (Dany Mabillard)
24 Oct-03:50:24 Re: collection of tobi (Tobi)
23 Oct-23:28:11 Re: help found with items (Philippe Durand)
23 Oct-15:08:16 Silvia's collection - classic portugese ferberite (Silvia)
23 Oct-13:59:45 Re: help found with items (Lou B.)
23 Oct-13:56:06 Help found with items (Lou B.)
23 Oct-13:45:52 Re: collection of tobi (Tobi)
23 Oct-12:29:41 Re: munich show (mineralientage) 2021 (Jordi Fabre)
23 Oct-12:23:20 Munich show (mineralientage) 2021 - like a new show (Jordi Fabre)
23 Oct-11:54:58 Re: collection of tobi (Tobi)
23 Oct-09:31:25 Re: munich show (mineralientage) 2021 (James Catmur)
23 Oct-04:07:52 Re: new generation for picture - ploum (Ploum)
23 Oct-01:16:33 Re: new generation for picture - ploum (Ploum)

For lists of newest topics and postings click here


RSS RSS

View unanswered posts

Why and how to register

Index Index
 FAQFAQ RegisterRegister  Log inLog in
 {Forgotten your password?}Forgotten your password?  

Like
85976


The time now is Oct 24, 2021 22:23

Search for a textSearch for a text   

A general guide for using the Forum with some rules and tips
The information provided within this Forum about localities is only given to allow reference to them. Any visit to any of the localities requires you to obtain full permission and relevant information prior to your visit. FMF is strictly against any illicit activities related to collecting minerals.
Just how fast can the crystals we love and collect grow?!
  
  Index -> Minerals and Mineralogy
Like
7


View previous topic :: View next topic  
Author Message

Peter Lemkin




Joined: 18 Nov 2016
Posts: 357
Location: Prague

View user's profile

Send private message

PostPosted: Sep 29, 2021 11:45    Post subject: Just how fast can the crystals we love and collect grow?!  

I can across this article which covered an area about mineralogy often neglected and very much of interest to me. I used to live in San Diego, and while I never got into the Stewart Mine [where this research article was done experimentally], I did collect often at a location geologically related, but South of the Pala area. This article is long, so I won't try to cut and paste it all here. It has a creative commons license, so I could. I hope other of you who find this research as amazing as I [that minerals like Quartz and other pegmatitic minerals may sometimes grow at a rate of 10m/day as the upper limit! [that is fast!!!]], will follow the url and read the full article with some interesting photos and diagrams. While perhaps the full text is quite technical for some, everyone will get something out of the article.

https://www.nature.com/articles/s41467-020-18806-w
(link normalized by FMF)

Episodes of fast crystal growth in pegmatites

Patrick R. Phelps, Cin-Ty A. Lee & Douglas M. Morton

Nature Communications volume 11, Article number: 4986 (2020)

Abstract

Pegmatites are shallow, coarse-grained magmatic intrusions with crystals occasionally approaching meters in length. Compared to their plutonic hosts, pegmatites are thought to have cooled rapidly, suggesting that these large crystals must have grown fast. Growth rates and conditions, however, remain poorly constrained. Here we investigate quartz crystals and their trace element compositions from miarolitic cavities in the Stewart pegmatite in southern California, USA, to quantify crystal growth rates. Trace element concentrations deviate considerably from equilibrium and are best explained by kinetic effects associated with rapid crystal growth. Kinetic crystal growth theory is used to show that crystals accelerated from an initial growth rate of 10−6–10−7 m s−1 to 10−5–10−4 m s−1 (10-100 mm day−1 to 1–10 m day−1), indicating meter sized crystals could have formed within days, if these rates are sustained throughout pegmatite formation. The rapid growth rates require that quartz crystals grew from thin (micron scale) chemical boundary layers at the fluid-crystal interfaces. A strong advective component is required to sustain such thin boundary layers. Turbulent conditions (high Reynolds number) in these miarolitic cavities are shown to exist during crystallization, suggesting that volatile exsolution, crystallization, and cavity generation occur together.

Introduction

Crystal size is traditionally linked to the cooling rate of a magmatic system1,2,3. In rapidly cooled systems (minutes to days), such as in a lava that erupts and quenches, the resulting rock is characterized by fine (nm–μm) grain sizes. In rocks that cool more slowly (10–100 ky), such as in magmas that intrude or stall deep in the crust to form plutons instead of erupting to the surface, grains are coarser (up to cm size). These patterns seem straightforward: crystals take time to grow, so rapid cooling limits grain size, while slow cooling permits crystals to grow large. Granite is a silicic magma that cooled slowly and generated a coarse-grained rock. That same magma, when cooled quickly, becomes a fine-grained volcanic rock called rhyolite.

There are, however, observations that suggest the above understanding of grain growth in nature is not so simple. For example, it is not uncommon to find igneous rocks with feldspar crystals (e.g., megacrysts) several cm to even 10 cm in size in shallow plutons (e.g., porphyries), which, presumably, cooled faster than their finer-grained, deeper counterparts4. More strikingly, decimeter to even meter-sized crystals typify pegmatitic systems, but the small sizes of pegmatite bodies suggest the systems cooled quickly5,6. Many of the largest crystals on Earth thus appear to be found in seemingly short-lived systems7. Are crystals larger because they have had more time to grow, or do they represent anomalously rapid crystal growth? If the latter, what conditions promote rapid crystal growth, and might these conditions be more widespread than currently thought?

In many of these systems, water may play a crucial role in augmenting grain growth. For granitic pegmatites, water is key for reducing melt viscosities and crystallization temperatures while enhancing melt transport, allowing for much lower temperature crystallization to occur8. Here, we show that late-stage centimeter-scale crystals grew within hours. If such growth rates can be sustained, decimeter to meter-sized crystals in pegmatitic systems could have grown in days, with growth occurring in highly dynamic rather than static conditions. We suggest that the crystals grew from a free fluid phase within a turbulent convective system.
Results
The Stewart pegmatite

To explore crystal growth rates in natural systems, we conducted a case study of quartz crystals from the Stewart pegmatite in the Pala district of southern California, USA. The pegmatites in this district are Cretaceous-aged and occur as dikes hosted in gabbroic plutons of the same age9,10. The pegmatites are thought to represent late-stage fluids fractionated from a hydrous gabbroic pluton that then form a tabular, compositionally zoned pegmatite body9. At its widest, the Stewart body is represented by a 50-m-wide dike with the outer layers composed of feldspar-rich granite (the upper and lower intermediate zones), a transitional perthite zone, and a coarse-grained internal core containing spodumene, lepidolite, heulandite, petalite, and albite (Fig. 1). Vertical structures (chimneys), which contain numerous miarolitic cavities, emanate from the internal core, penetrating upwards into the perthite zone. The chimneys host abundant course-grained tourmalines and albite, with tourmaline crystals in some cases approaching tens of centimeters in length. In these chimneys, tourmalines are set within a finer-grained matrix of albite, quartz, and lepidolite, with the base of the chimneys being lepidolite-rich. Quartz crystals in the chimneys, particularly within the miarolitic cavities, are often euhedral. The contacts between the chimneys and the pegmatite body are sharp, suggesting the chimneys represent late-stage fluids that diked into a mostly solidified upper zone of the pegmatite body. These chimneys make up ~50% by volume of the perthite zone and thus represent an important phase in the petrogenetic evolution of the pegmatite9. Based on conductive cooling models, Webber et al.6 suggest that the pegmatitic body cooled from 650 °C to below 550 °C within ~9 years after emplacement. This is an upper bound on cooling timescales (and lower bound on crystal growth rates) because any convective cooling, which was not considered by Webber et al.6, would decrease estimated cooling timescales. Crystal growth rates inferred from pegmatite thermal lifespans would also be under-estimated because crystal growth occurs at different stages within the pegmatite’s thermal evolution.
Fig. 1: Cartoon of the Stewart pegmatite with an image of a miarolitic cavity.
figure1

a Schematic diagram of the Stewart pegmatite modified from Morton et al.9. A 2 m tall person for scale on the lower left. Outer and intermediate zones are granitic in composition. b An example of a small miarolitic cavity (chimney) that propagated from the core zone into the solidified upper zone. The miarolitic cavity contains tourmaline and quartz crystals set within a fine-grained matrix of lepidolite. Lighter margins of the cavity are defined by white albite.

Full size image

We investigated euhedral quartz crystals from the miarolitic cavities to constrain crystallization rates. The quartzes in these cavities are transparent, euhedral (doubly terminated), and range in size from a few cm to 10 cm. Their euhedral shapes indicate that they nucleated and crystallized from within the fluid making up the miarolitic cavities, not by nucleating along the walls. They are peculiar in that they have flattened crystal habits with the <1000> axis being short (the angled braces referring to the family of axes). Tourmaline appears to penetrate some quartz crystals perhaps as a result of quartz overgrowths.

Quartz composition

Quartz crystals were imaged at Rice University using cold cathodoluminescence (CL) microscopy, a type of visible wavelength luminescence induced by high-energy electrons emitted from a cold cathode11. Defects in the crystal lattice cause luminescence in some minerals. Since defects are often correlated with trace element substitutions, changes in CL color may indirectly indicate a change in trace element concentration12. Thus, changes in CL intensity were used to guide our geochemical study. The crystals were mounted in epoxy, cut, and polished. For CL imaging, an accelerating voltage of 12 kV was used with a vacuum current of 0.4–0.5 mA. Camera exposure was set to 4 s to maximize image quality. The resulting image for one of the crystals is shown in Fig. 2a (see Supplementary Fig. 1 for another CL quartz crystal).
Fig. 2: Cathodoluminescence image of a quartz crystal and corresponding trace element transect.
figure2

a Cathodoluminescence (CL) image of a quartz crystal from the Stewart pegmatite. Yellow line represents trace of geochemical transect. The plane of the crystal is part of {1000}. b Geochemical transect A–A′ for Ge, Li, Al, and Ti in p.p.m. by weight. Dots represent data and red lines represent moving averages. Red envelopes represent standard error of the mean based on 3σ of point analyses and a three-point moving average. Gray envelope in Ti represents 1 standard deviation to the data points. The error values are 26%, 28%, 8.5%, and 57% for Ge, Li, Al, and Ti, respectively. Background colors correspond with CL colors in a. Ti is scattered due to measurements being near detection limits.
Full size image

Trace element compositions were quantified using a 213 nm New Wave laser ablation system coupled to a ThermoFinnigan Element 2 magnetic sector inductively coupled plasma mass spectrometry (ICP-MS) at Rice University. Laser spot analyses and line transects were made. Transects were performed using a 110 μm spot diameter with a ~13 J m−2 fluence, a 5 Hz repetition rate, and a lateral laser velocity of 9 μm s−1. The large spot size limits the spatial resolution of our transects to 100–200 μm. Analyses were done in medium mass resolution mode (m/Δm = 3000). The following masses were measured: Li7, Na23, Mg25, Al27, Si30, P31, K39, Ti49, Zn66, Ga69, Ge73, and Ge74. Signals were background-corrected, normalized to Si30, and converted to concentrations using synthetic glass NIST 61213 and two natural rhyolitic obsidians M3–33 and M3–8614 as external standards and forcing the sum of all measured metals (in oxide form) to equal 100%.

Three quartz crystals were mapped using CL. The CL of one crystal is shown in Fig. 2a and is representative of the color zones seen in CL from the pegmatite. In CL, the quartzes consist of three sharply demarcated zones: a white core zone, an orange middle zone, and a purple outer rim. The white core consists of two subtle internal zones: a dimmer inner zone surrounded by a thin, slightly brighter mantle. We consider these two zones to be internal parts of the core zone because their CL colors are only subtly different. Throughout the grain, transient blue luminescence mainly in the orange and purple zones, which may be related to electron beam-induced rearrangements of alkali trace elements at the nanoscale15,16, was observed.

A representative compositional transect (A–A′) across all three CL zones is shown in Fig. 2b (see Supplementary Figs. 2–7 for other transects). Transitions in CL color or intensity correlate with broad changes in trace element concentrations12. In the white core region, most elements show minimal variation with distance (no compositional change occurs across the internal CL transition within the white core zone discussed above). For example, in the white core zone, Li and Al are relatively constant, while Ge increases outwards from 10 to ~20 p.p.m. before decreasing back to 10 p.p.m. At the contact between the white core and orange middle zone, marked decreases in Ge and increases in Al and Li are observed. In the orange middle zone, Ge, Al, and Li show gradual outward increase in concentrations, punctuated by spikes in Al and Li due to inclusions (confirmed optically but of unknown mineralogy due to their small size). At the contact between the orange middle and purple outer zone, Al and Li increase markedly, but no sharp changes are seen in Ge. In the purple outer zone, Al and Li flatten quickly and even decrease slightly, while Ge continues to increase.

Ti concentrations within the core zone of the crystal are near detection limits and are to within error constant between 10 and 15 p.p.m. We can use these Ti contents to obtain a maximum bound on temperature using Ti in quartz thermometry. Assuming TiO2 activity of 0.6 (no Ti-bearing phases) and 2 kbar, we arrive at ~600 °C using the thermometer of Wark and Watson17. The thermometers of Huang and Audétat18 and Thomas et al.19 allow for pressure dependence, lowering the temperature to between 46019 and 540 °C18. Adopting a lower TiO2 activity would yield higher temperatures, but if TiO2 contents are high due to rapid crystal growth (see below), lower temperatures would be estimated. While we cannot derive an exact temperature for the formation of these quartzes, it would seem that the quartzes formed at temperatures no >540 °C.
Kinetic versus external control on compositional zoning

We interpret the CL transitions as reflecting original differences in the density of defects or trace element impurities associated with sudden changes in the crystal growth rate. The broader Li, Ge, and Al concentration profiles across the sharp CL zones are larger than the spatial resolution of our laser spot size and may suggest diffusive relaxation of these elements. However, the width of each of these concentration profiles is similar, which is reflected in the strong coupling of Li and Al concentrations (Fig. 3). The similarity in the length scales of these transitions, if related to diffusion, is unexpected because these elements have markedly different diffusivities. For example, at temperatures of ~600 °C, tracer diffusivities in quartz are 5 × 10−12 m2 s−1 (following Verhoogen20) for Li, 3 × 10−23 m2 s−1 for Al (following Tailby et al.21), and ~1 × 10−32 m2 s−1 for Ge (assuming Ge behaves similarly to Si self-diffusion22). In particular, the extremely low diffusivity of Al would require nearly 2 My to generate the Al profiles by diffusion assuming an initial step function profile—far longer than the lifespan of the pegmatite. No diffusive relaxation is thus predicted for Ge and Al.
Fig. 3: Al versus Li in quartz.
figure3

Line represents 1:1 correspondence of Li and Al. Correlation of Li with Al indicates possible charge coupling. Other monovalent cations, such as Na, may account for the remaining Al. Error bars are 3σ based on internal measurement precision.
Full size image

The presumably high diffusivity of Li, however, should predict measurable diffusive relaxation. For example, over the lifespan of the pegmatite (~9 years), Li should have diffused over a length scale of 7.5 cm—larger than the crystal sizes in this study. However, the fact that Li correlates strongly with Al suggests that Li+ in quartz is charge-coupled with Al3+, suggesting that published Li tracer diffusivities may not be applicable. Instead, Li might only diffuse as fast as Al due to charge coupling. The possibility that Li diffuses very slowly is at odds with some studies that have inferred short residence times of crystals in magmatic systems from Li concentration profiles23. Our findings are similar to the observation of very slow diffusion of Li in zircons due to coupling with slow-diffusing rare-earth elements24. The strong coupling of Li to slowly diffusing Al and Ge inferred here suggests that there may have been limited subsolidus diffusive relaxation of Li, Al, and Ge across these boundaries.

The origin of compositional differences between and within each CL zone are thus of interest. Possibilities include changes in fluid composition (external forcing) or changes in the nature of crystal growth (internal processes). For the following reasons, we do not favor changes in fluid composition as the primary driver. Although quartz Ti concentrations (varying from 10 to >50 p.p.m., depending on location within the crystals) in this study are similar to those in quartz from granitic plutons, Al concentrations (>1000 p.p.m.) are significantly higher than those in quartz from granitic plutons and most pegmatites (Fig. 4). Given that Al is highly incompatible in quartz25, these remarkably high Al concentrations would require fluids or melts with unreasonably high Al contents. We also note that the trace element composition of magmas or fluids undergoing differentiation should evolve in predictable ways: progressive crystallization leads to enrichment of incompatible elements and depletion of compatible elements. Li and Al are incompatible in quartz25,26. Ge is similar to Si27,28,29 and should thus be compatible in quartz, as confirmed by He et al.30. Li and Al should thus be systematically anticorrelated with Ge, if their concentrations are controlled by a fluid or magma undergoing crystal fractionation. Instead, the relationship of Li and Al with Ge across the entire quartz grain shows more complexity.
Fig. 4: Ti versus Al concentrations in quartz from granites and pegmatites.
figure4

Granitic data compiled from Ackerson et al.60 and pegmatite data from Garate-Olave et al.61 and Götze et al.15. Error bars represent 3σ. Data for this study are from point analyses.
Full size image

While the possibility of external changes in fluid composition cannot be completely excluded, adopting such a hypothesis would not only require highly unusual fluid compositions, but very large and sudden changes in fluid composition. For these reasons, we also explore the possibility that the trace element compositions reflect kinetic processes, specifically the effects of rapid crystal growth.
Discussion

In this section, we explore how changes in crystal growth rate can affect how trace elements are incorporated into crystals31,32,33. A trace element will partition into a mineral according to its equilibrium partition coefficient (k), which is the concentration in the crystal divided by the concentration in the growth medium (fluid or liquid) immediately adjacent to the growing crystal surface at equilibrium (k = Cs/Cl). However, if the crystal grows at rates comparable to or faster than the element can diffuse in the fluid, a chemical boundary layer in the fluid develops near the advancing crystal-fluid interface (Fig. 5). For an incompatible element (k <8201> 1), wherein preferential incorporation of the element into the crystal leads to a depleted chemical boundary layer in the fluid. Qualitatively, the width of the boundary layer is controlled by the relative rates of crystal growth and diffusive or advective homogenization in the fluid. For example, fast crystal growth compared to diffusion in the fluid will result in extremely enriched and narrow chemical boundary layers for incompatible elements (and extremely depleted and wide boundary layers for compatibles). If diffusive timescales in the crystal are longer than crystal growth timescales, crystal zoning will reflect changes in the composition of the boundary layer, not the composition of the far-field fluid or magmatic growth medium.
Fig. 5: Trace element response due to changes in growth rate.
figure5

Shown is a conceptual diagram of how crystal growth rate (R) affects trace element distributions associated with a sudden change in growth rate. The figure shows the transient case in which growth rate suddenly increases from a previous steady state. The concentration (C) in the fluid near the crystal rises or falls, depending on the element’s compatibility (k, the equilibrium ratio of the solid concentration to the liquid concentration), which is recorded in the crystal as an enrichment or depletion in the crystal, respectively. If the growth rate stays constant, the fluid concentration will eventually reach a new steady state, which will be reflected in a constant C in the crystal as a function of distance, as distance represents time. The concentration profile in the fluid at steady state will represent a static chemical boundary layer.
Full size image

These concepts can be placed in a quantitative framework following an advection-diffusion model with the crystal–liquid interface as a moving reference frame31,33,34,35
∂𝐶l∂𝑡=𝐷∂2𝐶l∂𝑥2+𝑅∂𝐶l∂𝑥,
(1)

where Cl is the concentration of the trace element in the liquid, D is the diffusivity of that element in the liquid, R is the growth rate, t is time, and x is the distance from the growth interface in the liquid. The boundary conditions are:
𝐷∂𝐶l∂𝑥∣∣𝑥=0=𝑅(𝑘−1)𝐶l|𝑥=0at𝑥=0,
(2a)
𝐶l=𝐶eqas𝑥→∞,
(2b)

where Ceq is the concentration in the far-field fluid. The interface boundary condition (2a) assumes flux continuity (Ch. 935). The diffusive flux of an element out of the fluid and into the crystal is governed by Fick’s law, 𝐽diff=−𝐷∂𝐶l∂𝑥

. This flux is balanced by the growth flux, which is given by Jrej = (Cl − Cs)R or Jrej = (1 − k)Cl R. The boundary condition at infinity (2b) assumes that the fluid composition is constant at distances greater than the diffusive length scale within the fluid. A fundamental assumption in this approach is that reaction kinetics at the crystal surface are fast compared to diffusion and can be neglected. If this is not the case, a modified partition coefficient, like that detailed in ref. 36, must be included. Because there are no constraints on rate dependent k for our system, we assume k to be constant from here on.

By setting the time derivative to zero in Eq. 1, we can explore the steady-state solution for the composition in the liquid31:
𝐶l=𝐶eq(1+1−𝑘𝑘exp(−𝑅𝑥𝐷)).
(3)

Equation 3 gives the form of the steady-state concentration profile in the fluid relative to the far-field fluid composition Ceq. The maximum deviation from the far-field concentration Ceq occurs at the interface and scales with 1/k, as seen by setting x = 0 in Eq. 3. The deviation (positive for k <8201> 1) of the liquid composition from equilibrium decreases exponentially with distance away from the growing crystal interface.

Importantly, at steady state, the concentration in the solid Cs remains constant as the crystal grows and is equal to the product of the partition coefficient k and the concentration in the liquid at x = 0, that is, k × (Ceq/k) = Ceq. This holds for a static growth medium. If the growth fluid flows, then a different steady state must be considered.

The concentration in the crystal at the onset of growth is Cs = kCeq and gradually evolves to the steady solution given by Burton et al.32,
𝐶s=𝐶eq(𝑘𝑘+(1−𝑘)exp(−𝑅δD)),
(4)

where δ is the steady-state width of the chemical boundary layer in the fluid, and Cs is the composition of the crystal at the interface. This is found by restating Eq. 2b with x → δ. The quantity Rδ/D is a Peclet number representing the relative strengths of advection versus diffusion. If the Peclet number is small, the composition of the crystal approaches equilibrium with the far-field fluid, that is, Cs = kCeq. For very large Peclet numbers, such as for large growth rates, the crystal composition Cs approaches the far-field liquid composition Ceq35. In a static fluid, the steady-state width of the chemical boundary layer δ scales with D/R and represents the distance over which the deviation from equilibrium with the far-field fluid decays by e-fold, as shown by Eq. 3. Thus, for constant diffusivity D in the fluid, an increase in the crystal growth rate R will translate to a smaller δ. If, however, the fluid itself is advecting, the chemical boundary layer can become physically thinned, such that δ becomes a complex function of D and fluid flow rate. A steady state can still exist, but δ is thinner due to advective erosion of the boundary layer in the fluid33.

The time-dependent solution of Eq. 1 in a static fluid can be derived for the case in which the system is initially at steady state under a constant initial growth rate Ri (Eq. 3), but then experiences an instantaneous change to a new and constant growth rate Rf34. How the system returns to steady state after a step change in R is given by Eq. 5 (see “Methods”). Qualitative solutions are shown in Fig. 5 to give examples of the behavior of this solution. The transient solution reveals that upon a sudden increase in growth rate, incompatible elements in the growing crystal will first increase before falling back to a steady state, whereas compatible elements first decrease before rising back to steady state. This behavior relates to the time required for δ to respond, translating to the peaks in concentration for the incompatible case and valleys for the compatible one. The primary controls on the behavior of the response are initial growth rate Ri, which controls initial steady-state concentration and hence the magnitude of the response, and the ratio of Ri/Rf and diffusivity, which together control the shape of the response after an increase in growth rate.

We first explore compositional variations in the quartz in the context of rapid but steady-state crystal growth. Exploring steady state gives insight into the direction systems should evolve towards. Equation 4 can be evaluated for different Peclet numbers. Doing so requires constraints on the mineral/fluid partition coefficient k, the diffusivity D of the element in the fluid, and the concentration of the far-field fluid Ceq. To better constrain results, we use two elements: Al and Ge. For Al, the mineral/fluid partition coefficient has not been measured, so we use published mineral/melt partition coefficients of ~0.0225. For Ge, following the approach of Hofmann et al.37, we can use the fact that Ge/SiO2 ratio changes only slightly over a range of SiO2 concentrations from andesitic to rhyolitic compositions (Supplementary Fig. 8) to infer that Ge partitioning into quartz is similar to that of Si30. Therefore, kGe can be estimated by constraining how Si partitions into quartz relative to a fluid, which can be constrained from Si solubility in the fluid. Because the quartz precipitated in miarolitic cavities, a water-dominated fluid is assumed rather than a hydrous silicate melt as the growth medium8,9,38,39. At 540 °C and 150–200 MPa6,10 (5–7 km emplacement depth) and using the Si solubility data of Fournier and Potter40 and Burnham et al.41, we arrive at a kSi of ~200. We assume that kGe for quartz is of similar magnitude. Recent work on natural quartz–silicate melt partitioning by He et al.30 confirms that Ge is compatible in quartz and only slightly less compatible than Si.

Diffusivity of Al and Ge in the fluid is determined following methods outlined in Nigrini42 and Yuan-Hui and Gregory43. Their parameterizations are applicable only to +1 to +3 valence states. The likely valence of Ge is +4, but because Ge4+ has a similar cation size to Al3+ of 0.5 pm27, we approximate Ge diffusivity with that of Al3+ for the above temperatures of interest and allow for an uncertainty of two orders of magnitude (~7 × 10−9 ± 1 m2 s−1). This uncertainty in diffusivity more than accounts for the uncertainty in temperature. Changes by ±100 °C do not significantly affect the diffusivities (<15% changes).

Finally, Ceq can be estimated by dividing the concentration at the center of the crystal (~15 p.p.m. for Ge and ~100 p.p.m. for Al) by k as inferred above. This approach is justified because the trace element concentration of a crystal at the onset of crystallization (core) will initially be in equilibrium with the far-field fluid, Cs = kCeq, because no boundary layer in the fluid should have yet developed.

In Fig. 6, we compare Ge and Al concentrations to modeled steady-state values of Ge and Al for different values of Rδ (Eq. 4), that is, the product of crystal growth rate and boundary layer thickness. We have assumed constant diffusivity so that variations in Rδ effectively represent changes in Peclet number. At steady state, δ is expected to be a constant for a given crystal growth rate and diffusivity (and fluid advection), but because we are considering changes in R, we do not a priori know what δ is, so Rδ is varied. In any case, increases in Rδ most likely reflect an increase in growth rate.
Fig. 6: Comparison of measured profiles with modeled steady-state concentrations.
figure6

Steady-state concentrations calculated using Eq. 4 for different Rδ and a constant D of 8 × 10−9 m2 s−1 for both Ge and Al. Data from both Ge and Al towards the core of the crystal correspond to similar Rδ, suggesting steady-state growth in the core. Beyond the core region, Ge and Al give different Rδ, suggesting a transient state. The lower Rδ values (10−10 to 10−12 m2 s−1) in the Al plot are nearly identical, causing them to plot on top of one another. Horizontal color bar corresponds to CL color regions. Error envelopes represent standard error of the mean based on 3σ of point analyses and a three-point moving average. The error values are 26% and 8.5% for Ge and Al, respectively.
Full size image

In the core region, Ge contents are high and Al contents are low, as expected because Ge is compatible and Al incompatible. We find that Ge and Al concentrations in the core can both be modeled by Rδ of 10−11 to 10−12 m2 s−1. In contrast, in the outer part of the crystal, both Ge and Al increase. If steady-state solutions are assumed, increasing Rδ is needed to match the Al increase, while decreases in Rδ are needed to match the Ge data.

In summary, our observations suggest that the crystal was initially growing at steady-state conditions (core white region), but then experienced a sudden change in growth rate (at the white to orange transition). Thus, while the core elemental concentrations reflect steady crystal growth, subsequent elemental profiles reflect transient responses to sudden increases in the crystal growth.

As noted above, Ge and Al profiles in the outer part of the crystal are inconsistent with steady-state growth, and therefore, transient conditions must be considered. Using a given element, the ratio of initial to final growth rate can be constrained but not absolute growth rates. However, for two or more elements, growth rates can be constrained with Eq. 5 (with laser smearing accounted for in Supplementary note 1, Supplementary Eqs. 1 and 2, and Supplementary Fig. 9), which describes the response to an instantaneous increase in growth rate from a crystal initially growing at steady state. We apply this equation to the transition from the white core zone to the orange middle zone because the constant profiles seen in the white core region suggest that steady state was attained as discussed above.

The two parameters we vary are the initial steady-state growth rate Ri and the final growth rate Rf. We adopt the same diffusivities discussed above. Ri and our assumed diffusivity control the steady-state concentration in the white core zone, while the ratio Rf/Ri controls the shape of the response after an increase in growth rate. By varying initial growth rate between 10−7 and 10−6 m s−1 and final growth between 10−5 and 10−4 m s−1, the behaviors of both Al and Ge can be captured (Fig. 7). Adopting slower growth rates (Ri = 10−7 m s−1, Rf = 10−6 m s−1) results in too slow of a response and thus will not fit the curvature of the data within the confines of the model and its assumptions. Adopting faster rates (Ri = 10−5 m s−1, Rf = 10−4 m s−1) also does not fit the data well, showing too short response times. Uncertainties in diffusivities can be accounted for by adopting higher or lower (by an order of magnitude) diffusivities, resulting in higher and lower growth rates, respectively. If the partition coefficient for Al is lower (more incompatible) than used here, we predict nearly the same trace element profile. A more compatible partition coefficient (i.e., 0.1) will produce a slightly lower concentration profile, but still within the array of the data.
Fig. 7: Comparison of measured profiles with modeled transient responses.
figure7

Solutions for transient responses for Ge and Al. Curves represent the response for a step function change in the crystal growth rate from Ri to Rf at ~1000 μm (the orange-white transition). A 2–2.5 order of magnitude change in growth rate is needed to match the data. Al should eventually reach a peak and fall back to steady state, yet the crystal did not grow long enough for this to occur. Color bar on top corresponds with CL color zones in crystal. Error envelopes represent standard error of the mean based on 3σ of point analyses and a three-point moving average. The error values are 26% and 8.5% for Ge and Al, respectively.
Full size image

Based on the above analysis, we estimate that growth rate increased by ~100 times from the core to the rim, with initial growth rates on the order of 10–100 mm day−1 (10−7–10−6 m s−1) in the white core region to growth rates approaching 1–10 m day−1 (10−5–10−4 m s−1) in the orange middle region (Fig. 8). We are unable to estimate the growth rate of the outer purple zone because the orange middle zone never reached steady state, invalidating the application of Eq. 5 to the orange-purple transition. However, the slight decrease in Al (incompatible) and increase in Ge (compatible) in the purple outer zone suggests a slight decrease in growth rate.
Fig. 8: Cartoon depicting how crystal growth rates may have changed from core to rim.
figure8

a A sketch of the CL zonation of the crystal. b Reconstructed growth history of the crystal with overlain CL colors.
Full size image

Based on these estimates of crystal growth rate, we find that the white core region grew in <3 h, while the orange middle zone grew in <4> ~2300 in pipe flow is considered turbulent57, our quartz crystals may have grown in a highly dynamic environment. For comparable conditions, but in silicate melt systems, Re for crystal growth is «1 due to the very high viscosities of silicate melts.

Pegmatites appear to form on much shorter timescales than typical plutonic rocks, but paradoxically, they have much larger crystals, which alone requires that crystal growth rates must be much higher in pegmatitic systems compared to typical plutonic systems. Exactly why is unclear, but it is widely thought that an abundance of a free fluid phase is important for the formation of pegmatites8,38,39. The reasoning is that diffusivities in fluids are higher than that in silicate melts, thereby permitting more rapid flux of nutrients to the crystal and enhancing growth rates58,59. However, an increase in diffusivity alone may not be enough to decrease growth times from thousands of years to hours as diffusivities of common cations in water at elevated temperature are on the order of 10−8 m2 s−142,43 and diffusivities in silicate melts are on the order of 10−10 m2 s−1. As discussed above, flux to the crystal can also be increased by decreasing the thickness of the chemical boundary layer in the fluid via advection. Our work shows that the boundary layer thickness of these crystals was ~104 times thinner than a purely diffusive boundary layer with no advective thinning, ultimately resulting in a 104-fold enhancement of growth rates. Such thin boundary layers, as concluded above, are maintained only under turbulent conditions.

These findings may shed light on the nature of mineralization in the Stewart pegmatite. We envision the pegmatite originating by dike emplacement of a hydrous silicate magma into the host gabbro rock. The outer zones (upper and lower) of the pegmatite body solidified first by cooling against already solidified host gabbro, while the central core zone solidified last as the crystallization front progressed inwards. This progressive crystallization concentrates volatiles and other incompatible elements/components in the residual silicate melt in the core zone. Eventual saturation in these volatiles leads to exsolution of a free fluid phase and a buildup of overpressure, triggering upwards fracturing and transport of the fluid phase into the solid upper zone and generating chimneys and miarolitic cavities.

As the fluid expands into the cavities, the resulting decompression leads to supersaturation of dissolved constituents in the fluid, driving rapid nucleation and crystal growth. Changes in pressure may also cause rapid changes in growth rates. Emplacement of these fluids into fractures appears be accompanied by turbulence, which significantly enhances crystal growth rates, allowing for pegmatitic crystals to form within hours. We thus envision miarolitic cavities to have formed under highly dynamic and ephemeral conditions.

Our crystal growth rates clearly only reflect processes operating in the formation of late-stage miarolitic cavities and do not apply to crystallization in other parts of the pegmatite. Similar studies using trace element chemistry of crystals in different parts of a pegmatite body should be measured in the future to better understand how crystal growth rates vary over the entire history of a pegmatite. In summary, our results suggest that it may be worth reconsidering whether large crystals in other magmatic systems, such as porphyries, are a consequence of rapid growth rather than by slow growth over long incubation times.
Methods

The transient response to a sudden change in the crystal growth after growing at steady state is described by Smith et al.34:
𝐶s=𝐶eq⎛⎝⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜1−12erfc(𝑅f𝑥′4D‾‾‾‾√)+(1−𝑘)(12−𝑅i𝑅f𝑘−𝑅i𝑅f)exp(−𝑅i𝑅f(1−𝑅i𝑅f)𝑅f𝐷𝑥′)erfc((𝑅i𝑅f−12)𝑅f𝑥′𝐷‾‾‾‾√)+2𝑘−12(1−𝑅i𝑅f𝑘−𝑅i𝑅f)exp(−𝑘(1−𝑘)𝑅f𝐷𝑥′)erfc((𝑘−12)𝑅f𝑥′𝐷‾‾‾‾√)⎞⎠⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟,
(5)

wherein time has been transformed to distance x′ within the crystal as that is the quantity that can be measured (x′ = Rf × t). Meanings of other symbols are outlined earlier in the text. This equation has three terms on the right-hand side (we consider 1−12erfc(𝑅f𝑥′4𝐷‾‾‾‾√)

to be one term). The first term is most important at long timescales as the other two terms decay exponentially with time. This first term describes the long-term diffusional response and only depends on the final growth rate Rf because on long timescales, the initial condition does not matter. The second and third terms are important on short timescales and are thus sensitive to the initial condition, which is controlled by Ri, D, and k. The second and third terms depend on Ri and Rf, but the first term only depends on Rf, which allows determination of crystal growth rates provided we have constraints from two elements, such as Ge and Al.
Data availability

All data generated in this study can be found in this published article or in Supplementary Data 1.
Code availability

All models were calculated using off the shelf software (MATLAB) and using the equations presented in this article and its Supplementary information.
References

1.

Winkler, H. G. F. Crystallization of basaltic magma as recorded by variation of crystal-size in dikes. Mineral. Mag. J. Mineral. Soc. 28, 557–574 (1949).

ADS

CAS

Article

Google Scholar

2.

Brandeis, G., Jaupart, C. & Allčgre, C. J. Nucleation, crystal growth and the thermal regime of cooling magmas. J. Geophys. Res. Solid Earth 89, 10161–10177 (1984).

CAS

Article

Google Scholar

3.

Cashman, K. V. Relationship between plagioclase crystallization and cooling rate in basaltic melts. Contrib. Mineral. Petrol. 113, 126–142 (1993).

ADS

CAS

Article

Google Scholar

4.

Vernon, R. H. & Paterson, S. R. How late are K-feldspar megacrysts in granites? Lithos 104, 327–336 (2008).

ADS

CAS

Article

Google Scholar

5.

Chakoumakos, B. C. & Lumpkin, G. R. Pressure–temperature constraints on the crystallization of the harding pegmatite, Taos County, New Mexico. Can. Mineral. 28, 287–298 (1990).

CAS

Google Scholar

6.

Webber, K. L., Simmons, W. B., Falster, A. U. & Foord, E. E. Cooling rates and crystallization dynamics of shallow level pegmatite-aplite dikes, San Diego County, California. Am. Mineral. 84, 708–717 (1999).

ADS

CAS

Article

Google Scholar

7.

London, D. & Morgan, G. B. The pegmatite puzzle. Elements 8, 263–268 (2012).

CAS

Article

Google Scholar

8.

Nabelek, P. I., Whittington, A. G. & Sirbescu, M.-L. C. The role of H2O in rapid emplacement and crystallization of granite pegmatites: resolving the paradox of large crystals in highly undercooled melts. Contrib. Mineral. Petrol. 160, 313–325 (2010).

ADS

CAS

Article

Google Scholar

9.

Morton, D. M., Sheppard, J. B., Miller, F. K. & Lee, C.-T. A. Petrogenesis of the cogenetic Stewart pegmatite-aplite, Pala, California: Regional implications. Lithosphere 11, 91–128 (2018).

ADS

Article

Google Scholar

10.

Jahns, R. H. & Wright, L. A. Gem-and Lithium-Bearing Pegmatites of the Pala District, San Diego County, California (State of California, Department of Natural Resources, Division of Mines, 1951).
11.

Marshall, D. J. & Mariano, A. N. Cathodoluminescence of Geological Materials (Taylor & Francis, 1988).
12.

Götze, J., Plötze, M., Graupner, T., Hallbauer, D. K. & Bray, C. J. Trace element incorporation into quartz: a combined study by ICP-MS, electron spin resonance, cathodoluminescence, capillary ion analysis, and gas chromatography. Geochim. Cosmochim. Acta 68, 3741–3759 (2004).

ADS

Article

CAS

Google Scholar

13.

Pearce, N. J. G. et al. A compilation of new and published major and trace element data for NIST SRM 610 and NIST SRM 612 glass reference materials. Geostand. Newsl. 21, 115–144 (1997).

CAS

Article

Google Scholar

14.

Savov, I. P., Leeman, W. P., Lee, C.-T. A. & Shirey, S. B. Boron isotopic variations in NW USA rhyolites: Yellowstone, Snake River Plain, Eastern Oregon. J. Volcanol. Geotherm. Res. 188, 162–172 (2009).

ADS

CAS

Article

Google Scholar

15.

Götze, J., Plötze, M. & Trautmann, T. Structure and luminescence characteristics of quartz from pegmatites. Am. Mineral. 90, 13–21 (2005).

ADS

Article

CAS

Google Scholar

16.

Perny, B., Eberhardt, P., Ramseyer, K., Mullis, J. & Pankrath, R. Microdistribution of Al, Li, and Na in $α$ quartz: Possible causes and correlation with short-lived cathodoluminescence. Am. Mineral. 77, 534–544 (1992).

CAS

Google Scholar

17.

Wark, D. A. & Watson, E. B. TitaniQ: a titanium-in-quartz geothermometer. Contrib. Mineral. Petrol. 152, 743–754 (2006).

ADS

CAS

Article

Google Scholar

18.

Huang, R. & Audétat, A. The titanium-in-quartz (TitaniQ) thermobarometer: a critical examination and re-calibration. Geochim. Cosmochim. Acta 84, 75–89 (2012).

ADS

CAS

Article

Google Scholar

19.

Thomas, J. B. et al. TitaniQ under pressure: the effect of pressure and temperature on the solubility of Ti in quartz. Contrib. Mineral. Petrol. 160, 743–759 (2010).

ADS

CAS

Article

Google Scholar

20.

Verhoogen, J. Ionic diffusion and electrical conductivity in quartz. Am. Mineral. J. Earth Planet. Mater. 37, 637–655 (1952).

CAS

Google Scholar

21.

Tailby, N. D., Cherniak, D. J. & Watson, E. B. Al diffusion in quartz. Am. Miner. 103, 839–847 (2018).

ADS

Article

Google Scholar

22.

Cherniak, D. J. Silicon self-diffusion in single-crystal natural quartz and feldspar. Earth Planet. Sci. Lett. 214, 655–668 (2003).

ADS

CAS

Article

Google Scholar

23.

Rubin, A. E. et al. Rapid cooling and cold storage in a silicic magma reservoir recorded in individual crystals. Science 356, 1154–1156 (2017).

CAS

PubMed

Article

PubMed Central

Google Scholar

24.

Tang, M., Rudnick, R. L., McDonough, W. F., Bose, M. & Goreva, Y. Multi-mode Li diffusion in natural zircons: evidence for diffusion in the presence of step-function concentration boundaries. Earth Planet. Sci. Lett. 474, 110–119 (2017).

ADS

CAS

Article

Google Scholar

25.

Nash, W. P. & Crecraft, H. R. Partition coefficients for trace elements in silicic magmas. Geochim. Cosmochim. Acta 49, 2309–2322 (1985).

ADS

CAS

Article

Google Scholar

26.

Neukampf, J., Ellis, B. S., Magna, T., Laurent, O. & Bachmann, O. Partitioning and isotopic fractionation of lithium in mineral phases of hot, dry rhyolites: the case of the Mesa Falls Tuff, Yellowstone. Chem. Geol. 506, 175–186 (2019).

ADS

CAS

Article

Google Scholar

27.

Shannon, R. D. & Prewitt, C. T. Revised values of effective ionic radii. Acta Crystallogr. Sect. B 26, 1046–1048 (1970).

CAS

Article

Google Scholar

28.

Wells, A. F. Structural Inorganic Chemistry (Oxford Univ. Press, 2012).
29.

Drewitt, J. W. E. et al. Structure of GeO2 glass at pressures up to 8.6 GPa. Phys. Rev. B 81, 14202 (2010).

ADS

Article

CAS

Google Scholar

30.

He, D., Lee, C.-T. A., Yu, X. & Farner, M. Ge/Si partitioning in igneous systems: constraints from laser ablation ICP-MS measurements on natural samples. Geochem. Geophys. Geosystems 20, 4472–4486 (2019).
31.

Tiller, W. A., Jackson, K. A., Rutter, J. W. & Chalmers, B. The redistribution of solute atoms during the solidification of metals. Acta Met. 1, 428–437 (1953).

CAS

Article

Google Scholar

32.

Burton, J. A., Prim, R. C. & Slichter, W. P. The distribution of solute in crystals grown from the melt. Part I. Theoretical. J. Chem. Phys. 21, 1987–1991 (1953).

ADS

CAS

Article

Google Scholar

33.

Watson, E. B. & Müller, T. Non-equilibrium isotopic and elemental fractionation during diffusion-controlled crystal growth under static and dynamic conditions. Chem. Geol. 267, 111–124 (2009).

ADS

CAS

Article

Google Scholar

34.

Smith, V. G., Tiller, W. A. & Rutter, J. A mathematical analysis of solute redistribution during solidification. Can. J. Phys. 33, 723–745 (1955).

ADS

CAS

Article

Google Scholar

35.

Glicksman, M. E. Principles of Solidification: An Introduction to Modern Casting and Crystal Growth Concepts (Springer Science & Business Media, 2010).
36.

DePaolo, D. J. Surface kinetic model for isotopic and trace element fractionation during precipitation of calcite from aqueous solutions. Geochim. Cosmochim. Acta 75, 1039–1056 (2011).

ADS

CAS

Article

Google Scholar

37.

Hofmann, A. W., Jochum, K. P., Seufert, M. & White, W. M. Nb and Pb in oceanic basalts: new constraints on mantle evolution. Earth Planet. Sci. Lett. 79, 33–45 (1986).

ADS

CAS

Article

Google Scholar

38.

Jahns, R. H. & Burnham, C. W. Experimental studies of pegmatite genesis; l, A model for the derivation and crystallization of granitic pegmatites. Econ. Geol. 64, 843–864 (1969).

CAS

Article

Google Scholar

39.

London, D. Ore-forming processes within granitic pegmatites. Ore Geol. Rev. 101, 349–383 (2018).

Article

Google Scholar

40.

Fournier, R. O. & Potter, R. W. II An equation correlating the solubility of quartz in water from 25 to 900 C at pressures up to 10,000 bars. Geochim. Cosmochim. Acta 46, 1969–1973 (1982).

ADS

CAS

Article

Google Scholar

41.

Burnham, C. W., Holloway, J. R. & Davis, N. F. The specific volume of water in the range 1000 to 8900 bars, 20 to 900 C. Am. J. Sci. 267, 70–95 (1969).

CAS

Google Scholar

42.

Nigrini, A. Diffusion in rock alteration systems; I, Prediction of limiting equivalent ionic conductances at elevated temperatures. Am. J. Sci. 269, 65–91 (1970).

ADS

CAS

Article

Google Scholar

43.

Yuan-Hui, L. & Gregory, S. Diffusion of ions in sea water and in deep-sea sediments. Geochim. Cosmochim. Acta 38, 703–714 (1974).

ADS

Article

Google Scholar

44.

Zhang, J. & Lee, C.-T. A. Disequilibrium crystallization and rapid crystal growth: a case study of orbicular granitoids of magmatic origin. Int. Geol. Rev. 1–17 (2020).
45.

Vance, D. & O’nions, R. K. Isotopic chronometry of zoned garnets: growth kinetics and metamorphic histories. Earth Planet. Sci. Lett. 97, 227–240 (1990).

ADS

CAS

Article

Google Scholar

46.

Dragovic, B., Baxter, E. F. & Caddick, M. J. Pulsed dehydration and garnet growth during subduction revealed by zoned garnet geochronology and thermodynamic modeling, Sifnos, Greece. Earth Planet. Sci. Lett. 413, 111–122 (2015).

ADS

CAS

Article

Google Scholar

47.

Swanson, S. E. Relation of nucleation and crystal-growth rate to the development of granitic textures. Am. Mineral. 62, 966–978 (1977).

CAS

Google Scholar

48.

Gualda, G. A. R. et al. Timescales of quartz crystallization and the longevity of the Bishop giant magma body. PLoS ONE 7, e37492 (2012).

ADS

CAS

PubMed

PubMed Central

Article

Google Scholar

49.

Pamukcu, A. S., Ghiorso, M. S. & Gualda, G. A. R. High-Ti, bright-CL rims in volcanic quartz: a result of very rapid growth. Contrib. Mineral. Petrol. 171, 105 (2016).

ADS

Article

CAS

Google Scholar

50.

DeMets, C., Gordon, R. G., Argus, D. F. & Stein, S. Current plate motions. Geophys. J. Int. 101, 425–478 (1990).

ADS

Article

Google Scholar

51.

Rowe, C. D. & Griffith, W. A. Do faults preserve a record of seismic slip: a second opinion. J. Struct. Geol. 78, 1–26 (2015).

ADS

Article

Google Scholar

52.

Wilson, L. O. On interpreting a quantity in the Burton, Prim and Slichter equation as a diffusion boundary layer thickness. J. Cryst. Growth 44, 247–250 (1978).

ADS

CAS

Article

Google Scholar

53.

Voloshin, A. E., Prostomolotov, A. I. & Verezub, N. A. On the accuracy of analytical models of impurity segregation during directional melt crystallization and their applicability for quantitative calculations. J. Cryst. Growth 453, 188–197 (2016).

ADS

CAS

Article

Google Scholar

54.

Turcotte, D. L. & Schubert, G. Geodynamics (Cambridge Univ. Press, 2002).
55.

Zhang, Y. & Xu, Z. Kinetics of convective crystal dissolution and melting, with applications to methane hydrate dissolution and dissociation in seawater. Earth Planet. Sci. Lett. 213, 133–148 (2003).

ADS

CAS

Article

Google Scholar

56.

Kerr, R. C. Convective crystal dissolution. Contrib. Mineral. Petrol. 121, 237–246 (1995).

ADS

CAS

Article

Google Scholar

57.

Kundu, P. K., Cohen, I. M. & Dowling, D. R. Fluid Mechanics (Elsevier Science, 2012).
58.

Behrens, H. & Zhang, Y. Ar diffusion in hydrous silicic melts: implications for volatile diffusion mechanisms and fractionation. Earth Planet. Sci. Lett. 192, 363–376 (2001).

ADS

CAS

Article

Google Scholar

59.

Watson, E. B. Diffusion of dissolved CO2 and Cl in hydrous silicic to intermediate magmas. Geochim. Cosmochim. Acta 55, 1897–1902 (1991).

ADS

CAS

Article

Google Scholar

60.

Ackerson, M. R., Tailby, N. D. & Watson, E. B. Trace elements in quartz shed light on sediment provenance. Geochem. Geophys. Geosyst. 16, 1894–1904 (2015).

ADS

CAS

Article

Google Scholar

61.

Garate-Olave, I., Müller, A., Roda-Robles, E., Gil-Crespo, P. P. & Pesquera, A. Extreme fractionation in a granite-pegmatite system documented by quartz chemistry: the case study of Tres Arroyos (Central Iberian Zone, Spain). Lithos 286, 162–174 (2017).

ADS

Article

CAS

Google Scholar

62.

Peacock, J. M. Temperature and leaf growth in four grass species. J. Appl. Ecol. 13, 225–232 (1976).

Download references
Acknowledgements

This work was supported by NSF-EAR 1753599 to C.-T.A.L. Discussions with Hehe Jiang, Ming Tang, Eytan Sharton-Bierig, Julin Zhang, and Blue Sheppard are appreciated.
Author information
Affiliations

Department of Earth, Environmental and Planetary Sciences, Rice University, Houston, TX, 77005, USA

Patrick R. Phelps & Cin-Ty A. Lee

Department of Earth and Planetary Sciences, University of California, Riverside, CA, 92521, USA

Douglas M. Morton

Contributions

P.R.P. analyzed samples, provided interpretation, and was involved in writing. C.-T.A.L. helped with interpretation and writing. D.M.M. supplied samples and geologic context.
Corresponding author

Correspondence to Patrick R. Phelps.
Ethics declarations
Competing interests

The authors declare no competing interests.
Additional information

Peer review information Nature Communications thanks Andrea Dini and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

One sentence summary: Trace element distributions of pegmatitic quartz crystals indicate rapid growth in highly dynamic environments, suggesting that large meter-scale crystals can be formed within days.
Supplementary information
Supplementary Information
Peer Review File
Description of Additional Supplementary Files
Supplementary Data 1


Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Phelps, P.R., Lee, CT.A. & Morton, D.M. Episodes of fast crystal growth in pegmatites. Nat Commun 11, 4986 (2020).

Download citation

Received21 January 2020

Accepted14 September 2020

Published05 October 2020
Back to top
Reply to topic Reply with quote
Like
   

James Catmur
Site Admin



Joined: 14 Sep 2006
Posts: 883
Location: Cambridge


Access to the FMF Gallery title=

View user's profile

Send private message

PostPosted: Sep 29, 2021 12:15    Post subject: Re: Just how fast can the crystals we love and collect grow?!  

Peter

For links to other web pages outside the forum you should use the same method when writing them as that used when adding email addresses.

For example: https://www.foro-minerales.com/forum/index.php, is a link within this forum so can be typed in directly, but had it been an external link you would have had to have added it as:

http(:)//www(.)foro(-)minerales(.)com/forum/index(.)php

How about:
https://www.nature.com/articles/s41467-020-18806-w
(link normalized by FMF)

James
____________________________________________

Recommended reading: Links within the message forum
Back to top
Reply to topic Reply with quote
Like
   

Peter Lemkin




Joined: 18 Nov 2016
Posts: 357
Location: Prague

View user's profile

Send private message

PostPosted: Sep 29, 2021 13:21    Post subject: Re: Just how fast can the crystals we love and collect grow?!  

Thanks James....I once knew that, but had forgotten. So all, if you want to read the article and the supplemental materials, diagrams, photos, etc. go to:

https://www.nature.com/articles/s41467-020-18806-w
(link normalized by FMF)
Back to top
Reply to topic Reply with quote
Like
1
   

alfredo
Site Admin



Joined: 30 Jan 2008
Posts: 871


Access to the FMF Gallery title=

View user's profile

Send private message

PostPosted: Oct 05, 2021 05:31    Post subject: Re: Just how fast can the crystals we love and collect grow?!  

Apart from pegmatite environments, I've personally seen specimens (miniatures to cabinet size) of other things growing in less than a month: Sulphur crystals in volcanic fumaroles. Gypsum crystals in acid mine wastewater.

Native copper specimens growing on wooden beams in abandoned mines can grow in a few decades.

At the opposite end of the time spectrum, one type of specimen that perhaps really does take millions of years to grow would be garnets in mica schist.
Back to top
Reply to topic Reply with quote
Like
3
   

Mathias




Joined: 01 Apr 2018
Posts: 132
Location: Vlodrop

View user's profile

Send private message

PostPosted: Oct 05, 2021 05:59    Post subject: Re: Just how fast can the crystals we love and collect grow?!  

....or Stalactites under a stairs platform........
Built in 1945, but perhaps a couple of times cleaned since then, and grown again.
Rainwater percolating through the concrete and creating stalactites.
Mathias.



001.jpg
 Description:
 Viewed:  370 Time(s)

001.jpg



002.jpg
 Description:
 Viewed:  362 Time(s)

002.jpg


Back to top
Reply to topic Reply with quote
Like
2
   

Peter Lemkin




Joined: 18 Nov 2016
Posts: 357
Location: Prague

View user's profile

Send private message

PostPosted: Oct 05, 2021 06:44    Post subject: Re: Just how fast can the crystals we love and collect grow?!  

For several years I lived in Aspen, Colorado - famous for its silver mines. Little known is that some of the old mines are being re-opened for specimens. One owner was kind enough to let me in. This mine was worked for the silver ore, although ocassionally in the mines they came across giant solid boulders of Ag. What I saw, and it was an amazing sight, were thousands of Ag wires of different shapes, lengths and thickness growing from cracks in the walls of the mine [which probably was closed around 1900]. So, these had grown in 90 years from when I saw them. I know there are also people who artificially grow Ag 'wires' using some solution and electricity. They grow quite large in a few days to weeks - a little more slowly under natural conditions.

I too used to grow crystals in jars and laboratory glassware when a child. Copper sulphate was one of my favorites and alum too - as they never failed to grow well, large and almost text-book shaped crystals [if one ignored the string in the center of some of them that the seed was hung from.... They grew in a few days, sometimes overnight depending on temperature and how saturated the solution was.
Back to top
Reply to topic Reply with quote
Like
1
   
Display posts from previous:   
   Index -> Minerals and Mineralogy   All times are GMT - 5 Hours
Page 1 of 1
    

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum
You cannot attach files in this forum
You can download files in this forum


All pictures, text, design © Forum FMF 2006-2021


Powered by FMF