Detection of a C4 Criegee Intermediate: Fourier-Transform Microwave Spectroscopy of Methacrolein Oxide.

Pure rotational transitions of methacrolein oxide (MACRO) were observed by Fourier-transform microwave spectroscopy. Among the four low-lying conformers existing within an energy window of 3 kcal/mol, only the lowest-energy conformer, the anti-trans conformer, was detected in a discharged jet of a 1,3-diiode-2-methylprop-1-ene and O2 mixture diluted in Ar. Nineteen pure rotational transitions, in the frequency range from 10 to 25 GHz, most of them showing A/E splitting due to the methyl-top internal rotation, were observed and analyzed by the XIAM program, yielding the internal rotation barrier of 559 cm-1, which very well agrees with a theoretically calculated value, 558 cm-1, at the CCSD(T)/cc-pVTZ level of theory.

Mechanism of proton-coupled electron transfer described with QM/MM implementation of coupled-perturbed density-functional tight-binding.

Coupled-perturbed equations for degenerate orbitals were implemented for third order density-functional tight binding, which allowed the use of Mulliken charges as reaction coordinates. The method was applied to proton-coupled electron transfer (PCET) reactions in a model system and thoroughly tested for QM and QM/MM setups (i.e., coupled quantum and molecular mechanics). The performed enhanced sampling simulations were stable, and the obtained potentials of the mean force were able to address the thermodynamic and kinetic features of the reactions by showing the expected topography and energy barriers. Hence, this method has the potential to distinguish between concerted and sequential mechanisms and could next be applied to proton-coupled electron transfer reactions in more complex systems like proteins.

Scanning pH-metry for Observing Reversibility in Protein Folding.

One of the main factors affecting protein structure in solution is pH. Traditionally, to study pH-dependent conformational changes in proteins, the concentration of the H+ ions is adjusted manually, complicating real-time analyses, hampering dynamic pH regulation, and consequently leading to a limited number of tested pH levels. Here, we present a programmable device, a scanning pH-meter, that can automatically generate different types of pH ramps and waveforms in a solution. A feedback loop algorithm calculates the required flow rates of the acid/base titrants, allowing one, for example, to generate periodic pH sine waveforms to study the reversibility of protein folding by fluorescence spectroscopy. Interestingly, for some proteins, the fluorescence intensity profiles recorded in such a periodically oscillating pH environment display hysteretic behavior indicating an asymmetry in the sequence of the protein unfolding/refolding events, which can most likely be attributed to their distinct kinetics. Another useful application of the scanning pH-meter concerns coupling it with an electrospray ionization mass spectrometer to observe pH-induced structural changes in proteins as revealed by their varying charge-state distributions. We anticipate a broad range of applications of the scanning pH-meter developed here, including protein folding studies, determination of the optimum pH for achieving maximum fluorescence intensity, and characterization of fluorescent dyes and other synthetic materials.

Density-Functional Tight-Binding Parameters for Bulk Zirconium: A Case Study for Repulsive Potentials.

Density-functional tight-binding (DFTB) parameters are presented for the simulation of the bulk phases of zirconium. Electronic parameters were obtained using a band structure fitting strategy, while two-center repulsive potentials were created by particle swarm optimization. As objective functions for the repulsive potential fitting, we employed the Birch-Murnaghan equations of state for hexagonal close-packed (HCP), body-centered cubic (BCC) and ω phases of Zr from density-functional theory (DFT). When fractional atomic coordinates are not allowed to change in the generation of the equation-of-state curves, long-range repulsive DFTB potentials are able to almost perfectly reproduce equilibrium structures, relative DFT energies of the bulk phases, and bulk moduli. However, the same potentials lead to artifacts in the DFTB potential energy surfaces when atom positions in the unit cell are allowed to fully relax during the change of unit cell parameters. Conventional short-range repulsive DFTB potentials, while inferior in their ability to reproduce DFT bulk energetics, are able to correctly reproduce the qualitative shape of the DFT potential energy surfaces, including the location of global minima, and can therefore be considered more transferable.

Structures of Pyridine-Water Clusters Studied with Infrared-Vacuum Ultraviolet Spectroscopy.

