Diffusion Time Scale Research Articles

Effects of elastic softening and helium accumulation kinetics on surface morphological evolution of plasma-facing tungsten

Based on a continuous-domain model, capable of accessing the spatiotemporal scales relevant to fuzz formation on the surface of plasma-facing component (PFC) tungsten, we report self-consistent simulation results that elucidate the effects of elastic softening and helium (He) accumulation kinetics on the surface morphological response of PFC tungsten. The model accounts for the softening of the elastic moduli in the near-surface region of PFC tungsten, including both thermal softening at high temperature and softening due to He accumulation upon He implantation. The dependence of the elastic moduli on the He content follows an exponential scaling relation predicted by molecular-dynamics simulations, while the He content in the near-surface region of PFC tungsten evolves according to a first-order saturation kinetics, consistent with experimental and simulation results reported in the literature. We establish that this elastic softening accelerates both nanotendril growth on the PFC surface and the onset of fuzz formation. We also explore the role of the rate of He accumulation to a saturation level in the near-surface region of irradiated tungsten in the onset of fuzz formation. For PFC tungsten surfaces such as W(110) where, under typical irradiation conditions, the characteristic time scale for stress-driven surface diffusion is comparable to the characteristic time scale for He accumulation, we find that accelerating the rate of He accumulation accelerates the growth rate of nanotendrils emanating from the surface. Additionally, we present a systematic parametric study of the PFC surface morphological response to explore its dependence on the He accumulation kinetics that is controlled by the irradiation conditions for low-energy implantation. Finally, we introduce an incubation time for nanotendril growth on the PFC surface, a concept equivalent to that of incubation fluence discussed in the literature, to predict and explain the minimum exposure time required to observe fuzz formation on PFC tungsten surfaces.

Correlated atomic dynamics in liquid seen in real space and time.

In liquids, the timescales for structure, diffusion, and phonon are all similar, of the order of a pico-second. This not only makes characterization of liquid dynamics difficult but also renders it highly questionable to describe liquids in these terms. In particular, the current definition of the structure of liquids by the instantaneous structure may need to be expanded because the liquid structure is inherently dynamic. Here, we advocate describing the liquid structure through the distinct-part of the Van Hove function, which can be determined by inelastic neutron and x-ray scattering measurements as well as by simulation. It depicts the dynamic correlation between atoms in space and time, starting with the instantaneous correlation function at t = 0. The observed Van Hove functions show that the atomic dynamics is strongly correlated in some liquids, such as water. The effect of atomic correlation on various transport properties of fluid, including viscosity and diffusivity, is discussed.

Quantification of fluctuations from fluorescence correlation spectroscopy experiments in reaction-diffusion systems.

Fluorescence correlation spectroscopy (FCS) is commonly used to estimate diffusion and reaction rates. In FCS the fluorescence coming from a small volume is recorded and the autocorrelation function (ACF) of the fluorescence fluctuations is computed. Scaling out the fluctuations due to the emission process, this ACF can be related to the ACF of the fluctuations in the number of observed fluorescent molecules. In this paper the ACF of the molecule number fluctuations is studied theoretically, with no approximations, for a reaction-diffusion system in which the fluorescence changes with binding and unbinding. Theoretical ACFs are usually derived assuming that fluctuations in the number of molecules of one species are instantaneously uncorrelated to those of the others and obey Poisson statistics. Under these assumptions, the ACF derived in this paper is characterized only by the diffusive timescale of the fluorescent species and its total weight is the inverse of the mean number of observed fluorescent molecules. The theory is then scrutinized in view of previous experimental results which, for a similar system, gave a different total weight and correct estimates of other diffusive timescales. The total weight mismatch is corrected by assuming that the variance of the number of fluorescent molecules depends on the variance of the particle numbers of the other species, as in the variance decomposition formula. Including the finite acquisition time in its computation, it is shown that the ACF depends on various timescales of the system and that its total weight coincides with the one obtained with the variance decomposition formula. This calculation implies that diffusion coefficients of nonobservable species can be estimated with FCS experiments performed in reaction-diffusion systems. Ways to proceed in future experiments are also discussed.

Gamma-Ray and Radio Background Constraints on Cosmic Rays in Milky Way Circumgalactic Medium

