Conductivity-dependent Flow Field-flow Fractionation of Fulvic and Humic Acid Aggregates

Fulvic (FAs) and humic acids (HAs) are chemically fascinating. In water, they have a strong propensity to aggregate, but this research reveals that tendency is regulated by ionic strength. In the environment, conductivity extremes occur naturally—freshwater to seawater—warranting consideration at low and high values. The flow field flow fractionation (flow FFF) of FAs and HAs is observed to be concentration dependent in low ionic strength solutions whereas the corresponding flow FFF fractograms in high ionic strength solutions are concentration independent. Dynamic light scattering (DLS) also reveals insight into the conductivity-dependent behavior of humic substances (HSs). Four particle size ranges for FAs and humic acid aggregates are examined: (1) <10 nm; (2) 10 nm–6 µm; (3) 6–100 µm; and (4) >100 µm. Representative components of the different size ranges are observed to dynamically coexist in solution. The character of the various aggregates observed—such as random-extended-coiled macromolecules, hydrogels, supramolecular, and micellar—as influenced by electrolytic conductivity, is discussed. The disaggregation/aggregation of HSs is proposed to be a dynamic equilibrium process for which the rate of aggregate formation is controlled by the electrolytic conductivity of the solution.


Introduction
Natural organic matter (NOM) is to the environmental sciences what proteins, nucleic acids, lipids, and carbohydrates are to the life sciences.All of these materials possess unique properties determined by their size, folding, and patterns at the nanoscale.The study of NOM is a very old, yet currently provocative, topic.Humic substances (HSs) comprise a major fraction of NOM, and reportedly represent up to 60%-70% of total organic carbon (TOC) in soils and 60%-90% of dissolved organic carbon (DOC) in natural water [1].Soil scientists have studied the chemical nature of HSs for two centuries, and their composition and conformation have been intensely researched.However, because of their complex nature, the primary, secondary, tertiary, and quaternary molecular arrangements of the fulvic (FAs) and humic (HAs) comprising TOC/DOC still remains ill defined.
This research is based on the hypothesis that multiple types of disaggregated/aggregated/ coiled/extended/micelles/pseudomicelles/micelle-like/vesicles/supramolecular/fractal/nanogel/microgel conformations of HSs exist, depending on concentration and solution physicochemistry, such as pH, ionic strength, and the nature of counterions in solution.Ultimately, increased understanding of the relationship between NOM composition, conformation, aggregation, and chemical reactivity will benefit society by improving potable water production and wastewater treatment of emerging contaminants.
This manuscript attempts to address a scientific conundrum: Why is the flow field-flow fractionation (FFF) of FAs and HAs concentration independent at high ionic strength and concentration dependent at low ionic strength?The aim of this paper is to knit together seemingly disparate puzzle pieces from experiments in flow FFF and dynamic light scattering (DLS) to postulate an answer to this quandary.
The answer to this question hinges upon our fundamental understanding of the aggregation of HSs in solution, and conversely, our understanding of the aggregation of HSs in solution informs the answer to this question.The purpose of this manuscript is to focus on HS aggregation and not to get into the terminology fray of whether to call the aggregates extended/coiled/micelles/pseudomicelles/ micelle-like/vesicles/supramolecular/fractal/nanogel/or microgel entities.In fact, variable types of these entities may coexist [2].The data presented here indicate the presence and size of aggregates and is not designed to absolutely define their character.Therefore, the terminology "aggregates" and "aggregation phenomena" will be preferred with two exceptions: (1) the terminology used by researchers in a cited reference will be restated and used without judgment; and (2) the term "micelle" will be retained in describing the physical measurement referred to as the "critical micelle concentration".

