Understanding the colloidal dispersion stability of 1D and 2D materials: Perspectives from molecular simulations and theoretical modeling

Insights into the Solution Processing of Non-van der Waals Boron Carbide

Exfoliation of non-layered materials is crucial to unleash their enormous potential in wide range of applications. However, the presence of strong non-van der Waals interactions in all three dimensions makes exfoliation challenging. Boron carbide (B4C), known for its high hardness, holds great potential for diverse applications. In this work, B4C is exfoliated in 29 solvents to determine its Hansen solubility parameters (HSP) and understand the dispersion stability. The experimentally determined HSP values are 19.6, 15.3 and 15.1 MPa1/2 for dispersive interactions (δD), polar interactions (δP) and hydrogen bonding interactions (δH), respectively. Subsequently, in situ dispersion stability is analyzed qualitatively and quantitatively from the space-time resolved extinction profiles employing an analytical centrifuge. The sedimentation kinetics of B4C dispersions are analyzed using instability index, integral extinction and sedimentation velocity. B4C dispersions in solvents like ϵ-caprolactone, isopropyl alcohol, propylene carbonate and formamide exhibited excellent stability. Additionally, B4C nanosheets are utilized to investigate their field emission properties. These nanosheets exhibit a turn-on electric field of 3.35 V/μm at an emission current density of 1 μA cm-2, and it delivered a large emission current density of ~375 μA cm-2 at an applied electric field of 5.7 V/μm.

Salinity-Driven Structural and Viscosity Modulation of Confined Polar Oil Phases by Carbonated Brine Films: Novel Insights from Molecular Dynamics

The structural and dynamic properties of fluids under confinement in a porous medium differ from their bulk properties. This study delves into the surface structuring and hydrodynamic characteristics of oil/thin film carbonated brine two-phase within a calcite channel upon salinity variation. To this end, both equilibrium and non-equilibrium molecular dynamics simulations are utilized to unveil the effect of the carboxylic acid component (benzoic acid) in a simple model oil (decane) confined between two thin films of carbonated brine on the oil-brine-calcite characteristics. The salinity effect was scrutinized under four saline carbonated waters, deionized carbonated water (DCW), carbonated low-salinity brine (CLSB, 30,000 ppm), carbonated seawater (CSW, 60,000 ppm), and carbonated high-salinity brine (CHSB, 180,000 ppm). An electrical double layer (EDL) is observed at varying salinities, comprising a Stern-like positive layer (formed by Na+ ions) followed by a negative one (formed by Cl- ions primarily residing on top of the adsorbed sodium cations). By lowering the salinity, the Na+ ions cover the interface regions (brine-calcite and brine-oil), depleting within the brine bulk region. The lowest positive surface charge on the rock surface was found in salinity corresponding to seawater. Two distinct Na+ peaks at the oleic phase interface have been observed in the carbonated high-salinity brine system, enhancing the adsorption of polar molecules at the thin brine film interfaces. There is a pronounced EDL formation at the oleic phase interface in the case of CSW, resulting in a strong interface region containing ions and functional fractions. Likewise, the oil region confined by CSW exhibited the lowest apparent viscosity, attributed to the optimized salinity distribution and inclination of benzoic acid fractions uniformly at the brine-oil interface, acting as a slippery surface. Moreover, the results reveal that the presence of polar fractions could increase the oil phase's apparent viscosity, and introducing ions to this system reduces the polar molecules' destructive effect on the apparent viscosity of the oil region. Therefore, the fluidity of confined systems is modulated by both composition of the brine and oil phases.

Surfactant-Aided Stabilization of Individual Carbon Nanotubes in Water around the Critical Micelle Concentration

Surfactants are widely used to disperse single-walled carbon nanotubes (SWCNTs) and other nanomaterials for liquid-phase processing and characterization. Traditional techniques, however, demand high surfactant concentrations, often in the range of 1-2 wt/v% of the solution. Here, we show that optimal dispersion efficiency can be attained at substantially lower surfactant concentrations of approximately 0.08 wt/v%, near the critical micelle concentration. This unexpected observation is achieved by introducing "bare" nanotubes into water containing the anionic surfactant sodium deoxycholate (DOC) through a superacid-surfactant exchange process that eliminates the need for ultrasonication. Among the diverse ionic surfactants and charged biopolymers explored, DOC exhibits the highest dispersion efficiency, outperforming sodium cholate, a structurally similar bile salt surfactant containing just one additional oxygen atom compared to DOC. Employing all-atomistic molecular dynamics simulations, we unravel that the greater stabilization by DOC arises from its higher binding affinity to nanotubes and a substantially larger free energy barrier that resists nanotube rebundling. Further, we find that this barrier is nonelectrostatic in nature and does not obey the classical Derjaguin-Landau-Verwey-Overbeek (DLVO) theory of colloidal stability, underscoring the important role of nonelectrostatic dispersion and hydration interactions at the nanoscale, even in the case of ionic surfactants like DOC. These molecular insights advance our understanding of surfactant chemistry at the bare nanotube limit and suggest low-energy, surfactant-efficient solution processing of SWCNTs and potentially other nanomaterials.

