Blandine Fontainea,
Youssef Benrkiaa,
Jean-François Blacha,
Christian Mathieua,
Pascal Rousselb,
Ahmad I. Ayeshc,
Adlane Sayedea and
Sébastien Saitzek*a
aUniv. Artois, CNRS, Centrale Lille, Univ. Lille, UMR 8181, Unité de Catalyse et Chimie du Solide (UCCS), F-62300 Lens, France. E-mail: sebastien.saitzek@univ-artois.fr; Fax: +33 321177955; Tel: +33 321791732
bUniv. Lille, CNRS, Centrale Lille, Univ. Artois, UMR 8181, Unité de Catalyse et Chimie du Solide (UCCS), F-59000 Lille, France
cPhysics Program, Department of Math. Stat. and Physics, College of Arts and Sciences, Qatar University, P. O. Box: 2713, Doha, Qatar
First published on 19th April 2023
Polymorphic phases of copper pyrovanadate (α- and β-Cu2V2O7) were synthesized by solid state reaction and the mechanisms governing the phase transitions have been highlighted by the ThermoGravimetric Analysis (TGA) and the Differential Scanning Calorimetry (DSC). The thermal evolution of the lattice parameters was determined by high temperature X-ray Diffraction revealing negative thermal expansion coefficients. The thermogravimetric analysis coupled with differential scanning calorimetry was also used to determine the optimal conditions to obtain a dense target in order to produce thin films by the Pulsed Laser Deposition (PLD) technique. Thin films elaborated under different oxygen pressures and temperatures exhibit a β-Cu2V2O7 polycrystalline phase and their band gap indicates absorption in the visible range. These oxides can be used as photoanodes and their photoelectrochemical properties were studied for both bulk (α-Cu2V2O7) and thin films (β-Cu2V2O7), as a function of the wavelength and/or intensity of the luminous flux. The best photocurrent efficiency was obtained under 450 nm illumination. Moreover, in the case of thin films, we have observed a linear evolution of the current density with the luminous flux. Finally, the photostability of thin films was measured and shows a reduction in the photocurrent of 8% after 1 h of measurement. This photocorrosion phenomenon was also highlighted by the elemental mapping performed on thin films by Scanning Electron Microscopy (SEM) coupled with Energy Dispersive X-ray Spectrometry (EDS).
Today, many ternary metal oxides are studied for applications as photoelectrodes because they have suitable bandgaps and bands edge positions for being used in water-splitting photoelectrochemical tandem cells under solar irradiation. Among these ternary metal oxides, there are main families such as spinels (CuFe2O4,8 CaFe2O4,9…), delafossite-type oxides (CuFeO2,10 AgRhO2,11…), perovskites (LaFeO3,12 BiFeO3,13…) and many others.14 Copper pyrovanadate (Cu2V2O7) oxides is also promising candidates for the design of photo-electrodes which can be used in PEC cells. For example, Guo et al.15 obtained a current density of 12 μA cm−2 (at 1.23 V vs. NHE) on thick films of β-Cu2V2O7 produced by the drop-casting technique in 0.1 M sodium borate buffer solution. However, these properties can be improved by controlling the microstructure of the oxide. Thus, Khan et al.16 highlighted a significant improvement in current density (0.70 mA cm−2) for thick films made up of micro-flakes arranged in channels; C. Gadiyar et al.17 indicated a current density of 0.23 mA cm−2 (always at 1.23 V vs. NHE) in 0.1 M sodium borate buffer solution (pH = 9.2) containing 0.1 M Na2SO3 as a hole scavenger. Their study shows an optimal crystallite size at 40 nm corresponding to the hole diffusion length, estimated between 20 and 40 nm.18 Copper pyrovanadates (α- or β-Cu2V2O7 oxides) have also been studied for other applications, in particular for: (i) magnetoelastic and ferromagnetic coupling properties; (ii) ferroelectric properties induced by a magnetic order;19 (iii) photo-degradation properties of azo-dye in water;20,21 (iv) solar water splitting;22,23 (v) catalytic SO3 decomposition;24 (vi) cathode materials in high voltage Li-batteries.25–27
Copper pyrovanadate oxides have three polymorphic phases: (i) α-Cu2V2O7 (blossite) has an orthorhombic structure with Fdd2 space group and lattice parameters: a = 20.68 Å, b = 8.411 Å and c = 6.448 Å;28 (ii) β-Cu2V2O7 (ziesite) crystallizes in a monoclinic structure with A2/a space group and lattice parameters: a = 10.094(15) Å, b = 8.020(13) Å, c = 7.711(10) Å and β = 110.43(9)°;29 (iii) γ-Cu2V2O7 (observed at high temperature in the range of 710 °C to 780 °C) has a triclinic structure with P space group and the following lattice parameters: a = 5.0873(10) Å, b = 5.8233(11) Å, c = 9.4020(18) Å, α = 99.780(3)°, β = 97.253(3)°, γ = 97.202(3)°.30 Under normal pressure and temperature conditions, α- or β-Cu2V2O7 phase can be synthesized. For example, the reaction between V2O5 and CuO, at low temperature (T < 550 °C), lead to the stabilization of β-Cu2V2O7. While for syntheses at higher temperature (T > 560 °C), α-Cu2V2O7 phase will be privileged.31 However, fairly long calcination times are required,32 to achieve a pure phase. The β-Cu2V2O7 → α-Cu2V2O7 phase transition takes place around 605 °C. However, the β-Cu2V2O7 phase can also be obtained from α-Cu2V2O7 phase by heating it up to 710 °C and then performing a rapid cooling. In this case, the phase transition process is as follows: α-Cu2V2O7 → γ-Cu2V2O7 → β-Cu2V2O7.31
In the literature, several synthetic routes are described, for example solid-state reaction,31 wet chemical method,21 Sol–gel synthesis,33 hydrothermal synthesis,34 thermal decomposition method35 and finally flux method in close crucible to obtain single crystals.36 In thin film form, studies are scarce and to the best of our knowledge only films produced by electrospraying,22 simplified Successive Ionic Layer Adsorption and Reaction (s-SILAR) method37 and drop casting method15 have been described.
In this present work, we studied the photoelectrochemical properties of Cu2V2O7 thin films synthesized by Pulsed Laser Deposition (PLD). The first part of the work, devoted to the synthesis of pure powder and the preparation of the PLD target will highlight the difficulties of obtaining a dense target. In the second part, and for the first time, we present the optimal conditions for the growth of thin films, in particular the influence of pressure and temperature. Finally, we will compare the photoelectrochemical properties of PLD-based thin films with bulk form.
The thin films were synthesized by PLD technique using a Compex Pro 102 Laser (KrF excimer laser with a wavelength of 248 nm). The synthesis conditions are: 1.5 J cm−2 for the laser density energy, 4.5 cm for the target-substrate distance and 4 Hz for the pulse repetition rate. Before the deposition, a pre-ablation of 15 min is carried out in order to clean the target's surface.
The substrate temperature during deposition and the oxygen pressure are variable parameters studied in this work.
The structural properties of both powders and thin films were performed by X-Ray Diffraction (XRD) and Raman Spectroscopy. For the bulk form, the XRD patterns were recorded in the 15–70° range with a step of 0.02° and a scan speed of 0.4° per min, using a Rigaku Ultima IV X-ray diffractometer equipped with Cu anticathode (λKα = 1.5418 Å). Nickel foil filter and Soller slits are used to attenuate Cu Kβ radiation and to limit the divergence of X-ray beam, respectively. While the thermo-diffraction patterns were recorded on a D8 Bruker device. For the thin films, the XRD patterns were recorded using a Rigaku SmartLab diffractometer attached to a 9 kW rotating anode X-ray generator (λKα1 = 1.54056 Å) and monochromatized with a double Ge (220) monochromator. Raman spectra were recorded in back scattering geometry using a Horiba Labram HR800 spectrometer coupled with a microscope which focused the 784 nm laser beam to a 1 μm spot. The spectra were recorded in the 100–1100 cm−1 range.
Thermal Analyses have been performed using a TGA/DSC3+ STARe System (Mettler Toledo). The sample was analyzed in the range 50 °C–750 °C with a rate of 10 °C min−1 under a flow of air of 50 mL min−1 in 70 μL Al2O3 crucible.
Infrared (IR) analysis was conducted using a Shimadzu IR Prestige 21 instrument with the addition of single-reflection diamond module (ATR). IR spectra were recorded in the range 500–4000 cm−1 with a 1 cm−1 resolution. A Nova made scanning electron microscope (SEM) model NanoSEM-450 was employed for investigation of morphology and grain size and the elementary maps were recorded with a Hitachi SU3800 SEM coupled with Energy Dispersive X-ray Spectrometry (EDS). The optical properties were characterized using a Shimadzu UV-visible spectrophotometer UV-2600 in the range of 190–800 nm with a resolution of 1 nm. This spectrophotometer was equipped with an integration sphere ISR-2600+, a photomultiplier tube and a InGaAs detector. The UV-visible diffuse reflectance spectra can be used to determine the optical gap band (Eg) using the Tauc plot defined by the following equation:
αhν = A(hν − Eg)n | (1) |
The photoelectrochemical measurements were performed using a three-electrode cell (PEC 15 mL Redox.me®). The cell is calibrated for an illumination of 1 cm2 from the backside of the ITO/glass substrate (Delta Technologies, USA). The Pt wire was used as counter electrode and Ag/AgCl in saturated KCl electrode was employed as the reference electrode. The working electrode consists of a Cu2V2O7 deposited either by drop casting using the powder directly (a suspension of nanopowders is prepared in N,N-dimethylformamide/distilled water with the 1:1 ratio, then annealed at 250 °C during 1 h in order to obtain a layer whose thickness is ∼50 μm) or by PLD technique.
The photoelectrochemical measurements were carried out using a PGSTAT204 potentiostat/galvanostat (Metrohm) expanded with Electrochemical Impedance Spectroscopy (EIS) module and coupled with optical bench consisting of highly focused LED source. The intensity of the light beam can be controlled as well as the excitation wavelength via the use of low spectral dispersion LEDs. The whole is controlled by NOVA 2.0 software. Finally, the electrolyte solution consists of 0.1 M sodium borate buffer solution also called NaBi Buffer (pH 9.2) with 0.1 M sodium sulfite as a hole scavenger. This electrolyte was used to preserve the chemical stability of Cu2V2O7 (ref. 38) and to improve the intensity of the collected photocurrent.39
Fig. 1 (a) XRD pattern of α-Cu2V2O7 as prepared and α-Cu2V2O7 reference; (b) SEM images of α-Cu2V2O7 powders; (c) absorbance spectrum of α-Cu2V2O7 (inset the calculated Tauc plot). |
The characteristic vibrational modes were also highlighted by IR spectroscopy using ATR configuration (Fig. S2 – ESI†). The bands located at 3330 and 1620 cm−1 are assigned to the traditional vibrations of OH stretching and bending vibration of adsorbed water molecules, respectively. The band centered to 993 cm−1 can be attributed to V–O stretching vibration and the band located around 700 cm−1 can be assigned to Cu–O bond. Moreover, the bands located between 750 cm−1 and 900 cm−1 can be assigned to the VO3 symmetric and asymmetric stretching vibrations.
The morphology of as prepared α-Cu2V2O7 is presented on SEM image (Fig. 1b). The powder has a granular morphology with a mostly spherical shape with heterogeneous size. The energy dispersive spectroscopy (EDS) confirms the correct stoichiometric ratio (% at. Cu/% at. V = 1) without heterogeneity.
The optical band gap energy was determined from the UV-visible absorption spectrum using Tauc method (Fig. 1c), knowing that α-Cu2V2O7 has an indirect band gap (n = 1/2).43 The determined value of the optical band gap is 2.35 eV (i.e. 528 nm). This value is particularly interesting because it is very close to the maximum intensity of the solar spectrum.
In order to prepare the PLD target, it is essential to study the Cu2V2O7 oxide stability by performing a thermal analysis by ThermoGravimetric Analysis (TGA) and Differential Scanning Calorimetry (DSC) measurements. For DSC, two heating/cooling thermal cycles are carried out between 50 °C to 750 °C with 10 °C min−1 and 20 °C min−1 rate, respectively. The results obtained are presented in Fig. 2a.
During the first cycle, an endothermic peak (I) is observed on heating, attributed to α → γ phase transition (677 °C).44 Then, on cooling, three exothermic peaks appear. The first peak (III) corresponds to γ → β phase transition (∼595–600 °C). The two other peaks (IV and V) correspond to crystallization at low temperature of β-Cu2V2O7 phase. Indeed, B. V. Slobodin et al.44,45 indicates that above 720 °C, a small amount of CuVO3 is formed which forms an eutectic with the α and β phases having a low crystallization temperature.44 For the second cycle, a difference is observed. Indeed, the appearance of a new endothermic peak (II) is observed, which is related to the formation of β-Cu2V2O7 phase during the cooling step of the first cycle. This peak is attributed to β → α transition (557 °C). The second cooling step took place in the same way as that of the first cycle (the difference between peaks IV and V is only due to the different cooling rates thus confirming that it is an enthalpy of crystallization). The analysis of the weight loss with temperature presented in Fig. 2b, highlights the reduction of copper observed during the cooling step of the first cycle. The weight loss after 715 °C thus highlights a copper reduction (oxygen release) and the start of peritectic transformation around 770 °C. Indeed, above 725 °C, Cu2V2O7 undergoes a congruent fusion followed by an incongruous transformation towards 770 °C (peritectic decomposition).45
In order to understand this low crystallization temperature (IV and V peaks) due to the presence of the eutectic, additional measurements were undertaken but this time with a temperature limited to 720 °C. At this temperature, Cu2V2O7 phase is not reduced and therefore should not form eutectics. The TGA/DSC analysis cycles limited to 720 °C are shown in Fig. 3. As already observed by B. V. Slobodin et al.,44 the first heating step leads to the transformation α → γ at 705 °C (the XRD reference pattern is shown in Fig. S3 – ESI†). Then, the first cooling is accompanied by a release of heat due to the transformation γ → α → β. The second cycle shows an endothermic peak due to the transition from β → α at 604 °C (this peak is due to β-Cu2V2O7 which is formed during the first cooling), then comes the α → γ phase transformation, as during the first heating. The following cycles are similar and the crystallization peaks at low temperature are no longer observed, thus confirming the presence of an eutectic with the phase formed above 725 °C. This decomposition at a temperature higher than 725 °C is problematic because it limits the sintering temperature of the PLD target and therefore its densification.
Fig. 3 (a) DSC and TGA measurements realized with two heating/cooling cycles; (b) fragment of the curves highlighting the observed transitions. |
To confirm this behavior, in situ high temperature X-ray diffraction was performed over the temperature range 50 °C to 725 °C (heating/cooling in the same cycle). Fig. 4a shows a partial 2D map (XRD patterns vs. temperature) of α-Cu2V2O7 powder (the complete 2D map is presented in Fig. S4 ESI†).
Initially, the presence of a pure α phase is observed. When the temperature increases, some diffraction peaks shift towards the large angles (as shown in Fig. S5 – ESI†), highlighting the presence of a negative Coefficient of Thermal Expansion (CTE).46 This shift is not identical according to the crystallographic planes. This reflects an anisotropic evolution of the unit cell as a function of temperature.47
Around 710 °C, a phase transition is observed (dashed lines in Fig. 4a). This transition corresponds to the transformation of the α phase towards a phase which can be attributed to the Cu1.98V1.96O6.92 phase with C2/c space group, previously described by S. A. Petrova et al.48 and called β′-Cu2V2O7 (Fig. S7 – ESI†). The presence of this phase also explains the mass loss observed in the TGA (oxygen release and Cu2+ reduction) measurements and we can note that a slight part seems to have melted in the sample holder, as observed at the end of the DSC thermal cycle. During cooling, a phase transition is also observed from the β′-Cu2V2O7 phase to the β-Cu2V2O7 phase at 600 °C. This phase is then preserved when returning to room temperature. From a structural point of view, the β′ → β phase transition causes a decrease in the V–O–V angle and a slight expansion of the VO4 tetrahedra.48 Here, a discrepancy with the TGA/DSC thermal analysis has to be noted. In the literature, α → γ transition is observed around 700 °C (as described in the DSC analysis part). However, in this present case, the XRD patterns indicate unambiguously that the obtained phase is β′-Cu2V2O7 phase (Cu1.98V1.96O6.92). Consequently, the α → γ phase transition previously claimed is in fact a α → β′ transition. Thus, all of the ATG data can be re-interpreted in the same way but this time by replacing the γ phase by the β′ phase (structural evolution of Cu2V2O7 oxide is presented in Fig. S8 – ESI†).
Concerning β-Cu2V2O7, X-ray thermodiffractograms show a negative expansion behavior (as depicted in Fig. S9 and S10 – ESI†) which correlates with the literature.49 In addition, the dotted circles in Fig. 4a reveal the formation of a small amount of β-Cu2V2O7 from 500 °C then disappears around 600 °C (corresponding to the β → α transition as observed on the DSC curve in Fig. 3). In summary, a good correlation between DSC and HT-XRD experiments is highlighted shading some light on the complex thermo-structural behavior of copper vanadate.
Those results also make possible to define the optimal conditions for developing the PLD target. For sintering, the temperature is limited to 720 °C (start of melting of a small part of Cu2V2O7 oxide) which allows solidification/densification of the target which is necessary for the PLD technique. The α-Cu2V2O7 powder was then shaped into a pellet via uniaxial press (15 kbar), then sintered at 720 °C for 1 hour. To avoid the diffusion of the elements into the alumina plate, a platinum sheet was used during the sintering step. The pellet obtained is dark brown and the XRD pattern indicates as expected the presence of β-Cu2V2O7 phase (Fig. 4b). It is important to note that after use, this target remains stable in air for a few weeks and then starts to disintegrate gradually. This phenomenon can be explained by the negative thermal expansion coefficient which induces strong mechanical stresses when the target is subjected to temperature variations caused by the laser impact.
The presence of a β phase with a (200) preferential orientation is highlighted for the temperatures of 230 and 300 °C. Note that a parasitic phase of CuV2O6 (legend with the Δ symbol) is also present at this pressure. At 150 °C and 400 °C, no reflection is highlighted except the reflections of Indium Tin Oxide (ITO). A temperature above 150 °C is necessary for the β-Cu2V2O7 crystallization. In addition, at high temperature also leads to the production of an amorphous thin film.
Thus, after setting the optimal temperature at 230 °C, we studied the influence of the dynamic oxygen pressure in the chamber. Fig. 5b shows the XRD patterns for different pressures. For thin film synthesized at 10−5 mbar, the pure β phase can be obtained but with the disappearance of the preferential orientation. At 10−3 mbar, a β phase is obtained but with the presence of a parasitic phase of CuV2O6 (Δ). For higher pressure (10−1 mbar), the (200) reflection of β-Cu2V2O7 goes out completely, indicating either the presence of an amorphous film or no film at all. The surface of the film synthesized at 10−1 mbar was thus analyzed by EDS at different points showing a homogeneous V/Cu ratio close to 1, in favor of the first hypothesis.
It should be noted that a freshly prepared target is necessary for the development of thin films and beyond ∼20 deposition cycles, the target crumbles, probably due to the negative expansion coefficient of the compound, under the effect of the thermal stresses induced by the laser impact.
As observed on Fig. S11 (ESI†), the lower oxygen pressure is, the lowest band gap is. This results in an increasingly dark appearance visible to the naked eye of the deposited film. Thus, the band gap evolves from 2.2 to 2.4 eV for 10−5 and 10−3 mbar, respectively. The photo-electrochemical measurements presented in next paragraph will be performed on the film elaborated at 10−5 mbar and at 230 °C. Indeed, these optimal parameters lead to a pure β-Cu2V2O7 phase (without CuV2O6 impurity).
Gpe = Bcnpe(npe + M) ∝ ϕ0 | (2) |
(3) |
(4) |
Note that for higher intensities of the light beam, saturation can be achieved. This model therefore explains the experimental observations and the behavior of the electrode under different light power values.
The photocurrent transient response and especially the instantaneous pike and the subsequent decay towards the steady state current is due to the build-up of carriers in the surfaces states.53 Likewise, the current overshoot observed when the light is switched off is due to the recombination of the remaining carriers in the surface states. So, the charge transfer and the recombination rate constants can be obtained by modeling the transient photocurrent response.53,54 Indeed, the time constant of the decay or the overshoot is given by the expression (5):55
(5) |
Moreover, the ratio of the steady state photocurrent to the instantaneous photocurrent observed when the illumination is switched on is given by the expression (6):55
(6) |
These two expressions make possible to calculate the values of ktr and krec. Considering 450 nm illumination (ϕ0 = 23 mW cm−2), the values of ktr and krec are 6.19 × 10−2 s−1 and 5.58 × 10−2 s−1, respectively. The close values show that only half of the holes arriving at the surface are transferred to the solution under steady state conditions.
In addition, the capacitance versus potential (C–V) measurement is used to determine the conduction type, the flat band (Efb) and carrier density (N), which can be determined from the Mott–Schottky (MS) plot (1/C2 – V). Fig. 6d shows the MS plot of the β-Cu2V2O7 film performed at 1 kHz. The positive slope observed on the depletion zone indicates a n-type semiconductor behaviour. The Efb is estimated at −0.08 V (vs. Ag/AgCl) which correspond to 0.66 V (vs. RHE). This value is consistent with the values described in the literature.15 The flat band potential reflects the position of the Fermi level.56 For n-type semi-conductors, the Fermi level lies close to the Conduction Band (CB). Therefore, knowing the band gap value, we can estimate the position of both valence and conduction bands (inset Fig. 6d). Moreover, in the depletion region, the MS plot follows the equation:57
(7) |
Regarding thin films, chronoamperometry measurements were also carried out on β-Cu2V2O7 films produced by PLD. The photocurrent response to different wavelengths and for several increasing luminous flux intensity is shown in Fig. S12 (ESI†). We can also see that the best response is obtained under 450 nm illumination. Indeed, that drops by 66% and 86% for an excitation of 505 nm and 590 nm, respectively (for the same light intensity). Fig. 7a shows the variation of the current density between illumination and black (Δj) as a function of the intensity of the luminous flux (ϕ0) for an applied potential of 0.4 V vs. Ag/AgCl. We can see a progressive and linear evolution of Δj vs. ϕ0 (Fig. 7b) as well as an improved current density response compared to the bulk of the α-Cu2V2O7 phase. The photocurrent scales linearly with the illumination power (Δj = 3.9 × 10−3ϕ0) indicating an ideal trap-free state. β-Cu2V2O7 can therefore be used as a photodetector. In addition, the instantaneous photocurrent spike is less intense (compared to the previous experiment performed on α-Cu2V2O7) indicating a low recombination constant (compared to ktr) in the case of β-Cu2V2O7 thin films. ktr = 9.58 × 10−2 s−1 et krec = 2.39 × 10−2 s−1. Considering these values, the charge transfer dominates over the recombination if ktr/krec > 1 for the thin films.
Finally, Cu2V2O7 is also known to exhibit photocorrosion phenomenon.8 The stability of the films was thus studied by chronoamperometry in 0.1 M NaBi buffer solution (pH = 9.2) and with a constant applied potential of 0.4 V vs. Ag/AgCl (Fig. 7c). For the β-Cu2V2O7 thin films, one can see that the anode current density is relatively stable, then undergoes a slight decrease with time. This drops to 8% over 1 hour under 450 nm illumination (ϕ0 = 95 mW cm−2). This decrease is mainly due to a gradual dissolution of the photoanode in the electrolyte. Indeed, this deterioration can be induced by the formation of O2 bubbles at the interface and/or a reaction with the electrolyte because the zone of stability (pH) of the oxide is very small.59 This phenomenon has also been observed in the case of Cu2V2O7 films, e.g. I. Khan et al.16 indicate a decrease of 2% and 14% in 1000 s for Cu2V2O7 produced at 500 °C and 250 °C, respectively. In order to check this point, we performed an image of the surface of the Cu2V2O7 thin film and maps for the V, In and Cu elements between a zone not subjected or subjected to the electrolyte. In Fig. 8, we can clearly observe the demarcation zone between a part subjected to the electrolyte showing circular corrosion phenomena and a part corresponding to the initial thin layer having a smooth appearance whose surface is constituted of granular conglomerates. The absence of copper and vanadium in these cracks confirms the corrosion of the film during prolonged use. The effect of the electrolyte on the stability of the film is demonstrated and will require work on the formulation of the electrolyte in order to improve this stability. Nevertheless, the results obtained on the films of β-Cu2V2O7 produced by PLD shows that it can be used for applications as a particularly active photoanode or as photodetector in the blue spectral region. Nevertheless, the photostability needs to be improved in the futures studies. We can consider two paths: (i) improved the chemical stability with respect to the electrolyte or (ii) the encapsulation in order to protect the photoanode. For this last point, this method has already been used successfully to protect hybrid halide perovskite by a double-layer graphite sheet/nickel (GS/Ni).60
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ra01509b |
This journal is © The Royal Society of Chemistry 2023 |