Experimental Section
This manuscript provides enhanced data analysis of previously published investigations; therefore, detailed experimental procedures exist in each of the publications cited [3,4].Briefly, the channel dimensions of the apparatus used for symmetrical flow FFF were 28.5 cm (tip-to-tip length) × 2 cm (breadth) × 0.0508 cm (width).The membrane consisted of polypropylene-backed polysulfone having a molecular weight cutoff of 10,000.
The low conductivity mobile phase was composed of HPLC-grade water containing a surfactant, FL-70, 0.05% (Fisher Scientific, Fair Lawn, NJ, USA) and sodium azide, 0.03% (pH 7) (Fisher Scientific).The channel flow was 0.6 mL•min −1 , while the cross flow was 3.1 mL•min −1 .The stop flow was initiated after 15 s and terminated after 2 min.The high conductivity mobile phase consisted of 0.095 M dibasic potassium phosphate, 0.005 M monobasic potassium phosphate, and 0.03% sodium azide in HPLC grade water (pH 8.1).Flow FFF was conducted with a cross flow rate of 0.4 mL•min −1 , a channel flow rate of 0.2 mL•min −1 , and a stop flow period of 7.5 min.The stop flow was initiated after 2.5 min and terminated after 10 min.
The hydrodynamic diameter of particles was calculated for normal mode flow FFF using Equation (1): where R is the retention ratio equal to the void time of the channel divided by the emergence time of the species, Vc is the channel volumetric flow rate and V is the void volume, w is the channel width, k is the Boltzmann's constant, T is the absolute temperature, d is the hydrodynamic diameter of the diffusing particle, and η is the solution viscosity.The hydrodynamic diameter is calculated for particles eluting in steric (lift) mode flow FFF using Equation (2): The diffusion coefficient (D) is related to the hydrodynamic diameter by the Stokes-Einstein equation: assuming the aggregate is spherical.

Results and Discussion
Concentration independence in an analytical procedure is a requirement for quantitative measurement.Consequently, concentration dependence is anathema to analytical chemists.Previously, the concentration dependence of HSs in flow FFF at low ionic strength [5,6] and the concentration independence at high ionic strength [4] was observed (Figure 1).Lead et al. [7] also demonstrated the trend toward a less broad and more Gaussian distribution of the flow FFF fractograms of Suwannee River FA as the ionic strength of the mobile phase increased.To investigate the conductivity-dependent behavior of HSs and postulate an answer to the question-Why is the flow field-flow fractionation of fulvic and humic acids concentration independent at high ionic strength and concentration dependent at low ionic strength?-various,seemingly unrelated, experiments in flow FFF and DLS at low and high ionic strength will be examined and compared.

Flow Field-Flow Fractionation (Flow FFF) of Humic Substances (HSs)
In these experiments, the flow FFF membrane consisted of polypropylene-backed polysulfone having a molecular weight cutoff of 10,000 Da.However, FAs and HAs ranging in weight-average molecular weight from 750 to 3400 were retained in the system [3,4].Monitoring the cross-flow revealed no detectable loss of analytes through the membrane.The negatively charged polysulfone membrane is postulated to repel the negatively charged polyelectrolytic humic macromolecules causing them to be retained in the channel.

Low Electrolytic Conductivity Solutions
In low ionic strength solution, the retention times of Aldrich HA (Aldrich Corporation, Milwaukee, WI, USA) and polystyrene sulfonate (PSS, MW = 5400) (Polysciences, Inc., Warrington, PA, USA) decreased as concentration increased (Figure 2).Additionally, at concentrations of approximately 5 mg/mL and greater, the retention order of the HA and the PSS is reversed relative to that observed at approximately 1 mg/mL.Rather than concluding that these data result from a shift to lower values of the measured particle size as concentration increases [6], it is alternatively here proposed that as concentration increased above 5 mg/mL, particle size increased.

