Novel Kraft Softwood Lignin-Derived Carbon Quantum Dots: Synthesis, Characterization, and In Vitro Cytocompatibility.

Carbon quantum dots (CQDs) have been investigated for biomedical applications in medical imaging due to their fluorescent properties, overall long-term stability, and excellent cytocompatibility and biocompatibility. Lignin is an organic polymer in the tissues of woody plants. It is also considered a byproduct of the wood and pulp industries. Hence, it presents as a renewable source of carbon nanoparticles. In this study, we report the synthesis and material and biological characterization of two colloidal suspensions of CQDs in water derived from lignin-based carbon. One was the native form of CQDs derived from lignin carbon, and the second was doped with nitrogen to evaluate material differences. Material characterization was carried out using various commonly used techniques, including Fourier transform infrared spectroscopy (FTIR), emission and absorbance spectra, zeta potential, transmission electron microscopy (TEM), and scanning electron microscopy (SEM). Thin films of CQDs were formed on glass and silicon substrates to assess the in vitro cytocompatibility with human mesenchymal stem cells (hMSCs). Observations suggest that the two forms of CQDs promote cell attachment within 24 h and sustain it for at least 7 days. The overall structure and shape of cells suggest a lack of any adverse or toxic effects of CQDs. The data lay down the novel foundation to support the use of lignin-derived CQDs in tissue engineering applications.

Nanoscale Stabilization Mechanism of Sodium Sulfate Decahydrate at Polyelectrolyte Interfaces.

Sodium sulfate decahydrate (SSD) is a low-cost phase-change material (PCM) for thermal energy storage applications that offers substantial melting enthalpy and a suitable temperature range for near-ambient applications. However, SSD's consistent phase separation with decreased melting enthalpy over repeated thermal cycles limits its application as a PCM. Sulfonated polyelectrolytes, such as dextran sulfate sodium (DSS), have shown great effectiveness in preventing phase separation in SSD. However, there is limited understanding of the stabilization mechanism of SSD by DSS at the atomic length and time scales. In this work, we investigate SSD stabilization via DSS using neutron scattering and molecular dynamics (MD) simulations. Neutron scattering and pair distribution function analysis revealed the structural evolution of the PCM samples below and above the phase change temperatures. MD simulations revealed that water from the hydrate structure migrates from the hydrate crystal to the SSD-DSS interfacial region upon melting. The water is stabilized at this interface by aggregation around the hydrophilic sulfonic acid groups attached to the backbone of the polyelectrolyte. This architecture retains water near the dehydrated sodium sulfate, preventing phase separation and, consequently, stabilizing SSD rehydration. This work provides atomistic insight into selecting and designing stable and high-performance PCMs for heating and cooling applications in building technologies.

Stabilization of low-cost phase change materials for thermal energy storage applications.

Sodium sulfate decahydrate (Na2SO4.10H2O, SSD), a low-cost phase change material (PCM), can store thermal energy. However, phase separation and unstable energy storage capacity (ESC) limit its use. To address these concerns, eight polymer additives-sodium polyacrylate (SPA), carboxymethyl cellulose (CMC), Fumed silica (SiO2), potassium polyacrylate (PPA), cellulose nanofiber (CNF), hydroxyethyl cellulose (HEC), dextran sulfate sodium (DSS), and poly(sodium 4-styrenesulfonate) (PSS)-were used to explore several stabilization mechanisms. The ESC of PCMs deteriorated when thickeners, SPA, PPA, and CNF, were added. DSS-modified PCMs exhibited greater stability up to 150 cycles. Rheology measurements indicated that DSS did not impact SSD viscosity significantly during stabilization. Dynamic light scattering showed that DSS reduces SSD particle size and electrostatically suspends salt particles in a stable homogeneous solution, avoiding phase separation. This study proposes a promising method to improve the thermal stability of salt hydrate PCMs by utilizing polyelectrolyte-salt hydrate mixture for thermal energy storage applications.

Exploration of Entropy Pair Functional Theory.

