Quantum Chemical Calculation and Evaluation of Partition Coefficients for Classical and Emerging Environmentally Relevant Organic Compounds

Temperature dependence of the rain-gas and snow-gas partition coefficients for nearly a thousand chemicals

Compared to the particle-gas partition coefficients (KPG), the rain-gas (KRG) and snow-gas (KSG) partition coefficients are also essential in studying the environmental behavior and fate of chemicals in the atmosphere. While the temperature dependence for the KPG have been extensively studied, the study for KRG and KSG are still lacking. Adsorption coefficients between water surface-air (KIA) and snow surface-air (KJA), as well as partition coefficients between water-air (KWA) and octanol-air (KOA) are vital in calculating KRG and KSG. These four basic adsorption and partition coefficients are also temperature-dependent, given by the well-known two-parameters Antoine equation logKXY = AXY+BXY/T, where KXY is the adsorption or partition coefficients, AXY and BXY are Antoine parameters (XY stand for IA, JA, WA, and OA), and T is the temperature in Kelvin. In this study, the parameters AXY and BXY are calculated for 943 chemicals, and logKXY can be estimated at any ambient temperature for these chemicals using these Antoine parameters. The results are evaluated by comparing these data with published experimental and modeled data, and the results show reasonable accuracy. Based on these coefficients, temperature-dependence of logKRG and logKSG is studied. It is found that both logKRG and logKSG are linearly related to 1/T, and Antoine parameters for logKRG and logKSG are also estimated. Distributions of the 943 chemicals in the atmospheric phases (gas, particle, and rain/snow), are illustrated in a Chemical Space Map. The findings reveal that, at environmental temperatures and precipitation days, the dominant state for the majority of chemicals is the gaseous phase. All the AXY and BXY values for logKSG, logKRG, and basic adsorption and partition coefficients, both modeled by this study and collected from published work, are systematically organized into an accessible dataset for public utilization.

Pesticide Identification Using Surface-Enhanced Raman Spectroscopy and Density Functional Theory Calculations: From Structural Insights to On-Site Detection

Pesticides play an important role in conventional agriculture. Yet, their harmful effects on the environment are becoming increasingly apparent. The occurrence of pesticides is hence being monitored worldwide. For fast, easy, yet sensitive identification, surface-enhanced Raman spectroscopy (SERS) is a powerful tool. In this study, a method is introduced that may be amended to in-field detection of pesticides. Gold and silver nanoparticles were synthesized, size-tailored, and characterized. The herbicide paraquat and the fungicide thiram served as model compounds. The preparation yielded reproducible SERS spectra. Using quantum chemical computation, Raman and SERS spectra were calculated and analyzed. The interpretation of vibrational modes in combination with SERS enhancement and attenuation allowed us to identify compound-specific bands. The assignment was interpreted in terms of the orientation of paraquat and thiram on the gold and silver nanoparticle surfaces. Paraquat preferred a co-planar arrangement parallel to the gold nanoparticle surface and a head-on orientation on the silver nanoparticle. For thiram, breaking of the disulfide bond was recognized, such that interaction with the surface occurred via the sulfur atoms. Successful detection of the pesticides after recollection from vegetable leaves demonstrated the method's applicability for pesticide identification.

QSPR modeling for the prediction of partitioning of VOCs and SVOCs to indoor fabrics: Integrating environmental factors

Porous fabrics have a significant impact on indoor air quality by adsorbing and emitting chemical substances, such as volatile organic compounds (VOCs) and semi-volatile organic compounds (SVOCs). Understanding the partition behavior between organic compound molecules and indoor fabrics is crucial for assessing their environmental fate and associated human exposure. The physicochemical properties of fabrics and compounds are fundamental in determining the free energy of partitioning. Moreover, environmental factors like temperature and humidity critically affect the partition process by modifying the thermal and moisture conditions of the fabric. However, existing methods for determining the fabric-air partition coefficient are limited to specific fabric-chemical combinations and lack a comprehensive consideration of indoor environmental factors. In this study, large amounts of experimental data on fabric-air partition coefficients (Kfa) of (S)VOCs were collected for silk, polyester, and cotton fabrics. Key molecular descriptors were identified, integrating the influences of physicochemical properties, temperature, and humidity. Subsequently, two typical quantitative structure-property relationship (QSPR) models were developed to correlate the Kfa values with the molecular descriptors. The fitting performance, robustness, and predictive ability of the two QSPR models were evaluated through statistical analysis and internal/external validation. This research provides insights for the high-throughput prediction of the environmental behaviors of indoor organic compounds.

