Fragmentation of polypropylene into microplastics promoted by photo-aging; release of metals, toxicity and inhibition of biodegradability.

The widespread presence of microplastics (MP) in water represents an environmental problem, not only because of the harmful effects of their size and potential to vector other pollutants, but also because of the release of additives, degradation products and residues contained in the polymer matrix. The latter includes metallic catalysts, which are often overlooked. This study focuses on the photo-aging of polypropylene (PP) and the resulting structural changes that promote its fragmentation microplastics (PP-MPs) and release of metals, as well as the resulting toxicity of leachates and their potential to inhibit biodegradation of organics in water. The pristine, photo-aged and waste PP are ground under the same regime to assess susceptibility to fragmentation. Obtained PP-MPs are submitted to leaching tests; the release of organics and metals is monitored by Total Organic Carbon (TOC) and Inductively Coupled Plasma Mass Spectrometry (ICP-MS) analysis, respectively. The leachates are assessed for their toxicity against Vibrio fischeri, Daphnia magna and Pseudokirchneriella subcapitata and their influence on the biodegradability of the glucose solution. Photo-aging induced changes in the crystallinity and morphology of the PP and manifested in the abundance of smaller MPs, as revealed by the particle size distribution. In the case of pristine PP, all particles were > 100 µm in size, while aged PP yielded significant mass fraction of MPs <100 µm. The toxicity of leachates from aged PP-MPs is higher than that of pristine and exhibits a positive correlation with portion of metals released. The biodegradability of glucose is strongly inhibited by PP-MPs leachates containing a mixture of metals in trace concentrations.

An Atomic-Level Perspective on the interactions between Organic Pollutants and PET particles: A Comprehensive Computational Investigation.

Microplastics (MPs) have recently attracted a lot of attention worldwide due to their abundance and potentially harmful effects on the environment and on human health. One of the factors of concern is their ability to adsorb and disperse other harmful organic pollutants in the environment. To properly assess the adsorption capacity of MP for organic pollutants in different environments, it is pivotal to understand the mechanisms of their interactions in detail at the atomic level. In this work, we studied interactions between polyethylene terephthalate (PET) MP and small organic pollutants containing different functional groups within the framework of density functional theory (DFT). Our computational outcomes show that organic pollutants mainly bind to the surface of a PET model via weak non-bonding interactions, mostly hydrogen bonds. The binding strength between pollutant molecules and PET particles strongly depends on the adsorption site while we have found that the particle size is of lesser importance. Specifically, carboxylic sites are able to form strong hydrogen bonds with pollutants containing hydrogen bond donor or acceptor groups. On the other hand, it is found that in such kind of systems π-π interactions play a minor role in adsorption on PET particles.

trans-Hydroalkynylation of Internal 1,3-Enynes Enabled by Cooperative Catalysis.

A cooperative catalyst system involving a Pd(0)/Senphos complex, tris(pentafluorophenyl)borane, copper bromide, and an amine base, is demonstrated to catalyze trans-hydroalkynylation of internal 1,3-enynes. For the first time, a Lewis acid catalyst is shown to promote the reaction involving the emerging outer-sphere oxidative reaction step. The resulting cross-conjugated dieneynes are versatile synthons for organic synthesis, and their characterization reveals distinct photophysical properties depending on the positioning of the donor/acceptor substituents along the conjugation path.

Structural Features Promoting Photocatalytic Degradation of Contaminants of Emerging Concern: Insights into Degradation Mechanism Employing QSA/PR Modeling.

Although heterogeneous photocatalysis has shown promising results in degradation of contaminants of emerging concern (CECs), the mechanistic implications related to structural diversity of chemicals, affecting oxidative (by HO•) or reductive (by O2•-) degradation pathways are still scarce. In this study, the degradation extents and rates of selected organics in the absence and presence of common scavengers for reactive oxygen species (ROS) generated during photocatalytic treatment were determined. The obtained values were then brought into correlation as K coefficients (MHO•/MO2•-), denoting the ratio of organics degraded by two occurring mechanisms: oxidation and reduction via HO• and O2•-. The compounds possessing K >> 1 favor oxidative degradation over HO•, and vice versa for reductive degradation (i.e., if K << 1 compounds undergo reductive reactions driven by O2•-). Such empirical values were brought into correlation with structural features of CECs, represented by molecular descriptors, employing a quantitative structure activity/property relationship (QSA/PR) modeling. The functional stability and predictive power of the resulting QSA/PR model was confirmed by internal and external cross-validation. The most influential descriptors were found to be the size of the molecule and presence/absence of particular molecular fragments such as C - O and C - Cl bonds; the latter favors HO•-driven reaction, while the former the reductive pathway. The developed QSA/PR models can be considered robust predictive tools for evaluating distribution between degradation mechanisms occurring in photocatalytic treatment.