The infrared (IR) spectra of the O-H stretching vibrations of pyridine-water clusters (Pyd)m(H2O)n, with m, n = 1-4, have been investigated with infrared-vacuum ultraviolet (VUV) spectroscopy under a jet-cooled condition. The time-of-flight mass spectrum of (Pyd)m(H2O)n+ by VUV ionization at ∼9 eV showed an unusual intensity pattern with very weak ion signals for m = 1 and 2 and stronger signals for m ≥ 3. This unusual mass pattern was explained by a drastic structural change of (Pyd)m(H2O)n upon the VUV ionization, which was followed by the elimination of water molecules. Among the recorded IR spectra, only one spectrum monitored, (Pyd)2+ cation, showed a well-resolved structure. The spectrum was analyzed by comparing with the simulated ones of possible stable isomers of (Pyd)2(H2O)n, which were obtained with quantum-chemical calculations. Most of the calculated (Pyd)2(H2O)n clusters had the characteristic structure in which H2O or (H2O)2 forms a hydrogen-bonded bridge between two pyridines to form the π-stacked (Pyd)2, and an additional H2O molecule(s) extends the H-bonded network. The π-stacked (Pyd)2(H2O)n moiety is very stable and is thought to exist as a local structure in a pyridine/water mixed solution. The Fermi resonance between the O-H stretch fundamentals and the overtones of the O-H bending vibrations in (Pyd)m(H2O)n was found to be less pronounced in the case of (Pyd)m(NH3)n studied previously.

Zhang-Zhang Polynomials of Multiple Zigzag Chains Revisited: A Connection with the John-Sachs Theorem.

Multiple zigzag chains Zm,n of length n and width m constitute an important class of regular graphene flakes of rectangular shape. The physical and chemical properties of these basic pericondensed benzenoids can be related to their various topological invariants, conveniently encoded as the coefficients of a combinatorial polynomial, usually referred to as the ZZ polynomial of multiple zigzag chains Zm,n. The current study reports a novel method for determination of these ZZ polynomials based on a hypothesized extension to John-Sachs theorem, used previously to enumerate Kekulé structures of various benzenoid hydrocarbons. We show that the ZZ polynomial of the Zm,n multiple zigzag chain can be conveniently expressed as a determinant of a Toeplitz (or almost Toeplitz) matrix of size m2×m2 consisting of simple hypergeometric polynomials. The presented analysis can be extended to generalized multiple zigzag chains Zkm,n, i.e., derivatives of Zm,n with a single attached polyacene chain of length k. All presented formulas are accompanied by formal proofs. The developed theoretical machinery is applied for predicting aromaticity distribution patterns in large and infinite multiple zigzag chains Zm,n and for computing the distribution of spin densities in biradical states of finite multiple zigzag chains Zm,n.

Sample Flow Rate Scan in Electrospray Ionization Mass Spectrometry Reveals Alterations in Protein Charge State Distribution.

Sample flow rate is one of the parameters that influence the sensitivity of electrospray ionization (ESI) mass spectrometry. By varying the sample flow rate, initial droplets of different sizes can be generated. Protein molecules in small droplets may form gas-phase ions earlier than the ones in large droplets. Here, we have systematically studied the influence of sample flow rate on the ESI charge state distributions (CSDs) of model proteins. A dedicated programmable sample flow rate scanner was used to infuse protein samples at different flow rates into a mass spectrometer. The synergistic influence of sample flow rate and various electrolytes (ammonium acetate, ammonium bicarbonate, ammonium formate, and piperidine) was studied. Significant alterations to the CSDs with increasing flow rate were observed. For example, in the presence of ammonium acetate, at low flow rates, lower charge states of proteins showed high intensities, while at high flow rates, ions related to higher charge states of proteins dominated the spectra. On the other hand, in the presence of piperidine, a significant reduction in the ion currents of all charge states was observed during the flow rate scan. Our observations suggest that at low flow rates the protein molecules follow a charged residue model of ionization mechanism, and at high flow rates-due to structural changes in protein molecules in large ESI droplets-the charged residue and chain ejection models can possibly coexist. We propose the use of sample flow rate scan as a way to reveal the influence of flow rate on the CSDs of the studied proteins.

Detection of a Criegee Intermediate with an Unsaturated Hydrocarbon Substituent: Fourier-Transform Microwave Spectroscopy of Methyl Vinyl Ketone Oxide.

