Accurate density functional calculations hinge on reliable approximations to the unknown exchange-correlation (xc) potential. The most popular approximations usually lack features of the exact xc potential that are important for an accurate prediction of the fundamental gap and the distribution of charge in complex systems. Two principal features in this regard are the spatially uniform shift in the potential, as the number of electrons infinitesimally surpasses an integer, and the spatial steps that form, for example, between the atoms of stretched molecules. Although both aforementioned concepts are well known, the exact relationship between them remained unclear. Here we establish this relationship via an analytical derivation. We support our result by numerically solving the many-electron Schrödinger equation to extract the exact Kohn-Sham potential and directly observe its features. Spatial steps in the exact xc potential of a full configuration-interaction (FCI) calculation of a molecule are presented in three dimensions.

Abstractorbkit is a toolbox for post-processing electronic structure calculations based on a highly modular and portable Python architecture. The program allows computing a multitude of electronic properties of molecular systems on arbitrary spatial grids from the basis set representation of its electronic wavefunction, as well as several gridindependent properties. The required data can be extracted directly from the standard output of a large number of quantum chemical programs. orbkit can be used as a standalone program to determine standard quantities, for example, the electron density, molecular orbitals, and derivatives thereof. The cornerstone of orbkit is its modular structure. The existing basic functions can be arranged in an individual way and can be easily extended by user-written modules to determine any other desired quantities. orbkit offers multiple output formats that can be processed by common visualization tools (VMD, Molden, etc.). Additionally, orbkit possesses routines to order molecular orbitals computed at different nuclear configurations according to their electronic character and to interpolate the wavefunction between these configurations. The program is open-source under GNU-LGPLv3 license and freely available at http://sourceforge.net/projects/orbkit/. This article provides an overview of orbkit with particular focus on its capabilities and applicability, and includes several example calculations.Keywords: quantum chemical calculation, electronic structure, molecular visualization, electron density, grid representation of one-electron quantities, molecular orbital ordering * Institut für Chemie und Biochemie, Freie Universität Berlin, Takustraße 3, 14195 Berlin, Germany † These authors contributed equally to this work. orbkit is an open-source toolbox for post-processing electronic structure calculations. Based on a highly modular and portable Python architecture, it comes both as a standalone program and a function library. The program allows computing electronic properties of molecular systems on arbitrary spatial grids from the output of standard quantum chemistry programs.2

Electronic dynamics in liquids are of fundamental importance, but time-resolved experiments have so far remained limited to the femtosecond time scale. We report the extension of attosecond spectroscopy to the liquid phase. We measured time delays of 50 to 70 attoseconds between the photoemission from liquid water and that from gaseous water at photon energies of 21.7 to 31.0 electron volts. These photoemission delays can be decomposed into a photoionization delay sensitive to the local environment and a delay originating from electron transport. In our experiments, the latter contribution is shown to be negligible. By referencing liquid water to gaseous water, we isolated the effect of solvation on the attosecond photoionization dynamics of water molecules. Our methods define an approach to separating bound and unbound electron dynamics from the structural response of the solvent.

Recent experimental investigations by the group of D. T. Anderson (Kettwich, S. C.; Raston, P. L.; Anderson, D. T. J. Phys. Chem. A 2009, 113, DOI 10.1021/jp811206a) show that the reaction Cl + H(2) --> HCl + H in the para-H(2) crystal can be induced by infrared (IR) + ultraviolet (UV) coirradiations causing vibrational pre-excitation of the molecular reactant, H(2)(v=1), and generation of the atomic reactant, Cl((2)P(3/2)), by near-resonant photodissociation of a matrix-isolated Cl(2) molecule in the C (1)Pi(u) state, respectively. The corresponding reaction probability P(v=1) for the reactants Cl + H(2)(v=1) is approximately 0.15; this is approximately 25 times larger than P(v=0) for Cl + H(2)(v=0) (as initiated by pure UV irradiation). We present a simple three-step quantum model which accounts for some important parts of the experimental results and allows predictions for other scenarios, for example, UV photodissociation of the Cl(2) molecule by a laser pulse. The first step, vibrational pre-excitation of H(2), yields the molecular initial state which is described using the Einstein model of the para-H(2) crystal. The second step, photodissociation of Cl(2), generates the Cl((2)P(3/2)) atom approaching H(2)(v=1). In the third step, Cl reacts with H(2)(v=1) much more efficiently than with H(2)(v=0) close to threshold. The ultrashort time domains (approximately 100 fs) of steps 2 plus 3 support one- and then two-dimensional models of photodissociation of Cl(2) by short laser pulses and of the subsequent reaction of the system Cl-H-H embedded in frozen environments. The widths of the corresponding wave function describing the translational motion of the reactants is revealed as a significant parameter which is determined not only by the duration of the laser pulse but, even more importantly, by the width of the Gaussian-type distribution of the center of mass of the H(2) molecule in its Einstein cell. As a consequence, the resulting P(v) are quite robust versus variations of the UV pulse durations, allowing extrapolations to continuous wave irradiation. Quantum dynamics simulations of the reaction reveal that the experimental results are due to energetic and dynamical effects.