Plasma Diagnostics Research Articles

Investigation of the poloidal magnetic flux at the PF-3 plasma focus within the framework of the program of laboratory simulation of astrophysical jets

Astrophysical jets are collimated plasma outflows observed in diverse astrophysical settings covering seven decades of spatial scale and twenty decades of power, which, nevertheless, share many common features. This similarity over wide range of scales indicates a common core of physics underlying this phenomenon, leading to considerable interest in observational, theoretical and numerical studies. Laboratory astrophysics experiments for simulating astrophysical jets are premised on this common core of physics responsible for multi-scale similarity of jets remaining valid down to laboratory spatial scales of millimeters. Jets formed after the disassembly of the non-cylindrical z-pinch formed in a plasma focus installation have recently been subjects of observational studies. They offer an important complementarity to the main lines of investigations in two respects. Firstly, the multi-faceted role of gravity, radiation, nuclear reactions and related astrophysics is eliminated retaining only a rapid implosion of a compact plasma object in a magnetohydrodynamic environment as a common feature. Secondly, observations can be made using techniques of laboratory plasma diagnostics. In this paper, we report preliminary results regarding presence of poloidal magnetic flux associated with the jets lasting long after the pinch disassembly. This is significant in the context of uncertainty regarding the origin of poloidal magnetic field postulated in several MHD models of astrophysical jet phenomena. Evidence indicating presence of a radial component of electric field suggests existence of plasma rotation as well. These results suggest that more refined experiments can provide insights into the astrophysical jetting phenomena not available from observational astronomy techniques.

A thorough experimental assessment of THz-TDS plasma diagnostic techniques for nuclear fusion applications.

In this paper, the study of a plasma diagnostic system based on the THz time domain spectroscopy technique is presented. Such a system could potentially probe a large part of the electromagnetic spectrum currently covered by several other diagnostics in a single measurement. This feature, keeping in mind the basic requirements for plasma diagnostics in nuclear fusion experiments, such as robustness and hard environment applicability, as well as durability and low maintenance, makes the diagnostic of great interest. A conceptual design of the THz-TDS diagnostic has been developed, starting from the well-established classical microwave and far infrared plasma diagnostics landscape. The physical constraints and required instrumental characteristics have been studied and are described in detail here, together with the solutions available for each type of plasma measurement. Specific experimental laboratory tests of the different experimental configurations have been carried out, evaluating the capacity and potential of the novel diagnostic, together with the instrumental constraint, within the diagnostic parameter space.

Plasma-chamber wall interaction and its impact on polymer deposition in inductively-coupled C4F8/Ar plasmas

Plasma processing chambers, key tools in semiconductor device fabrication, must be meticulously maintained to ensure tool-to-tool matching. Processing chamber walls are liable to contamination with fluorocarbon (CxFy) film, which acts as a source and/or sink for active species during plasma processes. This contamination potentially affects processing quality and alters the plasma properties, increasing deviations between processing chambers. However, quantitative studies on the impact of wall contamination on the substrate deposition outcomes have been limited. Additionally, spatially resolved investigations are crucial for understanding the plasma-wall interaction. Therefore, this study evaluated the influence of polymer (CxFy) film wall contamination on plasma deposition via various plasma diagnostics and surface analyses. When the chamber inner wall was coated with polymer film, the thickness, surface roughness, and chemical composition of the deposited film on the bottom substrate were more uniform. These results were linked to alterations in the uniformity of active species in the plasma, such as CF2 and F, resulting from interaction with the polymer film on the wall's surface. The C–Fx and C–CFx composition ratios in the wall-coated films dramatically reduced, indicating that this film supplied the depositing precursors to the plasma. These results suggest that controlling the chamber wall conditions is important to maintain tool-to-tool matching and modulate the processing performance.

Using the Cartwheel CME to Predict Off-limb Observations of CMEs for New and Upcoming UV and EUV Spectrometers

