Modelling of PFAS-surface interactions: Effect of surface charge and solution ions

PER-: and poly-fluoroalkyl substances (PFAS) are a class of substances of increasing concern as environmental contaminants. The interactions between PFAS and surfaces play an important role in PFAS transport and remediation. Previous studies have found PFAS adsorption to be dependent upon properties including pH, organic matter and particle size, along with PFAS functional group and carbon chain length. It is hypothesised that a theoretical examination of PFAS-surface interactions, via Monte Carlo molecular simulation, would show differences resulting from changes in surface charge, H⁺, OH^-, Ca²⁺ concentrations and PFAS carbon chain length. Monte Carlo molecular simulations of perfluorooctane and perfluorobutane sulfonic acids interacting with a graphite surface in an aqueous medium were performed. Variations in surface charge, H⁺, OH^- and Ca²⁺ concentrations were made. The distance-dependent density of molecules from the surface was analysed as a proxy for PFAS adsorption to the surface. Simulation results showed differences in surface behaviour that depended on surface charge, H⁺, OH^- and Ca²⁺ concentrations, along with carbon chain length, with surface charge playing the most prominent role in controlling PFAS adsorption. For negatively charged surfaces, adsorption due to divalent cation bridging was observed in Ca²⁺ solutions. Modelling, such as in this study, of the thermodynamic equilibrium behaviour of low concentrations of molecules, in scenarios where both adsorption and mobility of PFAS occur, can aid in the design and testing of sorptive surfaces for amendment-based PFAS remediation.

Microstructural and Dynamical Heterogeneities in Ionic Liquids

Ionic liquids (ILs) are a special category of molten salts solely composed of ions with varied molecular symmetry and charge delocalization. The versatility in combining varied cation-anion moieties and in functionalizing ions with different atoms and molecular groups contributes to their peculiar interactions ranging from weak isotropic associations to strong, specific, and anisotropic forces. A delicate interplay among intra- and intermolecular interactions facilitates the formation of heterogeneous microstructures and liquid morphologies, which further contributes to their striking dynamical properties. Microstructural and dynamical heterogeneities of ILs lead to their multifaceted properties described by an inherent designer feature, which makes ILs important candidates for novel solvents, electrolytes, and functional materials in academia and industrial applications. Due to a massive number of combinations of ion pairs with ion species having distinct molecular structures and IL mixtures containing varied molecular solvents, a comprehensive understanding of their hierarchical structural and dynamical quantities is of great significance for a rational selection of ILs with appropriate properties and thereafter advancing their macroscopic functionalities in applications. In this review, we comprehensively trace recent advances in understanding delicate interplay of strong and weak interactions that underpin their complex phase behaviors with a particular emphasis on understanding heterogeneous microstructures and dynamics of ILs in bulk liquids, in mixtures with cosolvents, and in interfacial regions.

Enhanced extraction of phenol from model oils using ionic liquids elucidated with neutron diffraction

Aromatic cation ionic liquids (ILs) based on alkylpyridiniums are shown to be good phenol extractants from model oils (hexane/toluene). ILs with hard basic anions are found to have best extraction efficiency consistent with tetraalkylammonium salts ([NR4]X). Key extraction interactions were analysed using small angle neutron diffraction. Trifluoromethanesulphonate ([OTf]- or triflate) anions provide the synergistic effects of reduced cation-phenol centre of mass (COM) distances and increased hydrogen bonding that are linked to the improved extraction efficiency. Increases in cation electron density (methylpyridinium ([Me-Py]+) vs. methylpicolinium ([Me-3-Pic]+)) also reduce cation-phenol COM interaction lengths consistent with small increases in extraction efficiency for the same ionic liquids.

Applying neutron diffraction with isotopic substitution to the structure and proton-transport pathways in protic imidazolium bis{(trifluoromethyl)sulfonyl}imide ionic liquids

Neutron diffraction with isotopic substitution has been applied to examine the potential for complex-ion formation in protic imidazolium bis{(trifluoromethyl)sulfonyl}imide ionic liquids.

Effect of cation structure on the oxygen solubility and diffusivity in a range of bis{(trifluoromethyl)sulfonyl}imide anion based ionic liquids for lithium–air battery electrolytes

Electrochemical measurements of oxygen mobility and solubility in a series of related [TFSI]⁻-based ILs relating to cation structure and physical properties of the ILs.

Molecular simulations of imidazolium-based tricyanomethanide ionic liquids using an optimized classical force field

An atomistic force field is optimized to accurately predict the equilibrium and transport properties of technologically important imidazolium-based tricyanomethanide ionic liquids.