Incorporating ion-specific van der Waals and soft repulsive interactions in the Poisson-Boltzmann theory of electrical double layers

Electrical double layers (EDLs) arise when an electrolyte is in contact with a charged surface, and are encountered in several application areas including batteries, supercapacitors, electrocatalytic reactors, and colloids. Over the last century, the development of Poisson-Boltzmann (PB) models and their modified versions have provided significant physical insight into the structure and dynamics of the EDL. Incorporation of physics such as finite-ion-size effects, dielectric decrement, and ion-ion correlations has made such models increasingly accurate when compared to more computationally expensive approaches such as molecular simulations and classical density functional theory. However, a prominent knowledge gap has been the exclusion of van der Waals (vdW) and soft repulsive interactions in modified PB models. Although short-ranged as compared to electrostatic interactions, we show here that vdW and soft repulsive interactions can play an important role in determining the structure of the EDL via the formation of a Stern layer and in modulating the differential capacitance of an electrode in an electrolyte. To this end, we incorporate ion-ion and wall-ion vdW attraction and soft repulsion via a 12-6 Lennard-Jones (LJ) potential, resulting in a modified PB-LJ approach. The wall-ion LJ interactions were found to have a significant effect on the electrical potential and concentration profiles, especially close to the wall. However, ion-ion LJ interactions do not affect the EDL structure at low bulk ion concentrations (<1 M). We also derive dimensionless numbers to quantify the impact of ion-ion and wall-ion LJ interactions on the EDL. Furthermore, in the pursuit of capturing ion-specific effects, we apply our model by considering various ions such as Na, K+, Mg2+, Cl-, and SO42-. We observe how varying parameters such as the electrolyte concentration and electrode potential affect the structure of the EDL due to the competition between ion-specific LJ and electrostatic interactions. Lastly, we show that the inclusion of vdW and soft repulsion interactions, as well as hydration effects, leads to a better qualitative agreement of the PB models with experimental double-layer differential capacitance data. Overall, the modified PB-LJ approach presented herein will lead to more accurate theoretical descriptions of EDLs in various application areas.

The elastic response of graphene oxide gels as a crumpling phenomenon

The broad spectrum of chemical and electronic properties of 2D nanomaterials makes them attractive in a wide range of applications, especially in the context of printed electronics. Therefore, understanding the rheological properties of nanosheet suspensions is crucial for many additive manufacturing techniques. Here, we study the viscoelastic properties of aqueous suspensions of graphene oxide nanosheets. We show that in the gel phase, the magnitude of the elastic response and its scaling with volume fraction is independent of the lateral size of the particles and the interaction strength between them. We explain this behavior by modelling the elasticity of these gels as a crumpling phenomenon where the magnitude of the response is determined by the bending stiffness and thickness of the sheets. Due to their low bending stiffness these nanosheets crumple upon deformation and may therefore be considered soft colloids. Furthermore, we provide an explanation why the yield strain decreases with packing fraction for these gels.

Impact of sulfhydryl ligands on the transformation of silver ions by molybdenum disulfide and their combined toxicity to freshwater algae

The transformation of silver ions (Ag⁺) mediated by engineered nanomaterials (ENMs) influences the biosafety of Ag-containing products in natural environments. Actually, modification of biomolecules to ENMs in aquatic ecosystems alters their interactions with Ag⁺. This study discovered that surface functionalization of glutathione (GSH, a sulfhydryl compound ubiquitous in natural waters) on molybdenum disulfide (MoS₂) nanoflakes suppressed the redox reaction between 1 T components and Ag⁺, inhibiting the MoS₂-mediated reduction of Ag⁺ to Ag nanoparticles (AgNPs) in aqueous phase in the dark. However, AgNPs formation (from 2.32 ± 0.35-3.25 ± 0.29 mg/L per day, pH 7.0) and oxidation of MoS₂ were remarkably accelerated after GSH binding under light conditions. The dominant electron donator of MoS₂ to Ag⁺ was transformed from the electron-hole pairs to surface ligands driven by the introduction of chromophoric groups was authenticated as the cause for the elevated Ag⁺ reduction. These processes also occurred between Ag⁺ and MoS₂ at low levels (50 μg/L). Additionally, the joint algal toxicity of GSH-modified MoS₂ with Ag⁺ was weaker than that of pristine MoS₂ due to increased retention of free Ag⁺ and AgNPs formation. Our findings improve the understanding of the interaction between ENMs and Ag⁺ in aquatic ecosystems.