Evaluation of the entropy from molecular dynamics (MD) simulation remains an outstanding challenge. The standard approach requires thermodynamic integration across a series of simulations. Recent work Nicholson et al. demonstrated the ability to construct a functional that returns excess entropy, based on the pair correlation function (PCF); it was capable of providing, with acceptable accuracy, the absolute excess entropy of iron simulated with a pair potential in both fluid and crystalline states. In this work, the general applicability of the Entropy Pair Functional Theory (EPFT) approach is explored by applying it to three many-body interaction potentials. These potentials are state of the art for large scale models for the three materials in this study: Fe modelled with a modified embedded atom method (MEAM) potential, Cu modelled with an MEAM and Si modelled with a Tersoff potential. We demonstrate the robust nature of EPFT in determining excess entropy for diverse systems with many-body interactions. These are steps toward a universal Entropy Pair Functional, EPF, that can be applied with confidence to determine the entropy associated with sophisticated optimized potentials and first principles simulations of liquids, crystals, engineered structures, and defects.

Local Structure Analysis and Modelling of Lignin-Based Carbon Composites through the Hierarchical Decomposition of the Radial Distribution Function.

Carbonized lignin has been proposed as a sustainable and domestic source of activated, amorphous, graphitic, and nanostructured carbon for many industrial applications as the structure can be tuned through processing conditions. However, the inherent variability of lignin and its complex physicochemical structure resulting from feedstock and pulping selection make the Process-Structure-Property-Performance (PSPP) relationships hard to define. In this work, radial distribution functions (RDFs) from synchrotron X-ray and neutron scattering of lignin-based carbon composites (LBCCs) are investigated using the Hierarchical Decomposition of the Radial Distribution Function (HDRDF) modelling method to characterize the local atomic environment and develop quantitative PSPP relationships. PSPP relationships for LBCCs defined by this work include crystallite size dependence on lignin feedstock as well as increasing crystalline volume fraction, nanoscale composite density, and crystallite size with increasing reduction temperature.

Lithium and Sodium Ion Binding Mechanisms and Diffusion Rates in Lignin-Based Hard Carbon Models.

Hard carbons are the primary candidate for the anode of next-generation sodium-ion batteries for large-scale energy storage, as they are sustainable and can possess high charge capacity and long cycle life. These properties along with diffusion rates and ion storage mechanisms are highly dependent on nanostructures. This work uses reactive molecular dynamics simulations to examine lithium and sodium ion storage mechanisms and diffusion in lignin-based hard carbon model systems with varying nanostructures. It was found that sodium will preferentially localize on the surface of curved graphene fragments, while lithium will preferentially bind to the hydrogen dense interfaces of crystalline and amorphous carbon domains. The ion storage mechanisms are explained through ion charge and energy distributions in coordination with snapshots of the simulated systems. It was also revealed that hard carbons with small crystalline volume fractions and moderately sized sheets of curved graphene will yield the highest sodium-ion diffusion rates at ∼10-7 cm2/s. Self-diffusion coefficients were determined by mean square displacement of ions in the models with extension through a confined random walk theory.

Hierarchical Lignin-Based Carbon Matrix and Carbon Dot Composite Electrodes for High-Performance Supercapacitors.

This work adopts an efficient chemical-wet method to build a three-dimensional (3D) carbon composite as an electrode material for high-performance supercapacitors (SCs). Carbon dots (CDs), prepared by thermal pyrolysis of citric acid and urea under microwaves at 280 °C, are homogeneously coated onto lignin-based activated carbons (ACs), thus forming the 3D composites possessing an interior surface decorated with CD binding sites. Benefiting from the hydrophilicity and ultrafine size of CDs, the affinity of the electrode surface toward aqueous electrolytes is significantly improved with the addition of CDs, leading to the enhanced effective surface area (i.e., abundant electroactive sites) and a decreased ionic diffusion path. The capacitance of the SCs is improved from 125.8 to 301.7 F g-1 with CD addition. The SC with CD addition possesses improved cycle stability with a coulombic efficiency around 100% after 3000 cycles. After cycling, the ion diffusion coefficient of the CD@AC-11 electrode is enhanced by 25.5 times as compared to that of the pristine AC one. This unique and robust carbon framework can be utilized for engineering the desired pore structure and micropore/mesopore fraction within the AC electrodes. This strategy of CD@AC electrodes demonstrates a promising route for using renewable porous carbon materials in advanced energy-storage devices.

