Exploring the synergistic effects of goethite intercalated coal in the presence of humic acids for enhanced growth of Sinapis alba

Background Iron oxide mineral–humic complexes serve as a reservoir of bioavailable Fe for plants, releasing metal ligands and providing Fe–humic complexes directly usable by plant Fe-uptake mechanisms. In this study, we synthesized and characterized goethite α-FeOOH (G) nanoparticles (NPs) intercalated in coal (GC) to estimate the bioactivity effect of humic acids (HA). The synthesized GC NPs were characterized by X-ray diffraction, scanning electron microscopy (SEM), Mössbauer spectroscopy, N 2 adsorption–desorption Brunauer–Emmett–Teller (BET) specific surface area, zeta potential, hydrodynamic particle diameter, iron ions release, and a phytoassay method of root elongation using the higher plant Sinapis alba . Results X-ray diffraction revealed that G was the primary phase in both GC and GC–HA complexes. Mössbauer spectroscopy analysis identified a goethite-doped Fe 2+ -in the GC samples. The intercalation of G into the coal matrix increased the specific surface area of GC, enhancing its HA sorption capacity. In addition, GC–HA demonstrated superior plant growth stimulation compared to HA and GC alone, indicating its role in colloidal stability. In contrast to GC, GC–HA exhibited a more consistent and time-dependent release of Fe 3+ and Fe 2+ . This sustained Fe release from GC– HA, coupled with the formation of Fe 3+ and more bioavailable (soluble) Fe 2+ humic complexes is a promising result in terms of iron nanofertilizers production. Conclusions The use of goethite nanoparticles intercalated within a coal matrix and subsequently complexed with HA contributes to prolonged phytoactivity by employing slowly released nutrient additives within the coal mesoporous matrix.


Introduction
Iron (Fe) is an indispensable element for the majority of living organisms, playing a crucial role in cellular respiration, energy conversion, and protein biosynthesis [1].In agriculture, the development of safe and effective formulations of biologically available iron is of paramount importance.
In the pH range between 7 and 9, the predominant iron mineral forms in the soil are represented by Fe(OH) +  2 , Fe(OH) 3 and Fe(OH) − 4 hydroxocomplexes.Their concentration in the soil solution is approximately 10 -10 M.However, for normal plant growth, the concentration of dissolved iron should be several orders of magnitude higher, ranging from 10 -6 to 10 -5 M. Iron is an essential micronutrient crucial for various physiological processes in plants, and Fe deficiencies are frequently observed in crops grown in calcareous soils [5].Iron deficiency causes a plant disease known as iron deficiency chlorosis.To address this issue, soluble iron salts are sometimes applied, but they prove ineffective at low doses and lack long-term efficacy due to Fe ion hydrolysis and subsequent removal from the soluble soil phase through hydroxide precipitation.Alternatively, more effective iron chelates with synthetic organic chelators are used in certain cases [6].However, when regularly applied to soil, they accumulate [7], leach into natural waters, and, due to their high complexing ability, contribute to an increase in the mobility and migration of heavy metals and radionuclides in the environment.Using biologically available forms of iron stabilized by humic substances (HS) represents an environmentally friendly method, because HS perform a number of important biospheric functions.These include soil structuring [8,9], accumulation of nutrients and microelements in a form available to plants [10,11], and regulation of geochemical fluxes of metals in water and soil ecosystems [12][13][14].HS show the ability to form stable complexes with metal ions [15] and to stabilize soil colloids containing oxide nanoparticles [16], while providing a protective effect to organisms under stress [17][18][19].Their ability to bind metal ions makes HS suitable for incorporation in the production of microfertilizers and feed and food additives containing trace elements [20].Furthermore, their redox properties allow the use of HS as reducing agents that limit the mobility of ecotoxicants, such as Cr(VI) [21].
The potential of iron-containing humic preparations as an environmentally safe alternative to synthetic iron chelates was explored in [22].The study found that iron in humic matrices predominantly exists as highly dispersed hydrated oxides.The effectiveness of iron compounds with HS in alleviating iron deficiency in plants has been consistently demonstrated [23][24][25][26].For example, in the study by [24], nutrient solutions containing iron-humic complexes were used as iron source in plant experiments.The findings indicated that iron-humic complexes serve as a readily absorbable iron source for plants, facilitated by easy uptake by plant roots [25,26].The impregnation of mineral oxides offers a novel approach to modifying virgin biochar [27] due to the abundance and low cost of iron oxides such as magnetite, maghemite, and goethite [28][29][30].
This study aimed to develop and characterize goethite nanoparticles intercalated within a coal matrix (GC) and to assess the impact of HA on the biological activity of higher plants, particularly Sinapis alba (S. alba).The effect of nanoparticles on plant growth is influenced by various factors, including the type, source, concentration, and size of NPs, the plant species (including developmental stage and growth rate), and the duration of NP exposure to the plants.The size and surface charge of NPs make them potential candidates for plant uptake [5].The incorporation of HA into formulations containing iron oxyhydroxide nanoparticles intercalated into a coal matrix contributes to a prolonged effect by employing slowly released nutrient additives from the mesoporous coal matrix [31][32][33][34].In addition, this approach appears to enhance bioavailability and environmental friendliness, as these substances naturally form co-precipitates with mineral components (e.g., Fe oxides and clays) or exist in solution, where they significantly contribute to dissolved organic matter.

