The Limitations of Certain Density Functionals in Modeling Neutral Water Clusters

Atmospheric microplastics and nanoplastics as vectors of primary air pollutants - A theoretical study on the polyethylene terephthalate (PET) case

Polyethylene terephthalate (PET) microplastics and nanoplastics are ubiquitously present in the atmosphere as atmospheric and airborne forms (PET-aMPs). Using first-principles calculations, we analyze the uptake of primary air pollutants onto PET-aMPs, focusing on their stabilities, adsorption mechanisms, and thermochemistry. The results show that PET-aMPs are selective for the spontaneous adsorption of CO, CO₂, NO, N₂O, NO₂, NH_3, and SO₂, reaching stable adsorption energies of 6-20 kcal/mol per molecule, with comparable uptake ability than carbon-based materials, metals/metalloids, and metal oxide surfaces. Then, PET-aMPs become a vector for coexisting air pollutants in the atmosphere, which adsorb by inner or outer adsorption depending on the molecular polarity (dipole moment) and atomic constitution (electronegativity) of gaseous molecules. Also, atmospheric H₂O and O₂ are not competitive molecules, and ozone could enhance adsorption due to surface oxidation and structure breakdown. The interplay of electrostatic (46-61%) and dispersion forces (21-58%) drives the adsorption mechanism, where low-polar pollutants display almost a balanced electrostatic vs. dispersion contribution, while high polar molecules display a higher electrostatic stabilization. The outer adsorption is reached by strong dispersion, hydrogen bonding, and dipole-dipole-induced pairs, while lone-pair-π interactions appear in the inner adsorption regime. These results expand the understanding of the hazards and risks of atmospheric and airborne microplastics/nanoplastics, their impacts, co-transport ability, and interaction with the environment.

Monomers of Glycine and Serine Have a Limited Ability to Hydrate in the Atmosphere

The role of atmospheric aerosols on climate change is one of the biggest uncertainties in most global climate models. Organic aerosols have been identified as potential cloud condensation nuclei (CCN), and amino acids are organic molecules that could serve as CCN. Amino acids make up a significant portion of the total organic material in the atmosphere, and herein we present a systematic study of hydration for two of the most common atmospheric amino acids, glycine and serine. We compute DLPNO/CCSD(T)//M08-HX/MG3S thermodynamic properties and atmospheric concentrations of Gly(H₂O)_n and Ser(H₂O)_n, where n = 1-5. We predict that serine-water clusters have higher concentrations at n = 1 and 5, while glycine-water clusters have higher concentrations at n = 2-4. However, both glycine and serine are inferred to exist primarily in their nonhydrated monomer forms in the absence of other species such as sulfuric acid.

Calculating Reliable Gibbs Free Energies for Formation of Gas-Phase Clusters that Are Critical for Atmospheric Chemistry: (H2SO4)3

The effects of atmospheric aerosols on our climate are one of the biggest uncertainties in global climate models. Calculating the pathway for the formation of pre-nucleation clusters that become aerosols is challenging, requiring a comprehensive analysis of configurational space and highly accurate Gibbs free energy calculations. We identified a large set of minimum energy configurations of (H₂SO₄)₃ using a sampling technique based on a genetic algorithm and a stepwise density functional theory (DFT) approach and computed the thermodynamics of formation of these configurations with more accurate wavefunction-based electronic energies computed on the DFT geometries. The DLPNO-CCSD(T) methods always return more positive energies compared to the DFT energies. Within the DLPNO-CCSD(T) methods, extrapolating to the complete basis set limit gives more positive free energies compared to explicitly correlated single-point energies. The CBS extrapolation was shown to be robust as both the 4-5 inverse polynomial and Riemann zeta function schemes were within chemical accuracy of one another.

Elucidating the co-transport of bisphenol A with polyethylene terephthalate (PET) nanoplastics: A theoretical study of the adsorption mechanism