Pure rotational transitions of methyl vinyl ketone oxide, or the so called methyl-vinyl Criegee intermediate, were observed by Fourier-transform microwave spectroscopy. Among the four possible conformers of this species, predicted by theory within an energy window of 3 kcal/mol, only the lowest-energy conformer, the syn-trans form, was detected in a discharged jet of a 1,3-diiode-but-2-ene (either in Z- or E-conformation) and O2 mixture diluted in Ar. Thirty pure rotational transitions with internal rotation splitting of the methyl top were observed. The observed frequencies were analyzed by the XIAM program, yielding an internal rotation barrier of 702.8(28) cm-1, which reasonably agrees with a theoretical value of 680 cm-1 determined with CCSD(T)/cc-pVTZ.

Correcting long-range electrostatics in DFTB.

We demonstrate that the atom-based charge model implemented in the current versions of the density functional tight binding (DFTB) method fails to reproduce the correct charge distribution of a range of systems, including homonuclear molecules, graphene, and nanotubes, resulting in serious distortions in the electrostatic interactions for such systems caused by the missing quadrupole moments. In particular, this failure seriously impacts the long- and medium-range interaction energies of the DFTB plus dispersion (DFTB-D) model, leading to incorrect predictions of translational or rotational barriers in such systems. We show explicitly on examples of H2 and N2 that correct quadrupole moments-and consequently correct electrostatic interactions-can be restored in such systems by adding additional bond (ghost) sites to the homonuclear molecules. Attempts to determine the point charges associated with the additional sites using the usual Mulliken population analysis lead to unphysical results. Instead, these charges can be determined using the actual DFTB densities used in the parameterization process. For homonuclear molecules, we propose an extension to the DFTB-D model by adding charges that reproduce the physically correct quadrupolar charge distribution. The resulting DFTB-D-Q model greatly improves the rotational barriers for interactions of molecular hydrogen and nitrogen with benzene.

Chemical clocks, oscillations, and other temporal effects in analytical chemistry: oddity or viable approach?

Most analytical methods are based on "analogue" inputs from sensors of light, electric potentials, or currents. The signals obtained by such sensors are processed using certain calibration functions to determine concentrations of the target analytes. The signal readouts are normally done after an optimised and fixed time period, during which an assay mixture is incubated. This minireview covers another-and somewhat unusual-analytical strategy, which relies on the measurement of time interval between the occurrences of two distinguishable states in the assay reaction. These states manifest themselves via abrupt changes in the properties of the assay mixture (e.g. change of colour, appearance or disappearance of luminescence, change in pH, variations in optical activity or mechanical properties). In some cases, a correlation between the time of appearance/disappearance of a given property and the analyte concentration can be also observed. An example of an assay based on time measurement is an oscillating reaction, in which the period of oscillations is linked to the concentration of the target analyte. A number of chemo-chronometric assays, relying on the existing (bio)transformations or artificially designed reactions, were disclosed in the past few years. They are very attractive from the fundamental point of view but-so far-only few of them have be validated and used to address real-world problems. Then, can chemo-chronometric assays become a practical tool for chemical analysis? Is there a need for further development of such assays? We are aiming to answer these questions.

Polarizability tensor invariants of H2, HD, and D2.

We report an exhaustive compilation of wavelength-dependent matrix elements over the mean polarizability (α¯) and polarizability anisotropy (γ) operators for the rovibrational states of the H2, HD, and D2 molecules together with an accompanying computer program for their evaluation. The matrix elements can be readily evaluated using the provided codes for rovibrational states with J = 0-15 and v = 0-4 and for any laser wavelengths in the interval 182.25-1320.6 nm corresponding to popular, commercially available lasers. The presented results substantially extend the scope of the data available in the literature, both in respect of the rovibrational transitions analyzed and the range of covered laser frequencies. The presented detailed tabulation of accurate polarizability tensor invariants is essential for successful realization of our main long-term goal: developing a universal standard for determining absolute Raman cross sections and absolute Raman intensities in experimental Rayleigh and Raman scattering studies of molecules.

Infrared absorption spectrum of the simplest deuterated Criegee intermediate CD2OO.

