PTEROSOR/sfBSE
10
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

reformating Be table

Browse Source
This commit is contained in:
Pierre-Francois Loos 2021-01-10 18:45:45 +01:00
parent 2da59d658e
commit 533c3d7818
2 changed files with 33 additions and 25 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
Manuscript/sfBSE.bib
View File

@ -1,13 +1,21 @@
%% This BibTeX bibliography file was created using BibDesk.
%% http://bibdesk.sourceforge.net/
%% Created for Pierre-Francois Loos at 2021-01-10 15:01:27 +0100
%% Created for Pierre-Francois Loos at 2021-01-10 18:27:42 +0100
%% Saved with string encoding Unicode (UTF-8)
@article{Sarkar_2021,
author = {R. Sarkar and M. Boggio-Pasqua and P. F. Loos and D. Jacquemin},
date-added = {2021-01-10 18:27:30 +0100},
date-modified = {2021-01-10 18:27:30 +0100},
journal = {J. Chem. Theory Comput.},
title = {Benchmark of TD-DFT and Wavefunction Methods for Oscillator Strengths and Excited-State Dipoles},
year = {submitted}}
@article{Krylov_2000b,
author = {Krylov,Anna I.},
date-added = {2021-01-10 15:01:14 +0100},
@ -1136,16 +1144,6 @@
volume = {94},
year = {2006}}
@misc{Sarkar_2020,
archiveprefix = {arXiv},
author = {Rudraditya Sarkar and Martial Boggio-Pasqua and Pierre-Fran{\c c}ois Loos and Denis Jacquemin},
date-added = {2020-12-09 09:59:26 +0100},
date-modified = {2020-12-09 09:59:26 +0100},
eprint = {2011.13233},
primaryclass = {physics.chem-ph},
title = {Benchmarking TD-DFT and Wave Function Methods for Oscillator Strengths and Excited-State Dipole Moments},
year = {2020}}
@article{Scemama_2018a,
author = {A. Scemama and Y. Garniron and M. Caffarel and P. F. Loos},
date-added = {2020-12-09 09:59:26 +0100},

38
Manuscript/sfBSE.tex
View File

@ -65,7 +65,7 @@ This issue can be alleviated by increasing the excitation order at a significant
Nowadays, spin-flip techniques are widely available for many types of methods such as equation-of-motion coupled cluster (EOM-CC), \cite{Krylov_2001a,Levchenko_2004,Manohar_2008,Casanova_2009a,Dutta_2013} configuration interaction (CI), \cite{Krylov_2001b,Krylov_2002,Mato_2018,Casanova_2008,Casanova_2009b} TD-DFT, \cite{Shao_2003,Wang_2004,Li_2011a,Bernard_2012,Zhang_2015} the algebraic-diagrammatic construction (ADC) scheme,\cite{Lefrancois_2015,Lefrancois_2016} and others \cite{Mayhall_2014a,Mayhall_2014b,Bell_2013,Mayhall_2014c} with successful applications in bond breaking processes, \cite{Golubeva_2007} radical chemistry, \cite{Slipchenko_2002,Wang_2005,Slipchenko_2003,Rinkevicius_2010,Ibeji_2015,Hossain_2017,Orms_2018,Luxon_2018} and photochemistry in general \cite{Casanova_2012,Gozem_2013,Nikiforov_2014,Lefrancois_2016} to mention a few.
Here we apply the spin-flip technique to the BSE formalism in order to access, in particular, double excitations. \cite{Authier_2020}
Here we apply the spin-flip technique to the BSE formalism in order to access, in particular, double excitations, \cite{Authier_2020} but not only.
The present BSE calculations are based on the spin unrestricted version of both $GW$ (Sec.~\ref{sec:UGW}) and BSE (Sec.~\ref{sec:UBSE}).
To the best of our knowledge, the present study is the first to apply the spin-flip formalism to the BSE method.
Moreover, we also go beyond the static approximation by taking into account dynamical effects (Sec.~\ref{sec:dBSE}) via an unrestricted generalization of our recently developed (renormalized) perturbative correction which builds on the seminal work of Strinati, \cite{Strinati_1982,Strinati_1984,Strinati_1988} Romaniello and collaborators, \cite{Romaniello_2009b,Sangalli_2011} and Rohlfing and coworkers. \cite{Rohlfing_2000,Ma_2009a,Ma_2009b,Baumeier_2012b,Lettmann_2019}
@ -526,7 +526,7 @@ where
(p_\sig|x|q_\sigp) = \int \MO{p_\sig}(\br) \, x \, \MO{q_\sigp}(\br) d\br
\end{equation}
are one-electron integrals in the orbital basis.
The total oscillator strength is given by
The total oscillator strength in the so-called length gauge \cite{Chrayteh_2021,Sarkar_2021} is given by
\begin{equation}
f_{m}^{\spc} = \frac{2}{3} \Om{m}{\spc} \qty[ \qty(\mu_{x,m}^{\spc})^2 + \qty(\mu_{x,m}^{\spc})^2 + \qty(\mu_{x,m}^{\spc})^2 ]
\end{equation}
@ -592,24 +592,34 @@ The TD-DFT calculations have been performed with Q-CHEM 5.2.1 \cite{qchem4} and
%%% TABLE I %%%
\begin{squeezetable}
\begin{table*}
\begin{table}
\caption{
Spin-flip excitations (in eV) of \ce{Be} obtained for various methods with the 6-31G basis.
The $GW$ calculations are performed with an UHF starting point.
Excitation energies [with respect to the $^1S(1s^2 2s^2)$ singlet ground state] of \ce{Be} obtained for various methods with the 6-31G basis.
\label{tab:Be}}
\begin{ruledtabular}
\begin{tabular}{lcccccccccccc}
State & TD-BLYP & TD-B3LYP & TD-BHHLYP & CIS
& BSE@{\GOWO} & BSE@{\evGW} & dBSE@{\GOWO} & dBSE@{\evGW}
& ADC(2) & ADC(2)-x & ADC(3) & FCI \\
\begin{tabular}{lcccc}
& \mc{4}{c}{Excitation energies (eV)} \\
\cline{2-5}
Method & $^3P(1s^22s2p)$ & $^1P(1s^22s2p)$
& $^3P(1s^22p^2)$ & $^1P(1s^22p^2)$ \\
\hline
$^3P(1s^22s2p)$ & 3.210 & 3.332 & 2.874 & 2.111 & 2.399& 2.407& 2.363& 2.369 & 2.433 & 2.866 & 2.863 & 2.862 \\
$^1P(1s^22s2p)$ & 3.210 & 4.275& 4.922& 6.036 & 6.191& 6.199 & 6.263 & 6.273 & 6.255 &6.581 & 6.579 & 6.577 \\
$^3P(1s^22p^2)$ & 6.691 & 6.864& 7.112 & 7.480 & 7.792 & 7.788 & 7.824 & 7.820 & 7.745 & 7.664 & 7.658 & 7.669 \\
$^1P(1s^22p^2)$ & 7.598 & 7.762& 8.188 & 8.945 &9.373 & 9.388 & 9.424 & 9.441 & 9.047 & 8.612 & 8.618 & 8.624 \\
SF-TD-BLYP & 3.210 & 3.210 & 6.691 & 7.598 \\
SF-TD-B3LYP & 3.332 & 4.275 & 6.864 & 7.762 \\
SF-TD-BHHLYP & 2.874 & 4.922 & 7.112 & 8.188 \\
SF-CIS & 2.111 & 6.036 & 7.480 & 8.945 \\
SF-BSE@{\GOWO}@UHF & 2.399 & 6.191 & 7.792 & 9.373 \\
SF-BSE@{\evGW}@UHF & 2.407 & 6.199 & 7.788 & 9.388 \\
SF-BSE@{\qsGW}@UHF & 2.376 & 6.241 & 7.668 & 9.417 \\
SF-dBSE@{\GOWO}@UHF & 2.363 & 6.263 & 7.824 & 9.424 \\
SF-dBSE@{\evGW}@UHF & 2.369 & 6.273 & 7.820 & 9.441 \\
SF-dBSE@{\qsGW}@UHF & 2.335 & 6.317 & 7.689 & 9.470 \\
SF-ADC(2)-s & 2.433 & 6.255 & 7.745 & 9.047 \\
SF-ADC(2)-x & 2.866 & 6.581 & 7.664 & 8.612 \\
SF-ADC(3) & 2.863 & 6.579 & 7.658 & 8.618 \\
FCI & 2.862 & 6.577 & 7.669 & 8.624 \\
\end{tabular}
\end{ruledtabular}
\end{table*}
\end{table}
\end{squeezetable}
%%% %%% %%% %%%