Impact of Early Coherences on the Control of Ultrafast Photodissociation Reactions.

By coherent control, the yield of photodissociation reactions can be maximized, starting in a suitable superposition of vibrational states. In ultrafast processes, the interfering pathways are born from the early vibrational coherences in the ground electronic potential. We interpret their effect from a purely classical picture, in which the correlation between the initial position and momentum helps to synchronize the vibrational dynamics at the Franck-Condon window when the pulse is at its maximum intensity. In the quantum domain, we show that this localization in time and space is mediated by dynamic squeezing of the wave packet.

Can we predict ambident regioselectivity using the chemical hardness?

The hard/soft acid/base (HSAB) principle is a cornerstone in our understanding of chemical reactivity preferences. Motivated by the success of the original ("global") version of this rule, a "local" counterpart was readily proposed to account for regioselectivity preferences, in particular, in ambident reactions. However, ample experimental evidence indicates that the local HSAB principle often fails to provide meaningful predictions. Here we examine the assumptions behind the standard proof of the local HSAB rule, showing that it is based on a flawed premise. By solving this issue, we show that it is critical to consider not only the charge transferred between the different reacting centers but also the charge reorganization within the non-reacting parts of the molecule. We propose different reorganization models and derive the corresponding regioselectivity rules for each.

Imaging the photodissociation dynamics of internally excited ethyl radicals from high Rydberg states.

The site-specific hydrogen-atom elimination mechanism previously reported for photoexcited ethyl radicals (CH3CH2) [D. V. Chicharro et al., Chem. Sci., 2019, 10, 6494] is interrogated in the photodissociation of the ethyl isotopologues CD3CD2, CH3CD2 and CD3CH2 through the velocity map imaging (VMI) detection of the produced hydrogen- and deuterium-atoms. The radicals, generated in situ from photolysis of a precursor using the same laser pulse employed in their excitation to Rydberg states, decompose along the Cα-H/D and Cß-H/D reaction coordinates through coexisting statistical and site-specific mechanisms. The experiments are carried out at two excitation wavelengths, 201 and 193 nm. The comparison between both sets of results provides accurate information regarding the primary role in the site-specific mechanism of the radical internal reservoir. Importantly, at 193 nm excitation, higher energy dissociation channels (not observed at 201 nm) producing low-recoil H/D-atoms become accessible. High-level ab initio calculations of potential energy curves and the corresponding non-adiabatic interactions allow us to rationalize the experimental results in terms of competitive non-adiabatic decomposition paths. Finally, the adiabatic behavior of the conical intersections in the face of several vibrational modes - the so-called vibrational promoting modes - is discussed.

Hardness of molecules and bandgap of solids from a generalized gradient approximation exchange energy functional.

The deviations from linearity of the energy as a function of the number of electrons that arise with current approximations to the exchange-correlation (XC) energy functional have important consequences for the frontier eigenvalues of molecules and the corresponding valence-band maxima for solids. In this work, we present an analysis of the exact theory that allows one to infer the effects of such approximations on the highest occupied and lowest unoccupied molecular orbital eigenvalues. Then, we show the importance of the asymptotic behavior of the XC potential in the generalized gradient approximation (GGA) in the case of the NCAPR functional (nearly correct asymptotic potential revised) for determining the shift of the frontier orbital eigenvalues toward the exact values. Thereby we establish a procedure at the GGA level of refinement that allows one to make a single calculation to determine the ionization potential, the electron affinity, and the hardness of molecules (and its solid counterpart, the bandgap) with an accuracy equivalent to that obtained for those properties through energy differences, a procedure that requires three calculations. For solids, the accuracy achieved for the bandgap lies rather close to that which is obtained through hybrid XC energy functionals, but those also demand much greater computational effort than what is required with the simple NCAPR GGA calculation.

Tryptophan regulates Drosophila zinc stores.