Abstract We study the interaction of cosmic rays (CRs) with the diffuse circumgalactic gas of the Milky Way (MW) galaxy that results in hadronic γ-ray emission and radio synchrotron emission. We aim to constrain the CR population in our circumgalactic medium (CGM) with the help of the observed isotropic γ-ray background (IGRB), its anisotropy, and radio continuum. We modify different models of CGM gas in hydrostatic equilibrium discussed in the literature by including a cosmic-ray population, parameterized by η ≡ P CR/P th. For the simplest isothermal model, while the IGRB intensity allows η ≲ 3, the anisotropy resulting from the solar system’s off-center position in MW rules out all values of η. For the precipitation model, in which the cooling of the CGM gas is regulated with an optimum ratio of cooling time to freefall time, while the observed IGRB intensity allows η ≲ 230, the observed anisotropy allows only very large values of η, of order ≳100. The radio continuum limits η ≲ 400 for the precipitation model and does not constrain the isothermal model; however, these constraints are mitigated by synchrotron loss time being comparable to CR diffusion timescales. These bounds are relevant for current numerical simulations that indicate a significant CR population in CGM of galaxies of MW mass.

Modeling and Experiments of Binary Electrolytes in the Presence of Diffusion, Migration, and Electro-Osmotic Flow

Combined diffusion, migration, and advection of ions in a binary electrolyte plays a role in various applications, including water electrolysis, electrodeposition, deionization, and electrophoresis. Here we analyze a dilute binary electrolyte with arbitrary ion valencies in a porous or nonporous medium using the one-dimensional Nernst-Planck equations. We examine how advection influences the limiting current, diffusion potential, and overall potential, deriving broadly useful analytical expressions. We provide experimental results for the electro-osmotic flow through a submicroporous separator in an alkaline water electrolysis setup. The time evolution of the potential is followed from the initial timescale of double-layer charging, followed by the diffusional timescale, to the time at which a limiting current is reached. For the longer timescales, a quasisteady model is shown to predict the time evolution of the advection-modified potential drop reasonably well. Additional interesting features arising due to electro-osmotic drag and unsteady diffusion are observed and explained.

An ultraviolet dyegraph for measuring the chemical disturbances of sinking particles and swimming plankton

AbstractSinking particles and swimming plankton may deform ambient chemical gradients, leaving long trails of chemical disturbance that dissipate on a diffusive timescale. Imaging this phenomenon in a laboratory setting requires a large field of view, a broad depth of field, and a passive chemical substance that does not influence the hydrodynamics. These requirements obviate the use of traditional schlieren optical systems, which rely on density stratification and their consequent buoyancy forces, and modern planar laser‐induced fluorescence, which has a narrow depth of field. Further, the light used to induce plankton movement must be kept separate from the recorded signal. Here, we describe a new ultraviolet “dyegraph” technique that quantifies dye disturbance caused by sinking spheres and swimming zooplankton. An experimental tank is filled with seawater and a gradient of an ultraviolet‐absorbing pigment. Using modified shadowgraph optics, ultraviolet light spreads from a pinhole to a collimating lens, passes through the tank to a telecentric lens, and is recorded by a camera fitted with an ultraviolet bandpass filter. The dye is transparent in visible light but appears dark in the imaging system. The 3D chemical disturbance is reconstructed from images using tomographic inversion. After outlining the method, we present examples of chemical deformation trails generated by sinking spheres and live, swimming copepods.

Diffusion of radial action in a galactic disc

Context. The stellar migration of the galactic disc stars has been invoked to explain the dispersion of stellar metallicity observed in the solar neighbourhood. Aims. We seek to identify the dynamical mechanisms underlying stellar migration in an isolated galaxy disc under the influence of a bar. Our approach is to analyse the diffusion of dynamical quantities. Methods. We extend our previous work by exploring Chirikov’s diffusion rate (and derived timescale) of the radial action JR in an idealised N-body simulation of an isolated disc galaxy. We limit our study to the evolution of the disc region well after the formation of the bar, in a regime of adiabatic evolution. Results. The JR diffusion timescale TD(JR) is less than 3 Gyr for roughly half the galaxy mass. It is always much shorter than the angular momentum diffusion timescale TD(Lz) outside the stellar bar. In the disc, ⟨TD(JR)⟩ ∼ 1 Gyr. All non-axisymmetric morphological structures that are characteristic of resonances and waves in the disc are associated to particles with TD(JR) < 3 Gyr and TD(Lz) > 10 Gyr. Short TD(JR) can be explained by the gradual de-circularisation of initially circular orbits (JR = 0) under the effect of intermittent. Inner Linblad resonance scattering by wave trains propagating in the disc, well beyond the outer Lindblad resonance of the bar (OLR). This leads to a moderate secular heating of the disc beyond the bar’s OLR for 7 Gyr, which is comparable to solar neighbourhood observations. The complex multi-wave structure, mixing permanent and intermittent modes, allows for multiple resonance overlaps.