The legalization of cannabis may result in increased indoor exposure to Δ9-tetrahydrocannabinol

Cannabis is a genus of plants in the Cannabaceae family that contains tetrahydrocannabinolic acid. When heated or burned, the acid decarboxylates to form tetrahydrocannabinol (THC). Its (-)-trans-Δ9-THC isomer is a psychoactive substance that has been used as a drug for centuries. In most countries, both the private sale of cannabis products and their use for non-medical purposes are still prohibited by law. However, for some time now there has been societal and political pressure to at least partially legalize cannabis products. It can be expected that such a measure will lead to a significant increase in the consumption of cannabis. However, this also increases the possibility of involuntary passive exposure to THC and contamination of the indoor environment. In indoor sciences, THC is still a largely unknown or underrepresented substance. In this perspective paper, THC will therefore first be presented on the basis of its physical properties. Then, the distribution of THC in different indoor compartments and potential routes of passive exposure are discussed. Finally, an assessment of the future importance of THC for indoor use is made. Previous experience has shown that early monitoring is always advantageous so that preventive and protective measures can be taken quickly if necessary.

Synthesis and absolute configuration of cyclic synthetic cathinones derived from α-tetralone

The emergence of new synthetic cathinones continues to be a matter of public health concern. In fact, already known products (drugs) are being rapidly replaced by new structurally related alternatives, whereby modifications in the basic cathinone structure are used by manufacturers to circumvent the legislation. On the other hand, some derivatives of synthetic cathinones represent important pharmaceuticals with antidepressant properties. In the search for pharmaceutically relevant analogs, the main goal of the present study was to design and characterize novel cyclic α-tetralone-based derivatives of synthetic cathinones. We synthesized a series of derivatives and verified their chemical structure. Subsequently, chiral separation has been accomplished by high-performance liquid chromatography (HPLC) equipped with a circular dichroism (CD) detector, which directly provided CD spectra of the enantiomers of the analyzed substances at 252 nm. Using density functional theory calculations, we have obtained stable conformers of selected enantiomers in solution and their relative abundances, which we used to simulate their spectra. The experimental and calculated data have been used to assign the absolute configuration of six as-yet unknown synthetic cathinones.

Building Chemical Intuition about Physicochemical Properties of C8-Per-/Polyfluoroalkyl Carboxylic Acids through Computational Means

We have predicted acid dissociation constants (pK a), octanol-water partition coefficients (K OW), and DMPC lipid membrane-water partition coefficients (K lipid-w) of 150 different eight-carbon-containing poly-/perfluoroalkyl carboxylic acids (C8-PFCAs) utilizing the COnductor-like Screening MOdel for Realistic Solvents (COSMO-RS) theory. Different trends associated with functionalization, degree of fluorination, degree of saturation, degree of chlorination, and branching are discussed on the basis of the predicted values for the partition coefficients. In general, functionalization closest to the carboxylic headgroup had the greatest impact on the value of the predicted physicochemical properties.

Cloth-Air Partitioning of Neutral Per- and Polyfluoroalkyl Substances (PFAS) in North Carolina Homes during the Indoor PFAS Assessment (IPA) Campaign

Partitioning of per- and polyfluoroalkyl substances (PFAS) to indoor materials, including clothing, may prolong the residence time of PFAS indoors and contribute to exposure. During the Indoor PFAS Assessment (IPA) Campaign, we measured concentrations of nine neutral PFAS in air and cotton cloth in 11 homes in North Carolina, for up to 9 months. Fluorotelomer alcohols (i.e., 6:2 FTOH, 8:2 FTOH, and 10:2 FTOH) are the dominant target species in indoor air, with concentrations ranging from 1.8 to 49 ng m-3, 1.2 to 53 ng m-3, and 0.21 to 5.7 ng m-3, respectively. In cloth, perfluorooctane sulfonamidoethanols (i.e., MeFOSE and EtFOSE) accumulated most significantly over time, reaching concentrations of up to 0.26 ng cm-2 and 0.24 ng cm-2, respectively. From paired measurements of neutral PFAS in air and suspended cloth, we derived cloth-air partition coefficients (Kca) for 6:2, 8:2, and 10:2 FTOH; ethylperfluorooctane sulfonamide (EtFOSA); MeFOSE; and EtFOSE. Mean log(Kca) values range from 4.7 to 6.6 and are positively correlated with the octanol-air partition coefficient. We investigated the effect of the cloth storage method on PFAS accumulation and the influence of home characteristics on air concentrations. Temperature had the overall greatest effect. This study provides valuable insights into PFAS distribution, fate, and exposure indoors.

