The 50 most-cited articles in The Journal of Physical Chemistry A, 2016

TiO 2 /g-C 3 N 4 photocatalysts with the ratio of TiO 2 to g-C 3 N 4 ranging from 0.3/1 to 2/1 were prepared by simple mechanical mixing of pure g-C 3 N 4 and commercial TiO 2 Evonik P25. All the nanocomposites were characterized by X-ray powder diffraction, UV–vis diffuse reflectance spectroscopy, photoluminescence, X-ray photoelectron spectroscopy, Raman spectroscopy, infrared spectroscopy, transmission electron microscopy, photoelectrochemical measurements, and nitrogen physisorption. The prepared mixtures along with pure TiO 2 and g-C 3 N 4 were tested for the photocatalytic reduction of carbon dioxide and photocatalytic decomposition of nitrous oxide. The pure g-C 3 N 4 exhibited the lowest photocatalytic activity in both cases, pointing to a very high recombination rate of charge carriers. On the other hand, the most active photocatalyst toward all the products was (0.3/1)TiO 2 /g-C 3 N 4 . The highest activity is achieved by combination of a number of factors: (i) specific surface area, (ii) adsorption edge energy, (iii) crystallite size, and (iv) efficient separation of the…

Careful analysis and comparison of the available acidity data of HCl, HBr, HI, HClO4, and CF3SO3H in water, dimethyl sulfoxide (DMSO), and gas-phase has been carried out. The data include experimental and computational pKa and gas-phase acidity data from the literature, as well as high-level computations using different approaches (including the W1 theory) carried out in this work. As a result of the analysis, for every acid in every medium, a recommended acidity value is presented. In some cases, the currently accepted pKa values were revised by more than 10 orders of magnitude.

The pKa's of substituted thiols are important for understanding their properties and reactivities in applications in chemistry, biochemistry, and material chemistry. For a collection of 175 different density functionals and the SMD implicit solvation model, the average errors in the calculated pKa's of methanethiol and ethanethiol are almost 10 pKa units higher than for imidazole. A test set of 45 substituted thiols with pKa's ranging from 4 to 12 has been used to assess the performance of 8 functionals with 3 different basis sets. As expected, the basis set needs to include polarization functions on the hydrogens and diffuse functions on the heavy atoms. Solvent cavity scaling was ineffective in correcting the errors in the calculated pKa's. Inclusion of an explicit water molecule that is hydrogen bonded with the H of the thiol group (in neutral) or S(-) (in thiolates) lowers error by an average of 3.5 pKa units. With one explicit water and the SMD solvation model, pKa's calculated with the M06-2X, PBEPBE, BP86, and LC-BLYP functionals are found to deviate from the experimental…

Anatase:rutile TiO2 junctions are often shown to be more photocatalytically active than anatase or rutile alone, but the underlying cause of this improvement is not fully understood. Herein, we employ transient absorption spectroscopy to study hole transfer across the anatase:rutile heterojunction in films as a function of phase composition. By exploiting the different signatures in the photoinduced absorption of trapped charges in anatase and rutile, we were able to separately track the yield and lifetime of holes in anatase and rutile sites within phase composites. Photogenerated holes transfer from rutile to anatase on submicrosecond time scales. This hole transfer can significantly increase the anatase hole yield, with a 20:80 anatase:rutile composite showing a 5-fold increase in anatase holes observed from the microsecond. Hole transfer does not result in an increase in charge-carrier lifetime, where an intermediate recombination dynamic between that of pure anatase (t1/2 ≈ 0.5 ms) and rutile (t1/2 ≈ 20 ms) is found in the anatase:rutile junction (t1/2 ≈ 4 ms). Irrespective of…

Although polymeric graphitic carbon nitride (g-C3N4) has been widely studied as metal-free photocatalyst, the description of its structure still remains a great challenge. Fourier transform infrared (FTIR) spectroscopy can provide complementary structural information. In this paper, we reconsider the representative crystalline melamine and develop a strategic approach to theoretically calculate the IR vibrations of this triazine-based nitrogen-rich system. IR calculations were based on three different models: a single molecule, a 4-molecule unit cell, and a 32-molecule cluster, respectively. By this comparative study the contribution of the intermolecular weak interactions were elucidated in detail. An accurate and visualized description on the experimental FTIR spectrum has been further presented by a combinatorial vibration-mode assignment based on the calculated potential energy distribution of the 32-molecule cluster. The theoretical approach reported in this study opens the way to the facile and accurate assignment for IR vibrational modes of other complex triazine-based…

Iodide-based chemical ionization mass spectrometry (CIMS) has been used to detect and measure concentrations of several atmospherically relevant organic and inorganic compounds. The significant electronegativity of iodide and the strong acidity of hydroiodic acid makes electron transfer and proton abstraction essentially negligible, and the soft nature of the adduct formation ionization technique reduces the chances of sample fragmentation. In addition, iodide has a large negative mass defect, which, when combined with the high resolving power of a high resolution time-of-flight chemical ionization mass spectrometer (HR-ToF-CIMS), provides good selectivity. In this work, we use quantum chemical methods to calculate the binding energies, enthalpies and free energies for clusters of an iodide ion with a number of atmospherically relevant organic and inorganic compounds. Systematic configurational sampling of the free molecules and clusters was carried out at the B3LYP/6-31G* level, followed by subsequent calculations at the PBE/SDD and DLPNO-CCSD(T)/def2-QZVPP//PBE/aug-cc-pVTZ-PP…

