Dense Melts Research Articles

Sulfide and Fe-Ti-P liquid immiscibility in the Ni-Cu-Co ovoid deposit of the Voisey’s Bay complex, Labrador, Canada

AbstractIn the Voisey’s Bay complex, sulfide-matrix breccias developed through the percolation of dense sulfide melt, leading to the displacement of the silicate melt within partially molten silicate-matrix breccias. In these sulfide matrix-breccias, hydrous silicate rims are commonly present at the interface between the sulfide matrix and the silicate framework. Multiple lines of evidence support a magmatic origin of these hornblende-biotite rims, which was largely coeval with the emplacement of the sulfide melt in the magmatic breccias. The formation of the hornblende-biotite rims required the addition of alkalis and water that could not have entirely been sourced from either the sulfide melt or the silicate framework. Through the integration of compositional maps with major and trace element analyses of the main accessory minerals, we propose that the critical components required for the development of the hydrous silicate rims in sulfide-matrix breccias originated from an immiscible Fe-Ti-P melt. Distinct textural and compositional features of apatite, hercynite, ilmenite and magnetite support the presence of small amounts of Fe-Ti-P melt in the sulfide melt. This Fe-Ti-P melt likely formed through melt immiscibility in the early stages of the development of the Voisey’s Bay complex, and was transported in the magma conduits together with the sulfide melt.

Unified Understanding of the Structure, Thermodynamics, and Diffusion of Single-Chain Nanoparticle Fluids.

Single-chain nanoparticles (SCNPs) are a fascinating class of soft nano-objects with promising properties and relevance to protein condensates, polymer nanocomposites, nanomedicine, bioimaging, catalysis, and drug delivery. We combine molecular dynamics simulations and equilibrium and time-dependent statistical mechanical theory to construct a unified understanding of how the internal conformational structure of SCNPs, of both a simple fractal globule-like form and more complex objects with multiple internal intermediate length scales, determines nm-scale intermolecular packing correlations, thermodynamic properties, and center-of-mass diffusion over a wide range of concentrations up to dense melts. The intermolecular pair correlations generically exhibit a distinctive deep correlation hole form due to SCNP internal connectivity structure and repulsive interparticle interactions associated with a globular-like conformation on the macromolecular scale, with concentration-dependent deviations at small separations. Unanticipated exponential-like dependences of the equation-of-state, osmotic compressibility, and center-of-mass diffusion constant on SCNP macromolecular packing fraction are theoretically predicted and confirmed via simulations. System-specific behaviors are found associated with SCNP internal structure, but overarching regularities are identified and understood based on a generalized effective globule conformation on macromolecular scales. Diffusivity slows down by 2-3 decades with increasing concentration and is understood as a consequence of a nonactivated excluded volume-driven weak-caging process associated with space-time correlated intermolecular forces experienced by the SCNP. Good agreement between the theory and simulations is established, testable predictions are made, and a quantitative comparison with viscosity measurements on a specific SCNP fluid is carried out. The basic theoretical approach can potentially be extended to treat the chemical and physical consequences of varying the structure of other classes of soft nanoparticles with distinctive internal nanoscale organization relevant in nanotechnology and nanomedicine, and the possible emergence of macromolecular kinetically arrested glasses.

Experimental Study on Calcination of Portland Cement Clinker Using Different Contents of Stainless Steel Slag.