Coronal mass ejections (CMEs) expel multithermal, magnetized plasma from the Sun, and when directed toward Earth, can cause extensive damage to space and ground-based electronics. To better understand the triggering, acceleration, and evolution of CMEs, it is critical to study CME plasma properties close to the Sun. High-resolution ultraviolet and extreme ultraviolet (UV-EUV) spectroscopy can give the most detailed plasma diagnostics of CMEs in the low solar corona. Unfortunately, very few spectrally resolved observations of CMEs in the low solar corona exist. However, with the recent launch of the Spectral Imaging of the Coronal Environment on board Solar Orbiter and the upcoming missions, including the EUV High-Throughput Solar Telescope (EUVST) on Solar-C and the Multi-slit Solar Explorer (MUSE), we will have the opportunity to obtain unprecedented, spectrally resolved CME observations. Using the only full EUV spectral observation of a CME by the Hinode/EUV Imaging Spectrometer, we predict the spectra that SPICE, EUVST, and MUSE are expected to observe during an off-limb CME eruption to investigate the diagnostic capabilities of each instrument. Finally, we provide a list of density-sensitive and temperature-sensitive ratios for CME plasma diagnostics along with the expected spectral atlas for each instrument to facilitate observing sequence planning.

Overview of the early campaign diagnostics for the SPARC tokamak (invited).

The SPARC tokamak is a high-field, Bt0 ∼12T, medium-sized, R0 = 1.85m, tokamak that is presently under construction in Devens, MA, led by Commonwealth Fusion Systems. It will be used to de-risk the high-field tokamak path to a fusion power plant and demonstrate the commercial viability of fusion energy. SPARC's first campaign plan is to achieve Qfus > 1 using an ICRF-heated, <10MW, high current, Ip ∼ 8.5 MA, L-mode fueled by D-T gas injection, and its second campaign will investigate H-mode operations in D-D. To facilitate plasma control and scientific learning, a targeted set of ∼50 plasma diagnostics are being designed and built for operation during these campaigns. While nearly all diagnostics are based on established techniques, the pace of deployment, relative to the first plasma, and the harshness of the thermal, electromagnetic, and radiation environment are unprecedented for medium-sized tokamaks. An overview of the SPARC diagnostic set is given, providing context to further details communicated by the SPARC team in companion publications that are system-specific. The system engineering philosophy for SPARC diagnostics is outlined, and the design and engineering verification process for components inside and outside the primary vacuum boundary are described. Diagnostics are mounted directly to the vacuum vessel as well as housed within a series of eight midplane and 24 off-midplane replaceable port plugs. With limited exceptions, signal conditioning, digitization electronics and cameras as well as lasers and microwave sources are localized to a series of five Diagnostic Lab spaces, totaling ∼350m2, located >15m from the center of the tokamak, on the other side of a 2.4m concrete shielding wall. A series of 31 large-scale penetrations have been included in the SPARC Tokamak Hall to facilitate integration of early campaign diagnostics and to provide upgradability.

High-resolution Imaging Spectroscopy of a Tiny Sigmoidal Minifilament Eruption

Minifilament eruptions producing small jets and microflares have mostly been studied based on coronal observations at extreme-ultraviolet and X-ray wavelengths. This study presents chromospheric plasma diagnostics of a quiet-Sun minifilament of size ∼ 2″ × 5″ with a sigmoidal shape and an associated microflare observed on 2021 August 7 17:00 UT using high temporal and spatial resolution spectroscopy from the Fast Imaging Solar Spectrograph (FISS) and high-resolution magnetograms from the Near InfraRed Imaging Spectropolarimeter (NIRIS) installed on the 1.6 m Goode Solar Telescope at Big Bear Solar Observatory. Using FISS Hα and Ca ii 8542 Å line spectra at the time of the minifilament activation we determined a temperature of 8600 K and a nonthermal speed of 7.9 km s−1. During the eruption, the minifilament was no longer visible in the Ca ii 8542 Å line, and only the Hα line spectra were used to find the temperature of the minifilament, which reached 1.2 × 104 K and decreased afterward. We estimated thermal energy of 3.6 × 1024 erg from the maximum temperature and kinetic energy of 2.6 × 1024 erg from the rising speed (18 km s−1) of the minifilament. From the NIRIS magnetograms we found small-scale flux emergence and cancellation coincident with the minifilament eruption, and the magnetic energy change across the conjugate footpoints reaches 7.2 × 1025 erg. Such spectroscopic diagnostics of the chromospheric minifilament complement earlier studies of minifilament eruptions made using coronal images.

Development of the neutral gas diagnostic system for neutral pressure measurements and gas analysis on the SPARC tokamak.

