Nanoscale Simulation of Plastic Contaminants Migration in Packaging Materials and Potential Leaching into Model Food Fluids.

Polymers are the most commonly used packaging materials for nutrition and consumer products. The ever-growing concern over pollution and potential environmental contamination generated from single-use packaging materials has raised safety questions. Polymers used in these materials often contain impurities, including unreacted monomers and small oligomers. The characterization of transport properties, including diffusion and leaching of these molecules, is largely hampered by the long timescales involved in shelf life experiments. In this work, we employ atomistic molecular simulation techniques to explore the main mechanisms involved in the bulk and interfacial transport of monomer molecules from three polymers commonly employed as packaging materials: polyamide-6, polycarbonate, and poly(methyl methacrylate). Our simulations showed that both hopping and continuous diffusion play important roles in inbound monomer diffusion and that solvent-polymer compatibility significantly affects monomer leaching. These results provide rationalization for monomer leaching in model food formulations as well as bulky industry-relevant molecules. Through this molecular-scale characterization, we offer insights to aid in the design of polymer/consumer product interfaces with reduced risk of contamination and longer shelf life.

Molecular-Scale Exploration of Mechanical Properties and Interactions of Poly(lactic acid) with Cellulose and Chitin.

Poly(lactic acid) (PLA), one of the pillars of the current overarching displacement trend switching from fossil- to natural-based polymers, is often used in association with polysaccharides to increase its mechanical properties. However, the use of PLA/polysaccharide composites is greatly hampered by their poor miscibility, whose underlying nature is still vastly unexplored. This work aims to shed light on the interactions of PLA and two representative polysaccharide molecules (cellulose and chitin) and reveal structure-property relationships from a fundamental perspective using atomistic molecular dynamics. Our computational strategy was able to reproduce key experimental mechanical properties of pure and/or composite materials, reveal a decrease in immiscibility in PLA/chitin compared to PLA/cellulose associations, assert PLA-oriented polysaccharide reorientations, and explore how less effective PLA-polysaccharide hydrogen bonds are related to the poor PLA/polysaccharide miscibility. The connection between the detailed chemical interactions and the composite behavior found in this work is beneficial to the discovery of new biodegradable and natural polymer composite mixtures that can provide needed performance characteristics.

Molecular-Level Examination of Amorphous Solid Dispersion Dissolution.

Amorphous solid dispersions (ASDs) are commonly used to orally deliver small-molecule drugs that are poorly water-soluble. ASDs consist of drug molecules in the amorphous form which are dispersed in a hydrophilic polymer matrix. Producing a high-performance ASD is critical for effective drug delivery and depends on many factors such as solubility of the drug in the matrix and the rate of drug release in aqueous medium (dissolution), which is linked to bioperformance. Often, researchers perform a large number of design iterations to achieve this objective. A detailed molecular-level understanding of the mechanisms behind ASD dissolution behavior would aid in the screening, designing, and optimization of ASD formulations and would minimize the need for testing a wide variety of prototype formulations. Molecular dynamics and related types of simulations, which model the collective behavior of molecules in condensed phase systems, can provide unique insights into these mechanisms. To study the effectiveness of these simulation techniques in ASD formulation dissolution, we carried out dissipative particle dynamics simulations, which are particularly an efficient form of molecular dynamics calculations. We studied two stages of the dissolution process: the early-stage of the dissolution process, which focuses on the dissolution at the ASD/water interface, and the late-stage of the dissolution process, where significant drug release would have occurred and there would be a mixture of drug and polymer molecules in a predominantly aqueous environment. Experimentally, we used Fourier transform infrared spectroscopy to study the interactions between drugs, polymers, and water in the dry and wet states and the chromatographic technique to study the rate of drug and polymer release. Both experiments and simulations provided evidence of polymer microstructures and drug-polymer interactions as important factors for the dissolution behavior of the investigated ASDs, consistent with previous work by Pudlas et al. (Eur. J. Pharm. Sci.2015, 67, 21-31). As experimental and simulation results are consistent and complementary, it is clear that there is significant potential for combined experimental and computational research for a detailed understanding of ASD formulations and, hence, formulation optimization.

Understanding the apparent fractional charge of protons in the aqueous electrochemical double layer.

A detailed atomic-scale description of the electrochemical interface is essential to the understanding of electrochemical energy transformations. In this work, we investigate the charge of solvated protons at the Pt(111) | H2O and Al(111) | H2O interfaces. Using semi-local density-functional theory as well as hybrid functionals and embedded correlated wavefunction methods as higher-level benchmarks, we show that the effective charge of a solvated proton in the electrochemical double layer or outer Helmholtz plane at all levels of theory is fractional, when the solvated proton and solvent band edges are aligned correctly with the Fermi level of the metal (EF). The observed fractional charge in the absence of frontier band misalignment arises from a significant overlap between the proton and the electron density from the metal surface, and results in an energetic difference between protons in bulk solution and those in the outer Helmholtz plane.

