# Simulant molten core-concrete interaction experiments in view of understanding Fukushima Daiichi Nuclear Power Station Cs-bearing particles generation mechanism

Hugo Laffolley, Christophe Journeau, Bernd Grambow

► To cite this version:

Hugo Laffolley, Christophe Journeau, Bernd Grambow. Simulant molten core-concrete interaction experiments in view of understanding Fukushima Daiichi Nuclear Power Station Cs-bearing particles generation mechanism. Scientific Reports, 2024, 14 (1), pp.6611. 10.1038/s41598-024-56972-9. hal-04515113

HAL Id: hal-04515113

https://hal.science/hal-04515113v1

Submitted on 14 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Simulant molten core-concrete interaction experiments in view of understanding Fukushima Daiichi Nuclear Power Station Cs-bearing particles generation mechanism

Hugo Laffolley<fnref n="1" />, Christophe Journeau<fnref n="1" />✉ & Bernd Grambow²

The Fukushima Daiichi accident resulted in the release of a novel form of radioactive Cs contamination into the environment, called Cs-bearing microparticles (CsMP). CsMPs constitute a substantial portion of the radioactive pollution near the nuclear power station and traveled beyond several hundred kilometers. Extensive characterization of the CsMPs revealed an amorphous silica matrix, along with Cs and other minor or trace elements such as Fe and Zn. This study explores the unclear generation mechanism of CsMPs by conducting experimental molten core concrete interactions (MCCI) as a source of Si and analyzing the resultant aerosols. The findings demonstrate that MCCI is in capacity to produce spherical submicronic and micronic particles, primarily composed of amorphous silica and incorporating elements akin to CsMPs. A humid atmosphere is found to favour an even closer chemical composition. Examination of the internal structure of the synthesized particles unveils pores and numerous crystalline nanoinclusions possibly serving as nucleation sites for CsMP formation through the condensation of Si-rich vapors.

**Keywords** CsMP, Fukushima Daiichi, Radioactive aerosols

On March 11, 2011, the Fukushima Daiichi Nuclear Power Station (1F) experienced three reactor core melt-downs, catastrophic events caused by the loss of power following the partial flooding of the 1F site<fnref n="1" />. This incident, ranked the highest level on the International Nuclear Event Scale, resulted in the release of radioactive materials into the environment a few days after the initiating events².

Subsequent airborne monitoring of radiocesium identified a novel form of radioactive particle, named Cs-bearing microparticles (CsMP), the first time at Tsukuba (170 km southwest of 1F) in the night of March 14th to 15th³. CsMPs are mostly made of an amorphous silica matrix containing various elements, including radioactive cesium, and are categorized into type-A (micrometric-sized particles, originating from Unit 2 or Unit 3, high specific radioactivity) and type-B (tens to hundreds of micrometers, originating from Unit 1, low specific radioactivity)⁴⁻⁸. Type-A CsMPs tend to be spherical particles and could travel on hundreds of kilometers due to their micrometric size⁹,¹⁰.

Determining the proportion of Cs radioactivity attributable to CsMPs in a given soil sample proves challenging. Nevertheless, it is proposed that in close proximity to 1F, CsMPs may locally contribute significantly to Cs radioactivity¹¹,¹². The cesium encapsulated within CsMPs has a limited transferability to the environment due to the low solubility of the silica matrix in water¹³⁻¹⁵. Combining their minute size and low solubility, aspirable CsMPs may cause prolonged radiation doses locally, in cases of deep penetration into an individual's respiratory system, surpassing the biological half-life of Cs in its ionic form. Moreover, the presence of this new type

<footnote n="1">¹CEA, DES, IRESNE, DTN, Severe Accident Experimental Laboratory, Cadarache, 13108 St-Paul-lez-Durance, France. ²SUBATECH (IMT Atlantique, CNRS-IN2P3, University of Nantes), 44307 Nantes, France. ³email: christophe.journeau@cea.fr</footnote>
of Cesium-bearing aerosols will affect the fission product chemistry in the containment while it is not currently modelled in severe accidents such as ASTEC¹⁶ or MAAP¹⁷.

Despite extensive characterization of type-A CsMPs in the years following their discovery, a consensus on the detailed mechanism of their generation remains elusive. It is evident that the source must entail a sufficient amount of silicon to form the particles' matrix, with other elements such as Fe, Zn, or Cs integrated into the particles.

One proposed origin of CsMPs suggests an interaction between Cs (as Cs₂O) and up to 1% silicon present in the stainless steel of the internals and the reactor pressure vessel (RPV) at high temperatures. Experimental and thermodynamic evidences support the formation of oxidic phases rich in Fe, Si, and Cs, which could later evolve into type-A CsMPs through remelting¹⁸. This hypothesis has the benefit of being valid for both in-vessel and ex-vessel phenomena. Another proposed mechanism suggests a gas-solid interaction between CsOH, likely formed in humid atmospheres¹⁹, and insulation materials such as calcium silicate at high temperatures²⁰. Experimental work demonstrated the formation of complex chemical species, such as CsAlSiO₄, potentially leading to CsMPs. Additionally, it has been theorized that during the Unit 3 hydrogen explosion, the safety gas treatment system, equipped with high-efficiency particulate air filters, may have atomized to form type-A CsMPs²¹,²². However, these hypotheses lack experimental validation as very similar particles could not be, to our best knowledge, synthesized.

Consequently, a fourth hypothesis, mentioned several times in the literature, suggests that type-A CsMPs may result from molten core-concrete interaction (MCCI) following RPV failure²³,²⁴. This theory relies on the potential interaction at high temperatures, between the zirconium fraction of sub-oxidized corium and the silica present in the concrete²⁵, particularly the basaltic concrete of 1F rich in SiO₂. This interaction could produce silica-rich gases that condense into micrometric spherical particles.

Preliminary investigations have demonstrated the feasibility of generating microparticles resembling CsMPs in shape and size, along with a comparable chemical composition, through small-scale MCCI experiments (about 20 g of reagents) in the VITI furnace at the CEA Cadarache²⁶. Simultaneously, another research team observed similar results while simulating MCCI using a distinctly different experimental setup²⁷. Consequently, additional experiments in the VITI furnace, coupled with more in-depth analyses, have been conducted to reinforce the theoretical framework and gain a nuanced understanding of the generation mechanism of type-A particles.

## Results

### Objective