Figure 2.
Flow FFF emergence time dependence on concentration for polystyrene sulfonate (MW = 5400 Da) at 20.2, 10.1, 5.0, 2.5, and 1.0 mg/mL and Aldrich HA at 20.0, 10.3, 5.2, 2.6, and 1.1 mg/mL in low electrolytic conductivity solution (0.05% surfactant/0.03%NaN3) [5].This behavior is postulated to result from increased aggregation, i.e., increased size, of the humic acid aggregate as concentration increased, resulting in a change in the separation mechanism from predominantly normal (Brownian) mode flow FFF to predominantly steric (hyperlayer or lift) mode flow FFF (Figure 3).The steric inversion point in flow FFF has been reported to occur at approximately 1 µm [8]; it is also influenced by shape.In normal (Brownian) mode flow FFF, submicrometer particles elute before larger particles, whereas in steric (hyperlayer, lift) mode, micrometer-sized particles elute before smaller particles [8].Because reversal of the elution order for Aldrich HA and PSS (Mw = 5400) was observed (Figures 2 and 3), the conclusion is that the predominant particle size became greater than 1 μm between a concentration of 2.6 and 5.2 mg•mL −1 .Peak narrowing also occurred as concentration increased indicating that the range of particle sizes in the sample reduced.
If normal flow FFF predominates during the low ionic strength fractionation of HSs at a concentration of 1 mg/mL, the hydrodynamic diameter of the HA is calculated Equation (1) to be 1.7 nm at the peak maximum (Figure 4), ranging from 0.5 nm to 3.3 nm for Suwannee River FA and from 0.7 nm to 4.4 nm for Suwannee River HA.These values correspond with a hydrodynamic diameter of 1.5 nm determined by Lead et al. [7] for Suwannee River FA at an ionic strength of 5 mM NaCl.The diffusion coefficients are reported [3] for 1 mg/mL solutions (normal mode) of 21 FAs/HSs-averaging 3.0 × 10 −10 m 2 •s −1 -that were measured in low ionic strength solution (0.05% surfactant/0.03%NaN3).Lead et al. [7] reported diffusion coefficients of 3.0 × 10 −10 m 2 •s −1 and 2.9 × 10 −10 m 2 •s −1 at pH 5.5 and 7.0, respectively, for Suwannee River FA measured by flow FFF at low ionic strength (5 mM•NaCl).
However, if steric flow FFF predominates in low ionic strength solution at a concentration of 20 mg•mL −1 , the hydrodynamic diameter of the HA is calculated Equation (2) to be 160 µm (Figure 4) at the peak maximum, which is approximately 30% of the width of the channel (i.e., 508 μm).At 10.32 mg•mL −1 the diameter is 125 μm, and 100 μm at 5.16 mg•mL −1 .

High Electrolytic Conductivity Solutions
In high ionic strength solution, the flow FFF of three types of humic substances exhibited concentration independence (Figure 1b) at 1-10 mg/L concentrations in 0.1 M potassium phosphate (pH 8.1) [4].The hydrodynamic diameter, calculated using Equation (1), averaged 2.2 nm at the peak maximum and ranged, on average, from 0.84 nm to 6.9 nm (Table 1).Comparably, Lead et al. [7] measured the hydrodynamic diameter of Suwannee River FA to be 2.1 nm at 50 mM, and Assemi et al. [13] reported hydrodynamic diameters of HSs ranging from 0.99 nm to 3.69 nm as measured by flow FFF in mobile phase ionic strength of 0.08 M.