We report a transient infrared (IR) absorption spectrum of the simplest deuterated Criegee intermediate CD2OO recorded using a step-scan Fourier-transform spectrometer coupled with a multipass absorption cell. CD2OO was produced from photolysis of flowing mixtures of CD2I2, N2, and O2 (13 or 87 Torr) with laser light at 308 nm. The recorded spectrum shows close structural similarity with the spectrum of CH2OO reported previously [Y.-T. Su et al., Science 340, 174 (2013)]. The four bands observed at 852, 1017, 1054, and 1318 cm(-1) are assigned to the OO stretching mode, two distinct in-plane OCD bending modes, and the CO stretching mode of CD2OO, respectively, according to vibrational wavenumbers, IR intensities, rotational contours, and deuterium-isotopic shifts predicted with extensive quantum-chemical calculations. The CO-stretching mode of CD2OO at 1318 cm(-1) is blue shifted from the corresponding band of CH2OO at 1286 cm(-1); this can be explained by a mechanism based on mode mixing and isotope substitution. A band near 936 cm(-1), observed only at higher pressure (87 Torr), is tentatively assigned to the CD2 wagging mode of CD2IOO.

Mechanism of back electron transfer in an intermolecular photoinduced electron transfer reaction: solvent as a charge mediator.

Back electron transfer (BET) is one of the important processes that govern the decay of generated ion pairs in intermolecular photoinduced electron transfer reactions. Unfortunately, a detailed mechanism of BET reactions remains largely unknown in spite of their importance for the development of molecular photovoltaic structures. Here, we examine the BET reaction of pyrene (Py) and 1,4-dicyanobenzene (DCB) in acetonitrile (ACN) by using time-resolved near- and mid-IR spectroscopy. The Py dimer radical cation (Py2(·+)) and DCB radical anion (DCB(·-)) generated after photoexcitation of Py show asynchronous decay kinetics. To account for this observation, we propose a reaction mechanism that involves electron transfer from DCB(·-) to the solvent and charge recombination between the resulting ACN dimer anion and Py2(·+). The unique role of ACN as a charge mediator revealed herein could have implications for strategies that retard charge recombination in dye-sensitized solar cells.

Critical interpretation of CH- and OH- stretching regions for infrared spectra of methanol clusters (CH3OH)n (n = 2-5) using self-consistent-charge density functional tight-binding molecular dynamics simulations.

Vibrational infrared (IR) spectra of gas-phase O-H⋅⋅⋅O methanol clusters up to pentamer are simulated using self-consistent-charge density functional tight-binding method using two distinct methodologies: standard normal mode analysis and Fourier transform of the dipole time-correlation function. The twofold simulations aim at the direct critical assignment of the C-H stretching region of the recently recorded experimental spectra [H.-L. Han, C. Camacho, H. A. Witek, and Y.-P. Lee, J. Chem. Phys. 134, 144309 (2011)]. Both approaches confirm the previous assignment (ibid.) of the C-H stretching bands based on the B3LYP/ANO1 harmonic frequencies, showing that ν3, ν9, and ν2 C-H stretching modes of the proton-accepting (PA) and proton-donating (PD) methanol monomers experience only small splittings upon the cluster formation. This finding is in sharp discord with the assignment based on anharmonic B3LYP/VPT2/ANO1 vibrational frequencies (ibid.), suggesting that some procedural faults, likely related to the breakdown of the perturbational vibrational treatment, led the anharmonic calculations astray. The IR spectra based on the Fourier transform of the dipole time-correlation function include new, previously unaccounted for physical factors such as non-zero temperature of the system and large amplitude motions of the clusters. The elevation of temperature results in a considerable non-homogeneous broadening of the observed IR signals, while the presence of large-amplitude motions (methyl group rotations and PA-PD flipping), somewhat surprisingly, does not introduce any new features in the spectrum.

Energy Decomposition Scheme for Rectangular Graphene Flakes.

We show-to our own surprise-that total electronic energies for a family of m × n rectangular graphene flakes can be very accurately represented by a simple function of the structural parameters m and n with errors not exceeding 1 kcal/mol. The energies of these flakes, usually referred to as multiple zigzag chains Z(m,n), are computed for m, n < 21 at their optimized geometries using the DFTB3 methodology. We have discovered that the structural parameters m and n (and their simple algebraic functions) provide a much better basis for the energy decomposition scheme than the various topological invariants usually used in this context. Most terms appearing in our energy decomposition scheme seem to have simple chemical interpretations. Our observation goes against the well-established knowledge stating that many-body energies are complicated functions of molecular parameters. Our observations might have far-reaching consequences for building accurate machine learning models.

Quantum chemical investigation of epoxide and ether groups in graphene oxide and their vibrational spectra.