Coal treatment
Brown coal from the Min-Kush deposit in the Kyrgyz Republic, designated as "C", was used for both HA extraction and goethite (G) nanoparticle intercalation.Following crushing and screening, the raw coal was sieved and the 1 mm size fraction was retained.To eliminate any impurities from the pores, the G nanoparticles were rinsed with deionized water (DI).They were then homogenized to a size of 50 µm in a high-energy ball mill (SPEX SamplePrep 8000, Mixer/Mill, Tungsten Carbide Vial, Metuchen, NJ, USA) at 1060 rotations per minute (rpm) for 30 min and dried in a centrifugal spray dryer (LPG-5, China) with an inlet air temperature of 210 °C and an outlet air temperature of 105 °C.This treated G was used for both HA extraction and the preparation of GC.

HA extraction procedure from coal
The methodology described by Lowe (1986) with modifications was followed for the extraction of HA from coal.The modifications included pre-homogenizing and grinding the coal to reduce HA extraction time, as well as increasing the extraction temperature as described in [35].In this study, 20 g of the pre-sieved coal fraction with a particle size of 50 µm was placed in 1 L bottles, and 200 mL of a 1% NaOH solution was added with continuous stirring at 600 rpm.The bottles were sealed with rubber stoppers and equilibrated at 70 °C for 3 h.Following extraction, the dark supernatant was separated by centrifugation at 1000 rpm for 15 min.Subsequently, 200 mL of 10% HCl solution (approximately pH 2) was added to the supernatant, and the HA was allowed to stand at room temperature for 24 h.The supernatant was then siphoned off from the acidified extracts.
To purify the HA samples, they were washed three times with DI until a nearly neutral pH was achieved.Each washing step involved centrifugation at 10,000 rpm for 10 min.This process effectively removed mineral matter.Following the washing procedure, the purified HA was dialyzed, dried in a centrifugal spray dryer (LPG-50, China) at an inlet air temperature of 105 °C and an outlet air temperature of 50 °C, and stored in a desiccator over P 2 O 5 .

Goethite intercalated coal preparation
G were synthesized following the method described by Hiemstra et al. [36].Briefly, 9.8 g of Mohr's salt ((NH 4 ) 2 Fe(SO 4 ) 2 • 6H 2 O) was dissolved in 250 mL DI under vigorous stirring at 600 rpm.Subsequently, 50 mL of a 25% NH 4 OH solution was added, and the resulting suspension was stirred at 70 °C for 24 h.After this period, the precipitate, G, was thoroughly washed with DI until a neutral pH was attained.Finally, the precipitate was air dried at room temperature for 24 h.
Next, 1.25 g of G and 5 g of C were dispersed in 90 mL of DI.The reaction mixture was stirred at 25 °C for 1 h.Subsequently, the precipitate GC was isolated by centrifugation at 1200 rpm for 30 min.Finally, the precipitate was placed in a desiccator with P 2 O 5 until it was completely dry.

GC-HA complex preparation
To prepare a 2 g L −1 HA solution, 1 g of HA was dissolved in 450 mL DI, followed by the addition of 50 mL of 10% NaOH.Subsequently, GC-HA complexes were formulated at 1:1 and 1:10 GC:HA (w/w) ratios.The suspensions were gently shaken on a shaker at 150 rpm for 24 h at room temperature without pH adjustment.After the 24-h shaking period, the suspension was centrifuged at 6000 rpm for 15 min.The resulting precipitates were washed three times with DI to remove unbound HA fractions.Following this iterative process, the wash waters were clear and nearly colorless, indicating the successful removal of unbound HA.
Subsequently, the washed precipitates were dried at 60 °C under vacuum conditions.The GC-HA complexes were designated GC-HA1 and GC-HA10, where the index indicates the initial HA concentration (in g L −1 ) in the complex.Elemental analysis for Fe and carbon for HA showed that the resulting Fe/HA precipitates contained 14% and 6% Fe and 12% and 37% HA, respectively.