Low Electrolytic Conductivity Solutions
Recently, using DLS, Esfahani et al. [2] demonstrated the abiotic reversible self-assembly of FA and HA aggregates in low electrolytic conductivity solutions, providing independent evidence for the existence of HSs aggregates.Particle size distributions in three ranges, i.e., 10-100 nm, 100-1000 nm, and >1 µm were observed in low electrolytic conductivity solutions with or without added sodium azide (0.03%).Sample concentrations ranged from 2 mg•L −1 to 330 mg•L −1 expressed as DOC.No particles were observed that were smaller than 10 nm.Although the submicrometer size distributions varied with time, they were always present in all samples analyzed, in contrast to a 5-µm-sized particle that formed within minutes, dissipated, and spontaneously re-formed over turbulent/quiescent cycles.The supramicrometer-sized aggregate self-assembled to a size ten-times greater than the nominal filter size (0.45 μm).The submicrometer aggregates are postulated to be precursors for the formation of the supramicrometer-sized aggregate.The supramicrometer particle was few-in-number because they were present in intensity-and volume-based DLS spectra but absent in number-based spectra.Because the volume-based particle size distributions are expressed as a percentage, DLS spectra collected at a specific time do not directly show the effects of changing concentration.However, volume-based spectra graphed over time indicate the increased stability of the supramicrometer-sized aggregate, and the diminished time at which formation of this aggregate was first observed, after onset of quiescence, as concentration increased.Larger or smaller aggregates may exist in these solutions, but were beyond the analytical range of the instrument, i.e., 0.6 nm-6 µm.

High Electrolytic Conductivity Solutions
The Verdugo research group [14][15][16] examined seawater (high electrolytic conductivity) amended with sodium azide by DLS and observed only freely dissolved HS polymers (~5-50 nm) during the first 2 h following filtration.Between 5 h and 10 h, nanogels (~100-200 nm) appeared that began to anneal to form larger nanogels.After 60 h, microgels (~3-6 µm) formed that were composed of nanogels annealed together.They also demonstrated reversible self-assembly of the microgel after filtration and verified their data using atomic force microscopy, flow cytometry, and environmental scanning electron microscopy.The same general types of particle size distributions, i.e., 10-100 nm, 100-1000 nm, and >1 µm, were observed using DLS by Esfahani et al. [2], although at low conductivity the microgels formed within minutes as opposed to within hours for the samples examined by the Verdugo research group at high conductivity [14][15][16].Significantly, in the DLS research, it should be recognized that the µm-sized-particle formed at the maximum operational range of the instrumentation-the existence of particles larger than 6 µm cannot be discounted.

Integrating the Data
The disaggregation/aggregation of HSs (Equation ( 4)) is here proposed to be a dynamic equilibrium process for which the rate of aggregate formation (at least for some types of aggregates) is controlled by the electrolytic conductivity of the solution: In low electrolytic conductivity solutions, k1 and k2 are larger (the disaggregation/aggregation of HSs occurs more rapidly-in minutes), whereas in high electrolytic conductivity solutions, k1 and k2 are smaller (the disaggregation/aggregation of HSs occurs more slowly-in hours).Both the forward and reverse processes appear to be kinetically favored in low electrolytic conductivity solutions.This explains why the fractograms of HSs are broad in low electrolytic conductivity solutions and are more narrow and Gaussian-shaped in high electrolytic conductivity solutions (Figure 1).In flow FFF conducted at high electrolytic conductivity, the time scale of the experiment is too short to allow for the formation of aggregates-therefore, the fractograms are concentration independent.
A substantial degree of variation in the hydrodynamic diameters of HSs has been reported in the data discussed here, overall representing four size ranges with one exception, i.e., (1) <10 nm; (2) 10 nm-6 μm; (3) 6-100 μm (no representatives in this data for this range); and (4) >100 μm.Wershaw and Hayes [17] noted that HAs form micelle-like aggregates above a critical concentration, but DOC molecules form pre-micellar aggregations at natural concentrations.NOM was considered to be a continuum of hydrodynamic diameters by Beckett and Ranville [18].Tarasevich et al. [12] reported that, at low concentrations, HAs exist as individual molecules; as the concentration is increased, they associate to form supramolecular structures at 5 mg•dm −3 (5 mg•L −1 ) and micelles at 8 g•dm −3 (8 g•L −1 ).Representative components of the different size ranges probably all dynamically coexist in solution.