A new method to analyze and understand molecular linear and nonlinear optical responses via field-induced functions: a straightforward alternative to sum-over-states (SOS) analysis.

The sum-over-states (SOS) method allows the computation of polarizabilities and hyperpolarizabilities additively from the contributions of different electronic excited states in a given molecule or cluster. Subsequent analysis of the main excited configurations contributing to the relevant excited states allows characterizing the orbitals involved in the linear and nonlinear optical response. Unfortunately, the chemically relevant information that can be obtained by SOS is hindered by a series of methodological and computational drawbacks. Among these drawbacks, we can highlight the high computational cost, problems arising from nonconvergent series and errors caused by the inaccurate description of excitation energies and/or higher excited state matrix elements. For this reason, coupled-perturbed schemes are currently widely used to determine the NLO potential of molecules and materials. However, such a choice limits the amount of intuitive chemical information that, on the other hand, can be retrieved by a successful SOS computation. In this work, we present and discuss a novel computational strategy that offers the means to extract the useful chemical insights from a coupled-perturbed calculation at almost negligible extra computational cost providing a transparent picture about orbital contributions to the properties of interest. The proposed method is based on the generation and further analysis of field-induced orbitals, FIOs, from the analytic or numerical derivatives of the dipole moment. Orbital symmetry rules are derived using group theory and the method is tested for a series of small and medium size systems.

Theoretical Investigation of the Infrared Spectrum of 5-Bromo-2,4-pentadiynenitrile from a CCSD(T)/B3LYP Anharmonic Potential.

On the grounds of a hybrid CCSD(T)/B3LYP/aug-cc-pVTZ anharmonic potential and the use of a variational and variational-perturbational methods, the IR spectra of 5-bromo-2,4-pentadiynenitrile (BrC5 N) is revisited in the mid-infrared region up to 4500 cm-1 . A position and intensity analysis of our theoretical results allow us to assign the fundamental bands together with their combinations and overtones, in the aforementioned range of frequencies. The main objective of this work is to give an "a priori" complete IR spectrum of BrC5 N, which can be used as a guide for the low-intensity bands in areas not completely studied so far.

Prediction of biodegradability of aromatics in water using QSAR modeling.

The study was aimed at developing models for predicting the biodegradability of aromatic water pollutants. For that purpose, 36 single-benzene ring compounds, with different type, number and position of substituents, were used. The biodegradability was estimated according to the ratio of the biochemical (BOD5) and chemical (COD) oxygen demand values determined for parent compounds ((BOD5/COD)0), as well as for their reaction mixtures in half-life achieved by UV-C/H2O2 process ((BOD5/COD)t1/2). The models correlating biodegradability and molecular structure characteristics of studied pollutants were derived using quantitative structure-activity relationship (QSAR) principles and tools. Upon derivation of the models and calibration on the training and subsequent testing on the test set, 3- and 5-variable models were selected as the most predictive for (BOD5/COD)0 and (BOD5/COD)t1/2, respectively, according to the values of statistical parameters R2 and Q2. Hence, 3-variable model predicting (BOD5/COD)0 possessed R2=0.863 and Q2=0.799 for training set, and R2=0.710 for test set, while 5-variable model predicting (BOD5/COD)1/2 possessed R2=0.886 and Q2=0.788 for training set, and R2=0.564 for test set. The selected models are interpretable and transparent, reflecting key structural features that influence targeted biodegradability and can be correlated with the degradation mechanisms of studied compounds by UV-C/H2O2.

Establishing the pivotal role of local aromaticity in the electronic properties of boron-nitride graphene lateral hybrids.