Entropy Pair Functional Theory: Direct Entropy Evaluation Spanning Phase Transitions.

We prove that, within the class of pair potential Hamiltonians, the excess entropy is a universal, temperature-independent functional of the density and pair correlation function. This result extends Henderson's theorem, which states that the free energy is a temperature dependent functional of the density and pair correlation. The stationarity and concavity of the excess entropy functional are discussed and related to the Gibbs-Bugoliubov inequality and to the free energy. We apply the Kirkwood approximation, which is commonly used for fluids, to both fluids and solids. Approximate excess entropy functionals are developed and compared to results from thermodynamic integration. The pair functional approach gives the absolute entropy and free energy based on simulation output at a single temperature without thermodynamic integration. We argue that a functional of the type, which is strictly applicable to pair potentials, is also suitable for first principles calculation of free energies from Born-Oppenheimer molecular dynamics performed at a single temperature. This advancement has the potential to reduce the evaluation the free energy to a simple modification to any procedure that evaluates the energy and the pair correlation function.

Local structure and distortions of mixed methane-carbon dioxide hydrates.

A vast source of methane is found in gas hydrate deposits, which form naturally dispersed throughout ocean sediments and arctic permafrost. Methane may be obtained from hydrates by exchange with hydrocarbon byproduct carbon dioxide. It is imperative for the development of safe methane extraction and carbon dioxide sequestration to understand how methane and carbon dioxide co-occupy the same hydrate structure. Pair distribution functions (PDFs) provide atomic-scale structural insight into intermolecular interactions in methane and carbon dioxide hydrates. We present experimental neutron PDFs of methane, carbon dioxide and mixed methane-carbon dioxide hydrates at 10 K analyzed with complementing classical molecular dynamics simulations and Reverse Monte Carlo fitting. Mixed hydrate, which forms during the exchange process, is more locally disordered than methane or carbon dioxide hydrates. The behavior of mixed gas species cannot be interpolated from properties of pure compounds, and PDF measurements provide important understanding of how the guest composition impacts overall order in the hydrate structure.

Li-Ion Localization and Energetics as a Function of Anode Structure.

In this work, we study the effect of carbon composite anode structure on the localization and energetics of Li-ions. A computational molecular dynamics study is combined with experimental results from neutron scattering experiments to understand the effect of composite density, crystallite size, volume fraction of crystalline carbon, and ion loading on the nature of ion storage in novel, lignin-derived composite materials. In a recent work, we demonstrated that these carbon composites display a fundamentally different mechanism for Li-ion storage than traditional graphitic anodes. The edges of the crystalline and amorphous fragments of aromatic carbon that exist in these composites are terminated by hydrogen atoms, which play a crucial role in adsorption. In this work, we demonstrate how differences in composite structure due to changes in the processing conditions alter the type and extent of the interface between the amorphous and crystalline domains, thus impacting the nature of Li-ion storage. The effects of structural properties are evaluated using a suite of pair distribution functions as well as an original technique to extract archetypal structures, in the form of three-dimensional atomic density distributions, from highly disordered systems. The energetics of Li-ion binding are understood by relating changes in the energy and charge distributions to changes in structural properties. The distribution of Li-ion energies reveals that some structures lead to greater chemisorption, while others have greater physisorption. Carbon composites with a high volume fraction of small crystallites demonstrate the highest ion storage capacity because of the high interfacial area between the crystalline and amorphous domains. At these interfaces, stable H atoms, terminating the graphitic crystallites, provide favorable sites for reversible Li adsorption.

Molecular Dynamics Simulations of Hydration Effects on Solvation, Diffusivity, and Permeability in Chitosan/Chitin Films.

