Simulating Attochemistry: Which Dynamics Method to Use?

Attochemistry aims to exploit the properties of coherent electronic wavepackets excited via attosecond pulses to control the formation of photoproducts. Such molecular processes can, in principle, be simulated with various nonadiabatic dynamics methods, yet the impact of the approximations underlying the methods is rarely assessed. The performances of widely used mixed quantum-classical approaches, Tully surface hopping, and classical Ehrenfest methods are evaluated against the high-accuracy DD-vMCG quantum dynamics. This comparison is conducted for the valence ionization of fluorobenzene. Analyzing the nuclear motion induced in the branching space of the nearby conical intersection, the results show that the mixed quantum-classical methods reproduce quantitatively the average motion of a quantum wavepacket when initiated on a single electronic state. However, they fail to properly capture the nuclear motion induced by an electronic wavepacket along the derivative coupling, the latter originating from the quantum electronic coherence property, key to attochemistry.

Unraveling the variational breakdown of core valence separation calculations: Diagnostic and cure to the over relaxation error of double core hole states.

The core valence separation (CVS) approximation is the most employed strategy to prevent the variational collapse of standard wave function optimization when attempting to compute electronic states bearing one or more electronic vacancies in core orbitals. Here, we explore the spurious consequences of this approximation on the properties of the computed core hole states. We especially focus on the less studied case of double core hole (DCH) states, whose spectroscopic interest has recently been rapidly growing. We show that the CVS error leads to a systematic underestimation of DCH energies, a property in stark contrast with the case of single core hole states. We highlight that the CVS error can then be interpreted as an over relaxation effect and design a new correction strategy adapted to these specificities.

ESIPT in the pyrrol pyridine molecule: mechanism, timescale and yield revealed using dynamics simulations.

Excited State Intramolecular Proton Transfer in pyrrol pyridine is theoretically investigated using non-adiabatic dynamics simulations. The photochemical process is completely characterised: the reaction time, the total yield and the accessibility of the conical intersection are evaluated. Finally, new mechanistic interpretation are extracted: the proton transfer reaction in this molecule is shown to be driven by two complementary mechanisms.

Alternative Pathway to Double-Core-Hole States.

Excited double-core-hole states of isolated water molecules resulting from the sequential absorption of two x-ray photons have been investigated. These states are formed through an alternative pathway, where the initial step of core ionization is accompanied by the shake-up of a valence electron, leading to the same final states as in the core-ionization followed by core-excitation pathway. The capability of the x-ray free-electron laser to deliver very intense, very short, and tunable light pulses is fully exploited to identify the two different pathways.

Specific chemical bond relaxation unraveled by analysis of shake-up satellites in the oxygen single site double core hole spectrum of CO2.

We recently developed [A. Ferté, et al., J. Phys. Chem. Lett., 2020, 11, 4359] a method to compute single site double core hole (ssDCH or K-2) spectra. We refer to that method as NOTA+CIPSI. In the present paper this method is applied to the O K-2 spectrum of the CO2 molecule, and we use this as an example to discuss in detail its convergence properties. Using this approach, theoretical spectra in excellent agreement with the experimental one are obtained. Thanks to a thorough interpretation of the shake-up states responsible for the main satellite peaks and through comparison with the O K-2 spectrum of CO, we can highlight the clear signature of the two non-equivalent carbon oxygen bonds in the oxygen ssDCH CO2 dication.

Advanced Computation Method for Double Core Hole Spectra: Insight into the Nature of Intense Shake-up Satellites.

Double core hole spectroscopy is an ideal framework for investigating photoionization shake-up satellites. Their important intensity in a single site double core hole (ssDCH) spectrum allows the exploration of the subtle mix of relaxation and correlation effects associated with the inherent multielectronic character of the shake-up process. We present a high-accuracy computation method for single photon double core-shell photoelectron spectra that combines a selected configuration interaction procedure with the use of non-orthogonal molecular orbitals to obtain unbiased binding energy and intensity. This strategy leads to the oxygen ssDCH spectrum of the CO molecule that is in excellent agreement with the experimental result. Through a combined wave function and density analysis, we highlight that the most intense shake-up satellites are characterized by an electronic reorganization that opposes the core hole-induced relaxation.