The primary objective of this study is to produce and gather aerosols resulting from MCCI, in conditions approaching those estimated during the accident of 1F, in order to examine their morphology and chemical composition. These particles will be compared with data accumulated over the last decade concerning type-A CsMP, aiming to enhance our understanding of the generation mechanism of type-A CsMP. To achieve this, a small-scale experimental setup has been devised to induce the reaction of a mixture of concrete and prototypic corium, at a temperature about 2000-2200 °C, defined by thermodynamic calculations²⁸, under atmospheric pressure, and within a controlled atmosphere (pure nitrogen or a mixture of steam and nitrogen). This experiment is based on the fact that at high temperature, Zr (from fuel cladding) is easily oxidised in an exothermic reaction, in particular by the silica contained in the basaltic concrete to form silicon-rich gas, and consequently particles²⁵. The aerosols generated are subsequently collected on a filter for post-experimental analyses.

### Size and morphology of the generated aerosols

The aerosols generated during our experiments were subjected to airborne particle mass size distribution analysis using a low-pressure impactor. The geometric diameter was determined based on the aerodynamic diameter, assuming spherical particles with a density similar to that of silica glass of low concentration in alkali species (2.25 kg m⁻³)²⁹. The 250 nm cutoff diameter stage collected the maximum mass of aerosols, corresponding to a geometric diameter of 170 nm (Figure S1). In comparison, the mass collected on the 1.6 µm and 2.5 µm cutoff stages, corresponding to geometric diameters of 1.1 µm and 1.7 µm respectively, was 15 times lower than that of the stage with the maximum mass.

A comparison of the aforementioned measurements was performed using Scanning Electron Microscopy (SEM) images. To count a large number of particles, the CELLPOSE algorithm, based on artificial intelligence (AI), was employed³⁰. Two SEM images taken at a magnification of approximately 18,000 and 10,000 allowed the identification of 354 and 821 regions of interest (ROI) respectively, which correspond to potential submicrometric-sized particles (Figure S2). The median diameter was determined by considering both the maximum Feret diameter and the mean of the major and minor axes of an ellipse that best fits the projected area of each particle without exceeding its boundaries. For the two images, the measured Feret diameters were found to be 208 nm and 211 nm, respectively. The ellipse method diameters were measured to be 165 nm and 166 nm, respectively. These measurements, particularly those related to the ellipse method diameter, are consistent with the size distribution obtained through impaction analysis.

SEM observations on the prepared samples revealed the presence of micrometric-sized particles resembling CsMPs. Their occurrence is significantly lower than that of submicrometric-sized particles, but their size makes them easily distinguishable even at low magnification levels on the SEM. Considering the cascade impactor stages 11 to 14, which have aerodynamic cutoff diameters of 1.6 µm, 2.5 µm, 3.6 µm, and 5.3 µm respectively, a rough estimation of the rate of generation of micrometric-sized particles can be achieved. Stage 15, with the highest cutoff diameter, is not considered, as large dust pollution can be collected. For stages 11 to 14, the corresponding geometric cutoff diameters are respectively 1.1 µm, 1.7 µm, 2.4 µm, and 3.5 µm, considering a particle material density of 2.25 kg m⁻³ and a spherical shape. By making the hypothesis that every particle collected on each
stage has a diameter equal to the geometric cutoff diameter, and considering the respective mass gain of these stages (Figure S1) during the 25 s collection time, the number of micrometric-sized particles collected on stages 11 to 14 can be estimated to be around $2 \times 10^8$ particles or $9 \times 10^6$ particles.s⁻¹.

The detailed observations of the micrometric-sized particles led to the initial identification of two distinct categories. Firstly, the majority of these particles have a nearly perfect spherical shape, with potential slight depositions on their surfaces, possibly resulting from the collection of other smaller aerosols (Fig. 1a,b). Within this category, the largest particle observed had a diameter of 6.4 µm (Fig. 1a). In general, most of the observed micrometric-sized particles have diameters ranging from about 2–3 µm. The second category consists of micro-metric-sized particles that appear to be aggregates of submicrometric-sized particles fused together (Fig. 1c). It is highly probable that these aggregates consist of micrometric-sized particles trapped beneath a layer of submicrometric-sized particles that remained attached during the water dispersion step of the sample preparation process, as some of them are only partly covered (Fig. 1d). This assumption is further supported by the fact that these particles typically appear larger than the spherical micrometric-sized particles found in the first category.

The efforts to optimize the experiments led to carry out experiments with a less confined reaction chamber, called open test section and using less fragile ceramic parts. The initial confined test section has a reaction volume of roughly 0.25 L (Figure S3a), while the whole volume of the vessel (70 L) can be considered as the reaction chamber for the open test section (Figure S3b). As a result, a new type of micrometric-sized particle variation was observed (Fig. 1e–g). These micrometric-sized particles will be termed “bumpy” particles hereinafter, owing to their distinct surface appearance. While predominantly having a global spherical shape, they possess a diverse range of surface roughness characteristics. In addition, the evolved experimental test section also resulted in the generation of particles closely resembling those depicted in the Fig. 1d, but entirely covered of submicrometric-sized particles (Fig. 1h). Upon closer examination of the surface, it appears that these particles consist

<figure><img src="image_3.png" /><figcaption>**Figure 1.** SEM images of synthesised microparticles. Images (a–d) have been captured at an acceleration voltage of 15 kV and are particles resulting from the test section with a closed reaction chamber. (a,b) are spherical particles with a smooth surface, (c) looks like an agglomeration of submicrometric-sized particles and (d) is a smooth surface particle half-covered with submicrometric-sized particles. Images (e–i) were captured at an acceleration voltage of 5 kV. These particles are resulting from the test section without reaction chamber and do not have a smooth surface. (e–g) are bumpy particles with a varying intensity and size of bumps. (h) is a particle with a distinctive layer stuck on the surface, and (i) a zoom on this particle showing with details the surface and the bonds between the submicrometric-sized particles.</figcaption></figure>
of micrometric-sized particles with attached submicrometric-sized particles on the surface, bound through
mechanical bonds. These bonds appear to be the result of the solidification of viscous wires (Fig. 1i).

Global chemical composition of synthesized CsMPs

In a previous study, we provided an analysis of the chemical composition of the particles synthesized via MCCI under an inert gas (N₂) atmosphere²⁶. The results demonstrated that the synthesized micrometer-sized particles exhibited a composition quite similar to CsMPs found in Fukushima. This conclusion was drawn through an analysis of their energy-dispersive X-ray (EDX) spectra, which revealed Si and O as the major elements, along with the presence of elements such as Na, Mg, Al, K, Ca, Fe, and Cs. It was observed, however, that the Fe and Cs concentrations in the synthesized particles appeared to be much lower than that found in real CsMPs. Additionally, a notable absence of Zn, an important element, was also identified in the synthesized particles. The source of Zn is yet to be determined as the addition of ZnO in a few experiments did not result in a detectable amount of Zn in the observed synthesized particles.