The effects of hydration on the solvation, diffusivity, solubility, and permeability of oxygen molecules in sustainable, biodegradable chitosan/chitin food packaging films were studied via molecular dynamics and confined random walk simulations. With increasing hydration, the membrane has a more homogeneous water distribution with the polymer chains being fully solvated. The diffusivity increased by a factor of 4 for oxygen molecules and by an order of magnitude for water with increasing the humidity. To calculate the Henry's constant and solubility of oxygen in the membranes with changing hydration, the excess chemical potential was calculated via free energy perturbation, thermodynamic integration and direct particle deletion methods. The simulations predicted a higher solubility and permeability for the lower humidity, in contradiction to experimental results. All three methods for calculating the solubility were in good agreement. It was found that the Coulombic interactions in the potential caused the oxygen to bind too strongly to the protonated amine group. Insight from this work will help guide molecular modeling of chitosan/chitin membranes, specifically permeability measurements for small solute molecules. Efforts to chemically tailor chitosan/chitin membranes to favor discrete as opposed to continuous aqueous domains could reduce oxygen permeability.

Ab Initio Molecular Dynamics Simulations of an Excess Proton in a Triethylene Glycol-Water Solution: Solvation Structure, Mechanism, and Kinetics.

We investigate the solvation shell structures, the distribution of protonic defects, mechanistic details, kinetics, and dynamics of proton transfer for an excess proton in bulk water and for an excess proton in an aqueous solution of triethylene glycol (TEG) via Car-Parrinello molecular dynamics simulations. The PW91, PBE, and PBE with the Tkatchenko-Scheffler (TS) density-dependent dispersion functionals were used and compared for bulk water and the TEG-water mixtures. The excess proton is found to reside predominantly on water molecules but also resides on hydroxyl groups of TEG. The lifetimes associated with structural diffusion time scales of the protonated water were found to be on the order of ∼1 ps. All three functionals studied support the presolvation requirement for structural diffusion. The highest level of theory shows a reduction in the free energy barrier for water-water proton transfer in TEG-water mixtures compared to bulk water. The effect of TEG shows no strong change in the kinetics for TEG-water mixtures compared to bulk water for this same level of theory. The excess proton displays burst-rest behavior in the presence of TEG, similar to that found in bulk water. We find that the TEG chain disrupts the hydrogen-bond network, causing the solvation shell around water to be populated by TEG chain groups instead of other waters, reducing the rigidity of the hydrogen-bond network. Methylene is a dominant hydrogen bond donor for the protonated water in hydrogen-bond networks associated with proton transfer and structural diffusion. This is consistent with previous studies that have found the hydronium ion to be amphiphilic in nature and to have higher proton mobility at oil-water interfaces.

Electron beam induced radiation damage in the catalyst layer of a proton exchange membrane fuel cell.

Electron microscopy is an essential tool for the evaluation of microstructure and properties of the catalyst layer (CL) of proton exchange membrane fuel cells (PEMFCs). However, electron microscopy has one unavoidable drawback, which is radiation damage. Samples suffer temporary or permanent change of the surface or bulk structure under radiation damage, which can cause ambiguity in the characterization of the sample. To better understand the mechanism of radiation damage of CL samples and to be able to separate the morphological features intrinsic to the material from the consequences of electron radiation damage, a series of experiments based on high-angle annular dark-field-scanning transmission scanning microscope (HAADF-STEM), energy filtering transmission scanning microscope (EFTEM), and electron energy loss spectrum (EELS) are conducted. It is observed that for thin samples (0.3-1 times λ), increasing the incident beam energy can mitigate the radiation damage. Platinum nanoparticles in the CL sample facilitate the radiation damage. The radiation damage of the catalyst sample starts from the interface of Pt/C or defective thin edge and primarily occurs in the form of mass loss accompanied by atomic displacement and edge curl. These results provide important insights on the mechanism of CL radiation damage. Possible strategies of mitigating the radiation damage are provided.

Rotating phenyl rings as a guest-dependent switch in two-dimensional metal-organic frameworks.