How To Identify Plasmons from the Optical Response of Nanostructures.

A promising trend in plasmonics involves shrinking the size of plasmon-supporting structures down to a few nanometers, thus enabling control over light-matter interaction at extreme-subwavelength scales. In this limit, quantum mechanical effects, such as nonlocal screening and size quantization, strongly affect the plasmonic response, rendering it substantially different from classical predictions. For very small clusters and molecules, collective plasmonic modes are hard to distinguish from other excitations such as single-electron transitions. Using rigorous quantum mechanical computational techniques for a wide variety of physical systems, we describe how an optical resonance of a nanostructure can be classified as either plasmonic or nonplasmonic. More precisely, we define a universal metric for such classification, the generalized plasmonicity index (GPI), which can be straightforwardly implemented in any computational electronic-structure method or classical electromagnetic approach to discriminate plasmons from single-particle excitations and photonic modes. Using the GPI, we investigate the plasmonicity of optical resonances in a wide range of systems including: the emergence of plasmonic behavior in small jellium spheres as the size and the number of electrons increase; atomic-scale metallic clusters as a function of the number of atoms; and nanostructured graphene as a function of size and doping down to the molecular plasmons in polycyclic aromatic hydrocarbons. Our study provides a rigorous foundation for the further development of ultrasmall nanostructures based on molecular plasmonics.

Heterometallic antenna-reactor complexes for photocatalysis.

Metallic nanoparticles with strong optically resonant properties behave as nanoscale optical antennas, and have recently shown extraordinary promise as light-driven catalysts. Traditionally, however, heterogeneous catalysis has relied upon weakly light-absorbing metals such as Pd, Pt, Ru, or Rh to lower the activation energy for chemical reactions. Here we show that coupling a plasmonic nanoantenna directly to catalytic nanoparticles enables the light-induced generation of hot carriers within the catalyst nanoparticles, transforming the entire complex into an efficient light-controlled reactive catalyst. In Pd-decorated Al nanocrystals, photocatalytic hydrogen desorption closely follows the antenna-induced local absorption cross-section of the Pd islands, and a supralinear power dependence strongly suggests that hot-carrier-induced desorption occurs at the Pd island surface. When acetylene is present along with hydrogen, the selectivity for photocatalytic ethylene production relative to ethane is strongly enhanced, approaching 40:1. These observations indicate that antenna-reactor complexes may greatly expand possibilities for developing designer photocatalytic substrates.

Aluminum Nanocrystals as a Plasmonic Photocatalyst for Hydrogen Dissociation.

Hydrogen dissociation is a critical step in many hydrogenation reactions central to industrial chemical production and pollutant removal. This step typically utilizes the favorable band structure of precious metal catalysts like platinum and palladium to achieve high efficiency under mild conditions. Here we demonstrate that aluminum nanocrystals (Al NCs), when illuminated, can be used as a photocatalyst for hydrogen dissociation at room temperature and atmospheric pressure, despite the high activation barrier toward hydrogen adsorption and dissociation. We show that hot electron transfer from Al NCs to the antibonding orbitals of hydrogen molecules facilitates their dissociation. Hot electrons generated from surface plasmon decay and from direct photoexcitation of the interband transitions of Al both contribute to this process. Our results pave the way for the use of aluminum, an earth-abundant, nonprecious metal, for photocatalysis.

Density Fitting and Cholesky Decomposition of the Two-Electron Integrals in Local Multireference Configuration Interaction Theory.

To treat large molecules with accurate ab initio quantum chemistry, reduced scaling correlated wave function methods are now commonly employed. Optimization of these wave functions in practice requires some approximation of the two-electron integrals. Both Cholesky decomposition (CD) and density fitting (DF) are widely used approaches to approximate these integrals. Here, we compare CD and DF for use in local multireference singles and doubles configuration interaction (LMRSDCI). DF-LMRSDCI provides less accurate total energies than CD-LMRSDCI, but both methods are accurate for energy differences. However, DF-LMRSDCI is significantly less computationally expensive than CD-LMRSDCI on the molecules tested, suggesting that DF-LMRSDCI is an efficient, often sufficiently accurate alternative to our previously reported CD-LMRSDCI method.

Ab Initio Kinetics of Hydrogen Abstraction from Methyl Acetate by Hydrogen, Methyl, Oxygen, Hydroxyl, and Hydroperoxy Radicals.