Zinc deficiency is commonly attributed to inadequate absorption of the metal. Instead, we show that body zinc stores in Drosophila melanogaster depend on tryptophan consumption. Hence, a dietary amino acid regulates zinc status of the whole insect­a finding consistent with the widespread requirement of zinc as a protein cofactor. Specifically, the tryptophan metabolite kynurenine is released from insect fat bodies and induces the formation of zinc storage granules in Malpighian tubules, where 3-hydroxykynurenine and xanthurenic acid act as endogenous zinc chelators. Kynurenine functions as a peripheral zinc-regulating hormone and is converted into a 3-hydroxykynurenine­zinc­chloride complex, precipitating within the storage granules. Thus, zinc and the kynurenine pathway­well-known modulators of immunity, blood pressure, aging, and neurodegeneration­are physiologically connected.

Structural and electronic analysis of the octarepeat region of prion protein with four Cu2+ by polarizable MD and QM/MM simulations.

Prions have been linked to neurodegenerative diseases that affect various species of mammals including humans. The prion protein, located mainly in neurons, is believed to play the role of metal ion transporter. High levels of copper ions have been related to structural changes. A 32-residue region of the N-terminal domain, known as octarepeat, can bind up to four copper ions. Different coordination modes have been observed and are strongly dependent on Cu2+ concentration. Many theoretical studies carried out so far have focused on studying the coordination modes of a single copper ion. In this work we investigate the octarepeat region coordinated with four copper ions. Molecular dynamics (MD) and hybrid quantum mechanics/molecular mechanics (QM/MM) simulations using the polarizable AMOEBA force field have been carried out. The polarizable MD simulations starting from a fully extended conformation indicate that the tetra-Cu2+/octarepeat complex forms a globular structure. The globular form is stabilized by interactions between Cu2+ and tryptophan residues resulting in some coordination sites observed to be in close proximity, in agreement with experimental results. Subsequent QM/MM simulations on several snapshots suggests the system is in a high-spin quintet state, with all Cu2+ bearing one single electron, and all unpaired electrons are ferromagnetically coupled. NMR simulations on selected structures provides insights on the chemical shifts of the first shell ligands around the metals with respect to inter-metal distances.

Threshold Photoelectron Spectroscopy of the CH2I, CHI, and CI Radicals.

VUV photoionization of the CHnI radicals (with n = 0, 1, and 2) is investigated by means of synchrotron radiation coupled with a double imaging photoion-photoelectron coincidence spectrometer. Photoionization efficiencies and threshold photoelectron spectra (TPES) for photon energies ranging between 9.2 and 12.0 eV are reported. An adiabatic ionization energy (AIE) of 8.334 ± 0.005 eV is obtained for CH2I, which is in good agreement with previous results [8.333 ± 0.015 eV, Sztáray J. Chem. Phys. 2017, 147, 013944], while for CI an AIE of 8.374 ± 0.005 eV is measured for the first time and a value of ∼8.8 eV is estimated for CHI. Ab initio calculations have been carried out for the ground state of the CH2I radical and for the ground state and excited states of the radical cation CH2I+, including potential energy curves along the C-I coordinate. Franck-Condon factors are calculated for transitions from the CH2I(X̃2B1) ground state of the neutral radical to the ground state and excited states of the radical cation. The TPES measured for the CH2I radical shows several structures that correspond to the photoionization into excited states of the radical cation and are fully assigned on the basis of the calculations. The TPES obtained for the CHI is characterized by a broad structure peaking at 9.335 eV, which could be due to the photoionization from both the singlet and the triplet states and into one or more electronic states of the cation. A vibrational progression is clearly observed in the TPES for the CI radical and a frequency for the C-I stretching mode of 760 ± 60 cm-1 characterizing the CI+ electronic ground state has been extracted.

Site-specific hydrogen-atom elimination in photoexcited alkyl radicals.