In order to increase the utilization rate of stainless steel slag, reduce storage needs, and mitigate environmental impacts, this study replaces a portion of limestone with varying amounts of stainless steel slag in the calcination of Portland cement clinker. The study primarily examines the influence of stainless steel slag on the phase composition, microstructure, compressive strength, and free calcium oxide (ƒ-CaO) content of Portland cement clinker. The results show the following: (1) Using stainless steel slag to calcine Portland cement clinker can lower the calcination temperature, reducing industrial production costs and energy consumption. (2) With an increase in the amount of stainless steel slag, the dicalcium silicate (C2S) and tricalcium silicate (C3S) phases in Portland cement clinker initially increase and then decrease; the C3S crystals gradually transform into continuous hexagonal plate-shaped distributions, while the tricalcium aluminate (C3A) and tetracalcium aluminoferrite (C4AF) crystal structures become denser. When the stainless steel slag content is 15%, the dicalcium silicate and tricalcium silicate phases are at their peak; the C3S crystals are continuously distributed with a relatively dense structure, and C3A and C4AF crystals melt and sinter together, becoming distributed around C3S. (3) As stainless steel slag content increases, the compressive strength of Portland cement clinker at 3 days, 7 days, and 28 days increases and then decreases, while ƒ-CaO content decreases and then increases. When the stainless steel slag content is 15%, the compressive strength at 28 days is at its highest, 64.4 MPa, with the lowest ƒ-CaO content, 0.78%. The test results provide a basis for the utilization of stainless steel slag in the calcination of Portland cement clinker.

Ultrafast x-ray detection of low-spin iron in molten silicate under deep planetary interior conditions.

The spin state of Fe can alter the key physical properties of silicate melts, affecting the early differentiation and the dynamic stability of the melts in the deep rocky planets. The low-spin state of Fe can increase the affinity of Fe for the melt over the solid phases and the electrical conductivity of melt at high pressures. However, the spin state of Fe has never been measured in dense silicate melts due to experimental challenges. We report detection of dominantly low-spin Fe in dynamically compressed olivine melt at 150 to 256 gigapascals and 3000 to 6000 kelvin using laser-driven shock wave compression combined with femtosecond x-ray diffraction and x-ray emission spectroscopy using an x-ray free electron laser. The observation of dominantly low-spin Fe supports gravitationally stable melt in the deep mantle and generation of a dynamo from the silicate melt portion of rocky planets.

Ultrahigh‐Pressure Acoustic Velocities of Aluminous Silicate Glass up to 155 GPa With Implications for the Structure and Dynamics of the Deep Terrestrial Magma Ocean

AbstractWe have carried out in situ high‐pressure acoustic velocity measurements of (Fe2+, Al)‐bearing MgSiO3 glass up to pressures of 155 GPa, which confirmed a distinct pressure‐induced trend change in the transverse acoustic velocity (VS) profile around 98 GPa, likely caused by the Si‐O coordination number (CN) change from 6 to 6+. Although it has been reported that the substitution of Fe2+ in MgSiO3 glass induces almost linear velocity reduction up to ∼160 GPa, we revealed that the VS profile of (Fe2+, Al)‐bearing MgSiO3 becomes anomalously steeper above ∼100 GPa and eventually came to be equivalent to MgSiO3 glass above ∼125 GPa. This implies the incorporation of Al into Fe‐bearing MgSiO3 glass significantly facilitates making it far elastically stiffer and thus the densification under pressures well within the Earth's lower mantle. Our results indicate the possible presence of stiff and highly dense silicate melts in deep MOs in the rocky terrestrial planets.

Quasi-stable structures in equilibrium dense bismuth melt: Experimental and first principles theoretical studies

Near the melting temperature, equilibrium bismuth melt is characterized by structural features that are absent in equilibrium monatomic simple liquids. In the present work, the structure of bismuth melt is studied by X-ray diffraction experiments and quantum chemical calculations. The presence of quasi-stable structures in the melt has been found, the lifetime of which exceeds the structural relaxation time of this melt. It is shown that these structures are characterized by a low degree of ordering and spatial localisation. It was found that up to 50% of the atoms in the melt can be involved in the formation of these structures. The elementary structural units of these structures are triplets of regular geometry with the characteristic lengths 3.25 Å and 4.7 Å as well as with the characteristic angles 45∘ and 90∘. The characteristic lengths of these triplets are fully consistent with correlation lengths associated with the short-range order in bismuth melt.