Aggregates >100 μm
Particle sizes greater than 100 μm dominated the flow FFF fractograms at concentrations of approximately 5 g•L −1 or greater in low conductivity solutions.Although it can be argued that these concentrations are not environmentally relevant in natural systems, some situations, such as in soils or in wastewaters (and definitely in laboratories), might approach such concentrations.The existence of these aggregates lends support to the micelle/vesicle conceptualization proposed by Wershaw [19].
Manning et al. [20] used laser diffraction to measure the particle size distribution by volume of a 10 ppm HA solution and reported particle diameters ranging from 100 µm to greater than 2500 µm, pointing out that a single 1000 µm particle occupies the same volume of HA as 4.6 billion particles that are 0.6 µm in size.Aggregates of this size/type were noted to form in low electrolytic conductivity solution above a specific concentration (Figure 1a); however, the data presented here do not establish whether or not aggregates of this size/type form in high electrolytic conductivity solutions.
Contrary to some comments in the literature, the critical micelle concentration is not the concentration at which micelles first begin to form.Speaking generally about micelles (not specifically HS micelles), Mukerjee and Mysels [21] stated that, depending on the sensitivity of the analysis, micelles may be undetectable in dilute solutions of the constituent monomers, and at low concentrations, the additional solute molecules which form micelles may be very low until a concentration is reached at which nearly all additional solute molecules form micelles.That concentration is recognized as the critical micelle concentration.Therefore, micelles may exist in solutions less concentrated than the critical micelle concentration, but be undetectable.

Aggregates 6-100 µm
Significantly, in the DLS research reported here, the maximum operational range of the instrumentation was 6 µm.In the steric mode flow FFF experiments at low electrolytic conductivity, the smallest particle measured was 100 µm.Therefore, in this data set, no aggregates were detected in this size category; however, the existence of such particles cannot be discounted.
3.3.3.Aggregates 10-100 nm, 100-1000 nm, and 1-6 µm Are aggregates in these size ranges densely packed associations of HS molecules, or are they primarily water-based hydrogels?Also, if they exist as weakly associated hydrogel structures, how will shear forces in FFF affect these types of aggregates?
By DLS, the formation of aggregates in each of these three size ranges was observed to occur in minutes in HS solutions of low ionic strength [2], but required hours to days to form in HS solutions of seawater at high ionic strength [14][15][16].The presence of an excess amount of counterions (high conductivity) could be hypothesized to neutralize the negatively charged HSs and promote hydrophobic association, but that is counterintuitive to the concept of the slow formation of the aggregates at high ionic strength.However, if the aggregates in the 10-100 nm, 100-1000 nm, and 1-6 µm particle size distributions are truly hydrogels (nanogels and microgels [14][15][16]), then the hydrophilic charged form of HSs would be more water-friendly, and would explain why they form more readily at low ionic strength.Chin et al. [16] described the DOC hydrogels as an ionic-bond-stabilized tangled network, rather than a covalently crosslinked network.The self-assembled marine microgel in seawater was envisioned to be composed of approximately 1% solids, which represents a higher local concentration of organic matter in the particle than in the surrounding water environment [16].If the concentration of organic matter in the hydrogel is as low as 1%, it may not be possible to identify these aggregates in flow FFF with ultraviolet or fluorescence detection.
The 1-6 µm-sized-particle was observed to reversibly self-assemble at very low DOC concentrations (≤2 ppm) in simulated and environmental water samples [2,16].This supramicrometer-sized particle was determined to be few-in-number because they were present in intensity-and volume-based DLS spectra but absent in number-based spectra [2].Manning et al. [20] used two instruments with a combined operational size range from <100 nm to 2000 µm to evaluate the effect of chemical matrix and measured the size of aggregates of Aldrich humic acid solutions ranging in ionic strength.Their light scattering research of humic acid also led to the conclusion that most of the volume existed in a few particles of large diameter.At pH 10 and 10 ppm (Day 4), they established that most particle diameters were in the 0.5 to 0.6 µm range, over 90% were smaller than 0.8 µm, and 99.9% were smaller than 1 µm.
Turbulence, in the form of pouring, filtering, or shaking, was demonstrated by DLS [2] to have a disruptive effect on the >1 µm-sized particle (microgel), after which it reversibly self-assembled upon re-establishment of quiescence.Conversely, by DLS, the presence of particles in the 10 nm to 1 µm-size range (nanogels) was observed before, and immediately after, inducing turbulence.The Verdugo group [14][15][16], posits that the microgels were formed by annealing of the nanogels, which could explain the different character observed in response to turbulence for these two entities.Therefore, the turbulence generated by shear forces of the cross and channel flows in flow FFF can be postulated to be more disruptive to particles in the 1-6 µm range.The application of flow FFF to the analysis of hydrogels of humic and other materials is an exciting research front to be pursued.