Within an attempt to unravel the conundrum of irregular bandgap variations in hybrids of white-graphene (hBN) and graphene (G) observed in both experiment and theory, strong proofs about the decisive role of aromaticity in their electronic properties are brought to light. Sound numerical experiments conducted on zero-, one- and two-dimensional hBNG hybrids demonstrate that upon structural and/or electronic perturbation caused by foreign doping agents, the uniformity in local cyclic electron delocalization of ideal graphene restructures locally creating carbon hexagons of contrasting cyclic electron delocalization (c.c. local aromatic patterns) which may dominate the bandgap size of the resulting systems. In addition, relying on the quantum chemical aspect of aromaticity in terms of quantitative computations of cyclic electron delocalization together with pictorial intrinsic polarizability density representations, this work provides a solid and handy rule-of-thumb to be used in qualitative and intuitive predictions. According to this empirical rule, the origin of any nonmonotonic bandgap variation observed in stoichiometric 0D (BN)n/graphene hybrids with increasing hBN segment lies in instabilities caused by partially substituted benzenoid rings formed locally at the hBNG interfaces. This relationship, established in 0D graphene flakes and extended to 1D periodic ribbons, can be used to understand and qualitatively predict conflicting bandgap variations of vacancy-free 2D periodic lattices, pointing at the property of aromaticity as the missing link needed to solve the puzzle of conflicting bandgap variations in hBNG hybrids observed in experiment.

A Computational Study of the Interaction and Polarization Effects of Complexes Involving Molecular Graphene and C60 or a Nucleobases.

A systematic analysis of the molecular structure, energetics, electronic (hyper)polarizabilities and their interaction-induced counterparts of C60 with a series of molecular graphene (MG) models, CmHn, where m = 24, 84, 114, 222, 366, 546 and n = 12, 24, 30, 42, 54, 66, was performed. All the reported data were computed by employing density functional theory and a series of basis sets. The main goal of the study is to investigate how alteration of the size of the MG model affects the strength of the interaction, charge rearrangement, and polarization and interaction-induced polarization of the complex, C60-MG. A Hirshfeld-based scheme has been employed in order to provide information on the intrinsic polarizability density representations of the reported complexes. It was found that the interaction energy increases approaching a limit of -26.98 kcal/mol for m = 366 and 546; the polarizability and second hyperpolarizability increase with increasing the size of MG. An opposite trend was observed for the dipole moment. Interestingly, the variation of the first hyperpolarizability is relatively small with m. Since polarizability is a key factor for the stability of molecular graphene with nucleobases (NB), a study of the magnitude of the interaction-induced polarizability of C84H24-NB complexes is also reported, aiming to reveal changes of its magnitude with the type of NB. The binding strength of C84H24-NB complexes is also computed and found to be in agreement with available theoretical and experimental data. The interaction involved in C60 B12N12H24-NB complexes has also been considered, featuring the effect of contamination on the binding strength between MG and NBs.

Hirshfeld-based intrinsic polarizability density representations as a tool to analyze molecular polarizability.

In this work, a general scheme to visualize polarizability density distributions is proposed and implemented in a Hirshfeld-based partitioning scheme. This allows us to obtain easy-to-interpret pictorial representations of both total and intrinsic polarizabilities where each point of the density is formed by the contribution of any atom or group of atoms in the molecule. In addition, the procedure used here permits the possibility of removing the size dependence of the electric-dipole polarizability. Such a development opens new horizons in exploring new applications for the analysis of the molecular polarizability tensor. For instance, this visualization shows which atoms or regions are more polarizable distinguishing, moreover, the fine structure of atoms affected by the vicinity, and might extend the dipole polarizability as a tool for aromaticity studies in polycyclic aromatic hydrocarbons. Additionally, this approach can serve us to assess the methods performance in describing the interaction of electric fields with a molecule and local electron correlation effects in intrinsic polarizabilities.

Unleashing the quadratic nonlinear optical responses of graphene by confining white-graphene (h-BN) sections in its framework.

