Pierre-François Loos

A theoretical study of Ru(II) polypyridyl DNA intercalators structure and electronic absorption spectroscopy of [Ru(phen)2(dppz)]2+ and [Ru(tap)2(dppz)]2+ complexes intercalated in guanine-cytosine base pairs.

The structural and spectroscopic properties of [Ru(phen)(2)(dppz)](2+) and [Ru(tap)(2)(dppz)](2+) (phen=1,10-phenanthroline; tap=1,4,5,8-tetraazaphenanthrene; dppz=dipyridophenazine ) have been investigated by means of density functional theory (DFT), time-dependent DFT (TD-DFT) within the polarized continuum model (IEF-PCM) and quantum mechanics/molecular mechanics (QM/MM) calculations. The model of the Delta and Lambda enantiomers of Ru(II) intercalated in DNA in the minor and major grooves is limited to the metal complexes intercalated in two guanine-cytosine base pairs. The main experimental spectral features of these complexes reported in DNA or synthetic polynucleotides are better reproduced by the theoretical absorption spectra of the Delta enantiomers regardless of intercalation mode (major or minor groove). This is especially true for [Ru(phen)(2)(dppz)](2+). The visible absorption of [Ru(tap)(2)(dppz)](2+) is governed by the MLCT(tap) transitions regardless of the environment (water, acetonitrile or bases pair), the visible absorption of [Ru(phen)(2)(dppz)](2+) is characterized by transitions to metal-to-ligand-charge-transfer MLCT(dppz) in water and acetonitrile and to MLCT(phen) when intercalated in DNA. The response of the IL(dppz) state to the environment is very sensitive. In vacuum, water and acetonitrile these transitions are characterized by significant oscillator strengths and their positions depend significantly on the medium with blue shifts of about 80 nm when going from vacuum to solvent. When the complex is intercalated in the guanine-cytosine base pairs the (1)IL(dppz) transition contributes mainly to the band at 370 nm observed in the spectrum of [Ru(phen)(2)(dppz)](2+) and to the band at 362 nm observed in the spectrum of [Ru(tap)(2)(dppz)](2+).

Hybrid stochastic-deterministic calculation of the second-order perturbative contribution of multireference perturbation theory

A hybrid stochastic-deterministic approach for computing the second-order perturbative contribution E(2) within multireference perturbation theory (MRPT) is presented. The idea at the heart of our hybrid scheme-based on a reformulation of E(2) as a sum of elementary contributions associated with each determinant of the MR wave function-is to split E(2) into a stochastic and a deterministic part. During the simulation, the stochastic part is gradually reduced by dynamically increasing the deterministic part until one reaches the desired accuracy. In sharp contrast with a purely stochastic Monte Carlo scheme where the error decreases indefinitely as t-1/2 (where t is the computational time), the statistical error in our hybrid algorithm displays a polynomial decay ∼t-n with n = 3-4 in the examples considered here. If desired, the calculation can be carried on until the stochastic part entirely vanishes. In that case, the exact result is obtained with no error bar and no noticeable computational overhead compared to the fully deterministic calculation. The method is illustrated on the F2 and Cr2 molecules. Even for the largest case corresponding to the Cr2 molecule treated with the cc-pVQZ basis set, very accurate results are obtained for E(2) for an active space of (28e, 176o) and a MR wave function including up to 2×107 determinants.

Deterministic construction of nodal surfaces within quantum Monte Carlo: the case of FeS

In diffusion Monte Carlo (DMC) methods, the nodes (or zeroes) of the trial wave function dictate the magnitude of the fixed-node (FN) error. In standard DMC implementations, the nodes are optimized by stochastically optimizing a short multideterminant expansion in the presence of an explicitly correlated Jastrow factor. Here, following a recent proposal, we pursue a different route and consider the nodes of selected configuration interaction (sCI) expansions built with the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) algorithm. By increasing the size of the sCI expansion, these nodes can be systematically and deterministically improved. The present methodology is used to investigate the properties of the transition metal sulfide molecule FeS. This apparently simple molecule has been shown to be particularly challenging for electronic structure theory methods due to the proximity of two low-energy quintet electronic states of different spatial symmetry and the difficulty to treat them on equal footing from a one-electron basis set point of view. In particular, we show that, at the triple-ζ basis set level, all sCI results-including those extrapolated at the full CI (FCI) limit-disagree with experiment, yielding an electronic ground state of 5Σ+ symmetry. Performing FN-DMC simulation with sCI nodes, we show that the correct 5Δ ground state is obtained if sufficiently large expansions are used. Moreover, we show that one can systematically get accurate potential energy surfaces and reproduce the experimental dissociation energy as well as other spectroscopic constants.

A Mountaineering Strategy to Excited States: Highly Accurate Reference Energies and Benchmarks