A prompt site-specific hydrogen-atom elimination from the α-carbon atom (Cα) has been recently reported to occur in the photodissociation of ethyl radicals following excitation at 201 nm [Chicharro et al., Chem. Sci., 2019, 10, 6494]. Such pathway was accessed by means of an initial ro-vibrational energy characterizing the radicals produced by in situ photolysis of a precursor. Here, we present experimental evidence of a similar dynamics in a series of alkyl radicals (C2H5, n-C3H7, n-C4H9, and i-C3H7) containing the same reaction coordinate, but different extended structures. The main requirements for the site-specific mechanism in the studied radicals, namely a rather high content of internal energy prior to dissociation and the participation of vibrational promoting modes, is discussed in terms of the chemical structure of the radicals. The methyl deformation mode in all alkyl radicals along with the CH bending motion in i-C3H7 appear to promote this fast H-atom elimination channel. The photodissociation dynamics of the simplest unsaturated alkyl radical, the vinyl radical (C2H3), is also discussed, showing no signal of site-specific fast H-atom elimination. The results are complemented with high-level ab initio electronic structure calculations of potential energy curves of the vinyl radical, which are compared with those previously reported for the ethyl radical.

Spectroscopic properties of open shell diatomic molecules using Piris natural orbital functionals.

Spectroscopic properties such as equilibrium distances, vibrational constants, rotational constants, dissociation energies, and excitation energies are calculated for nine heteronuclear diatomic molecules (PH, NF, NH, NO, CS, AlF, ClF, BeO and CF) using an interactive pair model (PNOF7s), that has been generalized for spin multiplet states, and its second order perturbation variant, NOF-MP2, which was also generalized for multiplets. The results obtained are compared with Complete Active Space (CASSCF) and Complete Active Space Perturbation Theory (CASPT2). It is shown that the potential energy curves provided by the PNOF functional for open shell diatomic molecules are in acceptable agreement with those from CASSCF and CASPT2. The spectroscopic constants depending at most on the second derivative of the potential energy are in good agreement with experiment, while those requiring the evaluation of the third and fourth derivatives show larger deviations from experiment and from those predicted by CASPT2. Thus, it is shown that the PNOF functional extension to multiplets is an alternative approach in predicting spectroscopic constants of molecules where static correlation plays an important role, like the open shell heteronuclear diatomic molecules studied in this work.

Temperature-Dependent Approach to Electronic Charge Transfer.

A charge transfer model is developed within the framework of the grand canonical ensemble through the analysis of the behavior of the fractional charge as a function of the chemical potential of the bath when the temperature and the external chemical potential are kept fixed. Departing from the fact that, before the interaction between two species, each one has a zero fractional charge, one can identify two situations after the interaction occurs where the fractional charge of at least one of the species is different from zero, indicating that there has been charge transference. One of them corresponds to the case when one of the species is immersed in a bath conformed by the other one, while the other is related to the case in which both species are present in equal amounts (stoichiometric proportion). Correlations between the fractional charges and average energies, thus obtained with experimental equilibrium constants, kinetic rate constants, hydration constants, and bond enthalpies, indicate that, although at the experimental temperatures, they are very small quantities, they have chemically meaningful information. Additionally, in the stoichiometric case, one also finds a rather good correlation between the equalized chemical potential and the one obtained from experimental information for a test set of diatomic and triatomic molecules.

Negative Electron Affinities and Derivative Discontinuity Contribution from a Generalized Gradient Approximation Exchange Functional.

Two methods to calculate negative electron affinities systematically from ground-state density functional methods are presented. One makes use of the lowest unoccupied molecular orbital energy shift provided by approximate inclusion of derivative discontinuity in the nearly correct asymptotic potential (NCAP) nonempirical, constraint-based generalized gradient approximation exchange functional. The other uses a second-order perturbation calculation of the derivative discontinuity based on the NCAP exchange-correlation potential. On a set of thirty-eight molecules, NCAP leads to a rather accurate description that is improved further through the perturbation correction. The results presented show the importance of the asymptotic behavior of the exchange-correlation potential in the calculation of negative electron affinities as well as demonstrating the versatility of the NCAP functional.