In this study, the experiments were conducted using the same prototypic corium and process than our previ-
ous study<sup>26</sup>. Some trials were partially conducted under a reactive atmosphere to investigate the potential impact
of steam on particle generation and composition, because it is estimated that the atmosphere inside the PCV
at the moment of the lower head failure of the concerned reactor was particularly rich in steam. The reactive
atmosphere was injected in the test section at the time of collection of the aerosols and the composition consisted
of 20 % (partial pressure) steam and 80 % nitrogen. Additionally, further experiments were performed in a pure
nitrogen atmosphere to expand the number of samples analyzed and explore the influence of slight variations
in the prototypic corium composition.

Seven particles from a prototypic corium of a higher stainless steel content than the previous study and gener-
ated in a pure nitrogen atmosphere were analyzed in various experiments (composition in Table S1 and Table S2)
(Fig. 2a). As expected, the overall composition appears quite similar to the previous findings, but significant
variations in specific element concentrations can be observed among the particles. Apart from Si and O from the
matrix, Na, Al, S, and Ca are consistently present in every particle. However, there are observable concentration
differences, particularly for Ca. Both Fe and Zn concentrations are very low or undetectable in all particles. The
concentrations of Mg, Cl, K, Sn, and Cs vary from one particle to another. In particular, spectra 10 and 24 of the
Fig. 2a show higher concentrations of K and Na, suggesting that some factor may have favored the enrichment
of alkali elements in these specific cases.

Ten particles from a single trial conducted under a mixed atmosphere were characterized using SEM-EDX (Fig. 2b). The elements found in these particles are generally similar to the experiments done under nitrogen atmosphere, except for Cl that is not detected. The presence of steam seems to have an influence on alkali elements as well, as Na, K, and Cs have been consistently measured in every particle and in significantly higher amounts. Fe was also detected in all particles. With the exception of spectrum 16 of the Fig. 2b, Ca was found to be present in very low concentrations in each particle. It is worth noting that the compositions of these particles might show more homogeneity compared to previous experiments, likely due to the fact that they originate from the same trial.

Secondary ion mass spectrometry (SIMS) analyses were conducted on particles generated in a pure nitrogen atmosphere to measure the concentration of the different minor elements and especially detect the presence of B and trace elements such as U. The identification of the particles was made considering Si and Cs clusters as being

<figure><img src="image_4.png" /><figcaption>**Figure 2.** (a) SEM-EDX spectra of 7 micrometric-sized particles synthesized over several experiments in dry atmosphere. (b) SEM-EDX spectra of 10 particles synthesized in one trial in humid atmosphere.</figcaption></figure>
particles (Fig. 3a). Using relative sensitivity factors, the relative composition of 108 particles was determined³¹. B was measured in significant concentrations, with a median B/Si atomic ratio value of 8 × 10⁻². A specific study³² reported that CsMPs contains 1518–6733 mg kg⁻¹ of ¹⁰⁺¹¹B. For a CsMP composed of 60 % SiO₂ by mass, the atomic ratio of B/Si thus varies from 1.4 × 10⁻² to 6.2 × 10⁻². It was found that the synthesized particles are slightly richer in B compared to real CsMPs. Among the 108 analyzed synthesized particles, 25 of them contained some U at low concentrations.

The synthesized particles were subjected to a comprehensive analysis of the relative atomic concentrations of Na, Mg, Al, K, Ca, Fe, Rb, Sn, Cs, Ba, and U relatively to Si. Figure 3b illustrate the statistical processing of the atomic concentration ratios for the 108 synthesized particles analyzed, juxtaposed with a parallel statistical analysis of the composition of authentic CsMPs data from the literature⁵,⁸,³³,³⁴. The investigation by Hagiwara et al.³⁴ presents the most exhaustive composition analysis over 37 particles, and Fig. 3b illustrates the commonly measured elements. Upon comparing the median concentration values between real CsMPs and synthesized particles, Fig. 3b reveals that Na, Mg, Al, Sn, Cs, and U show a Si relative concentration quite similar, with a difference of less than one order of magnitude. However, for K, Ca, Fe, and Ba, the median concentration diverges more significantly, exceeding one order of magnitude. Focusing specifically on the alkaline species concentration, the synthesized particles show a comparable concentration for Na. The concentration in Cs is almost identical, but marginally lower than that of the real particles from Okumura et al.⁴⁰, given the utilization of composition

<figure><img src="image_5.png" /><figcaption>**Figure 3.** (a) SIMS ionic images of several clusters of ¹¹B, ²⁹Si, ¹³³Cs and ²³⁸U identified as micrometric-sized particles. Potential CsMPs were identified using only Si and Cs ionic images. ²⁹Si has been chosen here due to the saturation of the detector caused by ²⁸Si. All of the identified particles contains B, and one particle contains U in this image. (b) Comparison of the atomic ratio values of Na, Mg, Al, K, Ca, Fe, Rb, Sn, Cs, Ba and U on Si, between 108 synthesized particles in a dry atmosphere and data from real CsMP. Hagiwara et al.³⁴ studied 37 particles exhaustively, Okumura et al.⁴⁰ provides alkaline species concentration of 5 particles and Kogure et al.⁵, Yamaguchi et al.⁸ and Imoto et al.³³ provide the concentration of major elements, all of them are not detailed here (7 particles). It should be noted that studies prior to Okumura et al.⁴⁰ did not identify Na due to overlap with Zn. Ratios without minimum value infer that some particles have a concentration equal to zero for this specific element.</figcaption></figure>
data from the particle centers for those having a Cs concentration gradient. For K, the concentration in the synthesized particles is higher in the same extent than observed in the previous study. Regarding Rb, the concentration in the synthesized particles is slightly higher; however, it is noteworthy that Rb is introduced as a distinct component in the simulant corium mixture, and its concentration can be readily adjusted. Considering the data from Kogure et al.⁵, Yamaguchi et al., and Imoto et al.³³, the concentration is consistent for Sn and Cs. Nonetheless, the findings align with the preceding data for K and Fe. It is important to note that these earlier data were collected prior to the confirmation of Na's presence in CsMPs. Consequently, the Na concentration was initially misinterpreted as Zn additional concentration due to peak overlap in EDX analysis. Overall, Na, Sn, and Cs show highly similar concentrations, while K and Ca concentrations are higher in synthesized particles, and Fe concentration is considerably lower.

## Internal structure of synthesized CsMP