Dilution Driven Self-assembly and Re-entrant Phase Transition in Molecular Hydrogels

TX-100 molecular hydrogels exhibited re-entrant melt-gel-sol phase transition driven solely by dilution, which sequentially altered the self-assembly of the micellar formations and their population was investigated through the monitoring of the physical parameters namely, solution viscosity, particle size histogram, ergodicity, and gel rigidity modulus. This phenomenon was noticed at 20°C in the TX-100 concentration region of 0.2 to 1.2 M much above the critical micellar concentration of 0.22 mM. The particle size histograms revealed the presence of spherical micelles (size »3 nm) in the solution ([TX-100] < 0.5 M) which formed entangled wormlike cylindrical micelles (apparent hydrodynamic radius » 50 nm) when (0.5 M< [TX-100] < 0.9 M) giving rise to a gel-like structure. Further increase in the TX-100 concentration increased the propensity of these wormlike cylindrical micelles that got randomly distributed creating a dense melt phase. Interestingly, we observed transition solely driven by dilution which defined complete re-entrant behavior at room temperature. These molecular gels could be created by dilution of the melt or concentration of the sol unlike in the polymer gels. Remarkably, this hitherto little known unique phenomenon was exhibited by a simple system of non-ionic surfactant solution. Thus, we have a hydration reversible gel at our disposal which has a special place in soft matter arena.

Pressure-induced structural evolution in boron-bearing model rhyolitic glasses under compression: Implications for boron isotope compositions and properties of deep melts in Earth’s interior

The pressure-induced structural evolutions of boron-bearing model rhyolitic melts under high pressures enable to infer the detailed geochemical processes (melting and fluid-rock-melt interactions) occurring in Earth interiors and to control the melt properties (viscosity and the boron isotope composition, δ11B) of complex magmatic melts, providing insights into the boron cycle toward the deeper part of the upper mantle (∼10 GPa). Despite the importance, the structures of multicomponent boron-bearing silicate melts above 3 GPa are currently unavailable. Here, we explore the structures, particularly, coordination transformation of constituent elements in boron-bearing nepheline and albite glasses – a model rhyolitic melts - upon compression to a depth of ∼270 km (∼9.2 GPa) in the mantle using multi-nuclear solid-state nuclear magnetic resonance (NMR) spectroscopy. The results showed that the conversion of [3]B into [4]B is prominent upon compression up to 6 GPa. In contrast, the formation of [5,6]Al is accompanied by the formation of oxygen tricluster above 6 GPa, where all the nonbridging oxygens are consumed. We quantify how the melt composition affects tendency to form highly coordinated B, Al, and Si upon compression. Particularly, the increase in the [4]B population tends to be larger for the glasses with low Si content as pressure increases to 9.2 GPa. We reveal the relationship between such structural adaptations of the compressed melts at high pressure and the melt properties, including viscosity and element partition coefficient in boron-bearing melts. The current NMR results also unravel the structural origins of 11B/10B ratios in rhyolitic melts at high pressure. Considering a preferential partitioning of 10B to [4]B, an increase in [4]B population in the melts leads to an pressure-induced enrichment of 10B. As the increase in Si/B ratio in the melts tends to decrease the pressure-induced increase in [4]B fraction, the contribution of boron coordination transformation on the 11B/10B ratios in silicate melt would be somewhat minor in deep mantle melts with increasing Si content. The detailed boron environments in rhyolitic melts at high pressure yield useful constraints for the isotope composition (11B/10B) of dense mantle melts, thereby enabling quantification of deep boron cycle.

Statistical potentials from the Gaussian scaling behaviour of chain fragments buried within protein globules.

