An "inverse" harpoon mechanism.

Electron-transfer reactions are ubiquitous in chemistry and biology. The electrons' quantum nature allows their transfer across long distances. For example, in the well-known harpoon mechanism, electron transfer results in Coulombic attraction between initially neutral reactants, leading to a marked increase in the reaction rate. Here, we present a different mechanism in which electron transfer from a neutral reactant to a multiply charged cation results in strong repulsion that encodes the electron-transfer distance in the kinetic energy release. Three-dimensional coincidence imaging allows to identify such "inverse" harpoon products, predicted by nonadiabatic molecular dynamics simulations to occur between H2 and HCOH2+ following double ionization of isolated methanol molecules. These dynamics are experimentally initiated by single-photon double ionization with ultrafast extreme ultraviolet pulses, produced by high-order harmonic generation. A detailed comparison of measured and simulated data indicates that while the relative probability of long-range electron-transfer events is correctly predicted, theory overestimates the electron-transfer distance.

Sequential and concerted C-C and C-O bond dissociation in the Coulomb explosion of 2-propanol.

We study the competing mechanisms involved in the Coulomb explosion of 2-propanol CH3 2CHOH2+ dication, formed by an ultrafast extreme ultraviolet pulse. Over 20 product channels are identified and characterized using 3D coincidence imaging of the ionic fragments. The momentum correlations in the three-body fragmentation channels provide evidence for a dominant sequential mechanism, starting with the cleavage of a C-C bond, ejecting CH3 + and CH3CHOH+ cations, followed by a secondary fragmentation of the hydroxyethyl cation that can be delayed for up to a microsecond after ionization. The C-O bond dissociation channels are less frequent, involving proton transfer and double proton transfer, forming H2O+ and H3O+ products, respectively, and exhibiting mixed sequential and concerted character. These results can be explained by the high potential barrier for the C-O bond dissociation seen in our ab initio quantum chemical calculations. We also observe coincident COH+ + C2Hn + ions, suggesting exotic structural rearrangements, starting from the Frank-Condon geometry of the neutral 2-propanol system. Remarkably, the relative yield of the H3 + product is suppressed compared with methanol and alkene dications. Ab initio potentials and ground state molecular dynamics simulations show that a rapid and direct C-C bond cleavage dominates the Coulomb explosion process, leaving no time for H2 roaming, which is a necessary precursor to the H3 + formation.

Absence of Triplets in Single-Photon Double Ionization of Methanol.

Despite the abundance of data concerning single-photon double ionization of methanol, the spin state of the emitted electron pair has never been determined. Here we present the first evidence that identifies the emitted electron pair spin as overwhelmingly singlet when the dication forms in low-energy configurations. The experimental data show that while the yield of the CH2O+ + H3+ Coulomb explosion channel is abundant, the metastable methanol dication is largely absent. According to high-level ab initio simulations, these facts indicate that photoionization promptly forms singlet dication states, where they quickly decompose through various channels, with significant H3+ yields on the low-lying states. In contrast, if we assume that the initial dication is formed in one of the low-lying triplet states, the ab initio simulations exhibit a metastable dication, contradicting the experimental findings. Comparing the average simulated branching ratios with the experimental data suggests a >3 order of magnitude enhancement of the singlet:triplet ratio compared with their respective 1:3 multiplicities.

Time-resolving the ultrafast H2 roaming chemistry and H3+ formation using extreme-ultraviolet pulses.

The time scales and formation mechanisms of tri-hydrogen cation products in organic molecule ionization processes are poorly understood, despite their cardinal role in the chemistry of the interstellar medium and in other chemical systems. Using an ultrafast extreme-ultraviolet pump and time-resolved near-IR probe, combined with high-level ab initio molecular dynamics calculations, here we report unambiguously that H3+ formation in double-ionization of methanol occurs on a sub 100 fs time scale, settling previous conflicting findings of strong-field Coulomb explosion experiments. Our combined experimental-computational studies suggest that ultrafast competition, between proton-transfer and long-range electron-transfer processes, determines whether the roaming neutral H2 dynamics on the dication result in [Formula: see text] or [Formula: see text] fragments respectively.

Making Sense of Coulomb Explosion Imaging.