Structural characterization of GC and GC-HA
Size and morphology of the samples were observed under a scanning electron microscope (Tescan Vega 3, Czech Republic), and conductive coating (Pt) was applied using magnetron sputtering (Jeol JFC-1600).The coating thickness is no more than 10 nm.The degree of polydispersity was calculated using the modified formula (2.1) derived from Carlos De La Vega's work [37]: where Kp-degree of polydispersity, ∂-weight average diameter (calculated from the fraction of each particle), ∆-arithmetic average diameter.
X-ray diffraction (XRD) analysis of the nanoparticles was performed using a Philips X-pert Cu-Kα diffractometer (λ = 1.5418Å) in the angular range 2θ = 15-85° at a scanning rate of 5 deg min −1 and at 25 °C.Diffraction pattern construction and XRD data analysis were performed using Match! and OriginPro software.
Mössbauer absorption spectra for 57 Fe were obtained using an Express Mössbauer Spectrometer MS1104EM (CJSC Kordon, Rostov-on-Don, Russia) at temperatures of 296 ± 3 K and 77.7 ± 0.3 K.The γ-radiation source was 57 Co in a matrix of metallic rhodium at room temperature.The noise-to-signal ratio for the spectra did not exceed 2%.The mathematical processing of the experimental Mössbauer spectra was performed for high-resolution spectra (1024 points) using the program SpectRelax 2.8 (Lomonosov Moscow State University, Russia).The values of isomer shifts are given relative to α-Fe.
Specific surface area and porous structure analysis of the samples were determined using a Sorptometer-M (Katakon, Russia) at liquid nitrogen temperature (77 K).The Brunauer-Emmett-Teller (BET) and Barrett-Joyner-Halenda (BJH) methods were used to calculate the specific surface area and porous structure characteristics based on the adsorption-desorption isotherms.
To eliminate absorbed gases and vapors from the surface, the samples underwent a "thermal training" process prior to testing.This involved heating them in a stationary nitrogen stream in a vacuum at a temperature of 150 °C.

Surface charging characterization of HA, GC, and GC-HA
The zeta potential of the nanoparticles was determined using electrophoretic light scattering with a NanoBrook Omni particle analyzer.Measurements were performed at a wavelength of 633 nm using a solid-state He-Ne (1) laser, a scattering angle of 173°, and a temperature of 25 °C.Each sample was appropriately diluted in DI to achieve a concentration of 0.1 g L −1 for analysis.Consistent time intervals were maintained for all measurements: 10 s of dispersion in an ultrasonic bath followed by an additional 100 s of equilibration.The experiments were performed at 25 ± 0.1 °C in BI-SCP disposable polystyrene plastic cuvettes.Small amounts of 0.1 M HCl or NaOH were added to the suspension to adjust the pH to a range between 2 and 10. pH adjustments were monitored using a pH electrode (Mettler Toledo InLab Expert Pro-ISMO) and a pH meter (SevenExcellence S400).

Elemental analyses (C, H, N, Fe)
Elemental analysis of the samples was performed using an elemental analyzer from Carlo Erba Strumentazione.
To determine the ash content, samples were subjected to heating at 800 °C for 4 h.The oxygen content was calculated as the difference between the total mass and the sum of the remaining elements.The H/C and O/C atomic ratios were derived from the elemental contents calculated on an ash-free and moisture-free basis.