We present a detailed analysis of the factors influencing the formation of epoxide and ether groups in graphene nanoflakes using conventional density functional theory (DFT), the density-functional tight-binding (DFTB) method, π-Hückel theory, and graph theoretical invariants. The relative thermodynamic stability associated with the chemisorption of oxygen atoms at various positions on hexagonal graphene flakes (HGFs) of D(6h)-symmetry is determined by two factors - viz. the disruption of the π-conjugation of the HGF and the geometrical deformation of the HGF structure. The thermodynamically most stable structure is achieved when the former factor is minimized, and the latter factor is simultaneously maximized. Infrared (IR) spectra computed using DFT and DFTB reveal a close correlation between the relative thermodynamic stabilities of the oxidized HGF structures and their IR spectral activities. The most stable oxidized structures exhibit significant IR activity between 600 and 1800 cm(-1), whereas less stable oxidized structures exhibit little to no activity in this region. In contrast, Raman spectra are found to be less informative in this respect.

Spectral Evidence of Bevel-Gear-Type Rotation of Benzene around Br in Solid p-H2: Infrared Spectrum of the C6H6Br Radical.

Whether the structure of C6H6X (X = halogen), an intermediate in the halogenation of benzene, is an open or a bridged form has been debated. We produced Br to react with C6H6 upon photolysis in situ of a Br2/C6H6/p-H2 matrix at 3.2 K. In contrast to the C6H6Cl σ-complex reported previously, the observed infrared spectrum indicates that C6H6Br is an open-form π-complex. Furthermore, lines of the two CH out-of-plane bending modes associated mainly with even- and odd-numbered carbons, predicted near 672 and 719 cm-1, merged into a broad line at 697.3 cm-1, indicating that these modes become nearly equivalent as Br migrates from one carbon atom to another. Quantum-chemical calculations support that the benzene ring performs a bevel-gear-type rotation with respect to Br. Observation of only trans-ortho- and trans-para-C6H6Br2 suggests that this gear-type motion allows the additional Br atom to attack C6H6Br only from the opposite side of the Br atom in C6H6Br.

Is our way of thinking about excited states correct? Time-resolved dispersive IR study on p-nitroaniline.

Low-lying excited electronic states of an important class of molecules known as push-pull chromophores are central to understanding their potential nonlinear optical properties. Here we report that a combination of high-sensitivity nanosecond time-resolved dispersive IR spectroscopy and DFT calculations on p-nitroaniline (PNA), a prototypical push-pull molecule, reveals that PNA in the lowest excited triplet state has a partial quinoid structure. In this structure, the quinoid configuration is restricted to a part of the phenyl ring adjacent to the NO(2) group. The partial quinoid structure of PNA cannot be explained by a commonly used hybrid of a neutral form and a zwitterionic charge-transfer form. Our findings not only cast doubt on the general applicability of the classical way of looking at excited states, based exclusively on characteristic resonance structures, but also provide deeper insights into excited-state structure of highly polarizable molecular systems.

Dramatic reduction of IR vibrational cross sections of molecules encapsulated in carbon nanotubes.

Combined temperature-programmed desorption and IR studies suggest that absorption cross sections of IR-active vibrations of molecules "strongly" bound to single-wall carbon nanotubes (SWCNTs) are reduced at least by a factor of 10. Quantum chemical simulations show that IR intensities of endohedrally encapsulated molecules are dramatically reduced, and identify dielectric screening by highly polarizable SWCNT sidewalls as the origin of such "screening". The observed intensity reduction originates from a sizable cancellation of adsorbate dipole moments by mirror charges dynamically induced on the nanotube sidewalls. For exohedrally adsorbed molecules, the dielectric screening is found to be orientation-dependent with a smaller magnitude for adsorption in groove and interstitial sites. The presented results clearly demonstrate and quantify the screening effect of SWCNTs and unequivocally show that IR spectroscopy cannot be applied in a straightforward manner to the study of peapod systems.

Relative isomer abundance of fullerenes and carbon nanotubes correlates with kinetic stability.

A methodology to evaluate the kinetic stability of carbon nanostructures is presented based on the assumption of the independent and random nature of thermal vibrations. The kinetic stability is directly correlated to the cleavage probability for the weakest bond of a given nanostructure. The application of the presented method to fullerenes and carbon nanotubes yields clear correlation to their experimentally observed relative isomer abundances. The general and simple formulation of the method ensures its applicability to other nanostructures for which formation is controlled by kinetic factors.