A semirigid bis(1,2,4-triazole) ligand binds in a syn conformation between copper(I) chains to form a series of two-dimensional metal-organic frameworks that display a topology of fused one-dimensional metal-organic nanotubes. These anisotropic frameworks undergo two different transformations in the solid state as a function of solvation. The 2D sheet layers can expand or contract, or, more remarkably, the phenyl rings can rotate between two distinct positions. Rotation of the phenyl rings allows for the adjustment of the tube size, depending on the guest molecules present. This "gate" effect along the 1D tubes has been characterized through single-crystal X-ray diffraction. The transformations can also be followed by powder X-ray diffraction (PXRD) and solid-state (13)C cross-polarization magic-angle-spinning (CP-MAS) NMR. Whereas PXRD cannot differentiate between transformations, solid-state (13)C CP-MAS NMR can be employed to directly monitor phenyl rotation as a function of solvation, suggesting that this spectroscopic method is a powerful approach for monitoring breathing in this novel class of frameworks. Finally, simulations show that rotation of the phenyl ring from a parallel orientation to a perpendicular orientation occurs at the cost of framework-framework energy and that this energetic cost is offset by stronger framework-solvent interactions.

Kinetics of methane hydrate decomposition studied via in situ low temperature X-ray powder diffraction.

Gas hydrate is known to have a slowed decomposition rate at ambient pressure and temperatures below the melting point of ice. As hydrate exothermically decomposes, gas is released and water of the clathrate cages transforms into ice. Based on results from the decomposition of three nominally similar methane hydrate samples, the kinetics of two regions, 180-200 and 230-260 K, within the overall decomposition range 140-260 K, were studied by in situ low temperature X-ray powder diffraction. The kinetic rate constants, k(a), and the reaction mechanisms, n, for ice formation from methane hydrate were determined by the Avrami model within each region, and activation energies, E(a), were determined by the Arrhenius plot. E(a) determined from the data for 180-200 K was 42 kJ/mol and for 230-260 K was 22 kJ/mol. The higher E(a) in the colder temperature range was attributed to a difference in the microstructure of ice between the two regions.

A coarse-grained model for polyethylene glycol polymer.

A coarse-grained (CG) model of polyethylene glycol (PEG) was developed and implemented in CG molecular dynamics (MD) simulations of PEG chains with degree of polymerization (DP) 20 and 40. In the model, two repeat units of PEG are grouped as one CG bead. Atomistic MD simulation of PEG chains with DP = 20 was first conducted to obtain the bonded structural probability distribution functions (PDFs) and nonbonded pair correlation function (PCF) of the CG beads. The bonded CG potentials are obtained by simple inversion of the corresponding PDFs. The CG nonbonded potential is parameterized to the PCF using both an inversion procedure based on the Ornstein-Zernike equation with the Percus-Yevick approximation (OZPY(-1)) and a combination of OZPY(-1) with the iterative Boltzmann inversion (IBI) method (OZPY(-1)+IBI). As a simple one step method, the OZPY(-1) method possesses an advantage in computational efficiency. Using the potential from OZPY(-1) as an initial guess, the IBI method shows fast convergence. The coarse-grained molecular dynamics (CGMD) simulations of PEG chains with DP = 20 using potentials from both methods satisfactorily reproduce the structural properties from atomistic MD simulation of the same systems. The OZPY(-1)+IBI method yields better agreement than the OZPY(-1) method alone. The new CG model and CG potentials from OZPY(-1)+IBI method was further tested through CGMD simulation of PEG with DP = 40 system. No significant changes are observed in the comparison of PCFs from CGMD simulations of PEG with DP = 20 and 40 systems indicating that the potential is independent of chain length.

Effective potentials between nanoparticles in suspension.

Results of molecular dynamics simulations are presented for the pair distribution function between nanoparticles in an explicit solvent as a function of nanoparticle diameter and interaction strength between the nanoparticle and solvent. The effect of including the solvent explicitly is demonstrated by comparing the pair distribution function of nanoparticles to that in an implicit solvent. The nanoparticles are modeled as a uniform distribution of Lennard-Jones particles, while the solvent is represented by standard Lennard-Jones particles. The diameter of the nanoparticle is varied from 10 to 25 times that of the solvent for a range of nanoparticle volume fractions. As the strength of the interactions between nanoparticles and the solvent increases, the solvent layer surrounding the nanoparticle is formed which increases the effective radii of the nanoparticles. The pair distribution functions are inverted using the Ornstein-Zernike integral equation to determine an effective pair potential between the nanoparticles mediated by the introduction of an explicit solvent.