A suite of plasma diagnostics will be installed on the SPARC tokamak to allow for real-time plasma control, an investigation of high-field tokamak physics, and to de-risk the design of ARC, a compact fusion power plant with the aim to supply electricity to the grid. Among these diagnostics is the neutral gas diagnostics system (NTGS), a set of pressure sensors and gas analyzers used to monitor neutral pressure and gas composition for plasma control, optimization of wall conditioning, and helium ash removal, among other measurement functions linked to operational and scientific research needs. While reliable measurements of neutral pressure and gas composition have been fielded on existing magnetic-confinement fusion devices, SPARC represents a step increase in challenge due to its larger power density, higher field, high vacuum vessel bake temperatures, and higher neutron flux environment, as well as a step decrease in the accessibility for maintenance of in-vessel sensors. Multiple sensor types will be employed to have defense-in-depth and mitigate common failure modes. The NTGS system is currently progressing through final design, working to close out decisions using prototyping and analysis, and then moving on to procuring sensors for assembly and installation on SPARC. This paper outlines the current status of the system design and the diagnostic requirements that motivate neutral gas measurements on SPARC, as well as highlights the planned prototyping activities.

Microwave diagnostics of flare plasma by the direct fitting method based on data from the Siberian Radioheliograph

In this paper, we analyze images and the frequency spectrum of microwave emission in the maximum of brightness distribution in the January 20, 2022 and July 16, 2023 flares recorded by the Siberian Radioheliograph in the 3–6 GHz and 6–12 GHz ranges. We use the obtained spectrum data for radio diagnostics of magnetic field strength and orientation, plasma density, and parameters of accelerated particles in a radio source. The radio diagnostics is carried out by a method based on minimizing the functional containing the intensities of theoretically calculated and observed frequency spectra of left-polarized and right-polarized emission. Since the form of such a multidimensional functional is quite complex, and it is not possible to minimize it by standard approaches, we employ a genetic minimization method. The radio diagnostics allows us to determine features of the dynamics of the magnetic field intensity and orientation, as well as the density and the energy spectral index of non-thermal electrons in the region of maximum brightness of the radio source. We have found that during the growth phase of the main radiation peaks the magnetic field decreases, whereas during the decay phase, on the contrary, it increases. The rate of these changes varies from a few G/s to 11 G/s for the January 20, 2022 flare and is about 1 G/s for the July 16, 2023 flare.

Collisional-Radiative modeling of unresolved transition array spectra near 200 Å from W17+-W25+ emissions for diagnostics of ITER edge plasma

Tungsten (W) is one of the major impurities in ITER and future DEMO reactors. However,diagnosing the ion density, temperature, and spatial distribution of tungsten ions in intermediate charge states, such as W8+-W27+, is difficult because of the lack of spectral line data. In this study, we observed a tungsten Unresolved Transition Array (UTA) spectrum, which has a pseudo-continuum structure, around W20+ in the Extreme Ultraviolet wavelength region in the Large Helical Device (LHD). We conducted Collisional-Radiative (CR) modeling for W17+-W25+. Two pseudo-continuum peaks corresponding to 5 s-5p transition and transition between doubly excited states of 5 s2-5s5p, are strongly emitted around 200 Å for each ion. The synthesized spectrum of W17+-W25+ ions reproduced the observed LHD spectrum near 200 Å. From the electron temperature dependence of the UTA spectral shape, the UTA shifted toward longer wavelength region with several pseudo-continuum peaks appearing for ions in lower charge states, with decreasing electron temperature from 0.4 keV to 0.2 keV. This result qualitatively explained the observed time evolution of the spectrum.

Plasma diagnostics of supernova remnant 3C 400.2 by Suzaku observations

Plasma diagnostics of supernova remnant 3C 400.2 by Suzaku observations

Electron transport barrier and high confinement in configurations with internal islands close to the plasma edge of W7-X