The internal structure of three synthesized CsMPs have been characterized through the preparation of thin sections using focused ion beam (FIB) techniques. It has been confirmed via diffraction analysis on a transmission electron microscope (TEM) that the silica matrix of the synthesized particles consists of an amorphous phase (Fig. 4d, 5c,f). Nanometer-sized crystalline inclusions and pores have been observed within the particles. The largest particle observed (OPT-1), having a diameter about 5 µm, does not contain any pore, but a high concentration of crystalline inclusions is present (Fig. 4a–c). The first type of crystalline inclusions mainly consists of round precipitates, each with a diameter smaller than 40 nm. These precipitates are composed of Fe, Cr, Mn, Si, and Zr and appear to be uniformly distributed throughout the particle (Fig. 4e,f). Another type of crystalline precipitates is distinguishable around the periphery of the particle’s thin section, at depths up to 200–300 nm from the particle’s surface. These inclusions have roughly the same size than the other inclusions but have an elongated projected shape and are predominantly composed of Zr without any visible impoverishment in oxygen content compared to the matrix (Fig. 4g).

The particle OPT-2 has a diameter of 1.4 µm and contains pores and precipitates (Fig. 5a and S4a). The porosity is mostly concentrated in the upper region of the thin section, occupying an area of an approximate diameter of 500 nm. The pores display numerous irregular shapes and varying sizes, with their maximum Feret size reaching around 50 nm. This particle also has an amorphous matrix and contains only a small number of nanometer-sized crystalline inclusions such as in OPT-1, but contains two significant crystalline inclusions occupying an area of approximately 400 × 400 nm (Fig. 5c).

The chemical composition analysis of these particles has proven challenging due to uncertainties arising from the use of a focused beam in scanning transmission electron microscopy (STEM) configuration during EDX measurements. This is known to alter the composition and structure of the materials being analyzed, particularly affecting alkali species, which have a propensity to migrate within the matrix. Nevertheless, the elemental map obtained from the upper part of the left precipitate suggests that this specific inclusion is abundant in chlorine (Cl) and deficient in alkali species, especially potassium (K) and sodium (Na) (Fig. 5b).

A third particle, designated as OPT-3, has been subjected to internal structure analysis. Its diameter is approximately 700 nm, which is smaller than the typical CsMPs size. Similar to the previous particle, OPT-3 also contains pores, with the largest pore having a diameter of about 100 nm (Fig. 5d and S4b). These pores are predominantly located in the upper left part of the particle. Additionally, OPT-3 contains spherical precipitates that are rich in Fe, Cr, Mn, Si, and Zr, similar to what was observed in OPT-1. Furthermore, some other nanometer-sized inclusions, enriched almost exclusively in Sn, have been identified near the surface of the particle (Fig. 5e,f). It is worth noting that the thin section of this specific particle was the thinnest among the three particles, measuring approximately 80 nm, while the others were 117 nm and 153 nm thick for OPT-2 and OPT-1 respectively (Figure S5a–c). This lower thickness probably increased this thin section structural modifications under the electron beam during the analysis process (Figure S6).

During the thin section preparation process, several submicrometer-sized particles were inadvertently trapped in the platinum (Pt) material. Due to the likelihood that the sections were not made near the equator of these particles, the observations obtained may not be fully representative from one particle to another. However, it has been observed that the majority of these trapped particles contain at least one nanometer-sized precipitate, similar in nature to those observed in OPT-1 and OPT-3.

## Discussion

It has been reported that aerosols resulting from MCCI typically have a diameter exceeding 0.4 µm, with a predominant size range of several micrometers³⁵. While the largest particles may exceed 10 µm, it is suggested that these particles primarily arise from mechanical generation processes, such as matter entrainment, rather than nucleation and condensation, as explored in our study³⁶,³⁷. Type-A CsMP have always been observed having a diameter smaller than 10 µm, which agrees with a formation by MCCI regarding their size distribution. A specific investigation has reported the existence of type-A CsMPs with a diameter less than 0.5 µm³⁸, implying that the isolation of CsMPs based on the amount of Cs they contain (i.e. radioactivity) remains the most reliable method for detecting the most harmful particles in terms of dose, as well as the largest particles. Given that the size distribution of the particles synthesized in this study is slightly smaller than the typical size distribution for MCCI aerosols, it is hypothesized that the rapid cooling induced by gas injection forced convection in the test section, coupled with a high prevalence of nucleation sites caused by the confined reaction volume, favored the creation of numerous smaller particles rather than promoting growth through continuous condensation, leading to the formation of larger particles.

The prevalence of the submicrometric-sized particles, together with the observation of the bumpy micrometric-sized particles has been beneficial in the understanding of the generation process. The observation of
<figure><img src="image_6.png" /><figcaption>**Figure 4.** (a) STEM dark-field image of the particle OPT-1 showing many crystalline inclusions (bright spots), more visible on (b,c) enlarged images. (d) is diffraction images of the area (c), showing a mostly amorphous phase and a few diffraction patterns from the crystalline nano-inclusions pointed out by red arrows. (e) shows the STEM-EDX elemental maps of the area (c) for interest elements constituting the crystalline inclusions of Fe, Cr, Si, Mn and Zr. (f) is the EDX spectrum of the same area, knowing that Cu and Ga have been brought by the preparation process or the sample holder, and that C has also mostly been brought by C coating. (g) is the elemental map of a smaller area near the top boundary of the particle showing inclusions rich in Zr.</figcaption></figure>

the inside of the particles revealed the existence of nanoinclusions such as inside numerous type-A CsMPs⁵,²³,³⁹. It also seems that among the submicrometric-sized particles that were unintentionally included in the thin sections, most of them contain at least a nanometric sized inclusion. It is therefore thought that these nano-metric inclusions might serve as nucleation sites for silicon rich aerosols generated by the MCCI to condense and create silica particles (Fig. 6). In the case of the experiments, many nucleation sites are available in a really confined volume, and the forced convection enhances this phenomenon by renewing the atmosphere above the corium. This favours the nucleation of many new particles by the condensable gases rich in Si. Particles with a short residence time will be extracted as submicrometric-sized particles, those with a longer residence time might grow by additional condensation and more likely coagulate with other submicrometric-sized particles to form micrometric particles. The conditions inside the containment at 1F were significantly different, in terms of volume, gas velocity and temperature gradient. It is supposed that nucleation sites (i.e. nanoparticles) were sparsely available, favouring a growth of microparticles by condensation mostly. It is worth noting that most SEM images of real CsMPs in the literature were captured at an acceleration voltage of 20 kV and with a back-scattered
<figure><img src="image_7.png" /><figcaption>**Figure 5.** (a) TEM image of the particle OPT-2 with a STEM dark-field enlargement of a crystalline inclusion area. (b) elemental map of the enlarged area of OPT-2 showing an enrichment in Cl and Na. (c) diffraction pattern of the same area showing a crystalline pattern. The detector was saturated rendering a poor contrast. (d) TEM image of the particle OPT-3. (e) Elemental map of an enlarged area of the particle OPT-3 showing inclusions similar to those observed in particle OPT-1, as well as Sn-rich inclusions near the boundary. (f) Diffraction pattern of the particle OPT-3 showing an amorphous matrix and diffraction patterns pointed out by the red arrows.</figcaption></figure>

