Effects of Natural Organic Matter on the Dissolution Kinetics and Bioavailability of Metal Oxide Nanoparticles

The rapid development of nanotechnology and wider applications of engineered nanomaterials (ENMs) in the last few decades have generated concerns regarding their environmental and health risks. After release into the environment, ENMs undergo aggregation, transformation, and, for metal-based nanomaterials, dissolution processes, which together determine their fate, bioavailability and toxicity to living organisms in the ecosystems. The rates of these processes are dependent on nanomaterial characteristics as well as complex environmental factors, including natural organic matter (NOM). As a ubiquitous component of aquatic systems, NOM plays a key role in the aggregation, dissolution and transformation of metal-based nanomaterials and colloids in aquatic environments.

The goal of this dissertation work is to investigate how NOM fractions with different chemical and molecular properties affect the dissolution kinetics of metal oxide ENMs, such as zinc oxide (ZnO) and copper oxide (CuO) nanoparticles (NPs), and consequently their bioavailability to aquatic vertebrate, with Gulf killifish (Fundulus grandis) embryos as model organisms.

ZnO NPs are known to dissolve at relatively fast rates, and the rate of dissolution is influenced by water chemistry, including the presence of Zn-chelating ligands. A challenge, however, remains in quantifying the dissolution of ZnO NPs, particularly for time scales that are short enough to determine rates. This dissertation assessed the application of anodic stripping voltammetry (ASV) with a hanging mercury drop electrode to directly measure the concentration of dissolved Zn in ZnO NP suspensions, without separation of the ZnO NPs from the aqueous phase. Dissolved zinc concentration measured by ASV ([Zn]ASV) was compared with that measured by inductively coupled plasma mass spectrometry (ICP-MS) after ultracentrifugation ([Zn]ICP-MS), for four types of ZnO NPs with different coatings and primary particle diameters. For small ZnO NPs (4-5 nm), [Zn]ASV was 20% higher than [Zn]ICP-MS, suggesting that these small NPs contributed to the voltammetric measurement. For larger ZnO NPs (approximately 20 nm), [Zn]ASV was (79±19)% of [Zn]ICP-MS, despite the high concentrations of ZnO NPs in suspension, suggesting that ASV can be used to accurately measure the dissolution kinetics of ZnO NPs of this primary particle size.

Using the ASV technique to directly measure dissolved zinc concentration, we examined the effects of 16 different NOM isolates on the dissolution kinetics of ZnO NPs in buffered potassium chloride solution. The observed dissolution rate constants (kobs) and dissolved zinc concentrations at equilibrium increased linearly with NOM concentration (from 0 to 40 mg-C L-1) for Suwannee River humic acid (SRHA), Suwannee River fulvic acid and Pony Lake fulvic acid. When dissolution rates were compared for the 16 NOM isolates, kobs was positively correlated with certain properties of NOM, including specific ultraviolet absorbance (SUVA), aromatic and carbonyl carbon contents, and molecular weight. Dissolution rate constants were negatively correlated to hydrogen/carbon ratio and aliphatic carbon content. The observed correlations indicate that aromatic carbon content is a key factor in determining the rate of NOM-promoted dissolution of ZnO NPs. NOM isolates with higher SUVA were also more effective at enhancing the colloidal stability of the NPs; however, the NOM-promoted dissolution was likely due to enhanced interactions between surface metal ions and NOM rather than smaller aggregate size.

Based on the above results, we designed experiments to quantitatively link the dissolution kinetics and bioavailability of CuO NPs to Gulf killifish embryos under the influence of NOM. The CuO NPs dissolved to varying degrees and at different rates in diluted 5‰ artificial seawater buffered to different pH (6.3-7.5), with or without selected NOM isolates at various concentrations (0.1-10 mg-C L-1). NOM isolates with higher SUVA and aromatic carbon content (such as SRHA) were more effective at promoting the dissolution of CuO NPs, as with ZnO NPs, especially at higher NOM concentrations. On the other hand, the presence of NOM decreased the bioavailability of dissolved Cu ions, with the uptake rate constant negatively correlated to dissolved organic carbon concentration ([DOC]) multiplied by SUVA, a combined parameter indicative of aromatic carbon concentration in the media. When the embryos were exposed to CuO NP suspension, changes in their Cu content were due to the uptake of both dissolved Cu ions and nanoparticulate CuO. The uptake rate constant of nanoparticulate CuO was also negatively correlated to [DOC]×SUVA, in a fashion roughly proportional to changes in dissolved Cu uptake rate constant. Thus, the ratio of uptake rate constants from dissolved Cu and nanoparticulate CuO (ranging from 12 to 22, on average 17±4) were insensitive to NOM type or concentration. Instead, the relative contributions of these two Cu forms were largely determined by the percentage of CuO NP that was dissolved.