Determination of total acidity of HA
Total acidity determination of purified HA was determined following the method described by Inbar et al. [38].
In this procedure, a 5-10 mL aliquot of HA solution containing 5-20 mg HA was transferred to a vial (~ 22 mL), followed by the addition of 10 mL of 0.03 M Ba(OH) 2 .The vial was tightly sealed, shaken thoroughly, and allowed to equilibrate for 24 h at room temperature.
After the equilibration, aliquots of the clear solution above the barium humate precipitate were transferred to a titration cell and titrated with a standard HCl solution (0.1 M) using phenolphthalein as an indicator.The total acidity (TA, mmol g −1 ) was calculated according to formula (2): where V 0 and V HA are the volumes of HCl consumed for the blank and sample titrations (mL), respectively, C HCl is the titrant concentration (mmol mL −1 ), and m is the mass (g) of HA in the aliquot.
A saturated Ba(OH) 2 solution was carefully prepared by dissolving BaO in CO 2 -free DI (boiled for 1 h) in a sealed volumetric flask with vigorous shaking.The solution was allowed to stand for 3-4 days until complete precipitation of BaCO 3 occurred.Working solutions were freshly prepared just before analysis by diluting an aliquot of the transparent supernatant and standardizing against HCl. (2

Determination of carboxylic acidity of HA
The Ca acetate method, as described by [39], was used to determine the content of strong acid groups in the humic samples.In this method, a 5-10 mL aliquot of HA solution containing 5-20 mg HA was transferred to a vial (22 mL), followed by the addition of 10 mL of 0.6 M Ca(CH 3 COO) 2 .The vial was tightly sealed, shaken thoroughly, and allowed to equilibrate for 24 h at room temperature.
After equilibration, aliquots of the clear solution above the precipitate of Ca-humates were transferred to a titration cell and titrated with a standard NaOH solution (0.05 M) using an autotitrator.The carboxyl acidity (CA, mmol g −1 ) was calculated according to formula (3): where V 0 and V HA are the volumes of NaOH consumed for the blank and sample titrations (mL), respectively, C NaOH is the titrant concentration (mmol mL −1 ), and m is the mass (g) of HA in the aliquot.

Size exclusion chromatography (SEC) of HA
SEC analysis was performed according to [40] using a column (25 mm × 20 cm) packed with Toyopearl HW-50S gel (Toso Haas, Japan).Polydextrans were used as markers for molecular weight calculations.Prior to analysis, HA solutions were adjusted to a concentration of 1-2 mg C L −1 by equilibration with the SEC mobile phase (0.028 M phosphate buffer, pH 6.8).The flow rate was set at 1 mL min −1 , and the absorbance of the eluate was monitored at 254 nm.Peak molecular weight (Mp) is the molecular weight at the peak of the distribution curve.
Number average molecular weight (Mn) is defined by formula (4): where n i is number of ith molecules with molecular weight M i .
The weight average molecular weight is calculated according to the following Eq.( 5):

Characterization of the HA
Information on the elemental composition, carboxylic and phenolic group content, and molecular weight distribution of HAs is summarized in Tables 1 and 2.

Detection of Fe ions release
The mass concentration of released Fe 3+ and Fe 2+ ions was determined following the methods with modification described in [41,42].GC and GC-HA powders were suspended in 0.1 M sodium acetate buffer (pH = 4.6, 0.6% NaCl) and subjected to centrifugation (6000 rpm, 5 min) at intervals of 0, 0.5, 1, 3, and 24 h.After separation, specific ion detection reagents were added to the supernatant.For Fe 3+ ions detection, 200 µL of potassium thiocyanate (KSCN) (50% solution) and 200 µL of HCl (18.25% solution) were added to 5 mL of the supernatant, initiating a reaction under strongly acidic conditions with a pH close to 2. The solution was allowed to equilibrate for 20 min, and the absorption spectrum was measured at 480 nm.
For Fe 2+ ions detection, 2 mL phenanthroline (C 12 H 8 N 2 -H 2 O) (2.5% solution) and 600 µL of sodium acetate were added to 5 mL of the supernatant.Similarly, the solution was allowed to equilibrate (with no observable color change) for 20 min and inspected at 510 nm.

Phytoassay method
The root elongation phytoassay method, widely used for quality assessment of metal-bearing samples, is based on various protocols (US EPA, 1996; ISO 11269-1, 2012) [43,44].This method is favored due to its cost effectiveness, simplicity, and short duration [43].
In this study, white mustard (Sinapis alba L.) was selected for the phytoassay to investigate the doseresponse relationship of aqueous suspensions of iron-based nanoparticles in a 96-h root elongation phytoassay test, in accordance with the recommendations of ISO 18763:2016 and the Russian standard method (Russian Federal Register FR 1.31.2012.11560)[44].
The experimental setup comprised specially designed boxes, each containing two compartments measuring 13.5 × 8.5 × 0.8 cm (length x width x height), corresponding to a volume of 92 cm 3 per compartment.The lower compartment, designated for root growth, was filled with filter paper soaked in an aqueous suspension of nanoparticles, while the upper compartment accommodated shoot growth.
The filter paper in the lower compartment was evenly moistened with 10 mL of the aqueous nanoparticle suspension.Subsequently, ten Sinapis alba L. seeds were placed on the filter paper.The test boxes were sealed with a transparent lid and initially incubated horizontally at 20 ± 2 °C in the dark for 24 h, followed by 72 h in a vertical position at 24 ± 2 °C under illumination for 16 h per day, maintaining a light intensity of 4000-7000 lx (light wavelength 400-700 nm, universal white).
At the conclusion of the incubation period, the length of the main root of the mustard seedlings was measured and the mean value was calculated and compared to the control values.The assay was performed in triplicate.

Statistics
The phytoassay effect of tested dilutions compared to the control was calculated as a percentage (%).Statistical processing of phytoassay results was performed using one-way ANOVA; p-values were calculated using the Statistica 10 statistical package (StatSoft Inc., USA).

X-ray diffractograms, SEM observation, and Mössbauer spectra of GC, GC-HA1, and GC-HA10
The X-ray diffractograms of GC, GC-HA1, and GC-HA10 exhibit narrow and symmetrical peaks, indicating the presence of crystalline material (Fig. 1).The noise-to-signal ratio was relatively high due to the carbon matrix.The primary structural parameters of the nanoparticles were determined and are summarized in Table 1.The phase states of the nanoparticles were determined using the Crystallography Open Database.
Table 3 presents a quantitative analysis of the XRD data, with a focus on the unit cell parameters determined by the Rietveld method.The primary phase identified in the synthesized NP samples matches the predicted goethite α-FeOOH phase, corroborating the available data [45].In addition, reflections corresponding to maghemite impurity were detected [46].The introduction of HA did not alter the goethite phase; maghemite also remained present.Slight variations in the average nanoparticle sizes were observed upon the addition of a high concentration of HA (10 wt.%).The sizes of the coherent scattering regions were 21 ± 4, 25 ± 3, and 37 ± 9 nm for GC, GC-HA1, and GC-HA10 NPs, respectively.
The size of the coherent scattering region was determined using the Scherrer method based on powder XRD data.Changes in nanoparticle size due to coating were investigated by SEM and XRD analysis, showing a correlation between the two techniques.Despite maintaining a  consistent spherical particle shape across all modification routes, subtle particle growth was observed (Fig. 2).
The bare GC particles had an average particle diameter of 21.0 nm and 65.5 nm according to XRD and SEM analyses, respectively.Following modification, the GC-HA1 particles exhibited an increase in size to 25.0 nm (XRD) and 66.7 nm (SEM) compared to bare GC.Notably, the GC-HA10 particles exhibited a more substantial increase in diameter to 37.0 nm (XRD) and 112.1 nm (SEM) compared to bare GC.
The larger crystallite size observed in HA-modified NPs compared to bare GC was attributed to the higher concentration of HA.Moreover, the addition of 10% HA induced particle structuring or an increased degree of crystallinity, characterized by a more pronounced contour and a strictly spherical shape, as evident in the SEM images.On the other side, supramolecular associates of HA, self-assembled due to hydrophobic interactions and hydrogen bonds of numerous low-molecular components [47,48], are also fractioned since the external forces from mineral surfaces on some small molecules exceed the intermolecular forces between the molecular components of HA [49].As a result, the well-ordered structures of crystalline are observed for GC-HA1 and GC-HA10.
Based on the coefficient of variation (CV) and polydispersity index (Kp) values obtained from the SEM analysis, it is noteworthy that GC-HA10 has a narrow size distribution (13.6%, 0.98) compared to GC (26.6%, 0.95) and GC-HA1 (18.3%, 0.97).Although all samples fall under the polydisperse category according to [37], the CV of 13.6% and Kp of 0.98 for GC-HA10 indicate a higher degree of homogeneity in the nanoparticles.This homogeneity is essential for achieving optimal performance in the final surface-activated material.
The addition of HA had a moderate effect on particle size but also influenced the lattice parameters of GS, as evidenced by XRD and Mössbauer spectroscopy.
Relying solely on XRD for identifying the crystal structure of iron (oxyhydr)oxides is insufficient because a significant portion of the material might be X-ray amorphous.Therefore, Mössbauer spectroscopy was used to complement the structural analysis of iron (oxyhydr)oxides.Given the remarkable similarity in XRD patterns for GC and GC-HA1, regardless of HA addition, Mössbauer spectral analysis was exclusively performed for GC and GC-HA10 samples.
The experimental Mössbauer spectra acquired at 296 K for GC and GC-HA10 samples exhibit a complex array of more than a dozen resonance lines with varying intensity, width, and profile (Fig. 3), indicating the intricate composition of the studied materials.The distortion in the profile of the resonance lines of the sextets towards the inner part of the spectrum reveals the small size of the corresponding magnetic domains, a characteristic feature of nanoscale materials.
These circumstances contribute to the complexity of describing the Mössbauer spectra model.A satisfactory mathematical description of the spectral profile could only be achieved by employing a superposition of two symmetric doublets and two nested sextets, each characterized by different probability distribution functions for hyperfine parameters (Table 2).Specifically, for the inner sextet, only the distribution probability of the hyperfine magnetic field was considered, while for the outer sextet, the probabilities of changes in both the magnetic field and the isomer shift were considered, and these were linked to the same distribution function (Table 2).
The probability distribution function for the inner sextet's magnetic field exhibits a unimodal, highly dispersed, and asymmetric profile due to a substantial broadening towards the low-field region (Fig. 3).The hyperfine parameters obtained for this sextet conclusively assign it to goethite [50,51].
Examination of the probability distribution function profiles for the outer sextet's isomer shift and magnetic field reveals two distinct modes.These modes correspond to iron atoms in octahedral (large isomer shift and smaller magnetic splitting-Table 4) and tetrahedral (smaller isomer shift and larger magnetic splitting) sites of non-stoichiometric magnetite-Fe 3-δ O 4 [52,53].The  2, subspectrum 3 and 4, respectively) in octahedral oxygen sites [54].
Comparing the high-temperature spectra of the samples, a noticeable broadening of the doublets responsible for the paramagnetic part of the spectrum is observed for GC-HA10 compared to GC.This broadening may indicate a degree of disorder in the corresponding crystalline fragments.In addition, the modes on the profile of the probability distribution function of the magnetite subspectrum for the GC-HA10 sample exhibit substantial dispersion and are poorly resolved.This could be attributed to the partial oxidation of the magnetite phase in the GC-HA10 sample by HA [55,56].
Upon cooling the samples to liquid nitrogen temperature, the Mössbauer spectra undergo a remarkable transformation: a sextet, albeit distorted, becomes clearly visible, and the intensity of the spectra increases approximately twofold (Fig. 4).For both samples, a satisfactory description of the experimental spectra is achieved by using a superposition of two symmetric doublets, similar to those described for the hightemperature spectra, and a sextet characterized by the probability distribution function of the quadrupole displacement and the hyperfine magnetic field (Table 2).This function comprises two modes with a regular symmetric profile, presumably corresponding to two ironcontaining phases (Fig. 4).
The more intense mode, characterized by substantial quadrupole shift and minimal hyperfine magnetic field, is attributed to the goethite phase [50,51].The second, less intense mode, exhibiting a quadrupole shift close to zero and a high hyperfine magnetic field, reflects the state of the iron atoms in non-stoichiometric magnetite.Based on the relative intensity of the corresponding mode, the magnetite content in the material is estimated to be 20-30%.
Therefore, Mössbauer spectroscopy revealed that the investigated samples are predominantly composed of goethite, with non-stoichiometric magnetite identified as the primary impurity.
The Mössbauer spectra of GC-HA10 corroborate that iron in the nodules primarily resides in the ferric oxyhydroxide form, constituting a goethite doped with Fe 2+ .HA are known to exhibit redox reactivity and the ability to chemically reduce metals, including Fe 3+ [57,58].Studies have demonstrated that both dissolved and solid-phase HA can accelerate the reduction of Fe(III) oxide in sediments [59,60] and facilitate the  bioreduction of Fe(III) minerals in soils [61] by acting as electron shuttles between bacteria and oxide surfaces.

Textural characteristics of GC, GC-HA1, and GC-HA10
Figure 5 presents the N 2 adsorption/desorption isotherms for nanoparticles.The samples exhibit type IV isotherms, characteristic of polymolecular adsorption and the presence of capillary condensation in mesopores.
The isotherms show H4 type hysteresis loops, characterized by a pronounced steep rise at low pressures, indicative of micropores.The results of the BET method for specific surface area (SSA) and the BJH method for average pore diameter are presented in Table 5.The data consistently show a decrease in SSA from 289.5 to 210.5 m 2 g −1 and a decrease in BET pore volume from 0.28 to 0.23 cm 3 g −1 .This suggests that higher HA concentrations during sorption lead to the blocking of more micropores.Interestingly, the pore diameter remains unchanged, probably due to the absence of HA adsorption on the pore walls, possibly attributed to the negative charge on the inner surface of the pores.

Zeta potential and hydrodynamic diameter of HA, GC, GC-HA1, and GC-HA10 as a function of pH
Zeta potential measurements provide valuable insights into the surface electrical properties of materials.The pH-dependent surface charging behavior of iron (oxyhydr)oxides, which significantly influences their affinity for anions and cations, is well documented in the literature [62].Figure 6 illustrates the pH-dependent variations in zeta potential and hydrodynamic particle diameter for HA and GC, considering GC both in its individual state and when combined with HA.
With increasing pH, the dissociation of surface OH groups induces surface recharging, resulting in a consequential increase in the negative charge of iron nanoparticles: The isoelectric point of GC in the absence of HA is observed at a pH (pH IEP ) of approximately 5.5, which is lower than most reported values for pristine goethite [62] owing to the influence of the carbon matrix.At pH values below 5.5, the GC surface is positively charged, transitioning to negative at higher pH values.
The speciation of HA in solution depends on ionic strength and pH.Initially, with increasing pH, carboxylic groups dissociate protons.Around neutral pH, a significant portion of these carboxylic groups have already dissociated, prompting the gradual dissociation of phenolic hydroxyls [63].
In the presence of HA, a noticeable shift in the pH isoelectric point is observed, from 5.3 for initial particles to 4 and 3.2 for 1 and 10 wt/wt GC-HA, respectively.This shift is consistent with classical specific adsorption of anions (in this case, HA polyanions), resulting in a lower pH IEP .
The reduction in negative charges on the G surface in the presence of 10 wt% HA polyanions can be attributed to the conformational structure of the polyelectrolytes [64].At higher concentrations of polyelectrolytes, there is a tendency towards a twisted structure, exposing fewer COOH/OH ions on the HA surface [64].
The decrease in zeta potential observed in the pH range 4.5-8 for the GC-HA10 sample suggests an aggregation process, further supporting the change in HA conformation with increasing concentration.
The hydrodynamic particle diameter of the samples was determined at pH 6.5 in DI, corresponding to phytoassay conditions.Under these conditions, electrostatic repulsive interactions between HA molecular segments contribute to an apparent increase in the average hydrodynamic diameter, reaching 460 nm at pH 6.5 [65].
The average hydrodynamic particle size for GC remained relatively unchanged after modification with 1 wt% HA, measuring approximately 300 nm for both samples.However, increasing the concentration of HA in GC to 10 wt% induced particle flocculation, resulting in a reduction in hydrodynamic diameter to about 150 nm.This phenomenon can be attributed to an enhanced Fig. 6 a Zeta potential (mV) vs pH and b hydrodynamic particle diameter (nm) at pH 6.5 for HA, and for GC with and without HA charge neutralization process arising from the excess fraction of protonated OH groups of HA (pKi 2 ~ 9 [66][67][68]), which do not directly participate in the coordination with G [69].

Phytoassay
The biological activity of iron nanoparticles was evaluated using the plant (S.alba) root elongation method, with the results presented in Fig. 7.The suspensions were maintained within the biologically and environmentally relevant pH range of 6-7.Statistical analysis using one-way ANOVA revealed a highly significant difference between the mean values of the groups (p-value < < 0.05).
The phytoassay revealed a subtle, wave-like stimulation of plant root elongation upon application of GC preparation at low concentrations (10-10 2 mg L −1 ).However, increasing the GC concentration to 104 mg L −1 resulted in a decrease length of plant root.These results are consistent with previous studies [70] that have demonstrated diverse bioeffects, including responses under stress conditions.In addition, our study highlighted the concentration-dependent toxic effects of iron oxide nanoparticles attributed to oxidative stress [70].
HA, known for their plant growth-promoting properties [19,71], consistently exhibited a stimulatory effect on root length growth, reaching up to 30% compared to the control.
The highest root length stimulation was observed for GC-HA1 and GC-HA10 compared to HA and GC alone.Notably, despite an order of magnitude increase in the concentration of HA in these samples, no change in bioactivity was observed.This suggests that the enhanced stimulation of mustard root growth by hybrid complexes GC-HA1 and GC-HA10 compared to HA and GC alone is likely due to the influence of Fe ions, including their presence as part of water-soluble complexes with HA.To validate this proposed mechanism, UV-Vis spectroscopy was used to investigate the kinetics and concentration of Fe 2+ and Fe 3+ released from the samples.

Release of Fe ions
UV-Vis spectroscopy was used to investigate the kinetics and concentration of Fe 2+ and Fe 3+ released from the samples (Fig. 8).Considering that the phytoassay results were evaluated after 96 h at pH ~ 6.5, the release of Fe 2+ and Fe 3+ within the same 96-h time frame at the specified pH was determined.Mössbauer spectroscopy data corroborated the presence of Fe 2+ in the samples.The concentrations of iron ions per gram of sample were recalculated.
Complexometric titrations with potassium rhodanide for Fe 3+ and potassium o-phenanthroline for Fe 2+ (Fig. 8) revealed that both GC and HA-modified complexes release both Fe 2+ and Fe 3+ .However, the concentration of ions released by the samples varies.GC nanoparticles released the highest concentration of Fe 3+ and Fe 2+ compared to those modified by HA.This dissolution of iron oxide nanoparticles is a well-documented phenomenon [72].
The kinetics of total ion release varied across all samples.Notably, the GC sample exhibited continuous dissolution over time, failing to reach equilibrium within 96 h.A rapid initial release of Fe 3+ and Fe 2+ from the GC sample was observed within the first hour, followed by a sharp decrease in concentration.This observation suggests a possible reverse sorption of ions or the formation of iron hydroxides [73,74].
The GC-HA1 and GC-HA10 complexes exhibited a sustained, low-level release of Fe 3+ and Fe 2+ over 96 h.The slightly reduced detectable release of Fe 3+ and Fe 2+ ions for the GC-HA10 sample is likely attributable to the formation of multiple surface bonds of COOH groups of HA and the OH groups of the GC surface via ligand exchange [75].
The dissolution rate of the studied samples also bears a proportional relationship to the particle surface area.This observation aligns with the Noyes-Whitney equation [76] and the Ostwald-Freundlich equation [77], both of which suggest that smaller nanoparticles dissolve more rapidly than larger ones.

Discussion
The XRD analysis results indicated that the crystal structure of GC-HA complexes did not exhibit significant alterations with increasing HA content compared to GC.However, the peak intensity of the complexes decreased, and the average nanoparticle size increased for GC-HA10.Mössbauer spectroscopy revealed the presence of Fe 2+ following the modification of GC with HA.The reduction of Fe 3+ occurred due to the well-known redox properties of HA [57,58].
In this study, Fe 3+ and Fe 2+ release was detected using the complexometric method to estimate Fe mobilization into the aqueous solution at neutral pH.The most abundant Fe 3+ species in the environment has low bioavailability [78] compared to the Fe 2+ form, which has better solubility under certain conditions and can be more easily absorbed by plants [79].
The observed reduction in hydrodynamic diameter and zeta potential in aqueous suspension of GC-HA10, determined by DLS and EDS methods, can be attributed to conformational changes in HA as previously documented in the literature [64].Importantly, this change in properties results in a bioactivity effect comparable to that observed with GC-HA1.
Moreover, the intercalation of G into the coal matrix resulted in a significant increase in the surface area of GC.The specific surface area of G, as measured by the BET N 2 method, is reported to be 94 m 2 g −1 [80].While the modification of GC with 10 wt% HA resulted in a decrease in surface area from 328.3 m 2 g −1 to 210.5 m 2 g −1 , it is noteworthy that the specific surface area of GC-HA remained higher than that of G-NPs alone.Weng et al.) [80] has previously calculated that the surface sites of G that form chemical bonds with HA are generally small.The adsorption of HA on these surface sites may physically block them, potentially limiting the formation of surface complexes between small ions and oxides [80].Crucially, the increased surface area resulting from the intercalation of G into the coal matrix was identified as a key factor contributing to the increased adsorption of HA.Consequently, the increase in specific surface area achieved by using coal as a matrix for G intercalation defines the extended surface sites and consequently the reactive behavior of GC-HA in this study.
The bioactivity of hybrid complexes containing iron oxides in a coal matrix, along with HA, is governed by intricate mechanisms that can vary depending on specific environmental conditions.In our model experiments, we evaluated the influence of Fe mobility on plant growth by measuring the release of Fe ions from the studied samples GC-HA1 and GC-HA10 using the complexometric method.However, the low percentages of Fe 3+ and Fe 2+ determined in our study likely resulted from the use of different extraction conditions compared to the real conditions encountered by seeds during germination.Specifically, we employed potassium thiocyanate for Fe 3+ and potassium o-phenanthroline for Fe 2+ in distilled water, whereas actual root exudates contain organic acids capable of desorbing iron from humic complexes.These organic acids possess high stability constants for iron complexes [81][82][83], enabling them to effectively displace the metal from the chelating centers of HS.Therefore, the extraction capacity of potassium thiocyanate and potassium o-phenanthroline may be significantly lower than that of the root exudates in natural environments.
Previous studies, including our own [84], have demonstrated that Fe ions at certain concentrations can trigger the Fenton reaction, leading to oxidative stress and plant growth inhibition.Both extracellular and intracellular reactive oxygen species (ROS) can be generated by iron-based NPs, with the former activating the latter [85].Elevated ROS levels induce oxidative damage and toxicity, as evidenced by increased malondialdehyde levels and reduced catalase and glutathione expression.The effects of ROS on superoxide dismutase activity vary across different studies [86][87][88].Recently, Chaithawiwat et al. [89] demonstrated that E. coli mutants lacking oxidative stress defense genes were more susceptible to nano-zero valent iron (nZVI) than the wild type, highlighting the crucial role of ROS-mediated oxidative stress in the toxicity mechanism of iron-based NPs at the molecular level.In addition to their negative effects, ROS also play a significant positive role in dormancy release, seed germination signaling, protection against pathogens, and regulation of internal cellular machinery in response to external environmental dynamics.There is a specific "oxidative window" that allows cellular events to unfold in sequential order for seed germination if ROS are maintained within a particular range [90].
Moreover, it has been shown [91] that nanomaterials promote seed germination by forming nanopores in seed coats, introducing reactive oxygen species (ROS), increasing enzyme activity at starch-degrading sites, and introducing ROS to the seed coat.
Therefore, as bioactivity of GC-HA1 did not differ from that of GC-HA10 though the latter released significantly more ferrous ions, an acquisition of iron by seedlings seems not to be important from this point of view.In addition, extra-data on iron beneficial effect on seed germination are needed.On the other hand, comparing EDS and phytoassay data, one can conclude on similarity of GC-HA1 and GC-HA10 in terms of zeta-potential dependence on pH.Data presented in Fig. 6a show higher colloidal stability of nanoparticles modified with HA at pH 6-7 as compared to GC. Besides, additional steric stabilization of GC-HA1 and GC-HA10 NPs as compared to GC can be expected.So, overall, studied GC-HA1 and GC-HA10 demonstrated similar bioactivity and colloidal stability.As for the other properties (XRD, SEM, Moessbauer, etc.) GC-HA1 and GC-HA10 differed one from another.Therefore, it is most logical to look for an explanation of the observed biological activity in the colloidal properties of the nanoparticles under study.Higher stability of nanoparticles with HA at pH 6-7 can be considered as a leading reason of their more pronounced biological activity.A further study of the synthetized NPs' aggregation and ROS generation and their relationship to biological activity should be proposed.