Enantio- and Diastereomerically Pure Titanocenes by Dynamic Conformational Locking

The synthesis of enantiomerically pure titanocenes is limited to cases with enantiomerically pure substituents at the cyclopentadienyl ligands and to ansa-titanocenes. For the latter complexes, the use of achiral ligands requires a resolution of enantiomers and frequently also a separation of the diastereoisomers obtained after metalation. Here, we introduce a new synthetic method that relies on the use of enantiomerically pure camphorsulfonate (CSA) ligands as control elements for the absolute and relative configuration of titanocene complexes. Starting from the conformationally flexible (RC5 H4 )2 TiCl2 , the desired conformationally locked and hence enantio- and diastereomerically pure complexes (RC5 H4 )2 Ti(CSA)2 are obtained in just two steps. According to X-ray crystallography the (RC5 H4 )2 Ti fragment is essentially C2 -symmetric and nuclear magnetic resonance displays overall C2 -symmetry. We applied density functional theory methods to unravel the dynamics of the complexes and the mechanisms and selectivities of their formation.

Extended Conductor-like Polarizable Continuum Solvation Model (CPCM-X) for Semiempirical Methods

We have developed a new method to accurately account for solvation effects in semiempirical quantum mechanics based on a polarizable continuum model (PCM). The extended conductor-like polarizable continuum model (CPCM-X) incorporates a computationally efficient domain decomposition conductor-like screening model (ddCOSMO) for extended tight binding (xTB) methods and uses a post-processing approach based on established solvation models, like the conductor-like screening model for real solvents (COSMO-RS) and the universal solvent model based on solute electron density (SMD). According to various benchmarks, the approach performs well across a broad range of systems and applications, including hydration free energies, non-aqueous solvation free energies, and large supramolecular association reactions of neutral and charged species. Our method for computing solvation free energies is much more accurate than the current methods in the xtb program package. It improves the accuracy of solvation free energies by up to 40% for larger supramolecular association reactions to match even the accuracy of higher-level DFT-based solvation models like COSMO-RS and SMD while being computationally more than 2 orders of magnitude faster. The proposed method and the underlying ddCOSMO model are readily available for a wide variety of solvents and are accessible in xtb for use in various computational applications.

Molecular Dynamics Simulation Prediction of the Partitioning Constants (KH, Kiw, Kia) of 82 Legacy and Emerging Organic Contaminants at the Water-Air Interface

The tendency of organic contaminants (OCs) to partition between different phases is a key set of properties that underlie their human and ecological health impacts and the success of remediation efforts. A significant challenge associated with these efforts is the need for accurate partitioning data for an ever-expanding list of OCs and breakdown products. All-atom molecular dynamics (MD) simulations have the potential to help generate these data, but existing studies have applied these techniques only to a limited variety of OCs. Here, we use established MD simulation approaches to examine the partitioning of 82 OCs, including many compounds of critical concern, at the water-air interface. Our predictions of the Henry's law constant (K_H) and interfacial adsorption coefficients (K_iw, K_ia) correlate strongly with experimental results, indicating that MD simulations can be used to predict K_H, K_iw, and K_ia values with mean absolute deviations of 1.1, 0.3, and 0.3 logarithmic units after correcting for systematic bias, respectively. A library of MD simulation input files for the examined OCs is provided to facilitate future investigations of the partitioning of these compounds in the presence of other phases.

Indazole-derived synthetic cannabinoids: Absolute configuration determination and structure characterization by circular dichroism and DFT calculations

An increasing number of products containing synthetic cannabinoids pose a growing crisis to public health worldwide. Recently, a rising number of cases of serious adverse health effects, intoxications, and death cases associated with synthetic cannabinoids were reported. The current study represents the comprehensive structural analysis of three new synthetic cannabinoids (AB-, ADB- and AMB-FUBINACA) in solution investigated by electronic and vibrational circular dichroism together with the conventional methods of infrared and ultraviolet absorption spectroscopy, all supported by the density functional theory (DFT) calculations. The best level of theory to reproduce the experimental wavenumbers and wavelengths was found to be the B3PW91 method with a 6-311++G(d,p) basis set including the implicit solvent effect simulation. Very good agreement between the experimental and simulated spectra allowed us to determine the absolute configuration and a detailed interpretation of the IR absorption, VCD, ECD and UV spectra of AB-, ADB- and AMB-FUBINACA. In addition, the HOMO and LUMO electronic transitions were calculated.