Applications of a general random-walk theory for confined diffusion.

A general random walk theory for diffusion in the presence of nanoscale confinement is developed and applied. The random-walk theory contains two parameters describing confinement: a cage size and a cage-to-cage hopping probability. The theory captures the correct nonlinear dependence of the mean square displacement (MSD) on observation time for intermediate times. Because of its simplicity, the theory also requires modest computational requirements and is thus able to simulate systems with very low diffusivities for sufficiently long time to reach the infinite-time-limit regime where the Einstein relation can be used to extract the self-diffusivity. The theory is applied to three practical cases in which the degree of order in confinement varies. The three systems include diffusion of (i) polyatomic molecules in metal organic frameworks, (ii) water in proton exchange membranes, and (iii) liquid and glassy iron. For all three cases, the comparison between theory and the results of molecular dynamics (MD) simulations indicates that the theory can describe the observed diffusion behavior with a small fraction of the computational expense. The confined-random-walk theory fit to the MSDs of very short MD simulations is capable of accurately reproducing the MSDs of much longer MD simulations. Furthermore, the values of the parameter for cage size correspond to the physical dimensions of the systems and the cage-to-cage hopping probability corresponds to the activation barrier for diffusion, indicating that the two parameters in the theory are not simply fitted values but correspond to real properties of the physical system.

Toward a predictive understanding of water and charge transport in proton exchange membranes.

An analytical model for water and charge transport in highly acidic and highly confined systems such as proton exchange membranes of fuel cells is developed and compared to available experimental data. The model is based on observations from both experiment and multiscale simulation. The model accounts for three factors in the system including acidity, confinement, and connectivity. This model has its basis in the molecular-level mechanisms of water transport but has been coarse-grained to the extent that it can be expressed in an analytical form. The model uses the concentration of H(3)O(+) ion to characterize acidity, interfacial surface area per water molecule to characterize confinement, and percolation theory to describe connectivity. Several important results are presented. First, an integrated multiscale simulation approach including both molecular dynamics simulation and confined random walk theory is capable of quantitatively reproducing experimentally measured self-diffusivities of water in the perfluorinated sulfonic acid proton exchange membrane material, Nafion. The simulations, across a range of hydration conditions from minimally hydrated to fully saturated, have an average error for the self-diffusivity of water of 16% relative to experiment. Second, accounting for three factors-acidity, confinement, and connectivity-is necessary and sufficient to understand the self-diffusivity of water in proton exchange membranes. Third, an analytical model based on percolation theory is capable of quantitatively reproducing experimentally measured self-diffusivities of both water and charge in Nafion across a full range of hydration.

Use of the Ornstein-Zernike Percus-Yevick equation to extract interaction potentials from pair correlation functions.

In this work, we test the ability of the Ornstein-Zernike equation in the Percus-Yevick approximation (OZPY) to generate interaction potentials from pair correlation functions (PCFs) of monatomic and diatomic Lennard-Jones fluids. The PCFs are generated by solving OZPY equation (monatomic fluid) and molecular-dynamics (MD) simulations (diatomic fluid). Since the interaction potentials are inputs in the OZPY method and the MD simulation, the extraction of the potential from the PCFs using OZPY is a test of self-consistency. This test is necessary if the procedure is to be used to generate coarse-grained (CG) potentials from PCFs. We find that the procedure is completely self-consistent for the monatomic fluid in the whole range of densities studied (reduced density up to 0.55, under reduced temperature of 2.0). In the diatomic case, we find that the procedure is generally self-consistent under both low and high densities, although there is a systematic deviation at high densities. The method is able to reproduce the two parameters (ε and σ) of the input Lennard-Jones potential model to within about 1%. This CG potential generating procedure can be straightforwardly extended to more complicated molecules.