CHO, we have developed a cost-effective approach to the practical implementation of multiconformer transition state theory for peroxy radical hydrogen shift reactions at atmospherically relevant temperatures. While conformer searching is crucial for accurate reaction rates, an energy cutoff can be used to significantly reduce the computational cost with little loss of accuracy. For the reaction barrier, high-level calculations are needed, but the highest level of electronic structure theory is not necessary for the relative energy between conformers. Improving the approach to both transition state theory and electronic structure theory decreases the calculated reaction rate significantly, so low-level calculations can be used to rule out slow reactions. Further computational time can be saved by approximating the tunneling coefficients for each transition state by only that of the lowest-energy transition state. Finally, we test and validate our approach using higher-level theoretical values for our test system and existing experimental results for additional peroxy radical hydrogen…

Intramolecular singlet fission (ISF) is a promising photophysical process to construct more efficient light energy conversion systems as one excited singlet state converts into two excited triplet states. Herein we synthesized and evaluated bent-shaped pentacene dimers as a prototype of ISF to reveal intrinsic characters of triplet states (e.g., lifetimes of triplet excited states). In this study, meta-phenylene-bridged TIPS-pentacene dimer (PcD-3Ph) and 2,2'-bipheynyl bridged TIPS-pentacene dimer (PcD-Biph) were newly synthesized as bent-shaped dimers. In the steady-state spectroscopy, absorption and emission bands of these dimers were fully characterized, suggesting the appropriate degree of electronic coupling between pentacene moieties in these dimers. In addition, the electrochemical measurements were also performed to check the electronic interaction between two pentacene moieties. Whereas the successive two oxidation peaks owing to the delocalization were observed in a directly linked-pentacene dimer (PcD) by a single bond, the cyclic voltammograms in PcD-Biph and PcD-3Ph…

This Feature Article presents an overview of the current status of ring polymer molecular dynamics (RPMD) rate theory. We first analyze the RPMD approach and its connection to quantum transition-state theory. We then focus on its practical applications to prototypical chemical reactions in the gas phase, which demonstrate how accurate and reliable RPMD is for calculating thermal chemical reaction rate coefficients in multifarious cases. This review serves as an important checkpoint in RPMD rate theory development, which shows that RPMD is shifting from being just one of recent novel ideas to a well-established and validated alternative to conventional techniques for calculating thermal chemical rate coefficients. We also hope it will motivate further applications of RPMD to various chemical reactions.

In this paper, stainless steel meshes with superhydrophobic and superoleophilic surfaces were fabricated by rapid and simple one-step immersion in a solution containing hydrochloric acid and stearic acid. The apparent contact angles were tested by a video contact angle measurement system (CA). Field emission scanning electron microscopy (FE-SEM), Fourier transform infrared spectroscopy (FTIR), and X-ray photoelectron spectroscopy (XPS) were conducted to characterize the surface topographies and chemical compositions. The SEM results showed that mesh surfaces were covered by ferric stearate (Fe[CH3(CH2)16COO]2) with low surface energy. The CA test results showed that the mesh had a maximum apparent contact angle of 160 ± 1.0° and a sliding angle of less than 5.0° for the water droplet, whereas the apparent contact angle for the oil droplet was zero. Ultrasound oscillation and exposure tests at atmospheric conditions and immersion tests in 3.5 wt % NaCl aqueous solution were conducted to confirm the mesh with excellent superhydrophobic and superoleophilic properties. On the basis of…

COSMO-RS (conductor-like screening model for real solvents) and three different group-contribution methods were used to compute saturation (subcooled) liquid vapor pressures for 16 possible products of ozone-initiated α-pinene autoxidation, with elemental compositions C10H16O4-10 and C20H30O10-12. The saturation vapor pressures predicted by the different methods varied widely. COSMO-RS predicted relatively high saturation vapor pressures values in the range of 10(-6) to 10(-10) bar for the C10H16O4-10 "monomers", and 10(-11) to 10(-16) bar for the C20H30O10-12 "dimers". The group-contribution methods predicted significantly (up to 8 order of magnitude) lower saturation vapor pressures for most of the more highly oxidized monomers. For the dimers, the COSMO-RS predictions were within the (wide) range spanned by the three group-contribution methods. The main reason for the discrepancies between the methods is likely that the group-contribution methods do not contain the necessary parameters to accurately treat autoxidation products containing multiple hydroperoxide, peroxy acid or…

Atmospheric aerosol acidity is an important characteristic of aqueous particles, which has been linked to the formation of secondary organic aerosol by catalyzing reactions of oxidized organic compounds that have partitioned to the particle phase. However, aerosol acidity is difficult to measure and traditionally estimated using indirect methods or assumptions based on composition. Ongoing disagreements between experiments and thermodynamic models of particle acidity necessitate improved fundamental understanding of pH and ion behavior in high ionic strength atmospheric particles. Herein, Raman microspectroscopy was used to determine the pH of individual particles (H2SO4+MgSO4) based on sulfate and bisulfate concentrations determined from νs(SO4(2-)) and νs(HSO4(-)), the acid dissociation constant, and activity coefficients from extended Debye-Hückel calculations. Shifts in pH and peak positions of νs(SO4(2-)) and νs(HSO4(-)) were observed as a function of relative humidity. These results indicate the potential for direct spectroscopic determination of pH in individual particles and…

This work provides new temperature-dependent mole fractions of elusive intermediates relevant to the low-temperature oxidation of dimethyl ether (DME). It extends the previous study of Moshammer et al. [ J. Phys. Chem. A 2015, 119, 7361–7374] in which a combination of a jet-stirred reactor and molecular beam mass spectrometry with single-photon ionization via tunable synchrotron-generated vacuum-ultraviolet radiation was used to identify (but not quantify) several highly oxygenated species. Here, temperature-dependent concentration profiles of 17 components were determined in the range of 450–1000 K and compared to up-to-date kinetic modeling results. Special emphasis is paid toward the validation and application of a theoretical method for predicting photoionization cross sections that are hard to obtain experimentally but essential to turn mass spectral data into mole fraction profiles. The presented approach enabled the quantification of the hydroperoxymethyl formate (HOOCH 2 OCH 2 O), which is a key intermediate in the low-temperature oxidation of DME. The quantification of this…