Gas Separations using Nanoporous Atomically Thin Membranes: Recent Theoretical, Simulation, and Experimental Advances

Porous graphene and other atomically thin 2D materials are regarded as highly promising membrane materials for high-performance gas separations due to their atomic thickness, large-scale synthesizability, excellent mechanical strength, and chemical stability. When these atomically thin materials contain a high areal density of gas-sieving nanoscale pores, they can exhibit both high gas permeances and high selectivities, which is beneficial for reducing the cost of gas-separation processes. Here, recent modeling and experimental advances in nanoporous atomically thin membranes for gas separations is discussed. The major challenges involved, including controlling pore size distributions, scaling up the membrane area, and matching theory with experimental results, are also highlighted. Finally, important future directions are proposed for real gas-separation applications of nanoporous atomically thin membranes.

Impact of algal extracellular polymeric substances on the environmental fate and risk of molybdenum disulfide in aqueous media

Molybdenum disulfide (MoS₂) poses great potential in water treatment as a popular transition metal dichalcogenide, arousing considerable concern regarding its fates and risk in aquatic environments. This study revealed that the interplay with extracellular polymeric substances (EPS) of freshwater algae significantly changed the properties and toxicity of MoS₂ to aquatic fish. The predominant binding of aromatic compounds, polysaccharides, and carboxyl-rich proteins in EPS on the 1T polymorph of MoS₂ via hydrophilic effects and the preferential adsorption of carboxylic groups contributed to morphological alterations, structural disorders (band gap and phase alterations), and the attenuated aggregation of MoS₂ in aqueous solutions. Electron charge transfer and n-π* interactions with EPS decreased the catalytic activity of MoS₂ by inhibiting its capability of generating reactive intermediates. The dissolution of MoS₂ slowed down after interacting with EPS (from 0.089 to 0.045 mg/L per day) owing to rapid initial oxidation (i.e., forming Mo-O bond) and carbon grafting. Notably, the morphological and structural alterations after EPS binding alleviated the toxicity (e.g., malformation and oxidative stress) of MoS₂ to infantile zebrafish. Our findings provide insights into the environmental fate and risk of MoS₂ by ubiquitous EPS in natural waters, serving as valuable information while developing water treatment processes accordingly.

Sulfur vacancies affect the environmental fate, corona formation, and microalgae toxicity of molybdenum disulfide nanoflakes

Sulfur vacancy (SV) defects have been engineered in two-dimensional (2D) transition metal dichalcogenides (TMDs) for high performance applications in various fields involving environmental protection. Understanding the influence of SVs on the environmental fate and toxicity of TMDs is critical for evaluating their risk. Our work discovered that SVs (with S/Mo ratios of 1.65 and 1.32) reduced the dispersibility and promoted aggregation of 2H phase molybdenum disulfide (2H-MoS₂, a hot TMD) in aqueous solution. The generation capability of •O₂^- and •OH was increased and the dissolution of 2H-MoS₂ was significantly accelerated after SVs formation. Different with pristine form, S-vacant 2H-MoS₂ preferentially harvested proteins (i.e., forming protein corona) involved in antioxidation, photosynthetic electron transport, and the cytoskeleton structure of microalgae. These proteins contain a higher relative number of thiol groups, which exhibited stronger affinity to S-vacant than pristine 2H-MoS₂, as elucidated by density functional theory calculations. Notably, SVs aggravated algal growth inhibition, oxidative damage, photosynthetic efficiency and cell membrane permeability reduction induced by 2H-MoS₂ due to increased free radical yield and the specific binding of functional proteins. Our findings provide insights into the roles of SVs on the risk of MoS₂ while highlighting the importance of rational design for TMDs application.

Electrochemical exfoliation of graphite in nanofibrillated kenaf cellulose (NFC)/surfactant mixture for the development of conductive paper

