Marco Caricato

Electronic Transition Energies: A Study of the Performance of a Large Range of Single Reference Density Functional and Wave Function Methods on Valence and Rydberg States Compared to Experiment

This work reports a comparison among wave function and DFT single reference methods for vertical electronic transition energy calculations toward singlet states, valence and Rydberg in nature. A series of 11 small organic molecules are used as test cases, where accurate experimental data in gas phase are available. We compared CIS, RPA, CIS(D), EOM-CCSD, and 28 multipurpose density functionals of the type LSDA, GGA, M-GGA, H-GGA, HM-GGA and with separated short and long-range exchange. The list of functionals is obviously not complete, but it spans more than 20 years of DFT development and includes functionals which are commonly used in the computation of a variety of molecular properties. Large differences in the results were found between the various functionals. The aim of this work is therefore to shed some light on the performance of the plethora of functionals available and compare them with some traditional wave function based methods on a molecular property of large interest as the transition energy.

Oscillator Strength: How Does TDDFT Compare to EOM-CCSD?

In this work, we compare a large variety of density functionals against the equation of motion coupled cluster singles and doubles (EOM-CCSD) method for the calculation of oscillator strengths. Valence and Rydberg states are considered for a test set composed of 11 small organic molecules. In our previous work, the same systems and methods were tested against experimental results for the excitation energies. The results from this investigation confirm our previous findings, i.e., that there is a large difference between the functionals. For the oscillator strength, the average best agreement with EOM-CCSD is provided by CAM-B3LYP followed by LC-ωPBE and, to a lesser extent, B3P86 and LC-BLYP.

A variational formulation of the polarizable continuum model

Continuum solvation models are widely used to accurately estimate solvent effects on energy, structural and spectroscopic properties of complex molecular systems. The polarizable continuum model (PCM) is one of the most versatile among the continuum models because of the variety of properties that can be computed and the diversity of methods that can be used to describe the solute from molecular mechanics (MM) to sophisticated quantum mechanical (QM) post-self-consistent field methods or even hybrid QM/MM methods. In this contribution, we present a new formulation of PCM in terms of a free energy functional whose variational parameters include the continuum polarization (represented by the apparent surface charges), the solutes atomic coordinates and-possibly-its electronic density. The problem of finding the optimized geometry of the (polarized) solute, with the corresponding self-consistent reaction field, is recast as the minimization of this free energy functional, simultaneously with respect to all its variables. The numerous potential applications of this variational formulation of PCM are discussed, including simultaneous optimization of solutes geometry and polarization charges and extended Lagrangian dynamics. In particular, we describe in details the simultaneous optimization procedure and we include several numerical examples.

A time-dependent polarizable continuum model: Theory and application

This work presents an extention of the polarizable continuum model to explicitly describe the time-dependent response of the solvent to a change in the solute charge distribution. Starting from an initial situation in which solute and solvent are in equilibrium, we are interested in modeling the time-dependent evolution of the solvent response, and consequently of the solute-solvent interaction, after a perturbation in this equilibrium situation has been switched on. The model introduces an explicit time-dependent treatment of the polarization by means of the linear-response theory. Two strategies are tested to account for this time dependence: the first one employs the Debye model for the dielectric relaxation, which assumes an exponential decay of the solvent polarization; the second one is based on a fitting of the experimental data of the solvent complex dielectric permittivity. The first approach is simpler and possibly less accurate but allows one to write an analytic expression of the equations. By contrast, the second approach is closer to the experimental evidence but it is limited to the availability of experimental data. The model is applied to the ionization process of N,N-dimethyl-aniline in both acetonitrile and water. The nonequilibrium free-energy profile is studied both as a function of the solvent relaxation coordinate and as a function of time. The solvent reorganization energy is evaluated as well.

Electronic excitation energies in solution at equation of motion CCSD level within a state specific polarizable continuum model approach.

We present a study of excitation energies in solution at the equation of motion coupled cluster singles and doubles (EOM-CCSD) level of theory. The solvent effect is introduced with a state specific polarizable continuum model (PCM), where the solute-solvent interaction is specific for the state of interest. Three definitions of the excited state one-particle density matrix (1PDM) are tested in order to gain information for the development of an integrated EOM-CCSD/PCM method. The calculations show the accuracy of this approach for the computation of such property in solution. Solvent shifts between nonpolar and polar solvents are in good agreement with experiment for the test cases. The completely unrelaxed 1PDM is shown to be a balanced choice between computational effort and accuracy for vertical excitation energies, whereas the response of the ground state CCSD amplitudes and of the molecular orbitals is important for other properties, as for instance the dipole moment.

A comparison between state-specific and linear-response formalisms for the calculation of vertical electronic transition energy in solution with the CCSD-PCM method.