Femtochemistry under scrutiny: Clocking state-resolved channels in the photodissociation of CH3I in the A-band.

Clocking of electronically and vibrationally state-resolved channels of the fast photodissociation of CH3I in the A-band is re-examined in a combined experimental and theoretical study. Experimentally, a femtosecond pump-probe scheme is employed in the modality of resonant probing by resonance enhanced multiphoton ionization (REMPI) of the methyl fragment in different vibrational states and detection through fragment velocity map ion (VMI) imaging as a function of the time delay. We revisit excitation to the center of the A-band at 268 nm and report new results for excitation to the blue of the band center at 243 nm. Theoretically, two approaches have been employed to shed light into the observations: first, a reduced dimensionality 4D nonadiabatic wavepacket calculation using the potential energy surfaces by Xie et al. [J. Phys. Chem. A 104, 1009 (2000)]; and second, a full dimension 9D trajectory surface-hopping calculation on the same potential energy surfaces, including the quantization of vibrational states of the methyl product. In addition, high level ab initio electronic structure calculations have been carried out to describe the CH3 3pz Rydberg state involved in the (2 + 1) REMPI probing process, as a function of the carbon-iodine (C-I) distance. A general qualitative agreement is obtained between experiment and theory, but the effect of methyl vibrational excitation in the umbrella mode on the clocking times is not well reproduced. The theoretical results reveal that no significant effect on the state-resolved appearance times is exerted by the nonadiabatic crossing through the conical intersection present in the first absorption band. The vibrationally state resolved clocking times observed experimentally can be rationalized when the (2 + 1) REMPI probing process is considered. None of the other probing methods applied thus far, i.e., multiphoton ionization photoelectron spectroscopy, soft X-ray inner-shell photoelectron spectroscopy, VUV single-photon ionization, and XUV core-to-valence transient absorption spectroscopy, have been able to provide quantum state-resolved (vibrational) clocking times. More experiments would be needed to disentangle the fine details in the clocking times and dissociation dynamics arising from the detection of specific quantum-states of the molecular fragments.

Introduction to celebrating recent chemical science in Mexico.

Gabriel Merino, María de Jesús Rosales and Alberto Vela introduce the PCCP, Dalton Transactions, New Journal of Chemistry and RSC Advances themed issue on celebrating recent chemical science in Mexico.

The 3s versus 3p Rydberg state photodissociation dynamics of the ethyl radical.

The photodissociation dynamics of the ethyl radical following excitation into the 3s and 3p Rydberg states are revisited in a joint experimental and theoretical study. Two different methods to produce the ethyl radical, pyrolysis and in situ photolysis, are employed in order to modify the initial ro-vibrational energy distribution characterizing the ethyl radical beam. H-atom velocity map images following excitation of the radical at 243 nm and at 201 nm are presented and discussed along with ab initio potential energy curves focussing on the bridged C2v geometry. The reported results show that the dynamics following excitation to the 3s Rydberg state is insensitive to the initial internal energy of the parent radical, in contrast to the dynamics on the 3p Rydberg state, which is strongly modified. The role of the bridged C2v geometry on both photodynamics is highlighted and discussed.

An Anionic Ring Locked into an Anionic Axle: A Metastable Rotaxane with Chemically Activated Electrostatic Stoppers.

The use of the electrostatic stoppers concept in the field of mechanically interlocked molecules is reported; these stoppers are chemically sensitive end groups on a linear guest molecule that allows for the conversion of a pseudo-rotaxane species into a rotaxane complex by a change in the medium acidity. The chemical stimulus causes the appearance of negative charges on both ends of the linear component, passing from cationic to anionic, and causing a significant ring-to-axle electrostatic repulsion. This phenomenon has two different and simultaneous effects: 1) destabilizes the complex as a consequence of confining an anionic ring into an anionic axle, and 2) increases the dissociation energy barrier, thus impeding ring extrusion. This newly formed metastable rotaxane species is resistant to solvent and temperature effects and performs as a two-state degenerated molecular shuttle in solution.