Abstract The low magnetic shear in the Wendelstein 7-X (W7-X) stellarator makes it feasible to shape the separatrix by the large islands constituting an island-divertor, and this can be exploited to access various magnetic configurations, including samples of different internal island sizes and locations. To investigate the configuration effects on the plasma confinement, a configuration scan was performed by changing the coil currents to vary the rotational transform between values 5 / 4 and 5 / 6 at the plasma boundary with different power levels (2, 4, 6 MW) of electron cyclotron resonance heating (ECRH) at a maximum plasma density of 8 × 10 19 m − 2 . neutral beam injection (NBI) heating was also applied during some configurations of the scan to create a density ramp and access high densities beyond the X2 ECRH cutoff. For the magnetic configurations, where the 5 / 5 and 5 / 6 island chains were moved closer to separatrix but remaining inside the last closed flux surface, the electron cyclotron emission shows that an electron temperature, T e , pedestal develops already during ECRH heated plasma buildup phase indicating a transport barrier, and the barrier sustains irrespective of changed plasma heating conditions such as NBI in the later part of discharge. The transport barrier is broken by subsequent fast crashes, observed with multiple plasma diagnostics with characteristics such as tokamak edge localized modes, and the corresponding crash amplitude and frequency vary with plasma pressure. The impact of the transport barrier on plasma confinement can be seen through the increased core T e profile, which could be responsible for the overall increase in the stored diamagnetic energy by approximately 10% for these configurations. After the plasma heating is terminated, a backwards transition to a degraded confinement state is also observed. These observations indicate a configuration triggered high confinement mode in low shear W7-X. This work focuses on the occurrence of this transport barrier for different magnetic configurations and its relation to internal magnetic islands.

Effect of radiofrequency bias power on transmission spectrum of flat-cutoff sensor in inductively coupled plasma

Real-time monitoring of plasma parameters at the wafer plane is important because it significantly affects the processing results, yield enhancement, and device integrity of plasma processing. Various plasma diagnostic sensors, including those embedded in a chamber wall and on-wafer sensors, such as flat-cutoff sensors, have been developed for plasma measurements. However, to measure the plasma density on the wafer surface in real-time when processing plasma with bias power, such as in the semiconductor etching process, one must analyze the transmission spectrum of the flat-cutoff sensor in an environment with bias power applied. In this study, the transmission-spectrum and measured plasma-density characteristics of an electrode-embedded flat-cutoff sensor are analyzed via electromagnetic simulations and experiments under applied bias power. Our findings indicate that the flat-cutoff sensor accurately measures the plasma density, which is equivalent to the input plasma density under low bias power. Conversely, under high bias power, the plasma density measured by the sensor is lower than the input plasma density. Also, a thick-sheath layer is formed owing to the high bias power, which may complicate the measurement of plasma parameters using the flat-cutoff sensor. Plasma diagnostics using a flat-cutoff sensor in thick-sheath environments can be achieved by optimizing the flat-cutoff sensor structure. Our findings can enhance the analysis of plasma parameters on-wafer surfaces in processing environments with bias power applied.

Characterization of the image plate multi-scan response to mono-energetic x-rays.

Image plates (IPs), or phosphor storage screens, are a technology employed frequently in inertial confinement fusion (ICF) and high energy density plasma (HEDP) diagnostics because of their sensitivity to many types of radiation, including, x rays, protons, alphas, beta particles, and neutrons. Prior studies characterizing IPs are predicated on the signal level remaining below the scanner saturation threshold. Since the scanning process removes some signal from the IP via photostimulated luminescence, repeatedly scanning an IP can bring the signal level below the scanner saturation threshold. This process, in turn, raises concerns about the signal response of IPs after an arbitrary number of scans and whether such a process yields, for example, a constant ratio of signal between the nth and n + 1st scan. Here, the sensitivity of IPs is investigated when scanned multiple times. It is demonstrated that the ratio of signal decay is not a constant with the number of scans and that the signal decay depends on the x-ray energy. As such, repeatedly scanning an IP with a mixture of signal types (e.g., x ray, neutron, and protons) enables ICF and HEDP diagnostics employing IPs to better isolate a particular signal type.

Frontier system-on-chip (SoC) technology for microwave diagnostics (invited).

The next generation of fusion reactors, exemplified by projects such as the Demonstration Power Plant following the International Thermonuclear Experimental Reactor, faces the monumental challenge of proving the viability of generating electricity through thermonuclear fusion. This pursuit introduces heightened complexities in diagnostic methodologies, particularly in microwave-based diagnostics. The increased neutron fluence necessitates significant reductions in vessel penetrations and the elimination of internal diagnostics, posing substantial challenges. SoC technology offers a promising solution by enabling the miniaturization, modularization, integration, and enhancing the reliability of microwave systems. After seven years of research, our team successfully pioneered the V- and W-band system-on-chip approach, leading to the development of active transmitters and passive receiver modules applied in practical settings, notably within the DIII-D tokamak project. Arrays of these modules have supported microwave imaging diagnostics. New physics measurement results from the Electron Cyclotron Emission Imaging system on DIII-D provide compelling evidence of improved diagnostics following the adoption of SoC technology. Furthermore, we achieved a breakthrough in developing an F-band SoC, advancing higher frequency capabilities for fusion devices. These achievements represent a significant leap forward in fusion diagnostic technology, marking substantial progress toward establishing reliable and efficient plasma diagnostics for future fusion reactors.