The effect of incorporating common dodecyl anionic and cationic surfactants such as dodecyltrimethylammonium bromide (DTAB), dodecylethyldimethylammonium bromide (DDAB), and sodium dodecylsulfate (SDS) in nanocomposites of reduced graphene oxide and nanocellulose are described. The stabilization and electrical properties of the nanocomoposites of reduced graphene oxide (RGO) and nanofibrillated kenaf cellulose (NFC) were characterized using four-point probe electrical conductivity measurements. Raman spectroscopy, field emission scanning electron microscopy, and high-resolution transmission electron microscopy were used to investigate dispersion morphology and the quality of RGO inside the NFC matrices. Small-angle neutron scattering (SANS) was used to study the aggregation behavior of the aqueous surfactant systems and RGO dispersions. The cationic surfactant DTAB proved to be the best choice for stabilization of RGO in NFC, giving enhanced electrical conductivity five orders of magnitude higher than the neat NFC. The results highlight the effects of hydrophilic surfactant moieties on the structure, stability and properties of RGO/NFC composites.

From GROMACS to LAMMPS: GRO2LAM : A converter for molecular dynamics software

Atomistic simulations have progressively attracted attention in the study of physical-chemical properties of innovative nanomaterials. GROMACS and LAMMPS are currently the most widespread open-source software for molecular dynamics simulations thanks to their good flexibility, numerous functionalities and responsive community support. Nevertheless, the very different formats adopted for input and output files are limiting the possibility to transfer GROMACS simulations to LAMMPS. In this article, we present GRO2LAM, a modular and open-source Python 2.7 code for rapidly translating input files and parameters from GROMACS to LAMMPS format. The robustness of the tool has been assessed by comparing the simulation results obtained by GROMACS and LAMMPS, after the format conversion by GRO2LAM. Specifically, three nanoscale configurations of interest in both engineering and biomedical fields are studied, namely a carbon nanotube, an iron oxide nanoparticle, and a protein immersed in water. In perspective, GRO2LAM may be the first step to achieve a full interoperability between molecular dynamics software. This would allow to easily exploit their complementary potentialities and post-processing functionalities. Moreover, GRO2LAM could facilitate the cross-check of simulation results, guaranteeing the reproducibility of molecular dynamics models and testing their robustness. Graphical Abstract GRO2LAM, a modular and open-source Python code for rapidly translating input files and parameters from GROMACS to LAMMPS format.

Liquids with Lower Wettability Can Exhibit Higher Friction on Hexagonal Boron Nitride: The Intriguing Role of Solid-Liquid Electrostatic Interactions

We investigate the wetting and frictional behavior of polar (water and ethylene glycol) and nonpolar (diiodomethane) liquids on the basal plane of hexagonal boron nitride (hBN) using molecular dynamics simulations. Our results for the wettability of water on the hBN basal plane (contact angle 81°) are in qualitative agreement with the experimentally deduced mild hydrophilicity of the hBN basal plane (contact angle 66°). We find that water exhibits the lowest wettability, as quantified by the highest contact angle, but the highest friction coefficient of (1.9 ± 0.4) × 10⁵ N-s/m³ on the hBN basal plane among the three liquids considered. This intriguing finding is explained in terms of the competition between dispersion and electrostatic interactions operating between the hBN basal plane and the three liquids. We find that electrostatic interactions do not affect the wetting behavior appreciably, as quantified by a less than 3° change in the respective contact angles of the three liquids considered. On the other hand, electrostatic interactions are found to increase the friction coefficients of the three liquids in contact with hBN to different extents, indicating that despite the increased friction of water on hBN, relative to that on graphene, nonpolar liquids may exhibit similar friction coefficients on hBN and graphene. Our findings reveal that the increase in the friction coefficient, upon incorporation of solid-liquid electrostatic interactions, is brought about by a greater increase in the solid-liquid mean-squared total lateral force, as compared to a smaller reduction in the decorrelation time of the solid-liquid force.

Exfoliation of Two-Dimensional Materials: The Role of Entropy