The kinetics of hydrogen abstraction by five radicals (H, O((3)P), OH, CH3, and HO2) from methyl acetate (MA) is investigated theoretically in order to gain further understanding of certain aspects of the combustion chemistry of biodiesels, such as the effect of the ester moiety. We employ ab initio quantum chemistry methods, coupled cluster singles and doubles with perturbative triples correction (CCSD(T)) and multireference averaged coupled pair functional theory (MRACPF2), to predict chemically accurate reaction energetics. Overall, MRACPF2 predicts slightly higher barrier heights than CCSD(T) for MA + H/CH3/O/OH, but slightly lower barrier heights for hydrogen abstraction by HO2. Based on the obtained reaction energies, we also report high-pressure-limit rate constants using transition state theory (TST) in conjunction with the separable-hindered-rotor approximation, the variable reaction coordinate TST, and the multi-structure all-structure approach. The fitted modified Arrhenius expressions are provided over a temperature range of 250 to 2000 K. The predictions are in good agreement with available experimental results. Abstractions from both of the methyl groups in MA are expected to contribute to consumption of the fuel as they exhibit similar rate coefficients. The reactions involving the OH radical are predicted to have the highest rates among the five abstracting radicals, while those initiated by HO2 are expected to be the lowest.

Plasmons in molecules: microscopic characterization based on orbital transitions and momentum conservation.

In solid state physics, electronic excitations are often classified as plasmons or single-particle excitations. The former class of states refers to collective oscillations of the electron density. The random-phase approximation allows for a quantum-theoretical treatment and a characterization on a microscopic level as a coherent superposition of a large number of particle-hole transitions with the same momentum transfer. However, small systems such as molecules or small nanoclusters lack the basic properties (momentum conservation and uniform exchange interaction) responsible for the formation of plasmons in the solid-state case. Despite an enhanced interest in plasmon-based technologies and an increasing number of studies regarding plasmons in molecules and small nanoclusters, their definition on a microscopic level of theory remains ambiguous. In this work, we analyze the microscopic properties of molecular plasmons in comparison with the homogeneous electron gas as a model system. Subsequently, the applicability of the derived characteristics is validated by analyzing the electronic excitation vectors with respect to orbital transitions for two linear polyenes within second order versions of the algebraic diagrammatic construction scheme for the polarization propagator.

Ultrafast branching in the excited state of coumarin and umbelliferone.

In the present work we have explored the ultrafast relaxation network of coumarin and umbelliferone (7-hydroxy-coumarin) using time-resolved femtosecond spectroscopy and quantum chemical calculations. Despite the importance of the photophysical properties of coumarin derivatives for applications in biomedicine, the low fluorescence quantum yield of coumarin itself has not been fully understood so far. On the basis of our combined experimental and theoretical results we suggest a model for the ultrafast decay after photoexcitation incorporating two parallel radiationless relaxation pathways: one within the initially excited state via ring opening and the other one by transition into a dark state along the carbonyl stretching mode. The fluorescence quantum yield is determined by the position of the branching point relative to the Franck-Condon region which is strongly influenced by interactions with the environment and the substitution pattern. This model is finally capable of giving a comprehensive account of the striking differences observed in the photophysical behavior of coumarin as opposed to umbelliferone.

Application of the scaled-opposite-spin approximation to algebraic diagrammatic construction schemes of second order.

With the concept of scaled-opposite-spin (SOS), a pragmatic semi-empirical approximation has been introduced to the extended algebraic diagrammatic construction scheme of second order (ADC(2)-x) that leads to a significant saving in computational effort. The parameters included were fitted with respect to a benchmark set of electronically excited states in standard organic molecules that include some doubly-excited states, as well. Like the original, unscaled ADC(2)-x scheme it can be used to identify electronically excited states with high double excitation character, however at reduced computational cost. At the same time, it is possible to reduce the overestimation of doubly-excited configurations that is inherent to ADC(2)-x. Additionally, a scheme for the strict variant (ADC(2)-s) was derived directly from SOS-MP2 by application of the intermediate state formalism and compared to an existing version of SOS-ADC(2)-s.

A configurational switch based on iridium-catalyzed allylic cyclization: application in asymmetric total syntheses of prosopis, dendrobate, and spruce alkaloids.

A method for the stereoselective synthesis of 2,6-disubstituted piperidines has been developed that is based on the use of an intramolecular iridium-catalyzed allylic substitution as a configurational switch. The procedure allows the preparation of 2-vinylpiperidines with enantiomeric excesses (ee) of greater than 99%. As applications, total syntheses of piperidine alkaloids have been elaborated, most often by using Ru-catalyzed cross-metatheses as a key step for introduction of a side chain. Asymmetric total syntheses of the prosopis alkaloids (+)-prosopinine, (+)-prosophylline, (+)-prosopine, and of the dendrobate alkaloid (+)-241D and its C6 epimer are described.

Stereoselective synthesis of 2,6-disubstituted piperidines using the iridium-catalyzed allylic cyclization as configurational switch: asymmetric total synthesis of (+)-241 d and related piperidine alkaloids.

Pressing the configurational switch: Use of enantiomeric Ir catalysts allows the vinylpiperidine building blocks 2 a and 2 b to be synthesized with high selectivity. Total syntheses of the dendrobate alkaloid (+)-241 D, its C6-epimer, and a spruce alkaloid are presented as applications.