Analysis of an induced Langmuir wave by ponderomotive forces and its applicability for plasma diagnostics

We present a numerical study on the electron and ion density perturbation in low-temperature plasmas driven by the frequency detuning of two intense laser beams. Our study is performed in the hydrodynamic regime, which becomes applicable when the plasma grating period induced by the beating of the laser beams is greater than the Debye length and collective processes such as plasma oscillations can be excited. Our findings show a resonance in electron density perturbation as the frequency detuning approaches a value consistent with the Bohm–Gross dispersion relation in low- and high-pressure plasmas. We discuss the potential of this resonance as a diagnostic tool for precisely measuring electron temperature and density in low-temperature plasmas through coherent scattering.

Bayesian modelling of multiple plasma diagnostics at Wendelstein 7-X

Abstract Inference of electron density and temperature has been performed using multiple, diverse sets of plasma diagnostic data at Wendelstein 7-X. Predictive models for the interferometer, Thomson scattering and helium beam emission spectroscopy (He-BES) systems have been developed within the Minerva framework and integrated into a unified model. Electron density and temperature profiles are modelled using Gaussian processes. Calibration factors for the Thomson scattering system and predictive uncertainties are considered as additional unknown parameters. The joint posterior probability distribution for the electron density and temperature profiles as well as Gaussian process hyperparameters and model parameters is explored through a Markov chain Monte Carlo algorithm. Samples from this distribution are numerically marginalised over the hyperparameters and model parameters to yield marginal posterior distributions for the electron density and temperature profiles. The profile inferences incorporate various data combinations from the interferometer and Thomson scattering as well as constraints at the limiter/divertor positions through virtual observations or edge data from He-BES. Additionally, the integration of x-ray imaging crystal spectrometer data into the model for ion temperature profiles is presented. All profiles presented in this study are inferred with optimally selected hyperparameters and model parameters by exploring the joint posterior distribution, inherently applying Bayesian Occam’s razor.

Laboratory Transition-rate Measurement of the Coronal Intercombination Line of Ar xv by Time-resolved Laser Spectroscopy

The extreme-ultraviolet emission line (424 Å) of the intercombination 1s 22s 2 1 S 0–1s 22s2p 3 P 1 transition of Ar xv can potentially characterize the electron temperature of astrophysical plasma. Various theoretical studies have investigated the intercombination transition rate, which is essential for the plasma diagnostics; however, experimental difficulties have prevented its measurement. We present here measurement of the lifetime of the 3 P 1 excited state of Ar xv, providing the experimental value of the intercombination transition rate. Employing time-resolved plasma-assisted laser spectroscopy, a method we recently demonstrated, enables us to measure this submicrosecond lifetime. The experimental result exhibits a 25%–43% higher transition rate than theoretical predictions.

Sustainable H2O2 production via solution plasma catalysis

Clean production of hydrogen peroxide (H2O2) with water, oxygen, and renewable energy is considered an important green synthesis route, offering a valuable substitute for the traditional anthraquinone method. Currently, renewable energy-driven production of H2O2 mostly relies on soluble additives, such as electrolytes and sacrificial agents, inevitably compromising the purity and sustainability of H2O2. Herein, we develop a solution plasma catalysis technique that eliminates the need for soluble additives, enabling eco-friendly production of concentrated H2O2 directly from water and O2. Screening over 40 catalysts demonstrates the superior catalytic performance of carbon nitride interacting with discharge plasma in water. High-throughput density functional theory calculations for 68 models, along with machine learning using 29 descriptors, identify cyano carbon nitride (CCN) as the most efficient catalyst. Solution plasma catalysis with the CCN achieves concentrated H2O2 of 20 mmol L-1, two orders of magnitude higher than photocatalysis by the same catalyst. Plasma diagnostics, isotope labeling, and COMSOL simulations collectively validate that the interplay of solution plasma and the CCN accounts for the significantly increased production of singlet oxygen and H2O2 thereafter. Our findings offer an efficient and sustainable pathway for H2O2 production, promising wide-ranging applications across the chemical industry, public health, and environmental remediation.