In this work, two different classes of polyaromatic hydrocarbon (PAH) systems have been investigated in order to characterize the amount of polyradical character and to localize the specific regions of chemical reactivity: (a) the non-Kekulé triangular structures phenalenyl, triangulene and a π-extended triangulene system with high-spin ground state and (b) PAHs based on zethrenes, p-quinodimethane-linked bisphenalenyl, and the Clar goblet containing varying polyradical character in their singlet ground state. The first class of structures already have open-shell character because of their high-spin ground state, which follows from the bonding pattern, whereas for the second class the open-shell character is generated either because of the competition between the closed-shell quinoid Kekulé and the open-shell singlet biradical resonance structures or the topology of the π-electron arrangement of the non-Kekulé form. High-level ab initio calculations based on multireference theory have been carried out to compute singlet-triplet splitting for the above-listed compounds and to provide…

Abstract Organic aerosol particles are known to often absorb/desorb water continuously with change in gas phase relative humidity (RH) without crystallization. Indeed, the prevalence of metastable ultraviscous liquid or amorphous phases in aerosol is well-established with solutes often far exceeding bulk phase solubility limits. Particles are expected to become increasingly viscous with drying, a consequence of the plasticizing effect of water. We report here measurements of the variation in aerosol particle viscosity with RH (equal to condensed phase water activity) for a range of organic solutes including alcohols (diols to hexols), saccharides (mono-, di-, and tri-), and carboxylic acids (di-, tri-, and mixtures). Particle viscosities are measured over a wide range (10–3 to 1010 Pa s) using aerosol optical tweezers, inferring the viscosity from the time scale for a composite particle to relax to a perfect sphere following the coalescence of two particles. Aerosol measurements compare well with bulk phase studies (well-within an order of magnitude deviation at worst) over ranges…

The prediction of electron ionization (EI) mass spectra (MS) from first principles has been a major challenge for quantum chemistry (QC). The unimolecular reaction space grows rapidly with increasing molecular size. On the one hand, statistical models like Eyring's quasi-equilibrium theory and Rice-Ramsperger-Kassel-Marcus theory have provided valuable insight, and some predictions and quantitative results can be obtained from such calculations. On the other hand, molecular dynamics-based methods are able to explore automatically the energetically available regions of phase space and thus yield reaction paths in an unbiased way. We describe in this feature article the status of both methodologies in relation to mass spectrometry for small to medium sized molecules. We further present results obtained with the QCEIMS program developed in our laboratory. Our method, which incorporates stochastic and dynamic elements, has been a significant step toward the reliable routine calculation of EI mass spectra.

In this study, a novel immobilized TiO 2 /Ti film with exposed {001} facets was prepared via a facile one-pot hydrothermal route for the degradation of norfloxacin from aqueous media. The effects of various hydrothermal conditions (i.e., solution pH, hydrothermal time ( H T ) and HF concentration) on the growth of {001} faceted TiO 2 /Ti film were investigated. The maximum photocatalytic performance of {001} faceted TiO 2 /Ti film was observed when prepared at pH 2.62, H T of 3 h and at HF concentration of 0.02 M. The as-prepared {001} faceted TiO 2 /Ti films were fully characterized by field-emission scanning electron microscope (FE-SEM), X-ray diffraction (XRD), high resolution transmission electron microscope (HR-TEM), and X-ray photoelectron spectroscopy (XPS). More importantly, the as-prepared {001} faceted TiO 2 /Ti film exhibited excellent photocatalytic performance toward degradation of norfloxacin in various water matrices (Milli-Q water, tap water, river water and synthetic wastewater). The individual influence of various anions (SO 4 2–, HCO 3 –, NO 3 –, Cl – ) and…

A theoretical study of the cooperativity in linear chains of (H3SiCN)n and (H3SiNC)n complexes connected by tetrel bonds has been carried out by means of MP2 and CCSD(T) computational methods. In all cases, a favorable cooperativity is observed, especially in some of the largest linear chains of (H3SiNC)n, where the effect is so large that the SiH3 group is almost equidistant to the two surrounding CN groups and it becomes planar. In addition, the combination of tetrel bonds with other weak interactions (halogen, chalcogen, pnicogen, triel, beryllium, lithium, and hydrogen bond) has been explored using ternary complexes, (H3SiCN)2:XY and (H3SiNC)2:XY. In all cases, positive cooperativity is obtained, especially in the (H3SiNC)2:ClF and (H3SiNC)2:SHF ternary complexes, where, respectively, halogen and chalcogen shared complexes are formed.

The Criegee intermediate acetone oxide, (CH 3 ) 2 COO, is formed by laser photolysis of 2,2-diiodopropane in the presence of O 2 and characterized by synchrotron photoionization mass spectrometry and by cavity ring-down ultraviolet absorption spectroscopy. The rate coefficient of the reaction of the Criegee intermediate with SO 2 was measured using photoionization mass spectrometry and pseudo-first-order methods to be (7.3 ± 0.5) × 10 –11 cm 3 s –1 at 298 K and 4 Torr and (1.5 ± 0.5) × 10 –10 cm 3 s –1 at 298 K and 10 Torr (He buffer). These values are similar to directly measured rate coefficients of anti -CH 3 CHOO with SO 2, and in good agreement with recent UV absorption measurements. The measurement of this reaction at 293 K and slightly higher pressures (between 10 and 100 Torr) in N 2 from cavity ring-down decay of the ultraviolet absorption of (CH 3 ) 2 COO yielded even larger rate coefficients, in the range (1.84 ± 0.12) × 10 –10 to (2.29 ± 0.08) × 10 –10 cm 3 s –1 . Photoionization mass spectrometry measurements with deuterated acetone oxide at 4 Torr show an inverse…