Conclusions
We present the synthesis and characterization of goethite (G) nanoparticles intercalated in coal (GC) and modified by humic acids (HA) to enhance their bioactivity towards Sinapis alba.The synthesized GC was identified and characterized by scanning electron microscopy, X-ray diffraction, Mössbauer spectroscopy, N 2 adsorption-desorption BET specific surface area, zeta potential, release of iron ions, and phytoassay with higher plants.
The primary phase of the goethite structure was identified by X-ray diffraction in both the GC and goethite-humic acids (GC-HA) complexes.Mössbauer spectroscopy analysis revealed a mixture of Fe 2+ -doped GC, possibly due to the redox potential of HA.The intercalation of GC in the coal matrix increased the specific surface area of GC, facilitating HA sorption.In addition, HA was found to promote plant growth stimulation due to colloidal stability of hybrid nanoparticles.
Compared to GC, the release of Fe 3+ and Fe 2+ from GC-HA was prolonged and more time-dependent.The revealed ferrous ion release from GC-HA10 can be considered as a promising result in terms of iron nanofertilizers production.

Fig. 3 4
Fig.3Experimental Mössbauer spectra obtained at 296 K, as well as models for their description and probability distribution functions of hyperfine magnetic fields for GC (a) and GC-HA10 (b) samples

Fig. 4 Fig. 5
Fig. 4 Experimental Mössbauer spectra obtained at 78 K, as well as models for their description and probability distribution functions of hyperfine magnetic fields for GC (a) and GC-HA10 (b) samples

Fig. 7 Fig. 8 a
Fig. 7 Dose-effect relationship of plants root length for different samples.All concentrations are nominal.Average values ± SD of in triplicates are shown

Table 1
Elemental, carboxylic, and phenolic group composition of HA

Table 2
Peak molecular weight Mp, number average molecular weight Mn, weight average molecular weight Mw, and polydispersity Mw/Mn of HA

Table 3
Quantitative evaluation of XRD data by the Rietveld method and SEM data for GC, GC-HA1, and GC-HA10