Striving to define very accurate vertical transition energies, we perform both high-level coupled cluster (CC) calculations (up to CCSDTQP) and selected configuration interaction (sCI) calculations (up to several millions of determinants) for 18 small compounds (water, hydrogen sulfide, ammonia, hydrogen chloride, dinitrogen, carbon monoxide, acetylene, ethylene, formaldehyde, methanimine, thioformaldehyde, acetaldehyde, cyclopropene, diazomethane, formamide, ketene, nitrosomethane, and the smallest streptocyanine). By systematically increasing the order of the CC expansion, the number of determinants in the CI expansion as well as the size of the one-electron basis set, we have been able to reach near full CI (FCI) quality transition energies. These calculations are carried out on CC3/ aug-cc-pVTZ geometries, using a series of increasingly large atomic basis sets systematically including diffuse functions. In this way, we define a list of 110 transition energies for states of various characters (valence, Rydberg, n → π*, π → π*, singlet, triplet, etc.) to be used as references for further calculations. Benchmark transition energies are provided at the aug-cc-pVTZ level as well as with additional basis set corrections, in order to obtain results close to the complete basis set limit. These reference data are used to benchmark a series of 12 excited-state wave function methods accounting for double and triple contributions, namely ADC(2), ADC(3), CIS(D), CIS(D∞), CC2, STEOM-CCSD, CCSD, CCSDR(3), CCSDT-3, CC3, CCSDT., and CCSDTQ. It turns out that CCSDTQ yields a negligible difference with the extrapolated CI values with a mean absolute error as small as 0.01 eV, whereas the coupled cluster approaches including iterative triples are also very accurate (mean absolute error of 0.03 eV). Consequently, CCSDT-3 and CC3 can be used to define reliable benchmarks. This observation does not hold for ADC(3) that delivers quite large errors for this set of small compounds, with a clear tendency to overcorrect its second-order version, ADC(2). Finally, we discuss the possibility to use basis set extrapolation approaches so as to tackle more easily larger compounds.

Excitation energies from diffusion Monte Carlo using selected configuration interaction nodes

Quantum Monte Carlo (QMC) is a stochastic method that has been particularly successful for ground-state electronic structure calculations but mostly unexplored for the computation of excited-state energies. Here, we show that within a Jastrow-free QMC protocol relying on a deterministic and systematic construction of nodal surfaces using selected configuration interaction (sCI) expansions, one is able to obtain accurate excitation energies at the fixed-node diffusion Monte Carlo (FN-DMC) level. This evidences that the fixed-node errors in the ground and excited states obtained with sCI wave functions cancel out to a large extent. Our procedure is tested on two small organic molecules (water and formaldehyde) for which we report all-electron FN-DMC calculations. For both the singlet and triplet manifolds, accurate vertical excitation energies are obtained with relatively compact multideterminant expansions built with small (typically double-ζ) basis sets.

Theoretical 0–0 Energies with Chemical Accuracy

Ab initio calculation of electronic excitation energies with chemical accuracy (ca. 1 kcal·mol-1 or 0.043 eV with respect to experiment) is a long-standing challenge in electronic structure theory. Indeed, the most advanced theories can, in practice, only be used to estimate vertical transition energies that cannot be measured experimentally, whereas the calculation of 0-0 energies requires excited-state structures and vibrations for both the ground and excited states, which drastically restrains the number of applicable methods. In this Letter, we present a composite computational protocol able to deliver chemically accurate theoretical 0-0 energies, with a mean absolute deviation of 0.018 eV for a set of 35 singlet valence states. Such accuracy, achievable for the valence states of small- and medium-sized molecules only, allows pinpointing questionable experimental assignments with very high confidence and constitutes a step toward quantitative prediction of excited-state properties.

Distributed Gaussian orbitals for the description of electrons in an external potential

In this work, we demonstrate the viability of using distributed Gaussian orbitals as a basis set for the calculation of the properties of electrons subjected to an external potential. We validate our method by studying one-electron systems for which we can compare to exact analytical results. We highlight numerical aspects that require particular care when using a distributedGaussian basis set. In particular, we discuss the optimal choice for the distance between two neighboring Gaussian orbitals. Finally, we show how our approach can be applied to many-electron problems.

Green Functions and Self-Consistency: Insights From the Spherium Model

We report an exhaustive study of the performance of different variants of Green function methods for the spherium model in which two electrons are confined to the surface of a sphere and interact via a genuine long-range Coulomb operator. We show that the spherium model provides a unique paradigm to study electronic correlation effects from the weakly correlated regime to the strongly correlated regime, since the mathematics are simple while the physics is rich. We compare perturbative GW, partially self-consistent GW and second-order Green function (GF2) methods for the computation of ionization potentials, electron affinities, energy gaps, correlation energies as well as singlet and triplet neutral excitations by solving the Bethe-Salpeter equation (BSE). We discuss the problem of self-screening in GW and show that it can be partially solved with a second-order screened exchange correction (SOSEX). We find that, in general, self-consistency deteriorates the results with respect to those obtained within perturbative approaches with a Hartree-Fock starting point. Finally, we unveil an important problem of partial self-consistency in GW: in the weakly correlated regime, it can produce artificial discontinuities in the self-energy caused by satellite resonances with large weights.

Selected configuration interaction dressed by perturbation

Selected configuration interaction (sCI) methods including second-order perturbative corrections provide near full CI (FCI) quality energies with only a small fraction of the determinants of the FCI space. Here, we introduce both a state-specific and a multi-state sCI method based on the configuration interaction using a perturbative selection made iteratively (CIPSI) algorithm. The present method revises the reference (internal) space under the effect of its interaction with the outer space via the construction of an effective Hamiltonian, following the shifted-Bk philosophy of Davidson and co-workers. In particular, the multi-state algorithm removes the storage bottleneck of the effective Hamiltonian via a low-rank factorization of the dressing matrix. Illustrative examples are reported for the state-specific and multi-state versions.