Liquid-phase exfoliation (LPE) is the best-known method for the synthesis of two-dimensional (2D) nanosheets. Compared to enthalpy, entropy is hardly considered to be a factor in choosing energy-efficient solvents and has not even been verified to be negligible. In this Letter, we explore the entropy contribution in LPE by performing molecular dynamics (MD) simulation of the structural flexibility effect in graphene, hexagonal boron nitride (hBN), and molybdenum disulfide (MoS₂). Our results show that surface vibration favors the exfoliation of graphene and hBN and destabilizes the reaggregation of nanosheets in water at 300 K, whereas the opposite is found for MoS₂. The entropy change is found to be 41%, 48%, and 4% of the enthalpy gain for graphene, hBN, and MoS₂ in LPE, respectively, and 64%, 32%, and 56% in reaggregation, which amounts to a step advancement for solvent screening in LPE of 2D materials.

Initial Studies Directed toward the Rational Design of Aqueous Graphene Dispersants

This study presents preliminary experimental data suggesting that sodium 4-(pyrene-1-yl)butane-1-sulfonate (PBSA), 5, an analogue of sodium pyrene-1-sulfonate (PSA), 1, enhances the stability of aqueous reduced graphene oxide (RGO) graphene dispersions. We find that RGO and exfoliated graphene dispersions prepared in the presence of 5 are approximately double the concentration of those made with commercially available PSA, 1. Quantum mechanical and molecular dynamics simulations provide key insights into the behavior of these molecules on the graphene surface. The seemingly obvious introduction of a polar sulfonate head group linked via an appropriate alkyl spacer to the aromatic core results in both more efficient binding of 5 to the graphene surface and more efficient solvation of the polar head group by bulk solvent (water). Overall, this improves the stabilization of the graphene flakes by disfavoring dissociation of the stabilizer from the graphene surface and inhibiting reaggregation by electrostatic and steric repulsion. These insights are currently the subject of further investigations in an attempt to develop a rational approach to the design of more effective dispersing agents for rGO and graphene in aqueous solution.

Stabilization of aqueous graphene dispersions utilizing a biocompatible dispersant: a molecular dynamics study

Investigation of the high efficiency of flavin mononucleotide sodium salt (FMNS) for the stabilization of aqueous graphene dispersions using all-atom molecular dynamics simulations.

Preparation of conductive cellulose paper through electrochemical exfoliation of graphite: The role of anionic surfactant ionic liquids as exfoliating and stabilizing agents

A facile electrochemical exfoliation method was established to efficiently prepare conductive paper containing reduced graphene oxide (RGO) with the help of single chain anionic surfactant ionic liquids (SAILs). The surfactant ionic liquids are synthesized from conventional organic surfactant anions and a 1-butyl-3-methyl-imidazolium cation. For the first time the combination of SAILs and cellulose was used to directly exfoliate graphite. The ionic liquid 1-butyl-3-methyl-imidazolium dodecylbenzenesulfonate (BMIM-DBS) was shown to have notable affinity for graphene, demonstrating improved electrical properties of the conductive cellulose paper. The presence of BMIM-DBS in the system promotes five orders of magnitude enhancement of the paper electrical conductivity (2.71 × 10^-5 S cm^-1) compared to the native cellulose (1.97 × 10^-10 S cm^-1). A thorough investigation using electron microscopy and Raman spectroscopy highlights the presence of uniform graphene incorporated inside the matrices. Studies into aqueous aggregation behavior using small-angle neutron scattering (SANS) point to the ability of this compound to act as a bridge between graphene and cellulose, and is responsible for the enhanced exfoliation level and stabilization of the resulting dispersion. The simple and feasible process for producing conductive paper described here is attractive for the possibility of scaling-up this technique for mass production of conductive composites containing graphene or other layered materials.

Ab Initio Molecular Dynamics and Lattice Dynamics-Based Force Field for Modeling Hexagonal Boron Nitride in Mechanical and Interfacial Applications

Hexagonal boron nitride (hBN) is an up-and-coming two-dimensional material, with applications in electronic devices, tribology, and separation membranes. Herein, we utilize density-functional-theory-based ab initio molecular dynamics (MD) simulations and lattice dynamics calculations to develop a classical force field (FF) for modeling hBN. The FF predicts the crystal structure, elastic constants, and phonon dispersion relation of hBN with good accuracy and exhibits remarkable agreement with the interlayer binding energy predicted by random phase approximation calculations. We demonstrate the importance of including Coulombic interactions but excluding 1-4 intrasheet interactions to obtain the correct phonon dispersion relation. We find that improper dihedrals do not modify the bulk mechanical properties and the extent of thermal vibrations in hBN, although they impact its flexural rigidity. Combining the FF with the accurate TIP4P/Ice water model yields excellent agreement with interaction energies predicted by quantum Monte Carlo calculations. Our FF should enable an accurate description of hBN interfaces in classical MD simulations.