electron detector, to proceed to the chemical analysis in the meantime. The few low voltage high-resolution images captured of spherical type-A CsMPs tend to show a not so regular surface⁸.

In the evolution of the test section design, particles are cooled down much faster due to the open reaction chamber that provides less thermal confinement into the outlet of the vessel (Figure S3). It is supposed that this is the reason for the existence of bumpy particles. The micrometric-sized particles are probably the result of the coagulation of submicrometric-sized particles at a viscous state, nucleated on the nanometric size mixture of stainless steel elements and zirconium inclusions. The bumpy particles are likely to be witnesses of this coagulation process, by showing the progressing coagulation, interrupted by a fast cooling during the extraction of the aerosols from the vessel. Experimentally, the aerosols are pushed out by a gas flow that forces their cooling by a high velocity rising, compared to what could have occurred in the containment of 1F. This forced cooling is enhanced during the experiments carried out with the open test section and increase the chance of generating "bumpy" particles.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Suggested mechanism of generation for the synthesised particles in relation to the observations on the particles size distribution and their internal structure.</figcaption></figure>

The chemical composition of the nano-inclusions is not as diverse and entirely identical to the observed inclusions in real CsMPs. In the actual accident scenario, nanoparticle sources have a much greater diversity than in our experimental setup. For instance, (U,Pu) oxides are likely the result of the mechanical deterioration of exposed fuel rods, Fe oxides might arise from the mechanical removal of the stainless steel oxide layer, and fission products inclusions may originate from the release of fuel-derived materials. However, the particle OPT-2 contains inclusions rich in Na and Cl (that could be considered as a sodium chloride precipitate), which is really similar to a precipitate of cesium chloride (by extension of alkaline properties) that has been observed in a real CsMP²³. In the same study, the authors also identify a zinc chloride inclusion that may have been misidentified and could, in reality, be a sodium chloride precipitate, given that Na has been confirmed to exist in CsMPs later and was misidentified as Zn initially⁴⁰. Furthermore, in another study, a particle contain an inclusion made of Fe and Cr, resembling those observed in synthesized particles and another particle an inclusion rich in Cl¹⁵.

The absence of Ni, a major constituent of stainless steel, in the nanoinclusions might seem surprising. The stainless steel (SS) utilized in these experiments is SS301, chosen due to a supply issue with SS304. Despite this substitution, the composition of SS301 closely resembles that of SS304, containing up to 2wt% of Mn and an equal amount of Si, the latter being slightly higher than in SS304. Mn and Si rank as the fourth most abundant elements in this alloy, following Fe, Cr, and Ni. Consequently, the presence of Mn and Si in these inclusions is unsurprising. Regarding Ni, studies on stainless steel welding have revealed that resulting fumes contain approximately 10 times less Ni than Fe or Cr, despite being only 5 to 6 times less concentrated than Fe and 2 times less concentrated than Cr⁴¹. Furthermore, SS301 contains a slightly lower amount of Ni in its composition compared to SS304.

The chemical analysis of the synthesized particles reveals differences, especially in the case of Ca, when compared to real CsMPs. Over the years of characterization, it has been observed that the composition of CsMPs may vary from one to another⁴. Ca has been detected in various particles from different locations or substrates, and it has been characterized using multiple techniques³⁴,⁴²,⁴³. In this study, experiments conducted in a humid atmosphere appear to generate particles with a low concentration of Ca, indicating that the atmosphere could play a crucial role in influencing the particle composition for certain elements.

SIMS analysis indicates a similarity in Cs concentration between the synthesized particles and the CsMPs, in contrast to observations from SEM-EDX analyses on synthesized particles. This discrepancy may be attributed to the particle identification process. Real CsMPs are isolated from samples, such as soil, utilizing an imaging plate to pinpoint highly radioactive spots. Consequently, the prioritized identification of the most radioactive CsMPs occurs, considering that Cs concentration can vary by approximately one order of magnitude from one particle to another, as reported by a previous study⁴⁰. In the SIMS analyses conducted for this study, particles were also identified by instructing the device to locate clusters rich in Si and Cs. Conversely, SEM-EDX relied solely on morphological features observed in images for particle identification, and the majority had significantly lower Cs concentrations.

It is suggested that the solubility of alkali species might have influenced their concentration for humid atmosphere experiments. In the presence of a highly humid atmosphere, Cs mainly forms CsOH, a highly reactive species¹⁹. The other alkali species tends to behave in the same way. This highly reactive chemical forms could have enhance the inclusions of alkali species into the particles. Thermodynamic calculations did not show a significative influence of the atmosphere in the concentration of these specific elements as well as Fe and Ca. However, it has been noticed that the data from the Thermodynamics of Advanced Fuels-International Database (version 13) might not depict all the real interactions that could occur.
These experiments represent a promising insight into understanding the potential of MCCI to easily generate particles closely resembling type-A CsMPs in terms of shape, size, structure and chemical composition. An exact synthesis of type-A CsMPs is hard due to their relative composition variation and due to the myriads of reaction conditions that occur simultaneously in a failed reactor, considering thermal and pressure gradients, gradients in reactive materials compositions, large compositional variations at different locations and only partial mixing of reactive fluids and with rapid variations over time. Regarding the chemical composition, refining certain parameters could strengthen our conclusions, particularly if more controlled experiments could be conducted in a humid atmosphere and if Zn origin could be identified. Some specific characteristics of the experimental setup are supposed to influence the particle formation, especially their size and shape. An upscaled experiment could lead to a slightly different observation due to a different aerodynamic condition and reaction volume. It is suggested that the particle result from condensation on nano-particles, then growth and agglomeration to form micrometric size particles. Further characterizations of environmental samples from the Fukushima area could concentrate on determining the size distribution of CsMPs compared to experimentally obtained particles.

These findings have shown that an experimental set up designed to simulate some aspects of MCCI in a laboratory has produced Cs-enriched silica microparticles which are in many aspects similar to CsMP observed in Fukushima. Although detecting and identifying this kind of particles during an ongoing accident pose inherent challenges, their recognition could serve as an early indication that MCCI has initiated, prompting a targeted response from crisis management centers. More research is necessary to confirm that the CsMP observed can be confirmed as an indicator of MCCI.

Moreover, considering CsMPs in severe accident codes could prove beneficial for accurately assessing the source term released into the environment. The Fukushima Daiichi incident has shown the substantial contribution of CsMPs to radioactive pollution.

Subsequently, the scientific community may seek to evaluate the impact of particulate-trapped Cs in comparison to ionic Cs, with regard to both the delivered dose to individuals and the transfer of radioactivity to the environment. In the context of designing future reactors, particularly in determining concrete composition, a critical decision must be made regarding whether CsMPs should be promoted or strictly avoided in the event of an inevitable release of radionuclides following the initiation of MCCI.