Quantum Package 2.0: An Open-Source Determinant-Driven Suite of Programs.

Quantum chemistry is a discipline which relies heavily on very expensive numerical computations. The scaling of correlated wave function methods lies, in their standard implementation, between O(N5) and O(eN) , where N is proportional to the system size. Therefore, performing accurate calculations on chemically meaningful systems requires (i) approximations that can lower the computational scaling and (ii) efficient implementations that take advantage of modern massively parallel architectures. Quantum Package is an open-source programming environment for quantum chemistry specially designed for wave function methods. Its main goal is the development of determinant-driven selected configuration interaction (sCI) methods and multireference second-order perturbation theory (PT2). The determinant-driven framework allows the programmer to include any arbitrary set of determinants in the reference space, hence providing greater methodological freedom. The sCI method implemented in Quantum Package is based on the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) algorithm which complements the variational sCI energy with a PT2 correction. Additional external plugins have been recently added to perform calculations with multireference coupled cluster theory and range-separated density-functional theory. All the programs are developed with the IRPF90 code generator, which simplifies collaborative work and the development of new features. Quantum Package strives to allow easy implementation and experimentation of new methods, while making parallel computation as simple and efficient as possible on modern supercomputer architectures. Currently, the code enables, routinely, to realize runs on roughly 2 000 CPU cores, with tens of millions of determinants in the reference space. Moreover, we have been able to push up to 12 288 cores in order to test its parallel efficiency. In the present manuscript, we also introduce some key new developments: (i) a renormalized second-order perturbative correction for efficient extrapolation to the full CI limit and (ii) a stochastic version of the CIPSI selection performed simultaneously to the PT2 calculation at no extra cost.

Range-separated multideterminant density-functional theory with a short-range correlation functional of the on-top pair density.

We introduce an approximation to the short-range correlation energy functional with multideterminantal reference involved in a variant of range-separated density-functional theory. This approximation is a local functional of the density, the density gradient, and the on-top pair density, which locally interpolates between the standard Perdew-Burke-Ernzerhof correlation functional at a vanishing range-separation parameter and the known exact asymptotic expansion at a large range-separation parameter. When combined with (selected) configuration-interaction calculations for the long-range wave function, this approximation gives accurate dissociation energy curves of the H2, Li2, and Be2 molecules and thus appears as a promising way to accurately account for static correlation in range-separated density-functional theory.

Curing basis-set convergence of wave-function theory using density-functional theory: A systematically improvable approach.

The present work proposes to use density-functional theory (DFT) to correct for the basis-set error of wave-function theory (WFT). One of the key ideas developed here is to define a range-separation parameter which automatically adapts to a given basis set. The derivation of the exact equations are based on the Levy-Lieb formulation of DFT, which helps us to define a complementary functional which corrects uniquely for the basis-set error of WFT. The coupling of DFT and WFT is done through the definition of a real-space representation of the electron-electron Coulomb operator projected on a one-particle basis set. Such an effective interaction has the particularity to coincide with the exact electron-electron interaction in the limit of a complete basis set, and to be finite at the electron-electron coalescence point when the basis set is incomplete. The non-diverging character of the effective interaction allows one to define a mapping with the long-range interaction used in the context of range-separated DFT and to design practical approximations for the unknown complementary functional. Here, a local-density approximation is proposed for both full-configuration-interaction (FCI) and selected configuration-interaction approaches. Our theory is numerically tested to compute total energies and ionization potentials for a series of atomic systems. The results clearly show that the DFT correction drastically improves the basis-set convergence of both the total energies and the energy differences. For instance, a sub kcal/mol accuracy is obtained from the aug-cc-pVTZ basis set with the method proposed here when an aug-cc-pV5Z basis set barely reaches such a level of accuracy at the near FCI level.