A multifaceted agreement between ab initio theoretical predictions and experimental measurements, including branching ratios, channel-specific kinetic energy release, and three-body momentum correlation spectra, leads to the identification of new mechanisms in Coulomb-explosion (CE) induced two- and three-body breakup processes in methanol. These identified mechanisms include direct nonadiabatic Coulomb explosion responsible for CO bond-breaking, a long-range " inverse harpooning" dominating the production of H2+ + HCOH+, a transient proton migration leading to surprising energy partitioning in three-body fragmentation and other complex dynamics forming products such as H2O+ and H3+. These mechanisms provide general concepts that should be useful for analyzing future time-resolved Coulomb explosion imaging of methanol as well as other molecular systems. These advances are enabled by a combination of recently developed experimental and computational techniques, using weak ultrafast EUV pulses to initiate the CE and a high-level quantum chemistry approach to follow the resulting field-free nonadiabatic molecular dynamics.

Computing Optimal Repairs for Functional Dependencies.

We investigate the complexity of computing an optimal repair of an inconsistent database, in the case where integrity constraints are Functional Dependencies (FDs). We focus on two types of repairs: an optimal subset repair (optimal S-repair) that is obtained by a minimum number of tuple deletions, and an optimal update repair (optimal U-repair) that is obtained by a minimum number of value (cell) up-dates. For computing an optimal S-repair, we present a polynomial-time algorithm that succeeds on certain sets of FDs and fails on others. We prove the following about the algorithm. When it succeeds, it can also incorporate weighted tuples and duplicate tuples. When it fails, the problem is NP-hard, and in fact, APX-complete (hence, cannot be approximated better than some constant). Thus, we establish a dichotomy in the complexity of computing an optimal S-repair. We present general analysis techniques for the complexity of computing an optimal U-repair, some based on the dichotomy for S-repairs. We also draw a connection to a past dichotomy in the complexity of finding a "most probable database" that satisfies a set of FDs with a single attribute on the left hand side; the case of general FDs was left open, and we show how our dichotomy provides the missing generalization and thereby settles the open problem.

Single-photon Coulomb explosion of methanol using broad bandwidth ultrafast EUV pulses.

Single-photon Coulomb explosion of methanol is instigated using the broad bandwidth pulse achieved through high-order harmonics generation. Using 3D coincidence fragment imaging of one molecule at a time, the kinetic energy release (KER) and angular distributions of the products are measured in different Coulomb explosion (CE) channels. Two-body CE channels breaking either the C-O or the C-H bonds are described as well as a proton migration channel forming H2O+, which is shown to exhibit higher KER. The results are compared to intense-field Coulomb explosion measurements in the literature. The interpretation of broad bandwidth single-photon CE data is discussed and supported by ab initio calculations of the predominant C-O bond breaking CE channel. We discuss the importance of these findings for achieving time resolved imaging of ultrafast dynamics.

A density functional theory for studying ionization processes in water clusters.

A generalized Kohn-Sham (GKS) approach to density functional theory (DFT), based on the Baer-Neuhauser-Livshits range-separated hybrid, combined with ab initio motivated range-parameter tuning is used to study properties of water dimer and pentamer cations. The water dimer is first used as a benchmark system to check the approach. The present brand of DFT localizes the positive charge (hole), stabilizing the proton transferred geometry in agreement with recent coupled-cluster calculations. Relative energies of various conformers of the water dimer cation compare well with previously published coupled cluster results. The GKS orbital energies are good approximations to the experimental ionization potentials of the system. Low-lying excitation energies calculated from time-dependent DFT based on the present method compare well with recently published high-level "equation of motion-coupled-cluster" calculations. The harmonic frequencies of the water dimer cation are in good agreement with experimental and wave function calculations where available. The method is applied to study the water pentamer cation. Three conformers are identified: two are Eigen type and one is a Zundel type. The structure and harmonic vibrational structure are analyzed. The ionization dynamics of a pentamer water cluster at 0 K shows a fast <50 fs transient for transferring a proton from one of the water molecules, releasing a hydroxyl radical and creating a protonated tetramer carrying the excess hole.

Tuned range-separated hybrids in density functional theory.

We review density functional theory (DFT) within the Kohn-Sham (KS) and the generalized KS (GKS) frameworks from a theoretical perspective for both time-independent and time-dependent problems. We focus on the use of range-separated hybrids within a GKS approach as a practical remedy for dealing with the deleterious long-range self-repulsion plaguing many approximate implementations of DFT. This technique enables DFT to be widely relevant in new realms such as charge transfer, radical cation dimers, and Rydberg excitations. Emphasis is put on a new concept of system-specific range-parameter tuning, which introduces predictive power in applications considered until recently too difficult for DFT.

Deleterious effects of long-range self-repulsion on the density functional description of O2 sticking on aluminum.