## Methods

### Experimental setup

The experimental setup comprises a furnace serving as a chemical reactor to induce MCCI and a sampling line for collecting and analyzing aerosols. The furnace, a water-cooled stainless steel vessel with a volume of approximately 70 L, contains the test section and a water-cooled inductor with four turns. A 50 kW electric generator (operating at about 110 kHz) supplies power to the inductor, establishing electromagnetic coupling with a conductive component (susceptor) to generate heat (Figures S7 and S8).

The test section is composed of an yttria fully stabilized zirconia (Y-FSZ) crucible, a tungsten cylinder (susceptor) for electromagnetic coupling, and a thermal shield made of porous graphite. To address the reactivity of tungsten and graphite with a humid atmosphere at high temperatures, a test section with a reaction chamber made of Y-FSZ has been developed to prevent the atmosphere from reaching the tungsten and graphite. This configuration, providing only moderate airtightness and limited mechanical resistance to thermal gradients, could not satisfy the requirements. An alternative version without the Y-FSZ reaction chamber has been developed, the entire vessel being the reaction chamber, only allows experiments in a dry atmosphere (Figure S3).

The sampling line extracts aerosols from the top of the vessel, dilutes the flux with nitrogen (at a 1:6 ratio) to reduce particle concentration, and splits the diluted flux into three streams. One stream passes through a particle counter (Pegasor*PPS-M), another through a HEPA filter (0.2 µm pores) for particle collection, and the third through a low-pressure impactor (Dekati DLPI+) to measure the inertial size distribution of particles. The sampling line is heated to over 100 °C to prevent water condensation.

Temperature measurements employ two bichromatic video pyrometers, one focused on the corium (from the top of the vessel) and another on the crucible or reaction chamber (depending on the test section version). However, measurements above 1700–1800 °C are deemed unreliable due to the substantial aerosol generation perturbing the optical path. Calibration tests with an empty crucible were conducted before real experiments to estimate crucible temperature based on power injection.

### Test protocol

The experiments involve the interaction between sub-oxidized corium and basaltic concrete resembling that of 1F. Corium compositions are determined using the results of the BSAF project⁴. The corium is prototypic as it utilizes depleted uranium oxide and does not incorporate any other radioactive isotopes (including PF), but only as stable isotopes or excludes them if non-existent in their stable form (Tables S1 and S2). This composition, rich in UO₂ and sub-oxidized Zr is typical of an early accidental phase.

The concrete has been made with materials from Japan close to the originals constituting 1F's concrete. The composition has been analyzed using inductively coupled plasma atomic emission spectroscopy (ICP-AES) and shows that the concrete is made from 58.7wt% of SiO₂, 9.9% of CaO, 8.4% of Al₂O₃, 5.8% of Fe₂O₃, 2.4% of Na₂O, 1.5% of K₂O and 1.3% of MgO. The concrete also contains many other minor or traces elements (as oxides) such as Ti, Pb, S, Ba, etc. The concentration of water and carbonates in the concrete could not be measured. Basaltic concretes usually have a low concentration in carbonates. The water content is estimated to be around 6–8wt%.

The mixture is loaded into the crucible as two configurations: a homogeneous mixture with all elements mixed as powders and a heterogeneous mixture where the sub-oxidized corium and concrete are vertically separated by a stainless steel sheet (no visible influence on the aerosols). After installing the test section in the vessel and
loading the mixture, the vessel is hermetically sealed and subjected to medium to high vacuum. It is then filled with nitrogen to atmospheric pressure. The heating process starts slowly to limit powder projection caused by rapid gas expansion (reaching 1200–1300 °C in about 30 min). Subsequently, the temperature is abruptly increased to minimize the interaction time until the target MCCI temperature of 2000–2200 °C is reached. Aerosol collection begins after reaching the target temperature by opening the sampling line and injecting a gas flow from the bottom of the test section to carry the aerosol to the top and prevent pressure decrease in the vessel. The injected gas flow may consist entirely of nitrogen or be a humid gas composed of mixed nitrogen and steam.

## SEM-EDX samples preparation and observations

Particles under study are collected on a HEPA filter. A portion of less than 1 cm² is taken from the filter with the deposit, transferred into 4 mL of distilled water, and vigorously agitated for about 1 min to disperse particles and dissolve soluble aerosols. The suspension is filtered through a 0.2 µm pore membrane, and the membrane is let to dry for a minimum of 24 h at room temperature. Subsequently, the membrane with the low-solubility particle deposit is transferred onto a stainless steel SEM stub, affixed with double-sided carbon tape, coated with carbon, and supplemented with a piece of aluminum tape to enhance electrical discharge.

Analyses have been conducted on a Zeiss Sigma 300 VP SEM. General observations use an Everhart-Thornley-type detector at 15 or 20 kV, simultaneously with EDX analyses. High-resolution images are captured at 5 kV using a secondary electron InLens detector.

## SIMS samples preparation and observations

SIMS samples are prepared identically to SEM samples until the coating step. Particles are then suctioned using a vacuum system onto a carbon stub previously coated with an adherence polymer.

Analyses have been carried out on a Cameca-Ametek IMS 1300HR³ LG-SIMS using an oxygen primary beam.

The intensity of the signal measured with SIMS is not directly related to the concentration of the isotopes in the sample. It is necessary to consider the ionization potential of the element of interest, as well as the matrix from which it is extracted. The best way to perform such kind of measurement is to use a reference material with a similar known composition. Unfortunately, no such kind of material was available, or any real CsMP. In consequence, the relative sensitivity factor (RSF) method has been used. Based on measurements carried out on a similar matrix (SiO₂) by other researchers, the RSF of each isotope can be determined. Considering the element X that has been measured through its isotope Nₓ, its relative atomic concentration to Si can be calculated as follow:

$$ \frac{X}{Si} = \frac{i_{N_X}}{i_{28Si}} \times \frac{a_{28Si}}{a_{N_X}} \times \frac{RSF_X}{\rho_{Si}} \quad (1) $$

with i being the ionization potential, a the abundance of the concerned isotope and ρ the matrix atomic density.

## TEM-EDX samples preparation and observations

Thin sections have been prepared using an FEI Helios 600 NanoLab. Particles are laid on a silicon wafer, affixed with silver coating. In the FIB-SEM, electronic and ionic Pt deposition secures the region of interest, and micro-machining is achieved by a Ga+ source starting from 30 kV.

Thin sections are then transferred onto a Cu holder for TEM observations, conducted using an FEI Talos F200S at 200 kV.

## Data availibility

The data will be made available upon request, and interested parties may obtain it by contacting the corresponding author (Christophe Journeau—christophe.journeau@cea.fr).