A kinetic study of the photolysis of sulfamethoxazole with special emphasis on the photoisomer

Abstract

The previously not studied photochemical degradation of sulfamethoxazole (SMX) to the isomer of SMX (ISO) was measured via a polychromatic (Xe) and a monochromatic (Hg) light source and accompanied by quantum chemical DFT calculations. In addition to the $$\mathrm{p}K_\mathrm{a} = \;7.0 \pm 0.1$$ p K a = 7.0 ± 0.1 of ISO, tautomer-dependent properties such as the $$K_\mathrm{OW}$$ K OW were measured and theoretically confirmed by DFT. The kinetics in solutions below and above the $$\mathrm{p}K_\mathrm{a} = 5.6$$ p K a = 5.6 of SMX were studied for the available and quantifiable products SMX, ISO, 3-amino-5-methylisoxazole (AMI), 2-amino-5-methyloxazole (AMO), and sulfanilic acid (SUA). The quantum yields of the neutral ($$\Phi _\mathrm{N}$$ Φ N ) and anionic $$\Phi _\mathrm{A}$$ Φ A ) forms of SMX ($$\Phi _\mathrm{A} = 0.03 \pm 0.001$$ Φ A = 0.03 ± 0.001 , $$\Phi _\mathrm{N} = 0.15 \pm 0.01$$ Φ N = 0.15 ± 0.01 ) and ISO ($$\Phi _\mathrm{A} = 0.05 \pm 0.01$$ Φ A = 0.05 ± 0.01 and $$\Phi _\mathrm{N} = 0.06 \pm 0.02$$ Φ N = 0.06 ± 0.02 ) were found to be wavelength-independent. In a competitive reaction to the formation of ISO from SMX, the degradation product TP271 is formed. Various proposed structures for TP271 described in the literature have been studied quantum mechanically and can be excluded for thermodynamic reasons. In real samples in a northern German surface water in summer 2021 mean concentrations of SMX were found in the range of 120 ng/L. In agreement with the pH-dependent yields, concentrations of ISO were low in the range of 8 ng/L.

Graphical abstract

Quantum chemical calculation of the vapor pressure of volatile and semi volatile organic compounds

The vapor pressure is a specific and temperature-dependent parameter that describes the volatility of a substance and thus its driving force for evaporation or sublimation into the gas phase. Depending on the magnitude of the vapor pressure, there are different methods for experimental determination. However, these are usually associated with a corresponding amount of effort and become less accurate as the vapor pressure decreases. For purposes of vapor pressure prediction, algorithms were developed that are usually based on quantitative structure-activity relationships (QSAR). The quantum mechanical (QM) approach followed here applies an alternative, much less empirical strategy, where the change in Gibbs free energy for the transition from the condensed to the gas phase is obtained from conformer ensembles computed for each phase separately. The results of this automatic, so-called CRENSO workflow are compared with experimentally determined vapor pressures for a large set of environmentally relevant compounds. In addition, comparisons are made with the single structure-based COSMO-RS QM approach, linear-free-energy relationships (LFER) as well as results from the SPARC program. We show that our CRENSO workflow is superior to conventional prediction models and provides reliable vapor pressures for liquids and sub-cooled liquids over a wide pressure range.

Decipher the molecular descriptors and mechanisms controlling sulfonamide adsorption onto mesoporous carbon: Density functional theory calculation and partial least-squares path modeling

Mesoporous carbons (MCs) exhibit excellent removal efficiencies to various organic chemicals. However, how the properties of chemicals influence the adsorption mechanisms and further determine their adsorption onto MCs are poorly understood. We investigated the adsorption of 22 sulfonamides (SAs) onto four MCs, and further uncovered the major molecular descriptors and adsorption mechanisms influencing the adsorption by density functional theory (DFT) and partial least-squares path modeling (PLS-PM). The results revealed that the excess molar refraction (E), McGowan's molar volume (V), energy of the highest occupied molecular orbital (E_HOMO), hardness (H), and most positive net charge on carbon atom (Q_c⁺) were identified as the indirect factors affecting the distribution coefficient (logK_D), by influencing the BE(π-π), BE(H), and logK_ow. BE(π-π) and logK_ow displayed significant direct impacts on logK_D (p < 0.05), while BE(H) showed insignificant direct influences on logK_D (p > 0.05). The PLS-PM results indicate the main driving forces for SAs adsorption including π-π interactions, hydrophobic effects, and hydrogen bonding. This study provides a new perspective on revealing the adsorption mechanisms, and the identified factors can be used to develop the quantitative model to further predict the adsorption of SAs onto MCs.