Six free base porphyrins bearing 0–4 meso substituents have been examined by steady-state and time-resolved absorption and fluorescence spectroscopy in both toluene and N, N -dimethylformamide (DMF). The lifetime of the lowest singlet excited state (S 1 ) decreases with an increase in the number of meso -phenyl groups; the values in toluene are H 2 P-0 (15.5 ns) > H 2 P-1 (14.9 ns) > H 2 P-2c (14.4 ns) > H 2 P-2t (13.8 ns) ∼ H 2 P-3 (13.8 ns) > H 2 P-4 (12.8 ns), where “H 2 P” refers to the core free base porphyrin, the numerical suffix indicates the number of meso -phenyl groups, and “c” and “t” refer to cis and trans, respectively. The opposite trend is found for the fluorescence quantum yield; the values in toluene are H 2 P-0 (0.049) < H 2 P-1 (0.063) ∼ H 2 P-2c (0.063) < H 2 P-2t (0.071) < H 2 P-3 (0.073) < H 2 P-4 (0.090). Similar trends occur in DMF. All radiative and nonradiative (internal conversion and intersystem crossing) rate constants for S 1 decay increase with the increasing number of meso -phenyl groups. The increase in the rate constant for fluorescence parallels…

We investigate the utilization of the domain local pair natural orbital coupled cluster (DLPNO-CCSD(T)) method for calculating binding energies of atmospherical molecular clusters. Applied to small complexes of atmospherical relevance we find that the DLPNO method significantly reduces the scatter in the binding energy, which is commonly present in DFT calculations. For medium sized clusters consisting of sulfuric acid and bases the DLPNO method yields a systematic underestimation of the binding energy compared to canonical coupled cluster results. The errors in the DFT binding energies appear to be more random, while the systematic nature of the DLPNO results allows the establishment of a scaling factor, to better mimic the canonical coupled cluster calculations. Based on the trends identified for the small and medium sized systems, we further extend the application of the DLPNO method to large acid - base clusters consisting of up to 10 molecules, which have previously been out of reach with accurate coupled cluster methods. Using the Atmospheric Cluster Dynamics Code (ACDC) we…