Received: 21 December 2023; Accepted: 13 March 2024

Published online: 19 March 2024

## References

1. IAEA (ed.) *The Fukushima Daiichi accident* (International Atomic Energy Agency, Vienna, 2015).

2. Stohl, A. et al. Xenon-133 and caesium-137 releases into the atmosphere from the Fukushima Dai-ichi nuclear power plant: Determination of the source term, atmospheric dispersion, and deposition. *Atmos. Chem. Phys.* **12**, 2313–2343. https://doi.org/10.5194/acp-12-2313-2012 (2012).

3. Adachi, K., Kajino, M., Zaizen, Y. & Igarashi, Y. Emission of spherical cesium-bearing particles from an early stage of the Fukushima nuclear accident. *Sci. Rep.* **3**, 2554. https://doi.org/10.1038/srep02554 (2013).

4. Igarashi, Y. et al. A review of Cs-bearing microparticles in the environment emitted by the Fukushima Dai-ichi nuclear power plant accident. *J. Environ. Radioact.* **205–206**, 101–118. https://doi.org/10.1016/j.jenvrad.2019.04.011 (2019).

5. Kogure, T. et al. Constituent elements and their distribution in the radioactive Cs-bearing silicate glass microparticles released from Fukushima nuclear plant. *Microscopy* **65**, 451–459. https://doi.org/10.1093/jmicro/dfw030 (2016).

6. Satou, Y. et al. Analysis of two forms of radioactive particles emitted during the early stages of the Fukushima Dai-ichi nuclear power station accident. *Geochem. J.* **52**, 137–143. https://doi.org/10.2343/geochemj.2.0514 (2018).

7. Miura, H. et al. Discovery of radiocesium-bearing microparticles in river water and their influence on the solid-water distribution coefficient (Kₐ) of radiocesium in the Kuchibuto River in Fukushima. *Geochem. J.* **52**, 145–154. https://doi.org/10.2343/geochemj.2.0517 (2018).

8. Yamaguchi, N. et al. Internal structure of cesium-bearing radioactive microparticles released from Fukushima nuclear power plant. *Sci. Rep.* **6**, 20548. https://doi.org/10.1038/srep20548 (2016).

9. Abe, Y. et al. Widespread distribution of radiocesium-bearing microparticles over the greater Kanto Region resulting from the Fukushima nuclear accident. *Progress Earth Planet. Sci.* **8**, 13. https://doi.org/10.1186/s40645-020-00403-6 (2021).
10. Utsunomiya, S. et al. Caesium fallout in Tokyo on 15th march, 2011 is dominated by highly radioactive, caesium-rich microparticles. https://doi.org/10.48550/arXiv.1906.00212 (2019).

11. Ikehara, R. et al. Abundance and distribution of radioactive cesium-rich microparticles released from the Fukushima Daiichi Nuclear Power Plant into the environment. *Chemosphere* https://doi.org/10.1016/j.chemosphere.2019.125019 (2020).

12. Ikehara, R. et al. Novel method of quantifying radioactive cesium-rich microparticles (CsMPs) in the environment from the Fukushima Daiichi Nuclear Power Plant. *Environ. Sci. Technol.* **52**, 6390–6398. https://doi.org/10.1021/acs.est.7b06693 (2018).

13. Okumura, T., Yamaguchi, N., Dohi, T., Iijima, K. & Kogure, T. Dissolution behaviour of radiocaesium-bearing microparticles released from the Fukushima nuclear plant. *Sci. Rep.* **9**, 3520. https://doi.org/10.1038/s41598-019-40423-x (2019).

14. Suetake, M. et al. Dissolution of radioactive, cesium-rich microparticles released from the Fukushima Daiichi Nuclear Power Plant in simulated lung fluid, pure-water, and seawater. *Chemosphere* https://doi.org/10.1016/j.chemosphere.2019.05.248 (2019).

15. Okumura, T., Yamaguchi, N. & Kogure, T. Dissolution behavior of radiocesium-bearing microparticles as a function of solution compositions. *Sci. Rep.* **13**, 4307. https://doi.org/10.1038/s41598-023-31519-6 (2023).

16. Cantrel, L., Cousin, F., Bosland, L., Chevalier-Jabet, K. & Marchetto, C. ASTEC v2 severe accident integral code: Fission product modelling and validation. *Nucl. Eng. Des.* **272**, 195–206. https://doi.org/10.1016/j.nucengdes.2014.01.011 (2014).

17. Maurice, A., Beuzet, E. & Maliverney, B. Use of phébus results for edf severe accidents studies. *Ann. Nucl. Energy* **61**, 199–205. https://doi.org/10.1016/j.anucene.2013.02.038 (2013).

18. Zheng, L. et al. A proposed formation mechanism of the Type-A radiocaesium-bearing microparticles released from Units 2/3 during the Fukushima Daiichi Nuclear Power Plant accident. *J. Nucl. Mater.* **563**, 153623. https://doi.org/10.1016/j.jnucmat.2022.153623 (2022).

19. Devell, L. & Johansson, K. Specific features of cesium chemistry and physics affecting reactor accident source term predictions. Tech. Rep. SKI Report 94:29, OECD/NEA and SKI (1994).

20. Rizaal, M., Nakajima, K., Saito, T., Osaka, M. & Okamoto, K. High-temperature gaseous reaction of cesium with siliceous thermal insulation: The potential implication to the provenance of enigmatic Fukushima cesium-bearing material. *ACS Omega* **7**, 29326–29336. https://doi.org/10.1021/acsomega.2c03525 (2022).

21. Hidaka, A. Identification of carbon in glassy cesium-bearing microparticles using electron microscopy and formation mechanisms of the microparticles. *Nucl. Technol.* https://doi.org/10.1080/00295450.2021.1929767 (2021).

22. Hagiwara, H., Kondo, K. & Hidaka, A. The formation mechanism of radiocesium-bearing microparticles derived from the Fukushima Daiichi nuclear power plant using electron microscopy. *J. Radioanal. Nucl. Chem.* **331**, 5905–5914. https://doi.org/10.1007/s10967-022-08434-7 (2022).

23. Furuki, G. et al. Caesium-rich micro-particles: A window into the meltdown events at the Fukushima Daiichi Nuclear Power Plant. *Sci. Rep.* https://doi.org/10.1038/srep42731 (2017).

24. Ohnuki, T., Satou, Y. & Utsunomiya, S. Formation of radioactive cesium microparticles originating from the Fukushima Daiichi Nuclear Power Plant accident: Characteristics and perspectives. *J. Nucl. Sci. Technol.* **56**, 790–800. https://doi.org/10.1080/00223131.2019.1595767 (2019).