Multiscale electronic and thermomechanical dynamics in ultrafast nanoscale laser structuring of bulk fused silica

We describe the evolution of ultrafast-laser-excited bulk fused silica over the entire relaxation range in one-dimensional geometries fixed by non-diffractive beams. Irradiation drives local embedded modifications of the refractive index in the form of index increase in densified glass or in the form of nanoscale voids. A dual spectroscopic and imaging investigation procedure is proposed, coupling electronic excitation and thermodynamic relaxation. Specific sub-ps and ns plasma decay times are respectively correlated to these index-related electronic and thermomechanical transformations. For the void formation stages, based on time-resolved spectral imaging, we first observe a dense transient plasma phase that departs from the case of a rarefied gas, and we indicate achievable temperatures in the excited matter in the 4,000–5,500 K range, extending for tens of ns. High-resolution speckle-free microscopy is then used to image optical signatures associated to structural transformations until the evolution stops. Multiscale imaging indicates characteristic timescales for plasma decay, heat diffusion, and void cavitation, pointing out key mechanisms of material transformation on the nanoscale in a range of processing conditions. If glass densification is driven by sub-ps electronic decay, for nanoscale structuring we advocate the passage through a long-living dense ionized phase that decomposes on tens of ns, triggering cavitation.

ECCD-induced sawtooth crashes at W7-X

The optimised superconducting stellarator W7-X generates its rotational transform by means of external coils, therefore no toroidal current is necessary for plasma confinement. Electron cyclotron current drive experiments were conducted for strikeline control and safe divertor operation. During current drive experiments periodic and repetitive crashes of the central electron temperature, similar to sawtooth crashes in tokamaks, were detected. Measurements from soft x-ray tomography and electron cyclotron emission show that the crashes are preceded by weak oscillating precursors and a displacement of the plasma core, consistent with a (m, n) = (1, 1) mode. The displacement occurs within , followed by expulsion and redistribution of the core into the external part of the plasma. Two types of crashes, with different frequencies and amplitudes are detected in the experimental program. For these non-stationary parameters a strong dependence on the toroidal current is found. A 1-D heuristic model for current diffusion is proposed as a first step to explain the characteristic crash time. Initial results show that the modelled current diffusion timescale is consistent with the initial crash frequency and that the toroidal current rise shifts the position where the instability is triggered, resulting in larger crash amplitudes.

Pre-eruptive storage, evolution, and ascent timescales of a high-Mg basaltic andesite in the southern Cascade Arc

Crustal storage and transit play a critical role in the compositional evolution of arc magmas; however, the enigmatic nature of lower crustal magma storage and early differentiation limit our understanding of the connections between the physical processes of subduction zones and the architecture of the arc crust. We present new geochemical data and applications of existing barometric and chronometric tools to interrogate the mantle source compositions, crustal storage depths, and ascent timescales of a primitive, high-Mg, clinopyroxene-bearing cinder cone, the basaltic andesite of Box Canyon, located in the vicinity of Lassen Volcanic National Park, CA, in the southern Cascades. Petrographic examination in addition to bulk and in situ geochemical analyses (XRF, LA-ICP-MS, and EPMA) of tephra and lava-derived samples reveals co-crystallization of clinopyroxene and olivine as phenocrysts and glomerocrysts with ~ Mg# 80 and > Fo85, respectively. Phase equilibria experiments of analogous Cascade Arc magma compositions estimate crustal storage of the observed phase assemblage at pressures in the lower crust > 700 MPa. Reverse zonation in olivine and clinopyroxene phenocryst interiors from core-rim analytical profiles record a lower crustal mafic mixing event. Results from one-dimensional, multi-elemental olivine and clinopyroxene diffusion models fit to these interior mixing zones provide an assessment of available trace element (Ti, La, Yb, Ce, and Nd) diffusion chronometers in clinopyroxene by considering multiple element profiles within two phases that experienced the same pre-eruptive conditions. While multi-elemental and multi-phase diffusion timescales span two to three orders of magnitude, Ni in olivine profiles provide the most-robust estimate of 19.1 ± 8.6 years from mixing-to-eruption. Our results provide new constraints on arc crustal differentiation processes, indicate rapid crustal transit timescales in agreement with a growing database of diffusion-based ascent timescales of basaltic magmas, and demonstrate significant, systematic deviation of diffusion timescales calculated from natural zonation of rare earth elements in clinopyroxene and Ti in clinopyroxene and olivine.