In an attempt to diversify the options in designing graphene-based systems bearing large second order nonlinear optical (NLO) responses of octupolar and/or dipolar character, the subject of the quadratic NLO properties of hybrid boron nitride (BN) graphene flakes is opened up. State of the art ab initio and density functional theory methods applied on a toolbox of book-text octupolar and arbitrary dipolar planar hybrid h-BN-graphene nanosized systems reveal that by confining finite h-BN sections in the internal network of graphene, the capacity of the π-electron network of graphene species in delivering giant second order NLO responses could be fully exploited. Configuration interaction (CIS) and time-dependent density functional (TD) computations, within the sum-overstate (SOS) perturbational approach, expose that the prevailing (hyper)polarization mechanism, lying under the sizable computed octupolar hyperpolarizabilities, is fueled by alternating positive and negative atomic charges located in the internal part of the hybrid flakes, and more precisely at the BN/graphene intersections. This type of charge transfer mechanism distinguishes, in fact, the elemental graphene dipoles/octupoles we report here from other conventional NLO dipoles or octupoles. More interestingly, it is shown that by controlling the shape, size, and covering area of the h-BN domain (or domains), one can effectively regulate "à volonté" both the magnitudes and types of the second order NLO responses switching from dipolar to octupolar and vice versa. Especially in the context of the latter class of NLO properties, this communication brings into surface novel, graphene-based, octupolar planar or quasiplanar motifs. The take home message of this communication is summarized as follows: When the right BN segment is incorporated in the right section of the right graphene flake, systems of giant quadratic NLO octupolar and/or dipolar responses may emerge.

Structural and static electric response properties of highly symmetric lithiated silicon cages: theoretical predictions.

It is shown by density functional theory calculations that high symmetry silicon cages can be designed by coating with Li atoms. The resulting highly symmetric lithiated silicon cages (up to D(5d) symmetry) are low-lying true minima of the energy hypersurface with binding energies of the order of 4.6 eV per Si atom and moderate highest occupied molecular orbital-lowest unoccupied molecular orbital gaps. Moreover, relying on a systematic study of the electric response properties obtained by ab initio (Hartree-Fock, MP2, and configuration interaction singles (CIS)) and density functional (B3LYP, B2PLYP, and CAM-B3LYP) methods, it is shown that lithium coating has a large impact on the magnitude of their second hyperpolarizabilities resulting to highly hyperpolarizable species. Such hyperpolarizable character is directly connected to the increase in the density of the low-lying excited states triggered by the interaction between the Si cage and the surrounding Li atoms.

Doping-enhanced hyperpolarizabilities of silicon clusters: a global ab initio and density functional theory study of Si10 (Li, Na, K)n (n=1, 2) clusters.

A global theoretical study of the (hyper)polarizabilities of alkali doped Si(10) is presented and discussed. First, a detailed picture about the low lying isomers of Si(10)Li, Si(10)Na, Si(10)K, Si(10)Li(2), Si(10)Na(2), and Si(10)K(2) has been obtained in a global manner. Then, the microscopic first (hyper)polarizabilities of the most stable configurations have been determined by means of ab initio methods of high predictive capability such as those based on the Møller-Plesset perturbation and coupled cluster theory, paying extra attention to the (hyper)polarizabilities of the open shell mono-doped systems Si(10)Li, Si(10)Na, Si(10)K, and the influence of spin contamination. These results were used to assess the performance of methods of low computational cost based on density functional theory (DFT) in the reliable computation of these properties in order to proceed with an in-depth study of their evolution as a function of the alkali metal, the cluster composition, and the cluster structure. The most interesting outcomes of the performed (hyper)polarizability study indicate that while alkali doping leaves the per atom polarizability practically unaffected, influences dramatically the hyperpolarizabilities of Si(10). The lowest energy structures of the mono-doped clusters are characterized by significantly enhanced hyperpolarizabilities as compared to the analogue neutral or charged bare silicon clusters Si(10) and Si(11), while, certain patterns governed by the type and the number of the doping agents are followed. The observed hyperpolarizability increase is found to be in close connection with specific cluster to alkali metal charge transfer excited states and to the cluster structures. Moreover, an interesting correlation between the anisotropy of the electron density, and the hyperpolarizabilities of these systems has been observed. Finally, it is important to note that the presented method assessment points out that among the various DFT functionals used in this work, (B3LYP, B3PW91, BhandHLYP, PBE0, CAM-B3LYP, LC-BLYP, LC-BPW91) only B3PW91 and PBE0 out of the seven provided a consistent quantitative performance for both polarizabilities and hyperpolarizabilities with respect to the ab initio methods utilized here. On the other hand, the long range corrected functionals LC-(U)BLYP and LC-(U)BPW91 (µ = 0.47) failed to supply quantitatively accurate hyperpolarizability results in all the studied clusters while the CAM-(U)B3LYP functional performs satisfactory only in the case of the Na and K doped systems.