Electrical conductivity and glass formation in nitrile-functionalized pyrrolidinium bis(trifluoromethylsulfonyl)imide ionic liquids: chain length and odd-even effects of the alkyl spacer between the pyrrolidinium ring and the nitrile group

The electrical conductivity of a series of pyrrolidinium bis(trifluoromethylsulfonyl)imide ionic liquids, functionalized with a nitrile (cyano) group at the end of an alkyl chain attached to the cation, was studied in the temperature range between 173 K and 393 K. The glass formation of the ionic liquids is influenced by the length of the alkyl spacer separating the nitrile function from the pyrrolidinium ring. The electrical conductivity and the viscosity do not show a monotonic dependence on the alkyl spacer length, but rather an odd-even effect. An explanation for this behavior is given, including the potential energy landscape picture for the glass transition.

Structure of ionic liquids with amino acid anions via neutron diffraction

Total neutron scattering studies on amino acid-based ionic liquids reveal significant directional interactions between anions, especially for serinate-derived ions.

Association and liquid structure of pyridine–acetic acid mixtures determined from neutron scattering using a ‘free proton’ EPSR simulation model

Hydrogen-bonded molecular acetic acid chains are observed in acid–base mixtures from small angle neutron diffraction.

4-Cyano-1-methyl-pyridinium perchlor-ate

The title salt, C7H7N2 (+)·ClO4 (-), crystallizes with alternating cations and anions in wavy sheets, which are formed by a number of C-H⋯O and C-H⋯N hydrogen bonds, lying approximately parallel to (001).

4-Cyano-1-methyl-pyridinium nitrate

The title mol-ecular salt, C7H7N2 (+)·NO3 (-), displays an inter-penetrating sheet structure parallel to a with each sheet containing nearly coplanar cations and anions, each ion being bis-ected by a crystallographic mirror plane. C-H⋯O hydrogen bonds involving both ring and methyl H atoms in addition to cation-cation C-H⋯N hydrogen bonds (ring H to cyano N) serve to link the sheets together. In each set of parallel layers, the cations and anions stack with short distances of 3.094 (2) (between aligned nitrate N and pyridine N atoms) and 3.057 (2) Å (between a nitrate O atom and the ring centroid). This motif is strikingly similar to the one that features in the isomeric salt 2-cyano-1-methyl-pyridinium nitrate.

Comparing reduced partial charge models with polarizable simulations of ionic liquids

Molecular ionic liquids are typically characterized by strong electrostatic interactions resulting in a charge ordering and retardation of their translational and rotational behaviour. Unfortunately, this effect is often overestimated in classical molecular dynamics simulations. This can be circumvented in a twofold way: the easiest way is to reduce the partial charges of the ions to sub-integer values of ±0.7-0.9 e. The more realistic model is to include polarizable forces, e.g. Drude-oscillators, but it comes along with an increasing computational effort. On the other hand, charge-scaled models are claimed to take an average polarizability into account. But do both models have the same impact on structure and dynamics of molecular ionic liquids? In the present study several molecular dynamics simulations of 1-ethyl-3-methylimidazolium trifluoromethanesulfonate are performed with different levels of polarization as well as with varying charge scaling factors of 0.74 to 0.90. The analysis of the structural and dynamical results are performed in different levels: from the atomic point of view over the molecular level to collective properties determined by the complete sample.

Molecular dynamics simulation of the interfacial structure of [C(n)mim][PF6] adsorbed on a graphite surface: effects of temperature and alkyl chain length

The structures and diffusion behaviors of a series of ionic liquids [C(n)mim][PF(6)] (n = 1, 4, 8 and 12) on a graphite surface have been investigated by means of molecular dynamics simulation. It was found that three or four ordering layers of ionic liquids were formed near the graphite surface, and this layering structure was stable over the temperature range investigated. At the liquid/vacuum interface, the ionic liquid with a butyl chain had a monolayer ordering surface, while [C(8)mim][PF(6)] and [C(12)mim][PF(6)] exhibited a bilayer ordering with a polar domain sandwiched between two orientational nonpolar domains. More impressively, the simulated results showed that for the ionic liquids with alkyl chains longer than C(4), the adjacent alkyl chains in the whole film tended to be parallel to each other, with the imidazolium rings packed closely together. This indicated that the ionic liquids have a better regulated short-range structure than was previously expected. It was also found that both in the bottom layer and in the bulk region, the diffusion of the alkyl chains was much faster than that of the polar groups. However, as the alkyl chain length increased, the charge delocalization in the cation and the enhanced van der Waals interaction between the nonpolar groups contributed by reducing this difference in the diffusivity of major groups.