Operando facet control of hydrophobic iron nitride nano-coral for surface protection via field-modulated plasma strategy

A facile facet reconfiguration has been considered as an essential issue for nitride coating in promoting surface protection behaviour. Compared to physical-chemical methods to achieve heterogeneous coating, plasma nitridation strategy is regarded as a favourable technique to in-situ deliver nitride protecting layer, but it is difficult to precisely control surface framework through plasma technique due to complicated interaction environment. To simply achieve nitride coating with favourable surface protection behaviour, we design a novel auxiliary insulator-confined plasma system to directly achieve hydrophobic iron nitride nano-coral (hFeNC) with desirable corrosion- and wear-resistant properties by controlling surface heating process during plasma nitridation. The resultant hFeNC nano-framework delivers corrosion rate of 1.2 mpy, 4- and 9-times lower than that of solid iron nitride film (sFeN) via normal plasma and initial Fe in NaCl condition, respectively. Operando plasma diagnostics along with numerical simulation further confirm the effect of surface heating on typical plasma parameters as well as the iron nitride nano-frameworks, indicating surface heating as the key factor responsible for the favourable surface protection coating.

Electrical and Optical Characterization of Atmospheric Pressure Plasmas for the Treatment of OPV Materials

In the last years, non-equilibrium cold plasma at atmospheric pressure have led to numerous opportunities for material treatment. Depending on the gas mixtures, the electrodes configurations, the electrical waveform used to generate the discharge, or the material, a wide range of surface properties can be created. Among the plethora of applications, the organic photovoltaic (OPV) community has recently explored the use of plasma [1]. Thanks to the versatility of the cold atmospheric pressure plasmas (CAPP), all the layers used in the OPV cells could be potentially treated and modified. The use of dry atmospheric pressure processes also opens the doors to fast treatments at reduced cost compared to low pressure [2].Possible CAPP configurations for the OPV include plasma jets and dielectric barrier discharges (DBD) in volume or surface [3]. Among the different use of plasma in OPV, surface cleaning has been largely studied. For example, the removal of organic contaminants from indium tin oxide (ITO) layer have successfully been achieved with all the aforementioned configurations as well as different gas mixtures such as N2, O2, or He [4–6]. Other interesting mechanisms, such as the passivation of the electron transport layer [7], the modification of the work function [8], or for the encapsulation of the whole OPV cells [9] have been also explored.Still, whatever the material, either the metallic electrodes, the semiconductive active layer, or the dielectric substrates, a fine control of the plasma properties is necessary to ensure the good reproducibility and the stability of the process. It is hence necessary to characterize the physical regime and clarify the chemical reactions suitable to assess the good operation of the process during the material treatment.This work focuses on the use of electrical measurements and optical emission spectroscopy (OES) to retrieve key parameters from the process, like the discharge regime (Townsend vs filamentary), the capacitances of the system, or the power dissipated during the treatment [10]. In this context, this work aims at characterizing different discharge regimes for the surface treatment of thin characteristic layers employed in OPV (i.e., ITO, fluorine doped tin oxide (FTO) and zinc oxide). The influence of different electrical signals (low frequency sinusoidal voltage vs nanopulse), as well as different configurations and dielectric materials are compared. In order to link and allow to use such measurement as monitoring tools for the treatment process of materials used in OPV cells, the treated surfaces are also analyzed by SEM and contact angle. The extracted quantities from the plasma diagnostics could hence be used as a predictive tool to forecast the final properties of the treated materials and improve the processes of material treatment in OPV.[1] Mariotti et al. https://doi.org/10.1002/ppap.201500187[2] Vida et al. https://doi.org/10.37904/nanocon.2019.8646[3] Homola et al. https://doi.org/10.1016/B978-0-323-89930-7.00001-7[4] Chiang et al. https://doi.org/10.1007/s11090-010-9237-4[5] Yi et al. https://doi.org/10.1016/j.surfcoat.2003.08.011[6] Hvojnik et al. https://doi.org/10.1016/j.mssp.2021.105850[7] Polydorou et al. https://doi.org/10.1039/C6TA03594A[8] Chaney et al. https://doi.org/10.1016/S0169-4332(01)00347-6[9] Juillard et al. https://doi.org/10.1002/admi.202000293[10] Pipa et al. https://doi.org/10.3390/atoms7010014