Knowledge-based approaches use the statistics collected from protein data-bank structures to estimate effective interaction potentials between amino acid pairs. Empirical relations are typically employed that are based on the crucial choice of a reference state associated to the null interaction case. Despite their significant effectiveness, the physical interpretation of knowledge-based potentials has been repeatedly questioned, with no consensus on the choice of the reference state. Here we use the fact that the Flory theorem, originally derived for chains in a dense polymer melt, holds also for chain fragments within the core of globular proteins, if the average over buried fragments collected from different non-redundant native structures is considered. After verifying that the ensuing Gaussian statistics, a hallmark of effectively non-interacting polymer chains, holds for a wide range of fragment lengths, although with significant deviations at short spatial scales, we use it to define a 'bona fide' reference state. Notably, despite the latter does depend on fragment length, deviations from it do not. This allows to estimate an effective interaction potential which is not biased by the presence of correlations due to the connectivity of the protein chain. We show how different sequence-independent effective statistical potentials can be derived using this approach by coarse-graining the protein representation at varying levels. The possibility of defining sequence-dependent potentials is explored.

Fractional crystallization of a basal lunar magma ocean: A dense melt-bearing garnetite layer above the core?

In the wake of its accretion, the Moon was likely partially or fully molten, forming a lunar magma ocean (LMO). A fully molten Moon may crystallize neutrally buoyant olivine, whose accumulation may form a barrier leading to a separation into two independently evolving melt reservoirs. This study investigates the crystallization of such a putative basal lunar magma ocean, ranging from the core mantle boundary (with radius r = 380 km) to the level of neutral olivine buoyancy (r = 600 km). Consecutive crystallization experiments determine liquidus temperatures and crystallize 10–30 wt% minerals, before conceptually segregating these according to their buoyancy, and then stepping to a new bulk composition that corresponds to the residual melt after removal of the cumulate minerals.The first olivine to crystallize from a Taylor Whole Moon composition (twm, XMg = Mg/(Mg + Fe2+) = 0.83) has XMg = 0.94 and is neutrally buoyant at 3.8 GPa according to the melt density model of Lange and Carmichael (1990). Crystallization begins at the core mantle boundary, but early olivine floats and re-dissolves until the magma ocean cools to the liquidus temperature at the depth of neutral buoyancy (1850 °C). At this point a > 500 km thick olivine-only layer could form, mainly fed from the upper magma shell. The crystallization sequence in the basal magma ocean is olivine-only at 1850–1675 °C, 0–26 pcs (wt% percent solidified) → olivine + opx (to 1600 °C, 42 pcs) → opx + cpx + garnet (to 1580 °C, 60 pcs) → cpx + garnet (to 1520 °C, 73 pcs) → cpx + garnet + olivine leaving 11 wt% residual melt at 1450 °C.Cooling this last melt to 1300 °C leads to 80% garnet + cpx + olivine + Ti-spinel, the residual liquid corresponding to 2.2% of the basal magma ocean. Olivine, opx and cpx remain buoyant over the entire crystallization interval and various pyroxenite layers are added to the olivine layer. Garnet and the final cotectic cpx + garnet + olivine + FeTi-oxide assemblage instead form a 70 km thick basal layer on the core-mantle boundary. Such a layer provides a gravitationally stable high-density lowermost lunar mantle, which would concur with a recent re-analysis of the lunar seismic data. The lowest temperature experiment at 1300 °C, not much above the 1250 °C proposed for the core-mantle boundary from inversion of geophysical data, had about 20% of a highly evolved, dense Fe-rich melt. It is feasible that this melt has remained in its liquid state to present day, providing an explanation for the proposed low vp and vs layer and high dissipation just above the core-mantle boundary.

Fully Automated Detection of Supraglacial Lake Area for Northeast Greenland Using Sentinel-2 Time-Series