A holistic modeling framework for estimating the influence of climate change on indoor air quality

The IPCC 2021 report predicts rising global temperatures and more frequent extreme weather events in the future, which will have different effects on the regional climate and concentrations of ambient air pollutants. Consequently, changes in heat and mass transfer between the inside and outside of buildings will also have an increasing impact on indoor air quality. It is therefore surprising that indoor spaces and occupant well-being still play a subordinate role in the studies of climate change. To increase awareness for this topic, the Indoor Air Quality Climate Change (IAQCC) model system was developed, which allows short and long-term predictions of the indoor climate with respect to outdoor conditions. The IAQCC is a holistic model that combines different scenarios in the form of submodels: building physics, indoor emissions, chemical-physical reaction and transformation, mold growth, and indoor exposure. IAQCC allows simulation of indoor gas and particle concentrations with outdoor influences, indoor materials and activity emissions, particle deposition and coagulation, gas reactions, and SVOC partitioning. These key processes are fundamentally linked to temperature and relative humidity. With the aid of the building physics model, the indoor temperature and humidity, and pollutant transport in building zones can be simulated. The exposure model refers to the calculated concentrations and provides evaluations of indoor thermal comfort and exposure to gaseous, particulate, and microbial pollutants.

Towards understanding solvation effects on the conformational entropy of non-rigid molecules

The absolute molecular entropy is a fundamental quantity for the accurate description of thermodynamic properties. For non-rigid molecules, a substantial part of the entropy can be attributed to a conformational contribution. Systems and properties where this is relevant, e.g., protein-ligand binding affinities or pK_a values refer usually to the liquid phase. In this work, the influence of solvation on the conformational entropy is investigated. A recently introduced state-of-the-art and automated computational protocol for the computation of conformational entropies [Pracht et al., Chem. Sci., 2021, 12, 6551-6568.] is applied in combination with fast and accurate semiempirical quantum-chemical methods and implicit solvation models for a set of 25 commercially available drug molecules and five transition metal compounds. Computed gas-phase conformational entropies are compared with values obtained in implicit n-hexane and water. It is found that implicit solvation can have a substantial effect of several cal mol^-1 K^-1 on the entropy as a result of large conformational changes in the different phases. We conclude that for flexible molecules chemical accuracy for free energies in solution can only be achieved if solvation effects on the conformational ensemble are considered.

Temperature and indoor environments

From the thermodynamic perspective, the term temperature is clearly defined for ideal physical systems: A unique temperature can be assigned to each black body via its radiation spectrum, and the temperature of an ideal gas is given by the velocity distribution of the molecules. While the indoor environment is not an ideal system, fundamental physical and chemical processes, such as diffusion, partitioning equilibria, and chemical reactions, are predictably temperature-dependent. For example, the logarithm of reaction rate and equilibria constants are proportional to the reciprocal of the absolute temperature. It is therefore possible to have non-linear, very steep changes in chemical phenomena over a relatively small temperature range. On the contrary, transport processes are more influenced by spatial temperature, momentum, and pressure gradients as well as by the density, porosity, and composition of indoor materials. Consequently, emergent phenomena, such as emission rates or dynamic air concentrations, can be the result of complex temperature-dependent relationships that require a more empirical approach. Indoor environmental conditions are further influenced by the thermal comfort needs of occupants. Not only do occupants have to create thermal conditions that serve to maintain their core body temperature, which is usually accomplished by wearing appropriate clothing, but also the surroundings must be adapted so that they feel comfortable. This includes the interaction of the living space with the ambient environment, which can vary greatly by region and season. Design of houses, apartments, commercial buildings, and schools is generally utility and comfort driven, requiring an appropriate energy balance, sometimes considering ventilation but rarely including the impact of temperature on indoor contaminant levels. In our article, we start with a review of fundamental thermodynamic variables and discuss their influence on typical indoor processes. Then, we describe the heat balance of people in their thermal environment. An extensive literature study is devoted to the thermal conditions in buildings, the temperature-dependent release of indoor pollutants from materials and their distribution in the various interior compartments as well as aspects of indoor chemistry. Finally, we assess the need to consider temperature holistically with regard to the changes to be expected as a result of global emergencies such as climate change.