Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBenchmarking TD-DFT and Wave Function Methods for Oscillator Strengths and Excited-State Dipole Moments
77
0
0.0
(
0
)
Added by
Pierre-Fran\\c{c}ois Loos Dr
Publication date
2020
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Using a set of oscillator strengths and excited-state dipole moments of near full configuration interaction (FCI) quality determined for small compounds, we benchmark the performances of several single-reference wave function methods (CC2, CCSD, CC3, CCSDT, ADC(2), and ADC(3/2)) and time-dependent density-functional theory (TD-DFT) with various functionals (B3LYP, PBE0, M06-2X, CAM-B3LYP, and $omega$B97X-D). We consider the impact of various gauges (length, velocity, and mixed) and formalisms: equation of motion (EOM) emph{vs} linear response (LR), relaxed emph{vs} unrelaxed orbitals, etc. Beyond the expected accuracy improvements and a neat decrease of formalism sensitivy when using higher-order wave function methods, the present contribution shows that, for both ADC(2) and CC2, the choice of gauge impacts more significantly the magnitude of the oscillator strengths than the choice of formalism, and that CCSD yields a notable improvement on this transition property as compared to CC2. For the excited-state dipole moments, switching on orbital relaxation appreciably improves the accuracy of both ADC(2) and CC2, but has a rather small effect at the CCSD level. Going from ground to excited states, the typical errors on dipole moments for a given method tend to roughly triple. Interestingly, the ADC(3/2) oscillator strengths and dipoles are significantly more accurate than their ADC(2) counterparts, whereas the two models do deliver rather similar absolute errors for transition energies. Concerning TD-DFT, one finds: i) a rather negligible impact of the gauge on oscillator strengths for all tested functionals (except for M06-2X); ii) deviations of ca.~0.10 D on ground-state dipoles for all functionals; iii) the better overall performance of CAM-B3LYP for the two considered excited-state properties.
rate research
Read More
The molecular dipole moment ($boldsymbol{mu}$) is a central quantity in chemistry. It is essential in predicting infrared and sum-frequency generation spectra, as well as induction and long-range electrostatic interactions. Furthermore, it can be extracted directly from high-level quantum mechanical calculations, making it an ideal target for machine learning (ML). In this work, we choose to represent this quantity with a physically inspired ML model that captures two distinct physical effects: local atomic polarization is captured within the symmetry-adapted Gaussian process regression (SA-GPR) framework, which assigns a (vector) dipole moment to each atom, while movement of charge across the entire molecule is captured by assigning a partial (scalar) charge to each atom. The resulting MuML models are fitted together to reproduce molecular $boldsymbol{mu}$ computed using high-level coupled-cluster theory (CCSD) and density functional theory (DFT) on the QM7b dataset. The combined model shows excellent transferability when applied to a showcase dataset of larger and more complex molecules, approaching the accuracy of DFT at a small fraction of the computational cost. We also demonstrate that the uncertainty in the predictions can be estimated reliably using a calibrated committee model. The ultimate performance of the models depends, however, on the details of the system at hand, with the scalar model being clearly superior when describing large molecules whose dipole is almost entirely generated by charge separation. These observations point to the importance of simultaneously accounting for the local and non-local effects that contribute to $boldsymbol{mu}$; further, they define a challenging task to benchmark future models, particularly those aimed at the description of condensed phases.
We present an improved approach for generating a set of optimized frontier orbitals (HOMO and LUMO) that minimizes the energy of one double configuration. We further benchmark the effect of including such a double within a CIS or TD-DFT configuration interaction Hamiltonian for a set of test cases. We find that, although we cannot achieve quantitative accuracy, the algorithm is quite robust and routinely delivers an enormous qualitative improvement to standard single-reference electronic structure calculations.
Excited state electron and hole transfer underpin fundamental steps in processes such as exciton dissociation at photovoltaic heterojunctions, photoinduced charge transfer at electrodes, and electron transfer in photosynthetic reaction centers. Diabatic states corresponding to charge or excitation localized species, such as locally excited and charge transfer states, provide a physically intuitive framework to simulate and understand these processes. However, obtaining accurate diabatic states and their couplings from adiabatic electronic states generally leads to inaccurate results when combined with low-tier electronic structure methods, such as time dependent density functional theory (TDDFT), and exorbitant computational cost when combined with high-level wavefunction-based methods. Here we introduce a DFT-based diabatization scheme, {Delta}-ALMO(MSDFT2), which directly constructs the diabatic states using absolutely localized molecular orbitals (ALMOs). We demonstrate that our method, which combines ALMO calculations with the {Delta}SCF technique to construct electronically excited diabatic states and obtains their couplings with charge-transfer states using our MSDFT2 scheme, gives accurate results for excited state electron and hole transfer in both charged and uncharged systems that underlie DNA repair, charge separation in donor-acceptor dyads, chromophore-to-solvent electron transfer, and singlet fission. This framework for the accurate and efficient construction of excited state diabats and evaluation of their couplings directly from DFT thus offers a route to simulate and elucidate photoinduced electron and hole transfer in large disordered systems, such as those encountered in the condensed phase.
Electronic structure of HCl+ and HBr+ molecular ions is calculated using the symmetry-adapted-cluster configuration interaction (SAC-CI) method. In this paper, we analyse dipole moments (DM) for a series of low-lying six 2Pi-states and transition dipole moments (TDM for transitions from the ground state X2Pi to the excited 2Pi-series. Behavior of DMs with change of interatomic distances is different for these ions for the excited 2Pi-states in correspondence with different dissociation paths. TDMs reveal the pronounced maxima at the beginning steps of dissociation.
We introduce natural transition geminals as a means to qualitatively understand a transition where double excitations are important. The first two $A_{1}$ singlet states of the CH cation are used as an initial example. We calculate these states with configuration interaction singles (CIS) and state-averaged Monte Carlo configuration interaction (SA-MCCI). For each method we compare the important natural transition geminals with the dominant natural transition orbitals. We then compare SA-MCCI and full configuration interaction (FCI) with regards to the natural transition geminals using the beryllium atom. We compare using the natural transition geminals with analyzing the important configurations in the CI expansion to give the dominant transition for the beryllium atom and the carbon dimer. Finally we calculate the natural transition geminals for two electronic excitations of formamide.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