The calculation of vertical electronic transition energies of molecular systems in solution with accurate quantum mechanical methods requires the use of approximate and yet reliable models to describe the effect of the solvent on the electronic structure of the solute. The polarizable continuum model (PCM) of solvation represents a computationally efficient way to describe this effect, especially when combined with coupled cluster (CC) methods. Two formalisms are available to compute transition energies within the PCM framework: State-Specific (SS) and Linear-Response (LR). The former provides a more complete account of the solute-solvent polarization in the excited states, while the latter is computationally very efficient (i.e., comparable to gas phase) and transition properties are well defined. In this work, I review the theory for the two formalisms within CC theory with a focus on their computational requirements, and present the first implementation of the LR-PCM formalism with the coupled cluster singles and doubles method (CCSD). Transition energies computed with LR- and SS-CCSD-PCM are presented, as well as a comparison between solvation models in the LR approach. The numerical results show that the two formalisms provide different absolute values of transition energy, but similar relative solvatochromic shifts (from nonpolar to polar solvents). The LR formalism may then be used to explore the solvent effect on multiple states and evaluate transition probabilities, while the SS formalism may be used to refine the description of specific states and for the exploration of excited state potential energy surfaces of solvated systems.

Absorption and Emission Spectra of Solvated Molecules with the EOM-CCSD-PCM Method.

The accurate calculation of transition energies and properties of isolated molecules is not enough for realistic simulations of their absorption and emission spectra in solution. In fact, the solvent influences the solute geometry, electronic structure, and response to external fields, and a proper description of the solvent effect is fundamental. However, the computational cost of including explicit solvent molecules around the solute becomes rather onerous when an accurate method such as the equation of motion coupled cluster singles and doubles (EOM-CCSD) is employed. The polarizable continuum model of solvation (PCM) may provide an efficient alternative to explicit models, since the sampling of solvent configurations is implicit and the solute-solvent mutual polarization is naturally accounted for. In this contribution, the absorption and emission spectra of molecules in solution are modeled through the EOM-CCSD-PCM method. The equilibrium solvation regime is employed for the geometry optimization of the solute molecule in the ground and excited states, while the nonequilibrium solvation regime is employed for vertical transitions. The theory, implementation, and prototypical applications of the method are presented. The numerical tests involve solvents that are particularly challenging for PCM: low-polar and protic polar solvents. Nonetheless, the experimental trends are well reproduced, and the overall agreement with the measured data is remarkable.

Exploring Potential Energy Surfaces of Electronic Excited States in Solution with the EOM-CCSD-PCM Method

The effect of the solvent on the structure of a molecule in an electronic excited state cannot be neglected. However, the computational cost of including explicit solvent molecules around the solute becomes rather onerous when an accurate method such as the equation of motion coupled cluster singles and doubles (EOM-CCSD) is employed. Solvation continuum models like the polarizable continuum model (PCM) provide an efficient alternative to explicit models, since the solvent conformational average is implicit and the solute-solvent mutual polarization is naturally accounted for. In this work, the coupling of EOM-CCSD and PCM in a state specific approach is presented for the evaluation of energy and analytic energy gradients. Also, various approximations are explored to maintain the computational cost comparable to gas phase EOM-CCSD. Numerical examples are used to test the different schemes.

Large Solvation Effect in the Optical Rotatory Dispersion of Norbornenone

The anomalously large chiroptical response of (1R,4R)-norbornenone has been probed under complementary vapor-phase and solution-phase conditions to assess the putative roles of environmental perturbations. Measurements of the specific rotation for isolated (gas-phase) molecules could not be reproduced quantitatively by comprehensive quantum-chemical calculations based on density-functional or coupled-cluster levels of linear-response theory, which suggests that higher-order treatments may be needed to accurately predict such intrinsic behavior. A substantial, yet unexpected, dependence of the dispersive optical activity on the nature (phase) of the surrounding medium has been uncovered, with the venerable Lorentz local-field correction reproducing solvent-mediated trends in rotatory dispersion surprisingly well, whereas more modern polarizable continuum models for implicit solvation performed less satisfactorily.

Energy-Specific Equation-of-Motion Coupled-Cluster Methods for High-Energy Excited States: Application to K-edge X-ray Absorption Spectroscopy

Single-reference techniques based on coupled-cluster (CC) theory, in the forms of linear response (LR) or equation of motion (EOM), are highly accurate and widely used approaches for modeling valence absorption spectra. Unfortunately, these equations with singles and doubles (LR-CCSD and EOM-CCSD) scale as O(N⁶), which may be prohibitively expensive for the study of high-energy excited states using a conventional eigensolver. In this paper, we present an energy-specific non-Hermitian eigensolver that is able to obtain high-energy excited states (e.g., XAS K-edge spectrum) at low computational cost. In addition, we also introduce an improved trial vector for iteratively solving the EOM-CCSD equation with a focus on high-energy eigenstates. The energy-specific EOM-CCSD approach and its low-scaling alternatives are applied to calculations of carbon, nitrogen, oxygen, and sulfur K-edge excitations. The results are compared to other implementations of CCSD for excited states, energy-specific linear response time-dependent density functional theory (TDDFT), and experimental results with multiple statistical metrics are presented and evaluated.