The usability of multispectral satellite data for detecting and monitoring supraglacial meltwater ponds has been demonstrated for western Greenland. For a multitemporal analysis of large regions or entire Greenland, largely automated processing routines are required. Here, we present a sequence of algorithms that allow for an automated Sentinel-2 data search, download, processing, and generation of a consistent and dense melt pond area time-series based on open-source software. We test our approach for a ~82,000 km2 area at the 79 °N Glacier (Nioghalvfjerdsbrae) in northeast Greenland, covering the years 2016, 2017, 2018 and 2019. Our lake detection is based on the ratio of the blue and red visible bands using a minimum threshold. To remove false classification caused by the similar spectra of shadow and water on ice, we implement a shadow model to mask out topographically induced artifacts. We identified 880 individual lakes, traceable over 479 time-steps throughout 2016–2019, with an average size of 64,212 m2. Of the four years, 2019 had the most extensive lake area coverage with a maximum of 333 km2 and a maximum individual lake size of 30 km2. With 1.5 days average observation interval, our time-series allows for a comparison with climate data of daily resolution, enabling a better understanding of short-term climate-glacier feedbacks.

Structural Transitions in MgSiO3 Glasses and Melts at the Core‐Mantle Boundary Observed via Inelastic X‐ray Scattering

AbstractThe structural adaptation in MgSiO3 melts under compression up to 130 GPa is the key to revealing the origins of the pronounced negative buoyancy of the melts at the core‐mantle boundary (CMB). A full understanding of the melt densification requires study of the pressure‐induced changes in the bonding configuration around oxygen at the CMB, which has proven to be difficult to measure. Here, the experimental breakthrough in O K‐edge inelastic X‐ray scattering enables collection of the spectra of MgSiO3 glasses up to ~130 GPa, along with ab initio molecular dynamics simulations, revealing the electronic bonding transitions around heavily compressed oxygen. The spectral results indicate the emergence of denser network structures around oxygen, stemming from contractions in the Mg‐O and O‐O distances associated with flexible topological and short‐range rearrangements around Si. The results unveil the electronic structure and thus the nature of densification in dense partial melts at the CMB.

Wetting Behavior of Iron–Carbon Melt in Silicates at Mid-Mantle Pressures With Implications for Earth’s Deep Carbon Cycle

Metallic melt containing iron (Fe) and carbon (C) may be present at depths greater than 250 km inside the Earth. Depending on its wetting behavior, such dense melt may be trapped locally or drain into deep mantle and core. Here, we report experimental data on the wetting behavior of Fe-C melt in silicates at the conditions of Earth’s mid-mantle between 10 and 23 GPa and 1600 and 1800 ˚C. The measured dihedral angles of Fe-C melt in olivine, ringwoodite or bridgmanite and ferropericlase matrixes are 117±14°, 120±14° and 107±16° respectively, well above the critical value of 60° for complete wetting. The estimated percolation thresholds are at least 7% in volume, far exceeding the amount of metal in the mantle. Consequently, slab-derived Fe-C melt in the mid-mantle is expected to occur as isolated pockets and would not percolate through its silicate matrix.

Partial Melting‐Induced Chemical Evolution in Shocked Crystalline and Amorphous Plagioclase From the Lunar Meteorite Mount DeWitt 12007

AbstractDetermining the formation mechanism of maskelynite is essential to understanding the shocked environments of meteorites on their parent bodies. Maskelynite has been accepted as a diaplectic glass for several decades, but there have been suggestions that it is a normal glass quenched from a dense melt. Morphological characteristics have been generally investigated to identify the formation mechanism of amorphous plagioclase in meteorites, but the chemical difference between crystalline and amorphous plagioclase has not been fully understood. In this study, we investigated the morphological, atomic‐scale structural, and chemical characteristics of amorphous plagioclase in the lunar meteorite DEW 12007 to constrain its formation mechanism via chemical analysis. The morphological characteristics showed that plagioclase was partially converted into amorphous phase through partial melting. Two‐dimensional Raman mapping confirmed the structural difference between amorphous and crystalline regions. Quantitative chemical analyses revealed that the amorphous regions were more albite‐rich than the crystalline regions, likely due to the partial melting of plagioclase. Under shocked conditions, the partial melting of plagioclase induced a chemical variation between amorphous and crystalline regions. The morphological and structural changes correspond well with the chemical variations, indicating that amorphization induced such variations. The chemical differences between amorphous and crystalline plagioclase in other meteorites also could be understood to be the results of partial melting. Thus, the chemical differences between amorphous and crystalline plagioclase in partially amorphized grains could elucidate the formation mechanism of amorphous plagioclase in many meteorites.