Density functional theory (DFT) with semilocal functionals such as the local-density and generalized gradients approximations predicts that the dissociative adsorption of oxygen on Al (111) goes through without a barrier in stark contradiction to experimental findings. This problem motivated our study of the reaction of oxygen colliding with a small aluminum cluster Al(5). We found semilocal functionals predict a minute barrier to sticking, associated with smeared long-range charge transfer from the metal to the oxygen. Hybrid B3LYP predicts a larger barrier while the range-separated the Baer-Neuhauser-Livshits (BNL, Phys. Chem. Chem. Phys. 2007, 9, 2932.) functional finds a more prominent barrier. BNL predicts short-ranged and more abrupt charge transfer from the surface to the oxygen. We conclude that spurious self-repulsion inherent in semilocal functionals causes early electron-transfer, long-range attraction toward the surface and low reaction barriers for these systems. The results indicate that the missing DFT barrier for O(2) sticking on Al (111) may be due to spurious self-repulsion.

A density functional theory for symmetric radical cations from bonding to dissociation.

It has been known for quite some time that approximate density functional (ADF) theories fail disastrously when describing the dissociative symmetric radical cations R(2)(+). By considering this dissociation limit, previous work has shown that Hartree-Fock (HF) theory favors the R(+1)-R(0) charge distribution, whereas DF approximations favor the R(+0.5)-R(+0.5). Yet, general quantum mechanical principles indicate that both these (as well as all intermediate) average charge distributions are asymptotically energy degenerate. Thus, HF and ADF theories mistakenly break the symmetry but in a contradicting way. In this letter, we show how to construct system-dependent long-range corrected (LC) density functionals that can successfully treat this class of molecules, avoiding the spurious symmetry breaking. Examples and comparisons to experimental data is given for R = H, He, and Ne, and it is shown that the new LC theory improves considerably the theoretical description of the R(2)(+) bond properties, the long-range form of the asymptotic potential curve, and the atomic polarizability. The broader impact of this finding is discussed as well, and it is argued that the widespread semiempirical approach which advocates treating the LC parameter as a system-independent parameter is in fact inappropriate under general circumstances.

A well-tempered density functional theory of electrons in molecules.

This Invited Article reports extensions of a recently developed approach to density functional theory with correct long-range behavior (R. Baer and D. Neuhauser, Phys. Rev. Lett., 2005, 94, 043002). The central quantities are a splitting functional gamma[n] and a complementary exchange-correlation functional E[n]. We give a practical method for determining the value of gamma in molecules, assuming an approximation for E is given. The resulting theory shows good ability to reproduce the ionization potentials for various molecules. However it is not of sufficient accuracy for forming a satisfactory framework for studying molecular properties. A somewhat different approach is then adopted, which depends on a density-independent gamma and an additional parameter w eliminating part of the local exchange functional. The values of these two parameters are obtained by best-fitting to experimental atomization energies and bond lengths of the molecules in the G2(1) database. The optimized values are gamma = 0.5 a and w = 0.1. We then examine the performance of this slightly semi-empirical functional for a variety of molecular properties, comparing to related works and experiment. We show that this approach can be used for describing in a satisfactory manner a broad range of molecular properties, be they static or dynamic. Most satisfactory is the ability to describe valence, Rydberg and inter-molecular charge-transfer excitations.

Time-dependent density-functional studies of the D2 coulomb explosion.

Real-time first principle simulations are presented of the D(2) Coulomb explosion dynamics detonated by exposure to very intense few-cycle laser pulse. Three approximate functionals within the time-dependent density functional theory (TDDFT) functionals are examined for describing the electron dynamics, including time-dependent Hartree-Fock theory. Nuclei are treated classically with quantum corrections. The calculated results are sensitive to the underlying electronic structure theory, showing too narrow kinetic energy distribution peaked at too high kinetic energy when compared with recent experimental results (Phys. Rev. Lett. 2003, 91, 093002). Experiment also shows a low energy peak which is not seen in the present calculation. We conclude that while Ehrenfest-adiabatic-TDDFT can qualitatively account for the dynamics, it requires further development, probably beyond the adiabatic approximation, to be quantitative.

Electrical or photocontrol of the rotary motion of a metallacarborane.

Rotary motion around a molecular axis has been controlled by simple electron transfer processes and by photoexcitation. The basis of the motion is intramolecular rotation of a carborane cage ligand (7,8-dicarbollide) around a nickel axle. The Ni(III) metallacarborane structure is a transoid sandwich with two pairs of carbon vertices reflected through a center of symmetry, but that of the Ni(IV) species is cisoid. The interconversion of the two provides the basis for controlled, rotational, oscillatory motion. The energies of the Ni(III) and Ni(IV) species are calculated as a function of the rotation angle.