Electric dipole (hyper)polarizabilities of selected X2Y2 and X3Y3 (X = Al, Ga, In and Y = P, As): III-V semiconductor clusters. An ab initio comparative study.

A systematic ab initio comparative study of the (hyper)polarizabilities of selected III-V stoichiometric semiconductor clusters has been carried out. Our investigation focuses on the ground state structures of the dimers and on two dissimilar trimer configurations of aluminum, gallium, indium phosphide and arsenide. The basis set effect on both the polarizabilities and hyperpolarizabilities of the studied systems has been explicitly taken into account relying on the augmented correlation consistent aug-cc-pVnZ (n = D, T, Q, and 5) basis sets series. In addition, a rough estimation of the effects of the relativistic effects on the investigated properties is provided by extension of the study to include calculations performed with relativistic electron core potentials (or pseudopotentials). Electron correlation effects have been estimated utilizing methods of increasing predictive reliability, e.g., the Møller-Plesset many body perturbation theory and the couple cluster approach. Our results reveal that in the considered semiconductor species the Group III elements (Al, Ga, In) play a vital role on the values of their relative (hyper)polarizability. At all levels of theory employed the most hyperpolarizable clusters are the indium derivatives while the aluminum arsenide clusters also exhibit high, comparable hyperpolarizabilities. The less hyperpolarizable species are those composed of gallium and this is associated with the strong influence of the nuclear charge on the valence electrons of Ga due to the poor shielding that is provided by the semicore d electrons. In addition, the analysis of the electronic structure and the hyperpolarizability magnitudes reveals that clusters, in which their bonding is characterized by strong electron transfer from the electropositive to the electronegative atoms, are less hyperpolarizable than species in which the corresponding electron transfer is weaker. Lastly, from the methodological point of view our results point out that the hyperpolarizabilities of those species converge when an augmented triple-zeta quality basis set is used and, also, that the second order Møller-Plesset approximation (MP2) overestimates considerably their second hyperpolarizabilities with respect to the highest level of coupled cluster theory applied in this study (CCSD(T)).

Polarizability evolution on natural and artificial low dimensional binary semiconductor systems: A case study of stoichiometric aluminum phosphide semiconductor clusters.

The dependences of the static dipole polarizabilities per atom (PPAs) on the bonding and shape of selected stoichiometric aluminum phosphide clusters (ground states and higher lying species) of small and medium sizes have been comprehensively studied at Hartree-Fock and the second order Moller-Plesset perturbation levels of theory. It is shown that the nonmonotonic size variations in the mean PPAs of AlP species which maintain closed cagelike structures, frequently observed in clusters, are directly related to covalent homoatomic bonds inside each cluster's framework. Accordingly, the PPAs of clusters which are characterized by one or more bonds between the Al and P atoms are larger than the PPAs of clusters with the uniform alternating Al-P bond matrix. This is caused by the electron transfer increase from the electropositive Al to the electronegative P atom with the cluster growth. This transfer is larger for the clusters characterized by alternating Al-P bonding. The later effect explains the decrease in the PPA of AlP species which maintain closed cage-like structures, with the cluster growth. However, this picture drastically changes for artificial metastable prolate species built up by the ground states of smaller clusters. It is demonstrated that for prolate binary AlP clusters of medium size, the shape dominates against any other structural or bonding factor, forcing the PPA to increase with the cluster size. Nonetheless, as the cluster size grows, it is predicted that the PPAs of the studied prolate clusters will saturate eventually with the cluster size. Also, it is verified that the theoretical predicted polarizabilities of AlP semiconductor clusters are larger than the bulk polarizability in accord with other theoretical predictions for similar systems. Lastly, it is pointed out that major bonding or structural changes should take place in order the convergence with the bulk polarizability to be accomplished since it is revealed that the size increase is a necessary but not a sufficient factor for the cluster to bulk transition.