Using a comparative evaporation kinetics approach, we describe a new and accurate method for determining the equilibrium hygroscopic growth of aerosol droplets. The time-evolving size of an aqueous droplet, as it evaporates to a steady size and composition that is in equilibrium with the gas phase relative humidity, is used to determine the time-dependent mass flux of water, yielding information on the vapor pressure of water above the droplet surface at every instant in time. Accurate characterization of the gas phase relative humidity is provided from a control measurement of the evaporation profile of a droplet of know equilibrium properties, either a pure water droplet or a sodium chloride droplet. In combination, and by comparison with simulations that account for both the heat and mass transport governing the droplet evaporation kinetics, these measurements allow accurate retrieval of the equilibrium properties of the solution droplet (i.e., the variations with water activity in the mass fraction of solute, diameter growth factor, osmotic coefficient or number of water…

We present an approach for parametrizing spectroscopic maps of carbonyl groups against experimental IR absorption spectra. The model correlates electric fields sampled from molecular dynamics simulations with vibrational frequencies and line shapes in different solvents. We perform an exhaustive search of parameter combinations and optimize the parameter values for the ester carbonyl stretching mode in ethyl acetate by comparing to experimental FTIR spectra of the small molecule in eight different solvents of varying polarities. Hydrogen-bonding solvents require that the peaks are fit independently for each hydrogen bond ensemble to compensate for improper sampling in molecular dynamics simulations. Spectra simulated using the optimized electrostatic map reproduce C═O IR absorption spectra of ethyl acetate with a line center RMSD error of 4.9 cm(-1) over 12 different solvents whose measured line centers span a 45 cm(-1) range. In combination with molecular dynamics simulations, this spectroscopic map will be useful in interpreting spectra of ester groups in heterogeneous…

Atomically deposited layers of SiO2 and Al2O3 have been recognized as promising coating materials to buffer the volumetric expansion and capacity retention upon the chemo-mechanical cycling of the nanostructured silicon- (Si-) based electrodes. Furthermore, silica (SiO2) is known as a promising candidate for the anode of next-generation lithium ion batteries (LIBs) due to its superior specific charge capacity and low discharge potential similar to Si anodes. In order to describe Li-transport in mixed silica/alumina/silicon systems we developed a ReaxFF potential for Li-Si-O-Al interactions. Using this potential, a series of hybrid grand canonical Monte Carlo (GCMC) and molecular dynamic (MD) simulations were carried out to probe the lithiation behavior of silica structures. The Li transport through both crystalline and amorphous silica was evaluated using the newly optimized force field. The anisotropic diffusivity of Li in crystalline silica cases is demonstrated. The ReaxFF diffusion study also verifies the transferability of the new force field from crystalline to amorphous…

One of the most significant aspects in the development of dye-sensitized solar cells is the exploration and design of high-efficiency and low-cost dyes. This paper reports the theoretical design of various triphenylamine analogues, wherein the central nitrogen moiety establishes an sp(2)-hybridization, which endows a significant participation in the charge-transfer properties. Density functional theory (DFT) and time-dependent DFT methodologies were utilized to investigate the geometry, electronic structure, photochemical properties, and electrochemical properties of these dyes. Different exchange-correlation functionals were initially evaluated to establish a proper methodology for calculating the excited-state energy of the reference dye, known as DIA3. Consequently, TD-LC-ωPBE with a damping parameter of 0.175 Bohr(-1) best correlates with the experimental value. Four new dyes, namely, Dhk1, Dhk2, Dhk3, and Dhk4, were designed by modifying the rigidity of the donor moiety. According to the results, altering the type and position of binding in the donor group leads to distinct…

We investigate the molecular interaction between methyl-substituted N,N,N',N'-ethylenediamines, propane-1,3-diamine, butane-1,4-diamine, and sulfuric acid using computational methods. Molecular structure of the diamines and their dimer clusters with sulfuric acid is studied using three density functional theory methods (PW91, M06-2X, and ωB97X-D) with the 6-31++G(d,p) basis set. A high level explicitly correlated CCSD(T)-F12a/VDZ-F12 method is used to obtain accurate binding energies. The reaction Gibbs free energies are evaluated and compared with values for reactions involving ammonia and atmospherically relevant monoamines (methylamine, dimethylamine, and trimethylamine). We find that the complex formation between sulfuric acid and the studied diamines provides similar or more favorable reaction free energies than dimethylamine. Diamines that contain one or more secondary amino groups are found to stabilize sulfuric acid complexes more efficiently. Elongating the carbon backbone from ethylenediamine to propane-1,3-diamine or butane-1,4-diamine further stabilizes the complex…

Abstract A new local property, the local electron attachment energy [E(r)], is introduced and is demonstrated to be a useful guide to predict intermolecular interactions and chemical reactivity. The E(r) is analogous to the average local ionization energy but indicates susceptibility toward interactions with nucleophiles rather than electrophiles. The functional form E(r) is motivated based on Janak’s theorem and the piecewise linear energy dependence of electron addition to atomic and molecular systems. Within the generalized Kohn–Sham method (GKS-DFT), only the virtual orbitals with negative eigenvalues contribute to E(r). In the present study, E(r) has been computed from orbitals obtained from GKS-DFT computations with a hybrid exchange-correlation functional. It is shown that E(r) computed on a molecular isodensity surface, ES(r), reflects the regioselectivity and relative reactivity for nucleophilic aromatic substitution, nucleophilic addition to activated double bonds, and formation of halogen bonds. Good to excellent correlations between experimental or theoretical measures…

High-valent metal-oxo oxidants are common reactive species in synthetic catalysts as well as heme and nonheme iron enzymes. In general, they efficiently react with substrates through oxygen atom transfer, and for a number of cases, experimental rate constants have been determined. However, because these rate constants are generally measured in a polar solution, it has been found difficult to find computational methodologies to reproduce experimental trends and reactivities. In this work, we present a detailed computational study into para-substituted thioanisole sulfoxidation by a nonheme iron(IV)-oxo complex. A range of density functional theory methods and basis sets has been tested for their suitability to describe the reaction mechanism and compared with experimentally obtained free energies of activation. It is found that the enthalpy of activation is reproduced well, but all methods overestimate the entropy of activation by about 50%, for which we recommend a correction factor. The effect of solvent and dispersion on the barrier heights is explored both at the single-point…

The structures and electronic properties of low-energy neutral and anionic Mg n ( n = 3–20) clusters have been studied by utilizing a widely adopted CALYPSO structure searching method coupled with density functional theory calculations. A large number of low-energy isomers are optimized at the B3PW91 functional with the 6-311+G(d) basis set. The optimized geometries clearly indicate that a structural transition from hollow three-dimensional configurations to filled-cage-like structures occurs at n = 16 for both neutral and anionic clusters. Based on the anionic ground state structures, photoelectron spectra are simulated using time-dependent density functional theory (TD-DFT) and compared with experimental results. The good agreement validates that the current ground state structures, obtained from the symmetry-unconstrained searches, are true global minima. A detailed chemical bonding analysis distinctly indicates that the Mg 17 cluster is the first neutral locally π-aromatic homonuclear all-metal cluster, which perfectly satisfies Hückel’s well-known 4 N + 2 rule.

Four mechanisms previously proposed for dioxygen activation catalyzed by α-keto acid dependent oxygenases (α-KAO) were studied with dispersion-corrected DFT methods employing B3LYP and TPSSh functionals in combination with triple-ζ basis set (cc-pVTZ). The aim of this study was to revisit mechanisms suggested in the past decade and resolve remaining issues related to dioxygen activation. Mechanism A, which runs on the quintet potential energy surface (PES) and includes formation of an Fe(III)-superoxide radical anion complex, subsequent oxidative decarboxylation, and O-O bond cleavage, was found to be most likely. However, mechanism B taking place on the septet PES involves a rate limiting barrier comparable to the one found for mechanism A, and thus it cannot be excluded, though two other mechanisms (C and D) were ruled out. Mechanism C is a minor variation of mechanism A, whereas mechanism D proceeds through formation of a triplet Fe(IV)-alkyl peroxo bridged intermediate. The study covered also full optimization of relevant minimum energy crossing points (MECPs). The relative…

We developed the ReaxFF force field for Pt/Ni/C/H/O interactions, specifically targeted for heterogeneous catalysis application of the Pt–Ni alloy. The force field is trained using the DFT data for equations of state of Pt 3 Ni, PtNi 3 and PtNi alloys, the surface energy of the Pt x Ni 1– x (111) ( x = 0.67–0.83), and binding energies of various atomic and molecular species (O, H, C, CH, CH 2, CH 3, CO, OH, and H 2 O) on these surfaces. The ReaxFF force field shows a Pt surface segregation at x ≥ 0.67 for the (111) surface and x ≥ 0.62 for the (100) surface in vacuum. In addition, from the investigation of the preferential alloy component of the adsorbates, it is expected that H and CH 3 on the alloy surface to induce a segregation of Pt whereas the oxidation intermediates and products such as C, O, OH, H 2 O, CO, CH, and CH 2 are found to induce Ni segregation. The relative order of binding strengths among adsorbates is a function of alloy composition and the force field is trained to describe the trend observed in DFT calculations, namely, H 2 < H 2 O < CH 3 ≈ O 2 ≈ CO < OH < CH 2…

A new allotrope of nitrogen in which the atoms are connected to form a novel N6 molecule is predicted to exist at ambient conditions. The N6 molecule is a charge-transfer complex with an open-chain structure containing both single and triple bonds. The charge transfer induces ionic characteristics in the intermolecular interactions and leads to a much higher cohesive energy for the predicted crystal compared to solid N2. The N6 solid is also more stable than a previously reported polymeric solid of nitrogen. Because of the kinetic stability of the molecules and strong intermolecular interactions, the N6 crystal is shown by metadynamics simulations to be dynamically stable around room temperature and to only dissociate to N2 molecules above 700 K. The N6 crystal can likely be synthesized under high-pressure high-temperature conditions, and the considerable metastability may allow for an ambient-pressure recovery of the crystal. Because of the large energy difference between the single and triple bonds, the dissociation of the N6 crystal is expected to release a large amount of…

The formation of new particles through condensation from the gas phase is an important source of atmospheric aerosols. The properties of the electrically neutral clusters formed in the very first steps of the condensation process are, however, not directly observable by experimental means. We present here electronic structure calculations on the hydrates of clusters of three molecules of sulfuric acid and three molecules of ammonia or dimethylamine. On the basis of the results of these new calculations together with previously published material we simulate the influence of hydration on the dynamic processes involved in particle formation. Most strongly affected by hydration and most important as a mediator for the effect on particle formation rates are the evaporation rates of clusters. The results give an estimate of the sensitivity of the atmospheric particle formation rate for humidity. The particle formation rate can change approximately two orders of magnitude in either direction due to hydration; the net effect, however, is highly dependent on the exact conditions.

The unimolecular decomposition of (CH3)2COO and (CD3)2COO was measured by direct detection of the Criegee intermediate at temperatures from 283 to 323 K using time-resolved UV absorption spectroscopy. The unimolecular rate coefficient kd for (CH3)2COO shows a strong temperature dependence, increasing from 269 ± 82 s(-1) at 283 K to 916 ± 56 s(-1) at 323 K with an Arrhenius activation energy of ∼6 kcal mol(-1). The bimolecular rate coefficient for the reaction of (CH3)2COO with SO2, kSO2, was also determined in the temperature range 283 to 303 K. Our temperature-dependent values for kd and kSO2 are consistent with previously reported relative rate coefficients kd/kSO2 of (CH3)2COO formed from ozonolysis of tetramethyl ethylene. Quantum chemical calculations of kd for (CH3)2COO are consistent with the experiment, and the combination of experiment and theory for (CD3)2COO indicates that tunneling plays a significant role in (CH3)2COO unimolecular decomposition. The fast rates of unimolecular decomposition for (CH3)2COO measured here, in light of the relatively slow rate for the…

We combine quantum chemical and molecular docking techniques to provide new insights into how piperine molecule in various forms of pepper enhances bioavailability of a number of drugs including curcumin in turmeric for which it increases its bioavailability by a 20-fold. We have carried out docking studies of quantum chemically optimized piperine structure binding to curcumin, CYP3A4 in cytochrome P450, p-Glycoprotein and UDP-glucuronosyltransferase (UGT), the enzyme responsible for glucuronosylation, which increases the solubility of curcumin. All of these studies establish that piperine binds to multiple sites on the enzymes and also intercalates with curcumin forming a hydrogen bonded complex with curcumin. The conjugated network of double bonds and the presence of multiple charge centers of piperine offer optimal binding sites for piperine to bind to enzymes such as UDP-GDH, UGT, and CYP3A4. Piperine competes for curcumin's intermolecular hydrogen bonding and its stacking propensity by hydrogen bonding with enolic proton of curcumin. This facilitates its metabolic transport,…

A wide variety of organic dyes form, under certain conditions, clusters know as J- and H-aggregates. Cyanine dyes are such a class of molecules where the spatial proximity of several dyes leads to overlapping electron orbitals and thus to the creation of a new energy landscape compared to that of the individual units. In this work, we create artificial H-aggregates of exactly two Cyanine 3 (Cy3) dyes by covalently linking them to a DNA molecule with controlled subnanometer distances. The absorption spectra of these coupled systems exhibit a blue-shifted peak, whose intensity varies depending on the distance between the dyes and the rigidity of the DNA template. Simulated vibrational resolved spectra, based on molecular orbital theory, excellently reproduce the experimentally observed features. Circular dichroism spectroscopy additionally reveals distinct signals, which indicates a chiral arrangement of the dye molecules. Molecular dynamic simulations of a Cy3-Cy3 construct including a 14-base pair DNA sequence verified chiral stacking of the dye molecules.

Atmospheric aerosols are usually complex mixtures of inorganic and organic compounds. The hygroscopicity of mixed particles is closely related to their chemical composition and interactions between components, which is still poorly understood. In this study, the hygroscopic properties of submicron particles composed of NaCl and dicarboxylic acids including oxalic acid (OA), malonic acid (MA), and succinic acid (SA) with various mass ratios are investigated with a hygroscopicity tandem differential mobility analyzer (HTDMA) system. Both the Zdanovskii-Stokes-Robinson (ZSR) method and extended aerosol inorganics model (E-AIM) are applied to predict the water uptake behaviors of sodium chloride/dicarboxylic acid mixtures. For NaCl/OA mixed particles, the measured growth factors were significantly lower than predictions from the model methods, indicating a change in particle composition caused by chloride depletion. The hygroscopic growth of NaCl/MA particles was well described by E-AIM, and that of NaCl/SA particles was dependent upon mixing ratio. Compared with model predictions, it…

Wide bandgap hybrid halide perovskites based on bromine and chlorine halide anions have emerged as potential candidates for various optoelectronic devices. However, these materials are relatively less explored than the iodine-based perovskites for microscopic details. We present experiment and first-principles calculations to understand the structural, optical, and electronic structure of wide bandgap CH3NH3Pb(Br1-xClx)3 (x = 0, 0.33, 0.66, and 1) 3D hybrid perovskite materials. We substituted Br(-) with Cl(-) to tune the bandgap from 2.4 eV (green emissive) to 3.2 eV (blue (UV) emissive) of these materials. We correlate our experimental results with first-principles theory and provide an insight into important parameters like lattice constants, electronic structure, excitonic binding energy (EX), dielectric constant, and reduced effective mass (μr) of charge carriers in these perovskite semiconductors. Electronic structure calculations reveal that electronic properties are mainly governed by Pb 6p and halide p orbitals. Our estimates of EX within a hydrogen model suggest that an…

Quantum chemical theory (DFT and MP2) and vibrational spectroscopy (ATR-IR and Raman) were employed to investigate the electronic structure and molecular interactions in the room-temperature ionic liquid 1-ethyl-3-methylimidazolium trifluoromethanesulfonate. Various possible conformers of a cation-anion pair based on their molecular interactions were simulated in the gas phase. All the different theoretical (MP2, B3LYP, and the dispersion-corrected wB97XD) methods assume the same ion-pair conformation for the lowest energy state. Basis set superimpose error (BSSE) correction was also introduced by using the counterpoise method. Strong C-H···O interactions between the most acidic hydrogen atom of the cation imidazole ring (C2H) and the oxygen atom of the anion were predicted where the anion is located at the top of (C2H). In this case, methyl and alkyl groups also interact with the anion in the form of a C-H···O hydrogen bond. Interestingly, the dispersion-corrected methodology neglects the C4/C5-H···O and C-H···F interaction in the ion-pair calculations. The theoretical results were…

Gaseous biogenic volatile organic compounds (BVOCs) are immediately oxidized by gaseous oxidants to form BVOC-acids that rapidly condense onto aqueous aerosol phase and thus contribute to the growth of atmospheric particles. Because BVOC-acids are highly hydrophobic and hence surface-active in nature, it seems critical to study the oxidation by gaseous hydroxyl radical (·OH(g)) at the air-water interface. Here we report on the fast (≤10 μs) oxidation of aqueous cis-pinonic acid (C10H16O3, CPA, cis-pinonate anion's m/z = 183), a representative BVOC-acid, by ·OH(g) at the air-water interface for the first time. We find that cis-pinonate anion is more enriched at the air-water interface by ∼4 and ∼14 times than n-octanoate anion at 10 and 100 μM, respectively, as revealed by an interface-specific mass spectrometry of the equimolar mixture of microjets. Exposure of aqueous CPA microjets to ·OH(g) pulses from the 266 nm laser photolysis of O3(g)/O2(g)/H2O(g)/N2(g) mixtures yields pinonic peroxyl radicals (m/z = 214) that lead to the functionalization products carbonyls (m/z = 197),…

Lignocellulosic biomass is a promising sustainable feedstock for the production of biofuels, biomaterials, and biospecialty chemicals. However, efficient utilization of biomass has been limited by our poor understanding of its molecular structure. Here, we report a dynamic nuclear polarization (DNP)-enhanced solid-state (SS)NMR study of the molecular structure of biomass, both pre- and postcatalytic treatment. This technique enables the measurement of 2D homonuclear 13 C– 13 C correlation SSNMR spectra under natural abundance, yielding, for the first time, an atomic-level picture of the structure of raw and catalytically treated biomass samples. We foresee that further such experiments could be used to determine structure–function relationships and facilitate the development of more efficient, and chemically targeted, biomass-conversion technologies.

Abstract The structures of gas-phase M+(CO2)n (M = Co, Rh, Ir; n = 2–15) ion–molecule complexes have been investigated using a combination of infrared resonance-enhanced photodissociation (IR-REPD) spectroscopy and density functional theory. The results provide insight into fundamental metal ion–CO2 interactions, highlighting the trends with increasing ligand number and with different group 9 ions. Spectra have been recorded in the region of the CO2 asymmetric stretch around 2350 cm–1 using the inert messenger technique and their interpretation has been aided by comparison with simulated infrared spectra of calculated low-energy isomeric structures. All vibrational bands in the smaller complexes are blue-shifted relative to the asymmetric stretch in free CO2, consistent with direct binding to the metal center dominated by charge−quadrupole interactions. For all three metal ions, a core [M+(CO2)2] structure is identified to which subsequent ligands are less strongly bound. No evidence is observed in this size regime for complete activation or insertion reactions.

Broadband optical pump and compressed white light continuum probe were used to measure the transient excited-state absorption, ground-state bleach, and stimulated emission signals of cresyl violet solution in methanol. Amplitude oscillations caused by wavepacket motion in the ground and excited electronic states were analyzed. It was found that vibrational coherences in the excited state persist for more than the experimental waiting time window of 6 ps, and the strongest mode had a dephasing time constant of 2.4 ps. We hypothesize the dephasing of the wavepacket in the excited state is predominantly caused by intramolecular vibrational relaxation (IVR). Slow IVR indicates weak mode-mode coupling and therefore weak anharmonicity of the potential of this vibration. Thus, the initially prepared vibrational wavepacket in the excited state is not significantly perturbed by nonadiabatic coupling to other electronic states, and hence the diabatic and adiabatic representations of the system are essentially identical within the Born-Oppenheimer approximation. The wavepacket therefore…

Enhanced spin-orbit coupling through external heavy-atom effect (EHE) has been routinely used to induce room-temperature phosphorescence (RTP) for purely organic molecular materials. Therefore, understanding the nature of EHE, i.e., the specific orbital interactions between the external heavy atom and the luminophore, is of essential importance in molecular design. For organic systems, halogens (e.g., Cl, Br, and I) are the most commonly seen heavy atoms serving to realize the EHE-related RTP. In this report, we conduct an investigation on how heavy-atom perturbers and aromatic luminophores interact on the basis of data obtained from crystallography. We synthesized two classes of molecular systems including N-haloalkyl-substituted carbazoles and quinolinium halides, where the luminescent molecules are considered as "base" or "acid" relative to the heavy-atom perturbers, respectively. We propose that electron donation from a π molecular orbital (MO) of the carbazole to the σ* MO of the C-X bond (π/σ*) and n electron donation to a π* MO of the quinolinium moiety (n/π*) are responsible…

Both proper, red-shifting and improper, blue-shifting hydrogen bonds have been well-recognized with enormous experimental and computational studies. The current consensus is that there is no difference in nature between these two kinds of hydrogen bonds, where the electrostatic interaction dominates. Since most if not all the computational studies are based on molecular orbital theory, it would be interesting to gain insight into the hydrogen bonds with modern valence bond (VB) theory. In this work, we performed ab initio VBSCF computations on a series of hydrogen-bonding systems, where the sole hydrogen bond donor CF3H interacts with ten hydrogen bond acceptors Y (═NH2CH3, NH3, NH2Cl, OH(-), H2O, CH3OH, (CH3)2O, F(-), HF, or CH3F). This series includes four red-shifting and six blue-shifting hydrogen bonds. Consistent with existing findings in literature, VB-based energy decomposition analyses show that electrostatic interaction plays the dominating role and polarization plays the secondary role in all these hydrogen-bonding systems, and the charge transfer interaction, which…

By employing the Dirac-Frenkel time-dependent variational principle, we study the dynamical properties of the Holstein molecular crystal model with diagonal and off-diagonal exciton-phonon coupling. A linear combination of the Davydov D1 (D2) ansatz, referred to as the "multi-D1 ansatz" ("multi-D2 ansatz"), is used as the trial state with enhanced accuracy but without sacrificing efficiency. The time evolution of the exciton probability is found to be in perfect agreement with that of the hierarchy equations of motion, demonstrating the promise the multiple Davydov trial states hold as an efficient, robust description of dynamics of complex quantum systems. In addition to the linear absorption spectra computed for both diagonal and off-diagonal cases, for the first time, 2D spectra have been calculated for systems with off-diagonal exciton-phonon coupling by employing the multiple D2 ansatz to compute the nonlinear response function, testifying to the great potential of the multiple D2 ansatz for fast, accurate implementation of multidimensional spectroscopy. It is found that the…

A computational study of the possibilities of contemporary theoretical chemistry as regards calculated thermodynamic properties for molecular crystals from first-principles is presented. The study is performed for a testing set of 22 low-temperature crystalline phases whose properties such as densities of phonon states, isobaric heat capacities, and densities are computed as functions of temperature within the quasi-harmonic approximation. Electronic structure and lattice dynamics are treated by plane-wave based calculations with optPBE-vdW functional. Comparison of calculated results with reliable critically assessed experimental data is especially emphasized.

This study investigates the possibility of detecting relativistic effects and electron correlation in single-crystal X-ray diffraction experiments using the examples of diphenyl mercury (HgPh2) and triphenyl bismuth (BiPh3). In detail, the importance of electron correlation (ECORR), relativistic effects (REL) [distinguishing between total, scalar and spin-orbit (SO) coupling relativistic effects] and picture change error (PCE) on the theoretical electron density, its topology and its Laplacian using infinite order two component (IOTC) wave functions is discussed. This is to develop an understanding of the order of magnitude and shape of these different effects as they manifest in the electron density. Subsequently, the same effects are considered for the theoretical structure factors. It becomes clear that SO and PCE are negligible, but ECORR and scalar REL are important in low- and medium-order reflections on absolute and relative scales-not in the high-order region. As a further step, Hirshfeld atom refinement (HAR) and subsequent X-ray constrained wavefunction (XCW) fitting have…

We carried out reactive molecular dynamics simulations by ReaxFF to study the initial events of an insensitive high explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) against various thermal stimuli including constant-temperature heating, programmed heating, and adiabatic heating to simulate TATB suffering from accidental heating in reality. Cluster evolution at the early stage of the thermal decomposition of condensed TATB was the main focus as cluster formation primarily occurs when TATB is heated. The results show that cluster formation is the balance of the competition of intermolecular collision and molecular decomposition of TATB, that is, an appropriate temperature and certain duration are required for cluster formation and preservation. The temperature in the range of 2000-3000 K was found to be optimum for fast formation and a period of preservation. Besides, the intra- and intermolecular H transfers are always favorable, whereas the C-NO2 partition was favorable at high temperature. The simulation results are helpful to deepen the insight into the thermal properties of…

The photophysics of a norbornyl-bridged covalent tetracene (Tc) dimer BT1 and a monomer analogue Tc-e were studied in room-temperature nonpolar solvents. Notably in BT1, a Davydov-split band is observed in UV absorption, heralding interchromophore electronic interactions. Emission spectra indicate an acene-like vibronic progression mirroring the lowest-energy visible absorption. For BT1, this argues against excited-state excimer formation. Evidence of intramolecular singlet fission (SF) comes from a comparison of time-resolved emission decay signals collected for BT1 versus Tc-e in toluene. In BT1, the multiexcitonic (1)TT state is produced in 70 ns in 6% yield. A ratio of fission versus fusion rate constants provides an experimental measure of the SF reaction free energy at 52 meV in good agreement with previous calculations. The low SF yield corroborates our expectations that orbital symmetry effects on diabatic coupling for SF are important for dimers that cannot rely on more favorable thermodynamics.