Site-specific hydrogen-atom elimination in photoexcited ethyl radical.

The photochemistry of the ethyl radical following excitation to the 3p Rydberg state is investigated in a joint experimental and theoretical study. Velocity map images for hydrogen atoms detected from photoexcited isotopologues CH3CH2, CH3CD2 and CD3CH2 at ∼201 nm, are discussed along with high-level ab initio electronic structure calculations of potential energy curves and non-adiabatic coupling matrix elements (NACME). A novel mechanism governed by a conical intersection allowing prompt site-specific hydrogen-atom elimination is presented and discussed. For this mechanism to occur, an initial ro-vibrational excitation is allocated to the radical permitting to access this reaction pathway and thus to control the ethyl photochemistry. While hydrogen-atom elimination from cold ethyl radicals occurs through internal conversion into lower electronic states followed by slow statistical dissociation, prompt site-specific Cα elimination into CH3CH + H, occurring through a fast non-adiabatic crossing to a valence bound state followed by dissociation through a conical intersection, is accessed by means of an initial ro-vibrational energy content into the radical. The role of a particularly effective vibrational promoting mode in this prompt photochemical reaction pathway is discussed.

Natural orbital functional for spin-polarized periodic systems.

Natural orbital functional theory is considered for systems with one or more unpaired electrons. An extension of the Piris natural orbital functional (PNOF) based on electron pairing approach is presented, specifically, we extend the independent pair model, PNOF5, and the interactive pair model PNOF7 to describe spin-uncompensated systems. An explicit form for the two-electron cumulant of high-spin cases is only taken into account, so that singly occupied orbitals with the same spin are solely considered. The rest of the electron pairs with opposite spins remain paired. The reconstructed two-particle reduced density matrix fulfills certain N-representability necessary conditions, as well as guarantees the conservation of the total spin. The theory is applied to model systems with strong non-dynamic (static) electron correlation, namely, the one-dimensional Hubbard model with periodic boundary conditions and hydrogen rings. For the latter, PNOF7 compares well with exact diagonalization results so the model presented here is able to provide a correct description of the strong-correlation effects.

Generalized Gradient Approximation Exchange Energy Functional with Near-Best Semilocal Performance.

We develop and validate a nonempirical generalized gradient approximation (GGA) exchange (X) density functional that performs as well as the SCAN (strongly constrained and appropriately normed) meta-GGA on standard thermochemistry tests. Additionally, the new functional (NCAP, nearly correct asymptotic potential) yields Kohn-Sham eigenvalues that are useful approximations of the density functional theory (DFT) ionization potential theorem values by inclusion of a systematic derivative discontinuity shift of the X potential. NCAP also enables time-dependent DFT (TD-DFT) calculations of good-quality polarizabilities, hyper-polarizabilities, and one-Fermion excited states without modification (calculated or ad hoc) of the long-range behavior of the exchange potential or other patches. NCAP is constructed by reconsidering the imposition of the asymptotic correctness of the X potential (-1/ r) as a constraint. Inclusion of derivative discontinuity and approximate integer self-interaction correction treatments along with first-principles determination of the effective second-order gradient expansion coefficient yields a major advance over our earlier correct asymptotic potential functional [CAP; J. Chem. Phys. 2015 , 142 , 054105 ]. The new functional reduces a spurious bump in the CAP atomic exchange potential and moves it to distances irrelevantly far from the nucleus (outside the tail of essentially all practical basis functions). It therefore has nearly correct atomic exchange-potential behavior out to rather large finite distances r from the nucleus but eventually goes as - c/ r with an estimated value for the constant c of around 0.3, so as to achieve other important properties of exact DFT exchange within the restrictions of the GGA form. We illustrate the results with the Ne atom optimized effective potentials and with standard molecular benchmark test data sets for thermochemical, structural, and response properties.