Correlations between bonding, size, and second hyperpolarizability (gamma) of small semiconductor clusters: ab initio study on Al(n)P(n) clusters with n=2, 3, 4, 6, and 9.

A comprehensive investigation of the correlation between the second hyperpolarizability and the bonding and structural characteristics of stoichiometric aluminum phosphide clusters up to 18 atoms is presented. Several aluminum phosphide species displaying different types of configurations and bonding have been studied. The obtained ab initio and density functional finite field results suggest that the ionic AlP clusters are considerably less hyperpolarizable than the covalent bonded species. Other structural features such as symmetry, atoms' arrangement, and shape also play an important role on the hyperpolarizabilities of those species. However, they are only noticeable among clusters characterized by the same bonding patterns. Furthermore, the results of this study demonstrate that the bonding which is determined by the atoms' arrangement of a cluster has a more profound effect on the second hyperpolarizability than the cluster's composition or size. In addition, the mean second hyperpolarizability increases with the increasing number of atoms, assuming that the bonding characteristics among the clusters of increasing size are similar. On the other hand, the hyperpolarizability per atom rapidly decreases with the increase of atoms' number in the cluster and converges to values of approximately 900e(4)a(0)(4) and approximately 1300e(4)a(0)(4)E(h) (-3) at the HF/cc-pVDZ and MP2/cc-pVDZ levels of theory respectively. Lastly, this work provides the first systematic study on the hyperpolarizabilities of small aluminum phosphide clusters which, in their covalent forms, exhibit larger second order hyperpolarizabilities than the well studied small gallium arsenide clusters.

Ab initio finite field (hyper)polarizability computations on stoichiometric gallium arsenide clusters GanAsn (n=2-9).

We report reliable ab initio finite field (hyper)polarizability values at Hartree-Fock and second order Moller-Plesset perturbation theory (MP2) levels of theory for different geometrical configurations of small gallium arsenide clusters Ga(n)As(n) with n=2-5. We relied on all-electron basis sets and pseudopotentials suitable for (hyper)polarizability calculations. In each case, we used structures that have been established in the literature after we optimized their geometries at B3LYP/cc-pVTZ-PP level of theory. Our results suggest that the first order hyperpolarizability (beta) is much more sensitive to the special geometric features than the second order hyperpolarizability (gamma). For the most stable configurations up to ten atoms the second order hyperpolarizability at MP2 level of theory varies between 15 x 10(4) and 32 x 10(4) e(4)a0 (4)Eh(-3). In addition, we examined the polarizability per atom evolution versus the cluster size for Ga(n)As(n) with n=2-9. Our work extends earlier theoretical studies which were limited to eight atoms and exposes that the polarizability/atom of the most stable stoichiometric configurations up to Ga(9)As(9) continues the monotonic downward trend with increasing size. Lastly, from the methodological point of view, our analysis shows that apart from polarizabilities, augmented pseudopotentials yield reliable first and second hyperpolarizability values as well.

Hyperpolarizability of GaAs dimer is not negative.

We present a systematic study of the static electric hyperpolarizability of Ga(2)As(2). The authors rely on finite-field high-level ab initio calculations with carefully optimized basis sets. Their best values for the mean and the anisotropy of the dipole polarizability are alpha=158.57 and Deltaalpha=130.33e(2)a(0) (2)E(h) (-1). For the hyperpolarizability we propose an estimate gamma=(155+/-15)x10(3)e(4)a(0) (4)E(h) (-3), which does not agree with the negative value predicted by Lan et al. [J. Chem. Phys. 124, 094302 (2006)]. Density functional theory based methods yield values close to those predicted by conventional ab initio methods. The (hyper)polarizability components are particularly enhanced along the direction defined by the Ga-Ga axis.

Molecular geometry and polarizability of small cadmium selenide clusters from all-electron ab initio and Density Functional Theory calculations.

We have calculated molecular geometries and electric polarizabilities for small cadmium selenide clusters. Our calculations were performed with conventional ab initio and density functional theory methods and Gaussian-type basis sets especially designed for (CdSe)(n). We find that the dipole polarizability per atom converges rapidly to the bulk value.