A simple random walk model explains the disruption process of hierarchical, Eccentric three-body systems

ABSTRACT We study the disruption process of hierarchical three-body systems with bodies of comparable mass. Such systems have long survival times that vary by orders of magnitude depending on the initial conditions. By comparing with three-body numerical integrations, we show that the evolution and disruption of such systems can be statistically described as a simple random walk process in the outer orbit’s energy, where the energy exchange per pericenter passage (step size) is calculated from the initial conditions. In our derivation of the step size, we use previous analytic results for parabolic encounters, and average over the (Kozai–Lidov) oscillations in orbital parameters, which are faster then the energy diffusion time-scale. While similar random walk models were studied before, this work differs in two manners: (a) this is the first time that the Kozai–Lidov averaged step size is derived from first principles and demonstrated to reproduce the statistical evolution of numerical ensembles without fitting parameters, and (b) it provides a characteristic lifetime, instead of answering the binary question (stable/unstable), set by case-specific criteria.

The Habitability of the Galactic Bulge.

We present a new investigation of the habitability of the Milky Way bulge, that expands previous studies on the Galactic Habitable Zone. We discuss existing knowledge on the abundance of planets in the bulge, metallicity and the possible frequency of rocky planets, orbital stability and encounters, and the possibility of planets around the central supermassive black hole. We focus on two aspects that can present substantial differences with respect to the environment in the disk: (i) the ionizing radiation environment, due to the presence of the central black hole and to the highest rate of supernovae explosions and (ii) the efficiency of putative lithopanspermia mechanism for the diffusion of life between stellar systems. We use analytical models of the star density in the bulge to provide estimates of the rate of catastrophic events and of the diffusion timescales for life over interstellar distances.

Competition between slow and fast regimes for extreme first passage times of diffusion * * The authors were supported by the National Science Foundation (Grant Nos. DMS-1944574, DMS-1814832, and DMS-1148230).

Many physical, chemical, and biological systems depend on the first passage time (FPT) of a diffusive searcher to a target. Typically, this FPT is much slower than the characteristic diffusion timescale. For example, this is the case if the target is small (the narrow escape problem) or if the searcher must escape a potential well. However, many systems depend on the first time a searcher finds the target out of a large group of searchers, which is the so-called extreme FPT. Since this extreme FPT vanishes in the limit of many searchers, the prohibitively slow FPTs of diffusive search can be negated by deploying enough searchers. However, the notion of ‘enough searchers’ is poorly understood. How can one determine if a system is in the slow regime (dominated by small targets or a deep potential, for example) or the fast regime (dominated by many searchers)? How can one estimate the extreme FPT in these different regimes? In this paper, we answer these questions by deriving conditions which ensure that a system is in either regime and finding approximations of the full distribution and all the moments of the extreme FPT in these regimes. Our analysis reveals the critical effect that initial searcher distribution and target reactivity can have on extreme FPTs.

Ions in an AC Electric Field: Strong Long-Range Repulsion between Oppositely Charged Surfaces.

Two oppositely charged surfaces separated by a dielectric medium attract each other. In contrast we observe a strong repulsion between two plates of a capacitor that is filled with an aqueous electrolyte upon application of an alternating potential difference between the plates. This long-range force increases with the ratio of diffusion coefficients of the ions in the medium and reaches a steady state after a few minutes, which is much larger than the millisecond timescale of diffusion across the narrow gap. The repulsive force, an order of magnitude stronger than the electrostatic attraction observed in the same setup in air, results from the increase in osmotic pressure as a consequence of the field-induced excess of cations and anions due to lateral transport from adjacent reservoirs.

Partial Viscosity Decoupling of Solute Solvation, Rotation, and Translation Dynamics in Lauric Acid/Menthol Deep Eutectic Solvent: Modulation of Dynamic Heterogeneity with Length Scale.