The degree of polymerization and structural disorder in (Mg,Fe)SiO3 glasses and melts: Insights from high-resolution 29Si and 17O solid-state NMR

Probing the degree of polymerization and structural disorder of (Mg,Fe)SiO3 glasses and melts remains a challenging problem, despite their strong implications for the properties of natural volcanic glasses and dense melts in Earth’s surface and interior. Here, we explore the iron-induced changes in silicon and the oxygen configurations of (Mg1−XFeX)SiO3 glasses (Fe3+/ƩFe = 0.22 ± 0.07) with varying XFe [=Fe/(Mg + Fe), containing up to 13.5 wt% FeO] using high-resolution 29Si and 17O MAS, and 2D 17O 3QMAS NMR, unveiling the effect of Fe2+ on the extent of disorder in these glasses. The 29Si MAS NMR spectra show the apparent peak broadening in less polymerized species (i.e., Q1) and an increase in the fraction of highly polymerized (i.e., Q3) species with increasing XFe. These changes imply the selective decrease in NMR signal intensity of Q1 and Q2 due to their spatial proximity to Fe2+ (i.e., paramagnetic effect), and slight increase in the degree of polymerization due to the presence of minor fraction of Fe3+. The 2D 17O NMR spectra revealed that the fraction of Si-O-Si apparently increases with XFe and more significant changes in peak width for Mg-O-Si are observed with varying XFe. The results indicate that Mg-O-Si is more strongly affected by Fe2+ than Si-O-Si, suggesting a moderate degree of preferential proximity between Fe2+ and non-bridging oxygen. The results allow us to conceptualize the scheme to account for paramagnetic interaction between the nuclear spins and Fe2+. The scheme provides a systematic protocol to interpret the NMR spectra of diverse iron-bearing silicate glasses. The current NMR results, in conjunction with the results in the previous study for Fe3+-dominant sodium silicate glasses, suggest that paramagnetic effect depends on the valence state of Fe, indicating the presence of a moderate degree of preferential interaction between Fe2+ and NBO and less polymerized Si units as well as that between Fe3+ and highly polymerized Si units and BO. Though qualitative, the current NMR results for (Mg,Fe)SiO3 glasses shed light on an opportunity to explore the detailed structures of iron-bearing natural volcanic glasses and can improve understanding of the transport properties of the basaltic melts.

Fingerprint of silicic magma degassing visualised through chlorine microscopy

Volatile-rich silicic magma erupts either explosively as a jet of a mixture of pyroclasts and high-temperature gas, or non-explosively to effuse lava. The bifurcation of the eruption style is widely recognised as being controlled by the efficiency of open-system gas loss from vesiculated magma during ascent. However, the fundamental question of how the gas escapes from highly viscous magma still remains unsolved because the pathways of gas flow are rarely preserved in dense lava. Here we show that such pathways are visualised in groundmass glass using high-resolution chlorine (Cl) mapping analysis on the rhyolitic lava of the Mukaiyama volcano, Japan. The results showed that the glass was highly heterogeneous in Cl content. A spatial distribution of the Cl content in the groundmass glass showed that volatiles diffused towards most bubbles, but the bubbles collapsed into the dense melt rather than growing. All observations, in combination with melt inclusion analysis, indicate that vesiculation, the formation of interconnected bubble channels, open-system gas loss via the channels, and channel collapse repeated within the period of a few days to two weeks during ascent. This cycle repeated individually in centimetre-sized portions of magma with different timing.

Osmotic pressure of confined square lattice self-avoiding walks