Polyethylene terephthalate (PET) is a possible key component of nanoplastics in water environments, which can migrate pollutants through co-transport. In this regard, the co-transport of endocrine disruptors (such as bisphenol A, BPA) by nanoplastics is of emergent concern because of its cytotoxicity/bioaccumulation effects in aquatic organisms. In this work, a computational study is performed to reveal the BPA adsorption mechanism onto PET nanoplastics (nanoPET). It is found that the outer surface of nanoPET has a nucleophilic nature, allowing to increase the mass transfer and intraparticle diffusion into the nanoplastic to form stable complexes by inner and outer surface adsorption. The maximum adsorption energy is similar (even higher) in magnitude with respect to nanostructured adsorbents such as graphene, carbon nanotubes, activated carbon, and inorganic surfaces, indicating the worrying adsorption properties of nanoPET. The adsorption mechanism is driven by the interplay of dispersion (38-49%) and electrostatics effects (43-50%); specifically, dispersion effects dominate the inner surface adsorption, while electrostatics energies dominate the outer surface adsorption. It is also determined that π-π stacking is not a reliable interaction mechanism for aromatics on nanoPET. The formed complexes are also highly soluble, and water molecules behave as non-competitive factors, establishing the high risk of nanoPET to adsorb and migrate pollutants in water ecosystems. Furthermore, the adsorption performance is decreased (but not inhibited) at high ionic strength in salt-containing waters. Finally, these results give relevant information for environmental risk assessment, such as quantitative data and interaction mechanisms for non-biodegradable nanoplastics that establish strong interactions with pollutants in water.

Assessing the accuracy of some popular DFT methods for computing harmonic vibrational frequencies of water clusters

A wide range of density functional theory (DFT) methods (37 altogether), including pure, hybrid, range-separated hybrid, double-hybrid, and dispersion-corrected functionals, have been employed to compute the harmonic vibrational frequencies of eight small water clusters ranging in size from the dimer to four different isomers of the hexamer. These computed harmonic frequencies have been carefully compared to recently published benchmark values that are expected to be very close to the CCSD(T) complete basis set limit. Of the DFT methods examined here, ωB97 and ωB97X are the most consistently accurate, deviating from the reference values by less than 20 cm(-1) on average and never more than 60 cm(-1). The performance of double-hybrid methods including B2PLYP and mPW2-PLYP is only slightly better than more economical approaches, such as the M06-L pure functional and the M06-2X hybrid functional. Additionally, dispersion corrections offer very little improvement in computed frequencies.

Solvation Energies of the Proton in Methanol

pKa's, proton affinities, and proton dissociation free energies characterize numerous properties of drugs and the antioxidant activity of some chemical compounds. Even with a higher computational level of theory, the uncertainty in the proton solvation free energy limits the accuracy of these parameters. We investigated the thermochemistry of the solvation of the proton in methanol within the cluster-continuum model. The scheme used involves up to nine explicit methanol molecules, using the IEF-PCM and the strategy based on thermodynamic cycles. All computations were performed at B3LYP/6-31++G(dp) and M062X/6-31++G(dp) levels of theory. It comes out from our calculations that the functional M062X is better than B3LYP, on the evaluation of gas phase clustering energies of protonated methanol clusters, per methanol stabilization of neutral methanol clusters and solvation energies of the proton in methanol. The solvation free energy and enthalpy of the proton in methanol were obtained after converging the partial solvation free energy of the proton in methanol and the clustering free energy of protonated methanol clusters, as the cluster size increases. Finally, the recommended values for the solvation free energy and enthalpy of the proton in methanol are -257 and -252 kcal/mol, respectively.

Accurate predictions of water cluster formation, (H₂O)(n=2-10)

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2-10 water molecules. This approach reveals new low energy conformers for (H(2)O)(n=7,9,10). The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.

Atmospheric implications for formation of clusters of ammonium and 1-10 water molecules

A mixed molecular dynamics/quantum mechanics model has been applied to the ammonium/water clustering system. The use of the high level MP2 calculation method and correlated basis sets, such as aug-cc-pVDZ and aug-cc-pVTZ, lends confidence in the accuracy of the extrapolated energies. These calculations provide electronic and free energies for the formation of clusters of ammonium and 1-10 water molecules at two different temperatures. Structures and thermodynamic values are in good agreement with previous experimental and theoretical results. The estimated concentration of these clusters in the troposphere was calculated using atmospheric amounts of ammonium and water. Results show the favorability of forming these clusters and implications for ion-induced nucleation in the atmosphere.