Deep eutectic solvents (DESs) are new-generation media that can be fine-tuned to have desired properties circumventing economic and environmental issues. Typically, these are ionic, and only recently, nonionic DESs, having interesting properties, are being explored. In this report, we examined the structure and dynamics of a nonionic lauric acid/menthol (LA/Men) DES through steady-state emission, solvation dynamics, time-resolved fluorescence anisotropy, and translational diffusion dynamics. The zero shift in the emission spectra of coumarin 153 (a solvatochromic dye) as a function of the excitation wavelength suggests that LA/Men DES is spatially homogenous. Decoupling (p = 0.63) of the average solvation time, ⟨τs⟩, from medium viscosity suggests the presence of mild dynamic heterogeneity in the system. Rotational time, ⟨τr⟩, which reflects the nature of the first solvation shell, shows little decoupling (p = 0.81), suggesting it to be fairly dynamically homogeneous at a shorter length scale. An Arrhenius-type analysis also proves that rotation is mainly controlled by medium viscosity. Translational diffusion time, ⟨τD⟩, which provides information at a larger length scale, is strongly decoupled from medium viscosity (p = 0.29). This indicates that at a larger length scale, the DES is quite dynamically heterogeneous. The slow component of solvation time, which is believed to originate at a larger length scale, correlates well with the translational diffusion timescale having similar activation energies. This suggests that their origin is same. Expectedly, for the long component of solvation time, the decoupling is quite strong (p = 0.30). Overall, our result demonstrates the structure and dynamics of the nonionic LA/Men DES, and the existence of length scale-dependent heterogeneity has been proposed.

Constraints on the pre-eruptive magmatic history of the Quaternary Laacher See volcano (Germany)

Ba zonation patterns in sanidine phenocrysts from mafic and intermediate phonolite and crystal-rich cumulates from the Laacher See volcano (12.9 ka) in western Germany document diffusion times suggestive of periodic recharge events throughout the magma reservoir’s entire lifespan of ~ 24 ky. Phenocrysts analysed from samples that formed late at the base of the compositionally zoned magma reservoir by mixing and mingling between a resident phonolite magma and recharging basanite show resorption and thin (2–10 μm) late-stage Ba-rich overgrowth. Short diffusion profiles across these boundaries give diffusion times of ~ 1.5–3 years at most, which are interpreted to be the maximum duration between the most recent recharge by the basanite and eruption. The lack of such late overgrowth in samples from other parts of the phonolite reservoir suggests that effect of this mixing and mingling was limited to the crystal-rich base. Sanidines in the cumulates, by contrast, are generally devoid of zoned crystals. Only rare cumulate crystals with resorbed outer boundaries and very thin overgrowths (a few microns) with very sharp compositional changes imply the remobilization of cumulates only months before eruption. Based on the diffusion timescales and storage temperatures obtained in a previous study, we present a genetic model for the conditions and timing of storage and (re-)activation of the magma system prior to the eruption of Laacher See, which is the largest volcanic event in Central Europe since the last glaciation.

On the prediction of hot spot induced ignition by the Livengood-Wu integral

The Livengood-Wu (L-W) integral is popularly used to predict the occurrence of homogeneous autoignition without heat and mass transport; yet it has not been used for non-homogeneous forced ignition which usually occurs in fire or explosion accidents. In this study, the forced ignition process caused by a hot spot in a flammable mixture is considered, and a method using the L-W integral is developed to predict the critical ignition temperature. In this method, the evolution of the temperature at the hot spot center is first obtained by solving the one-dimensional unsteady heat conduction equation. Then the L-W integral is evaluated based on the central temperature evolution and the corresponding homogenous ignition delay time. The critical ignition temperature is determined at the condition when the L-W integral reaches unity. The present method based on the L-W integral can reduce the computational cost by two to three orders compared to the transient simulation considering detailed chemistry and transport. Different fuels including methane, hydrogen, n-heptane, and dimethyl ether are considered. The predicted critical ignition temperature based on the L-W integral is compared with that predicted by one-dimensional transient simulations considering detailed chemistry and transport. For hydrogen, the present method based on the L-W integral can accurately predict the critical ignition temperature. However, for the other three fuels, the predicted critical temperature is lower than that from detailed simulation. Such under-prediction is interpreted by the ratio between the time scales for ignition and mass diffusion. Nevertheless, the method based on the L-W integral can give a conservative prediction of the critical ignition temperature, which is important for fire safety considerations. Besides, the effects of low temperature chemistry (LTC) on hot spot induced ignition are examined. It is found that the hot spot induced ignition is mainly controlled by high temperature chemistry.