Overall, this dissertation elucidated the important role that dissolved NOM plays in affecting the environmental fate and bioavailability of soluble metal-based nanomaterials. This dissertation work identified aromatic carbon content and its indicator SUVA as key NOM properties that influence the dissolution, aggregation and biouptake kinetics of metal oxide NPs and highlighted dissolution rate as a useful functional assay for assessing the relative contributions of dissolved and nanoparticulate forms to metal bioavailability. Findings of this dissertation work will be helpful for predicting the environmental risks of engineered nanomaterials.

Physics of Hexagonal Limit-Periodic Phases: Thermodynamics, Formation and Vibrational Modes

Limit-periodic (LP) structures exhibit a type of nonperiodic order yet to be found in a natural material. A recent result in tiling theory, however, has shown that LP order can spontaneously emerge in a two-dimensional (2D) lattice model with nearest-and next-nearest-neighbor interactions. In this dissertation, we explore the question of what types of interactions can lead to a LP state and address the issue of whether the formation of a LP structure in experiments is possible. We study emergence of LP order in three-dimensional (3D) tiling models and bring the subject into the physical realm by investigating systems with realistic Hamiltonians and low energy LP states. Finally, we present studies of the vibrational modes of a simple LP ball and spring model whose results indicate that LP materials would exhibit novel physical properties.

A 2D lattice model defined on a triangular lattice with nearest- and next-nearest-neighbor interactions based on the Taylor-Socolar (TS) monotile is known to have a LP ground state. The system reaches that state during a slow quench through an infinite sequence of phase transitions. Surprisingly, even when the strength of the next-nearest-neighbor interactions is zero, in which case there is a large degenerate class of both crystalline and LP ground states, a slow quench yields the LP state. The first study in this dissertation introduces 3D models closely related to the 2D models that exhibit LP phases. The particular 3D models were designed such that next-nearest-neighbor interactions of the TS type are implemented using only nearest-neighbor interactions. For one of the 3D models, we show that the phase transitions are first order, with equilibrium structures that can be more complex than in the 2D case.

In the second study, we investigate systems with physical Hamiltonians based on one of the 2D tiling models with the goal of stimulating attempts to create a LP structure in experiments. We explore physically realizable particle designs while being mindful of particular features that may make the assembly of a LP structure in an experimental system difficult. Through Monte Carlo (MC) simulations, we have found that one particle design in particular is a promising template for a physical particle; a 2D system of identical disks with embedded dipoles is observed to undergo the series of phase transitions which leads to the LP state.

LP structures are well ordered but nonperiodic, and hence have nontrivial vibrational modes. In the third section of this dissertation, we study a ball and spring model with a LP pattern of spring stiffnesses and identify a set of extended modes with arbitrarily low participation ratios, a situation that appears to be unique to LP systems. The balls that oscillate with large amplitude in these modes live on periodic nets with arbitrarily large lattice constants. By studying periodic approximants to the LP structure, we present numerical evidence for the existence of such modes, and we give a heuristic explanation of their structure.

Out-of-equilibrium dynamical fluctuations in glassy systems

"In this paper we extend the earlier treatment of out-of-equilibrium mesoscopic fluctuations in glassy systems in several significant ways. First, via extensive simulations, we demonstrate that models of glassy behavior without quenched disorder display scalings of the probability of local two-time correlators that are qualitatively similar to that of models with short-ranged quenched interactions. The key ingredient for such scaling properties is shown to be the development of a criticallike dynamical correlation length, and not other microscopic details. This robust data collapse may be described in terms of a time-evolving "extreme value" distribution. We develop a theory to describe both the form and evolution of these distributions based on a effective sigma model approach."

Phase behavior and far-from-equilibrium gelation in charged attractive colloids

In this Rapid Communication we demonstrate the applicability of an augmented Gibbs ensemble Monte Carlo approach for the phase behavior determination of model colloidal systems with short-ranged depletion attraction and long-ranged repulsion. This technique allows for a quantitative determination of the phase boundaries and ground states in such systems. We demonstrate that gelation may occur in systems of this type as the result of arrested microphase separation, even when the equilibrium state of the system is characterized by compact microphase structures.

Phase coexistence of cluster crystals: beyond the Gibbs phase rule

We report a study of the phase behavior of multiple-occupancy crystals through simulation. We argue that in order to reproduce the equilibrium behavior of such crystals it is essential to treat the number of lattice sites as a constraining thermodynamic variable. The resulting free-energy calculations thus differ considerably from schemes used for single-occupancy lattices. Using our approach, we obtain the phase diagram and the bulk modulus for a generalized exponential model that forms cluster crystals at high densities. We compare the simulation results with existing theoretical predictions. We also identify two types of density fluctuations that can lead to two sound modes and evaluate the corresponding elastic constants.