25. Alsmeyer, H. et al. Molten corium/concrete interaction and corium coolability: a state of the art report. *EUR(Luxembourg)* (1995).

26. Laffolley, H., Journeau, C., Delacroix, J., Grambow, B. & Suteau, C. Synthesis of Fukushima Daiichi Cs-bearing microparticles through molten core-concrete interaction in nitrogen atmosphere. *Nucl. Mater. Energy* **33**, 101253. https://doi.org/10.1016/j.nme.2022.101253 (2022).

27. Ohnuki, T. et al. Chemical species of cesium and iodine in condensed vaporized microparticles formed by melting nuclear fuel components with concrete materials. *Environ. Sci. Processes & Impacts* **25**, 1204–1212. https://doi.org/10.1039/D3EM00074E (2023).

28. Laffolley, H., Journeau, C., Thilliez, S. & Grambow, B. Thermodynamics of aerosols during a molten core-concrete interaction at Fukushima Daiichi Unit 2 estimated conditions. *Ann. Nucl. Energy* **186**, 109770. https://doi.org/10.1016/j.anucene.2023.109770 (2023).

29. Tischendorf, B. et al. The density of alkali silicate glasses over wide compositional ranges. *J. Non-Cryst. Solids* **239**, 197–202 (1998). https://doi.org/10.1016/S0022-3093(98)00740-6

30. Stringer, C., Wang, T., Michaelos, M. & Pachitariu, M. Cellpose: A generalist algorithm for cellular segmentation. *Nat. Methods* **18**, 100–106. https://doi.org/10.1038/s41592-020-01018-x (2021).

31. Wilson, R. G. & Novak, S. W. Systematics of secondary-ion-mass spectrometry relative sensitivity factors versus electron affinity and ionization potential for a variety of matrices determined from implanted standards of more than 70 elements. *J. Appl. Phys.* **69**, 466–474. https://doi.org/10.1063/1.347687 (1991).

32. Fueda, K. et al. Volatilization of B₄C control rods in Fukushima Daiichi nuclear reactors during meltdown: B-Li isotopic signatures in cesium-rich microparticles. *J. Hazard. Mater.* **428**, 128214. https://doi.org/10.1016/j.jhazmat.2022.128214 (2022).

33. Imoto, J. et al. Isotopic signature and nano-texture of cesium-rich micro-particles: Release of uranium and fission products from the Fukushima Daiichi Nuclear Power Plant. *Sci. Rep.* **7**, 5409. https://doi.org/10.1038/s41598-017-05910-z (2017).

34. Hagiwara, H., Funaki, H., Shiribiki, N., Kanno, M. & Sanada, Y. Characterization of radiocesium-bearing microparticles with different morphologies in soil around the Fukushima Daiichi Nuclear Power Plant. *J. Radioanal. Nucl. Chem.* https://doi.org/10.1007/s10967-021-08061-8 (2021).

35. OECD/NEA. State-of-the-art report on nuclear aerosols. Tech. Rep., OECD Publishing (2009).

36. Powers, D. A., Brockmann, J. E. & Shiver, A. W. VANESA: A mechanistic model of radionuclide release and aerosol generation during core debris interactions with concrete. Tech. Rep. NUREG/CR-4308, SAND-85-1370, 5150148, Sandia National Laboratories (1986). https://doi.org/10.2172/5150148.

37. Ginsberg, T. (1985) Aerosol Generation by Liquid Breakup Resulting from Sparging of Molten Pools of Corium by Gases Released During Core/Concrete Interactions. *Nucl. Sci. Eng.* **89**, 36–48. https://doi.org/10.13182/NSE85-A17881.

38. Okumura, T., Yamaguchi, N. & Kogure, T. Finding Radiocesium-bearing microparticles more minute than previously reported, emitted by the Fukushima Nuclear Accident. *Chem. Lett.* **48**, 1336–1338. https://doi.org/10.1246/cl.190581 (2019).

39. Okumura, T., Yamaguchi, N., Dohi, T., Iijima, K. & Kogure, T. Inner structure and inclusions in radiocesium-bearing microparticles emitted in the Fukushima Daiichi Nuclear Power Plant accident. *Microscopy* **68**, 234–242. https://doi.org/10.1093/jmicro/dfz004 (2019).

40. Okumura, T. et al. Reactor environment during the Fukushima nuclear accident inferred from radiocaesium-bearing microparticles. *Sci. Rep.* **10**, 1352. https://doi.org/10.1038/s41598-020-58464-y (2020).

41. Karlsen, J. T., Farrants, G., Torgrimsen, T. & Reith, A. Chemical composition and morphology of welding fume particles and grinding dusts. *Am. Ind. Hyg. Assoc. J.* **53**, 290–297. https://doi.org/10.1080/15298669291359681 (1992).

42. Higaki, S., Kurihara, Y., Yoshida, H., Takahashi, Y. & Shinohara, N. Discovery of non-spherical heterogeneous radiocesium-bearing particles not derived from unit 1 of the Fukushima Dai-ichi Nuclear Power Plant, in residences 5 years after the accident. *J. Environ. Radioact.* **177**, 65–70. https://doi.org/10.1016/j.jenvrad.2017.06.006 (2017).

43. Yamaguchi, N. et al. Structures of radioactive Cs-bearing microparticles in non-spherical forms collected in Fukushima. *Geochem. J.* **52**, 123–136. https://doi.org/10.2343/geochemj.2.0483 (2018).

44. Pellegrini, M. et al. Benchmark Study of the Accident at the Fukushima Daiichi Nuclear Power Plant (BSAF Project)—Phase I Summary Report (Tech. Rep, OECD/NEA, 2016).
## Acknowledgements

The authors would like to express their deep appreciation to Jules DELACROIX (CEA Cadarache) for the management of the VITI facility and his contributions to the experiments. The authors express their sincere gratitude to Laurent BRISSONNEAU, Sid VIDALLER and Nicolas CILLIERES (CEA Cadarache) for their contributions to the SEM analysis. Special appreciation is extended to Claire ONOFRI, Doris DROUAN, Catherine SABATHIER (CEA Cadarache) and Martiane CABIE (Aix-Marseille University) for their expertise and assistance in conducting the TEM analysis and the FIB preparation. The authors also acknowledge the outstanding support received from Anne-Laure FAURE and Manon CORNATON (CEA) for their significant role in SIMS analysis.

## Author contributions

H.L.: Writing—original draft, Conceptualization, Software, Data curation, Investigation. C.J.: Writing—review & editing, Conceptualization, Investigation, Supervision. B.G.: Writing—review & editing, Supervision.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-024-56972-9.

**Correspondence** and requests for materials should be addressed to C.J.

**Reprints and permissions information** is available at www.nature.com/reprints.

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

<figure><img src="image_9.png" /></figure>

**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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024