Rapid assembly of high-Mg andesites and dacites by magma mixing at a continental arc stratovolcano

Abstract Studies of pre-eruptive processes at active volcanoes require precise petrochronological constraints if they are to contribute to hazard assessment during future eruption events. We present petrological and geochemical data and orthopyroxene diffusion time scales for samples from Late Pleistocene high-Mg andesite-dacite lavas (Mg# 53–69) at Ruapehu volcano, New Zealand, as a case study of rapid magma genesis and eruption at a continental arc stratovolcano. Assembly of Ruapehu high-Mg magmas involved the mixing of primitive magmas plus entrained mantle-equilibrated olivines with mid-crustal felsic mush bodies, yielding hybridized magmas with ubiquitous pyroxene reverse-zoning patterns. Orthopyroxene Fe-Mg interdiffusion time scales linked to quantitative crystal orientation data show that most lavas erupted <10 days after resumption of crystal growth following magma mixing events. The eruption of lavas within days of mixing events implies that pre-eruptive warnings may be correspondingly short.

Mass measurements during lymphocytic leukemia cell polyploidization decouple cell cycle- and cell size-dependent growth

Cell size is believed to influence cell growth and metabolism. Consistently, several studies have revealed that large cells have lower mass accumulation rates per unit mass (i.e., growth efficiency) than intermediate-sized cells in the same population. Size-dependent growth is commonly attributed to transport limitations, such as increased diffusion timescales and decreased surface-to-volume ratio. However, separating cell size- and cell cycle-dependent growth is challenging. To address this, we monitored growth efficiency of pseudodiploid mouse lymphocytic leukemia cells during normal proliferation and polyploidization. This was enabled by the development of large-channel suspended microchannel resonators that allow us to monitor buoyant mass of single cells ranging from 40 pg (small pseudodiploid cell) to over 4,000 pg, with a resolution ranging from ∼1% to ∼0.05%. We find that cell growth efficiency increases, plateaus, and then decreases as cell cycle proceeds. This growth behavior repeats with every endomitotic cycle as cells grow into polyploidy. Overall, growth efficiency changes 33% throughout the cell cycle. In contrast, increasing cell mass by over 100-fold during polyploidization did not change growth efficiency, indicating exponential growth. Consistently, growth efficiency remained constant when cell cycle was arrested in G2 Thus, cell cycle is a primary determinant of growth efficiency. As growth remains exponential over large size scales, our work finds no evidence for transport limitations that would decrease growth efficiency.

In situ XANES study of the influence of varying temperature and oxygen fugacity on iron oxidation state and coordination in a phonolitic melt

Iron oxidation state and environment in magmas affect their phase diagram and their properties, including viscosity and density, which determine magma mobility and eruptive potential. In turn, magma composition, pressure, temperature and oxygen fugacity affect iron oxidation state and coordination, potentially leading to complex feedbacks associated with magma ascent, degassing and eruption. While equilibrium experiments and models have led to a deep understanding of the role of iron in melts, our knowledge of the effects of disequilibrium processes on iron oxidation state and its structural role in lavas and magmas remains limited. Accordingly, we performed a series of dynamic disequilibrium experiments on a natural melt composition (a phonolite lava from Erebus volcano, Antarctica) at atmospheric pressure, in which oxygen fugacity and temperature were controlled and varied. During the experiments, we continuously measured iron oxidation and coordination using Fe K-edge dispersive X-ray Absorption Spectroscopy (XAS). We found that iron oxidation state changes in the phonolite melt are reversible and well reproduced by existing models. Changes in iron oxidation state are driven by joint diffusion of alkali cations and oxygen anions at magmatic temperatures (~ 1000 °C for Erebus phonolite). However, redox diffusion timescales are too slow for any significant oxygen exchange with the atmosphere at the lava/air interface or via air entrainment. Turning to iron coordination, while Fe2+ and Fe3+ are present mostly in an average five-fold coordination, complex coordination variations decoupled from redox changes were detected. The data suggest transitions between Fe3+ in four-fold and six-fold coordination prior to reduction or as a consequence of oxidation. This questions the possible implication of Fe coordination changes in triggering crystallisation of magnetite nanolites upon magma ascent, and, through such crystallisation events, in promoting magma explosivity.