Flory–Huggins theory (Flory 1942 J. Chem. Phys. 10 51–61; Huggins 1942 J. Am. Chem. Soc. 64 2716–8) is a mean field theory for modelling the free energy of dense polymer solutions and polymer melts. In this paper we use Flory–Huggins theory as a model of a dense two-dimensional self-avoiding walk compressed in a square in the square lattice. The theory describes the free energy of the walk well, and we estimate the Flory interaction parameter of the walk to be .

Effects of Copolymer Sequence on Adsorption and Dynamics Near Nanoparticle Surfaces in Simulated Polymer Nanocomposites

ABSTRACT We study a simple nanocomposite, consisting of a single spherical nanoparticle in a dense melt of coarse-grained copolymer chains, using molecular dynamics simulations. The polymers contain equal amounts of two monomer types, which differ only in their monomer-nanoparticle interaction (adsorption) strengths, and are placed in a random sequence or in a sequence of alternating blocks with various block lengths. This model is motivated by a need to understand copolymer sequence effects relevant to designing tire tread compounds, given the synthetic ability to adjust copolymer architecture (e.g., amount of blockiness) in styrene-butadiene rubbers and the ability to choose fillers and covering agents that functionalize the filler particles to give them different affinities for the styrene and butadiene monomers (it is also possible to use coupling agents, which form covalent bonds between filler particles and polymer chains). Our simple, generalized model considers linear polymers of uniform length that are not cross-linked, and thus we do not attempt to capture the overall mechanical properties of tire tread materials. Instead, we focus our analysis on how copolymer sequence impacts polymer adsorption on the nanoparticle, known to be a significant factor in determining polymer nanocomposite properties. We find that copolymer sequence impacts both the range and magnitude of the interphase of slowed dynamics surrounding the nanoparticle, with longer block lengths yielding a greater reduction in mobility over a wider region.

Simulation of ring polymer melts with GPU acceleration

We implemented the elastic lattice polymer model on the GPU (Graphics Processing Unit), and show that the GPU is very efficient for polymer simulations of dense polymer melts. The implementation is able to perform up to 4.1⋅109 Monte Carlo moves per second. Compared to our standard CPU implementation, we find an effective speed-up of a factor 92. Using this GPU implementation we studied the equilibrium properties and the dynamics of non-concatenated ring polymers in a melt of such polymers, using Rouse modes. With increasing polymer length, we found a very slow transition to compactness with a growth exponent ν≈1/3. Numerically we find that the longest internal time scale of the polymer scales as N3.1, with N the molecular weight of the ring polymer.

Percolative core formation in planetesimals enabled by hysteresis in metal connectivity

The segregation of dense core-forming melts by porous flow is a natural mechanism for core formation in early planetesimals. However, experimental observations show that texturally equilibrated metallic melt does not wet the silicate grain boundaries and tends to reside in isolated pockets that prevent percolation. Here we use pore-scale simulations to determine the minimum melt fraction required to induce porous flow, the percolation threshold. The composition of terrestrial planets suggests that typical planetesimals contain enough metal to overcome this threshold. Nevertheless, it is currently thought that melt segregation is prevented by a pinch-off at melt fractions slightly below the percolation threshold. In contrast to previous work, our simulations on irregular grain geometries reveal that a texturally equilibrated melt network remains connected down to melt fractions of only 1 to 2%. This hysteresis in melt connectivity allows percolative core formation in planetesimals that contain enough metal to exceed the percolation threshold. Evidence for the percolation of metallic melt is provided by X-ray microtomography of primitive achondrite Northwest Africa (NWA) 2993. Microstructural analysis shows that the metal-silicate interface has characteristics expected for a texturally equilibrated pore network with a dihedral angle of ∼85°. The melt network therefore remained close to textural equilibrium despite a complex history. This suggests that the hysteresis in melt connectivity is a viable process for percolative core formation in the parent bodies of primitive achondrites.