Aggregates <10 nm
Aggregates smaller than 10 nm were not observed in either of the DLS studies reviewed here.Verdugo [14] illustrated the presence of extended free DOC polymers in angstrom-sized molecular dimensions by atomic force microscopy and environmental scanning electron microscopy.Conversely, aggregates in this size range were measured by flow FFF at both low and high ionic strength mobile phases.
The smallest unit size reported here is 0.5 nm or 5 Angstroms, which approaches the size expected for individual molecules.The largest size in this range determined by flow FFF, i.e., 7.8 nm, can be considered to consist of several humic subunits joined together in extended or linear random-coiled chains-and most nearly deserve to be called macromolecules of all the aggregates studied-different in character compared to the members of the larger aggregate size ranges.

Conclusions
FAs and HAs are fascinating chemicals.The more they are studied, the greater is our respect for their complexity.The question that began as-Why is the flow field-flow fractionation of fulvic and humic acids concentration independent at high ionic strength and concentration dependent at low ionic strength?-led to investigation of the fundamental differences in the character of HS aggregates as ionic strength varies.The data presented here are interpreted to indicate that the kinetics of the association-dissociation equilibrium for HSs favor the increased rate of aggregation in low ionic strength solutions in preference to high ionic strength solutions.In flow FFF conducted at high electrolytic conductivity, the time scale of generating the fractogram is too short to allow for the formation of aggregates >10 nm-therefore, the fractograms are concentration independent and Gaussian.In flow FFF conducted at low electrolytic conductivity, the disaggregation/aggregation processes are fast enough to be observed within the time scale of generating the fractogram-therefore, the fractograms are concentration dependent and non-Gaussian.
Eventually, equilibrium will be reached in both low and high ionic strength solutions; however, the equilibrium state itself does not reveal anything about how kinetically fast it will be reached.The kinetics of the disaggregation/aggregation equilibrium best supports conceptualization that the HS aggregates between 10 nm and 6 µm are hydrogels rather than densely packed associations of HS molecules.If the 1-6 µm particle size distribution is truly a hydrogel, and if the >100 µm particle is truly a micelle or vesicle (resulting from hydrophobic association) they will have different properties that can be examined by further research.To build upon this research, future examination of the fractionation of FAs and HAs as ionic strength varies should be pursued by coupling FFF with light scattering detectors and with more sophisticated inductively coupled plasma-mass spectrometry (ICP-MS) or ICP-atomic emission spectroscopy (AES) detectors to investigate the influence of different counterions on the formation of aggregates.

Table 1 .
Calculation of hydrodynamic diameters at the lower range of the peak (dmin), at the upper range of the peak (dmax), and at the peak maximum (dp) measured in high electrolytic conductivity solution by ultraviolet (UV) or fluorescence (FL) detection.SRFA: Suwannee River fulvic acid.

Table 2 .
Calculation of diffusion coefficients at the peak maximum (Dp) measured in high electrolytic conductivity solution by ultraviolet (UV) or fluorescence (FL) detection.