Manuscript
This commit is contained in:
parent
bd8fa5b9af
commit
17318c887a
@ -1,12 +1,16 @@
|
||||
\documentclass[aip,jcp,reprint,noshowkeys,superscriptaddress]{revtex4-1}
|
||||
\usepackage{graphicx,dcolumn,bm,xcolor,microtype,multirow,amscd,amsmath,amssymb,amsfonts,physics,longtable,wrapfig,txfonts}
|
||||
\usepackage[version=4]{mhchem}
|
||||
\usepackage{natbib}
|
||||
\bibliographystyle{achemso}
|
||||
|
||||
%\usepackage{natbib}
|
||||
%\bibliographystyle{achemso}
|
||||
\usepackage[utf8]{inputenc}
|
||||
\usepackage[T1]{fontenc}
|
||||
\usepackage{txfonts}
|
||||
\usepackage{siunitx}
|
||||
\DeclareSIUnit[number-unit-product = {\,}]
|
||||
\cal{cal}
|
||||
\DeclareSIUnit\kcal{\kilo\cal}
|
||||
\newcommand{\kcalmol}{\si{\kcal\per\mole}}
|
||||
|
||||
\usepackage[
|
||||
colorlinks=true,
|
||||
@ -28,7 +32,7 @@
|
||||
\newcommand{\fnm}{\footnotemark}
|
||||
\newcommand{\fnt}{\footnotetext}
|
||||
\newcommand{\tabc}[1]{\multicolumn{1}{c}{#1}}
|
||||
\newcommand{\SI}{\textcolor{blue}{supporting information}}
|
||||
%\newcommand{\SI}{\textcolor{blue}{supporting information}}
|
||||
\newcommand{\QP}{\textsc{quantum package}}
|
||||
\newcommand{\T}[1]{#1^{\intercal}}
|
||||
|
||||
@ -158,7 +162,7 @@
|
||||
\newcommand{\IneV}[1]{#1 eV}
|
||||
\newcommand{\InAU}[1]{#1 a.u.}
|
||||
\newcommand{\InAA}[1]{#1 \AA}
|
||||
\newcommand{\kcal}{kcal/mol}
|
||||
%\newcommand{\kcal}{kcal/mol}
|
||||
|
||||
% orbitals, gaps, etc
|
||||
\newcommand{\eps}{\varepsilon}
|
||||
@ -217,7 +221,22 @@ Write an abstract
|
||||
\label{sec:intro}
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
|
||||
Despite the fact that excited states are involved in ubiquitious processes such as photochemistry, catalysis or in solar cell technology, none of the many methods existing is the reference in providing accurate excited states energies. Indeed, each method has its own flaws and there are so many chemical scenario that can occur, so it is still one of the biggest challenge in theoretical chemistry. Speaking of difficult task, cyclobutadiene (CBD) molecule has been a real challenge for experimental and theoretical chemists for many decades \cite{bally_cyclobutadiene_1980}. Due to his antiaromaticity \cite{noauthor_aromaticity_nodate} and his large angular strain \cite{baeyer_ueber_1885} the CBD molecule presents a high reactivity which made the synthesis of this molecule a particularly difficult exercise. Hückel molecular orbital theory gives a triplet state with square ($D_{4h}$) geometry for the ground state of the CBD,with the two singly occupied frontier orbitals that are degenerated by symmetry. This degeneracy is lifted by the Jahn-Teller effect, meaning by distortion of the molecule (lowering symmetry), and gives a singlet state with rectangular ($D_{2h}$) geometry for the ground state. Indeed, synthetic work from Pettis and co-workers \cite{reeves_further_1969} gives a rectangular geometry to the singlet ground state of CBD and then was confirmed by experimental works \cite{irngartinger_bonding_1983,ermer_three_1983,kreile_uv_1986}. At the ground state structrure ($D_{2h}$), the ${}^1A_g$ state has a weak multi-configurational character because of the well separated frontier orbitals and can be described by single-reference methods. But at the square ($D_{4h}$) geometry, the singlet state (${}^1B_{1g}$) has two singly occupied frontier orbitals that are degenerated so has a two-configurational character and single-reference methods are unreliable to describe it. The singlet ($D_{4h}$) is a transition state in the automerization reaction between the two rectangular structures. The energy barrier for the automerization of the CBD was predicted, experimentally, in the range of 1.6-10 kcal.mol$^{-1}$ \cite{whitman_limits_1982} and multi-reference calculations gave an energy barrier in the range of 6-7 kcal.mol$^{-1}$ \cite{eckert-maksic_automerization_2006}. All the specificities of the CBD molecule make it a real playground for excited-states methods. In the present work we investigate excited states represented in Fig. \ref{fig:CBD} but not only, we also investigate higher states that are not present in Fig. \ref{fig:CBD} due to energy scaling. In those states we have doubly excited states that are known to be challenging to represent for adiabatic time-dependent density functional theory (TD-DFT) and even for state-of-the-art methods like the approximate third-order coupled-cluster (CC3) or equation-of-motion coupled-cluster with singles, doubles and triples (EOM-CCSDT). In order to tackle the problems of multi-configurational character and double excitations we use multi-reference perturbation theory methods like complete active space perturbation theory (CASPT2) and N-electron valence state perturbation theory (NEVPT2), we also use spin-flip formalism established by Krylov in 2001. To briefly introduce the spin-flip idea we can present it like: instead of taking the singlet ground state as reference, the reference configuration is taken as the lowest triplet state. So one can access the singlet ground state and the singlet doubly-excited state via a spin-flip deexcitation and excitation (respectively), the difference of these two excitation energies providing an estimate of the double excitation. Both excited states and automerization barrier of the CBD are studied using a collection of methods.
|
||||
Despite the fact that excited states are involved in ubiquitious processes such as photochemistry, catalysis or in solar cell technology, none of the many methods existing is the reference in providing accurate excitation energies. Indeed, each method has its own flaws and there are so many chemical scenario that can occur, so it is still one of the biggest challenge in theoretical chemistry. Speaking of difficult task, cyclobutadiene (CBD) molecule has been a real challenge for experimental and theoretical chemists for many decades \cite{bally_cyclobutadiene_1980}. Due to his antiaromaticity \cite{noauthor_aromaticity_nodate} and his large angular strain \cite{baeyer_ueber_1885} the CBD molecule presents a high reactivity which made the synthesis of this molecule a particularly difficult exercise. Hückel molecular orbital theory gives a triplet state with square ($D_{4h}$) geometry for the ground state of the CBD,with the two singly occupied frontier orbitals that are degenerated by symmetry. This degeneracy is lifted by the Jahn-Teller effect, meaning by distortion of the molecule (lowering symmetry), and gives a singlet state with rectangular ($D_{2h}$) geometry for the ground state.
|
||||
Indeed, synthetic work from Pettis and co-workers \cite{reeves_further_1969} gives a rectangular geometry to the singlet ground state of CBD and then was confirmed by experimental works \cite{irngartinger_bonding_1983,ermer_three_1983,kreile_uv_1986}.
|
||||
|
||||
At the ground state structrure ($D_{2h}$), the ${}^1A_g$ state has a weak multi-configurational character because of the well separated frontier orbitals and can be described by single-reference methods. But at the square ($D_{4h}$) geometry, the singlet state ${}^1B_{1g}$ has two singly occupied frontier orbitals that are degenerated so has a two-configurational character and single-reference methods are unreliable to describe it. The singlet ($D_{4h}$) is a transition state in the automerization reaction between the two rectangular structures (see Fig.\ref{fig:CBD}). The energy barrier for the automerization of the CBD was predicted, experimentally, in the range of 1.6-10 kcal.mol$^{-1}$ \cite{whitman_limits_1982} and multi-reference calculations gave an energy barrier in the range of 6-7 kcal.mol$^{-1}$ \cite{eckert-maksic_automerization_2006}. All the specificities of the CBD molecule make it a real playground for excited-states methods.
|
||||
|
||||
Excited states of the CBD molecule in both geometries are represented in Fig.\ref{fig:CBD}. Are represented ${}^1A_g$ and $1{}^3B_{1g}$ states for the rectangular geometry and ${}^1B_{1g}$and $1{}^3A_{2g}$ for the square one. Due to energy scaling doubly excited states $1{}^1B_{1g}$ and $2{}^1A_{1g}$ for the $D_{2h}$ and $D_{4h}$ structures, respectively, are not drawn. Doubly excited states are known to be challenging to represent for adiabatic time-dependent density functional theory (TD-DFT) and even for state-of-the-art methods like the approximate third-order coupled-cluster (CC3) or equation-of-motion coupled-cluster with singles, doubles and triples (EOM-CCSDT). (\textcolor{red}{cite papier spin-flip}).
|
||||
|
||||
In order to tackle the problems of multi-configurational character and double excitations several ways are explored. The most evident way that one can think about to describe multiconfigurational and double excitations are multiconfigurational methods. Among these methods, one can find complete active space self-consistent field (CASSCF) (\textcolor{red}{voir papier reference energies for double excitations}), the second perturbation-corrected variant (CASPT2) (\textcolor{red}{voir papier reference energies for double excitations}) and the second-order $n$-electron valence state perturbation theory (NEVPT2). (\textcolor{red}{voir papier reference energies for double excitations}) The exponential scaling of these methods with the size of the active space is the limitation to the application of these ones to big molecules.
|
||||
|
||||
Another way to deal with double excitations is to use high level truncation of the equation-of-motion (EOM) formalism of coupled-cluster (CC) theory. However, to provide a correct description of doubly excited states one have to take into account contributions from the triple excitations in the CC expansion. Again, due to the scaling of CC methods with the number of basis functions the applicability of these methods is limited to small molecules.
|
||||
|
||||
An alternative to multiconfigurational and CC methods is the use of selected CI (SCI) methods for the computation of transition energies for singly and doubly excited states that are known to reach near full CI energies for small molecules. These methods allow to avoid an exponential increase of the size of the CI expansion by retaining the most energetically relevant determinants only, using a second-order energetic criterion to select perturbatively determinants in the FCI space.
|
||||
|
||||
Finally, to describe doubly excited states, one can think of spin-flip formalism established by Krylov in 2001. To briefly introduce the spin-flip idea we can present it like: instead of taking the singlet ground state as reference, the reference configuration is taken as the lowest triplet state. So one can access the singlet ground state and the singlet doubly-excited state via a spin-flip deexcitation and excitation (respectively), the difference of these two excitation energies providing an estimate of the double excitation. Obviously spin-flip methods have their own flaws, especially the spin contamination (\textcolor{red}{voir papier spin-flip}) (i.e., an artificial mixing of electronic states of differents spin multiplicities) due to spin incompleteness of the spin-flip expansion and by spin contamination of the reference configuration. One can adress part of this problem by expansion of the excitation order but with an increase of the computational cost or by complementing the spin-incomplete configuration set with the missing configurations.
|
||||
|
||||
In the present work we investigate ${}^1A_g$, $1{}^3B_{1g}$, $1{}^1B_{1g}$, $2{}^1A_{g}$ and ${}^1B_{1g}$, $1{}^3A_{2g}$, $2{}^1A_{1g}$,$1{}^1B_{2g}$ excited states for the $D_{2h}$ and $D_{4h}$, respectively, geometries. Computational details are reported in \ref{sec:compmet} for SCI (\ref{sec:SCI}), EOM-CC (\ref{sec:CC}), multiconfigurational (\ref{sec:Multi}) and spin-flip (\ref{sec:sf}) methods. Section \ref{sec:res} is devoted to the discussion of our results for excited states \ref{sec:states} and autoisomerization barrier \ref{sec:auto} of the CBD molecule.
|
||||
|
||||
\begin{figure}
|
||||
\includegraphics[width=0.6\linewidth]{figure2.png}
|
||||
@ -226,8 +245,9 @@ Despite the fact that excited states are involved in ubiquitious processes such
|
||||
\end{figure}
|
||||
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
\section{Computational methods}
|
||||
\section{Computational Details}
|
||||
\label{sec:compmet}
|
||||
The system under investigation in this work is the cyclobutadiene (CBD) molecule, rectangular ($D_{2h}$) and square ($D_{4h}$) geometries are considered. The ($D_{2h}$) geometry is obtained at the CC3 level without frozen core using the aug-cc-pVTZ and the ($D_{4h}$) geometry is obtained at the RO-CCSD(T) level using aug-cc-pVTZ again without frozen core. All the calculations are performed using four basis set, the 6-31+G(d) basis and the aug-cc-pVXZ with X$=$D, T, Q. In the following we will use the notation AVXZ for the aug-cc-pVXZ basis, again with X$=$D, T, Q. The $\%T_1$ metric that gives the percentage of single excitation calculated at the CC3 level in \textcolor{red}{DALTON} allows to characterize the various states.Throughout all this work, spin-flip and spin-conserved calculations are performed with a UHF reference.
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
@ -239,26 +259,24 @@ States energies and excitations energies calculations in the SCI framework are p
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
\subsection{Coupled-Cluster}
|
||||
\label{sec:CC}
|
||||
Different flavours of coupled-cluster (CC) calculations are performed using different codes. Indeed, CC theory provides a hierarchy of methods that provide increasingly accurate energies via the increase of the maximum excitation degree of the cluster operator. However, due to the computational exponential scaling with system size we have to use truncated CC methods. The CC with singles and doubles (CCSD), CC with singles, doubles and triples (CCSDT) calculations are achieved with \textcolor{red}{CFOUR}. The calculations in the context of CC response theory or ``approximate'' series (CC3,CC4) are performed with \textcolor{red}{DALTON}. The CC with singles, doubles, triples and quadruples (CCSDTQ) are done with the \textcolor{red}{CFOUR} code.
|
||||
Different flavours of coupled-cluster (CC) calculations are performed using different codes. Indeed, CC theory provides a hierarchy of methods that provide increasingly accurate energies via the increase of the maximum excitation degree of the cluster operator. Without any truncation of the cluster operator one has the full CC (FCC) that is equivalent to the full configuration interaction (FCI) giving the exact energy and wave function of the system for a fixed atomic basis set. However, due to the computational exponential scaling with system size we have to use truncated CC methods. The CC with singles and doubles (CCSD), CC with singles, doubles and triples (CCSDT) calculations are achieved with \textcolor{red}{CFOUR}. The calculations in the context of CC response theory or ``approximate'' series (CC3,CC4) are performed with \textcolor{red}{DALTON}. The CC with singles, doubles, triples and quadruples (CCSDTQ) are done with the \textcolor{red}{CFOUR} code. The CC2, CC3 and CC4 methods can be seen as cheaper approximations of CCSD, CCSDT and CCSDTQ by skipping the most expensive terms and avoiding the storage of higher-excitations amplitudes.
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
\subsection{Multiconfigurational methods}
|
||||
\label{sec:Multi}
|
||||
State-averaged complete-active-space self-consistent field (SA-CASSCF) calculations are performed with \textcolor{red}{MOLPRO}. On top of those, NEVPT2 calculations, both partially contracted (PC) and strongly contracted (SC) scheme are considered. The PC-NEVPT2 is theoretically more accurate to the strongly contracted version due to the larger number of pertubers and greater flexibility. CASPT2 and extended multistate (XMS) CASPT2 are also performed.
|
||||
State-averaged complete-active-space self-consistent field (SA-CASSCF) calculations are performed with \textcolor{red}{MOLPRO}. On top of those, NEVPT2 calculations, both partially contracted (PC) and strongly contracted (SC) scheme are considered. The PC-NEVPT2 is theoretically more accurate to the strongly contracted version due to the larger number of pertubers and greater flexibility. CASPT2 is performed and extended multistate (XMS) CASPT2 for strong mixing between states with same spin and spatial symmetries is also performed.
|
||||
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
\subsection{Spin-Flip}
|
||||
\label{sec:sf}
|
||||
In both structures the CBD has a singlet ground state, for the spin-flip calculations we consider the lowest triplet state as reference. Spin-flip techniques are broadly accessible and here, among them, we explore equation-of-motion coupled-cluster singles and doubles (EOM-CCSD), configuration interaction singles (CIS), algebraic-diagrammatic construction (ADC) scheme and TD-DFT. The standard and extended spin-flip ADC(2) (SF-ADC(2)-s and SF-ADC(2)-x respectively) and SF-ADC(3) are performed using Q-CHEM 5.2.1. Spin-flip TD-DFT calculations are also performed using Q-CHEM 5.2.1. The BLYP, B3LYP, PBE0 and BH\&HLYP functionals are considered, they contain 0\%, 20\%, 25\%, 50\% of exact exchange and are labeled, respectively, as SF-BLYP, SF-B3LYP, SF-PBE0, SF-BH\&HLYP. We also have done spin-flip TD-DFT calculations using RSH functionals as: CAM-B3LYP, LC-$\omega$PBE08 and $\omega$B97X-V. The main difference here between these RSHs functionals is the amount of exact-exchange at long-range: 75$\%$ for CAM-B3LYP and 100$\%$ for LC-$\omega$PBE08 and $\omega$B97X-V. To complete the use of spin-flip TD-DFT we also considered GH meta-GGA functional M06-2X, RSH meta-GGA functional M11. EOM-SF-CCSD and EOM-SF-CC(2,3) are also performed with Q-CHEM 5.2.1.
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
\section{Computational Details}
|
||||
\label{sec:compdet}
|
||||
|
||||
The system under investigation in this work is the cyclobutadiene (CBD) molecule, rectangular ($D_{2h}$) and square ($D_{4h}$) geometries are considered. The ($D_{2h}$) geometry is obtained at the CC3 level without frozen core using the aug-cc-pVTZ and the ($D_{4h}$) geometry is obtained at the RO-CCSD(T) level using aug-cc-pVTZ again without frozen core. In both structures the CBD has a singlet ground state, for the spin-flip calculations we consider the lowest triplet state as reference. Spin-flip techniques are broadly accessible and here, among them, we explore equation-of-motion coupled-cluster singles and doubles (EOM-CCSD), configuration interaction singles (CIS), algebraic-diagrammatic construction (ADC) scheme and TD-DFT. The standard and extended spin-flip ADC(2) (SF-ADC(2)-s and SF-ADC(2)-x respectively) and SF-ADC(3) are performed using Q-CHEM 5.2.1. Spin-flip TD-DFT calculations are also performed using Q-CHEM 5.2.1. The BLYP, B3LYP, PBE0 and BH\&HLYP functionals are considered, they contain 0\%, 20\%, 25$\%$, 50\% of exact exchange and are labeled, respectively, as SF-BLYP, SF-B3LYP, SF-PBE0, SF-BH\&HLYP. We also have done spin-flip TD-DFT calculations using RSH functionals as: CAM-B3LYP, LC-$\omega$PBE08 and $\omega$B97X-V. The main difference here between these RSHs functionals is the amount of exact-exchange at long-range: 75$\%$ for CAM-B3LYP and 100$\%$ for LC-$\omega$PBE08 and $\omega$B97X-V. To complete the use of spin-flip TD-DFT we also considered GH meta-GGA functional M06-2X, RSH meta-GGA functional M11 and DH functionals B2PLYP and B2GPPLYP. EOM-SF-CCSD and EOM-SF-CC(2,3) are also performed with Q-CHEM 5.2.1. All the calculations are performed using four basis set, the 6-31+G(d) basis and the aug-cc-pVXZ with X$=$D, T, Q. The $\%T_1$ metric that gives the percentage of single excitation calculated at the CC3 level in \textcolor{red}{DALTON} allows to characterize the various states.Throughout all this work, spin-flip and spin-conserved calculations are performed with a UHF reference.
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
@ -270,7 +288,7 @@ As said in \ref{sec:intro} we study both excited states and automerization barri
|
||||
|
||||
%================================================
|
||||
\subsection{Excited States}
|
||||
|
||||
\label{sec:states}
|
||||
|
||||
%%% TABLE I %%%
|
||||
\begin{squeezetable}
|
||||
@ -285,46 +303,41 @@ As said in \ref{sec:intro} we study both excited states and automerization barri
|
||||
Method & Basis & $1\,{}^3B_{1g}$ & $1\,{}^1B_{1g}$ & $2\,{}^1A_{g}$ \\
|
||||
\hline
|
||||
SF-TD-BLYP & 6-31+G(d) & $1.829$ & $1.926$ & $3.755$ \\
|
||||
& aug-cc-pVDZ & $1.828$ & $1.927$ & $3.586$ \\
|
||||
& aug-cc-pVTZ & $1.825$ & $1.927$ & $3.546$ \\
|
||||
& aug-cc-pVQZ & $1.825$ & $1.927$ & $3.528$ \\[0.1cm]
|
||||
& AVDZ & $1.828$ & $1.927$ & $3.586$ \\
|
||||
& AVTZ & $1.825$ & $1.927$ & $3.546$ \\
|
||||
& AVQZ & $1.825$ & $1.927$ & $3.528$ \\[0.1cm]
|
||||
SF-TD-B3LYP & 6-31+G(d) & $1.706$ & $2.211$ & $3.993$ \\
|
||||
& aug-cc-pVDZ & $1.706$ & $2.204$ & $3.992$ \\
|
||||
& aug-cc-pVTZ & $1.703$ & $2.199$ & $3.988$ \\
|
||||
& aug-cc-pVQZ & $1.703$ & $2.199$ & $3.989$\\[0.1cm]
|
||||
& AVDZ & $1.706$ & $2.204$ & $3.992$ \\
|
||||
& AVTZ & $1.703$ & $2.199$ & $3.988$ \\
|
||||
& AVQZ & $1.703$ & $2.199$ & $3.989$\\[0.1cm]
|
||||
SF-TD-PBE0 & 6-31+G(d) & $1.687$ & $2.314$ & $4.089$ \\
|
||||
& aug-cc-pVDZ & $1.684$ & $2.301$ & $4.085$ \\
|
||||
& aug-cc-pVTZ & $1.682$ & $2.296$ & $4.081$ \\
|
||||
& aug-cc-pVQZ & $1.682$ & $2.296$ & $4.079$\\[0.1cm]
|
||||
& AVDZ & $1.684$ & $2.301$ & $4.085$ \\
|
||||
& AVTZ & $1.682$ & $2.296$ & $4.081$ \\
|
||||
& AVQZ & $1.682$ & $2.296$ & $4.079$\\[0.1cm]
|
||||
SF-TD-BH\&HLYP & 6-31+G(d) & $1.552$ & $2.779$ & $4.428$ \\
|
||||
& aug-cc-pVDZ & $1.546$ & $2.744$ & $4.422$ \\
|
||||
& aug-cc-pVTZ & $1.540$ & $2.732$ & $4.492$ \\
|
||||
& aug-cc-pVQZ & $1.540$ & $2.732$ & $4.415$ \\[0.1cm]
|
||||
& AVDZ & $1.546$ & $2.744$ & $4.422$ \\
|
||||
& AVTZ & $1.540$ & $2.732$ & $4.492$ \\
|
||||
& AVQZ & $1.540$ & $2.732$ & $4.415$ \\[0.1cm]
|
||||
SF-TD-M06-2X & 6-31+G(d) & $1.477$ & $2.835$ & $4.378$ \\
|
||||
& aug-cc-pVDZ & $1.467$ & $2.785$ & $4.360$ \\
|
||||
& aug-cc-pVTZ & $1.462$ & $2.771$ & $4.357$ \\
|
||||
& aug-cc-pVQZ & $1.458$ & $2.771$ & $4.352$ \\[0.1cm]
|
||||
& AVDZ & $1.467$ & $2.785$ & $4.360$ \\
|
||||
& AVTZ & $1.462$ & $2.771$ & $4.357$ \\
|
||||
& AVQZ & $1.458$ & $2.771$ & $4.352$ \\[0.1cm]
|
||||
SF-TD-CAM-B3LYP & 6-31+G(d) & $1.750$ & $2.337$ & $3.315$ \\
|
||||
& aug-cc-pVDZ & $1.745$ & $2.323$ & $4.140$ \\
|
||||
& aug-cc-pVTZ & $1.742$ & $2.318$ & $4.138$ \\
|
||||
& aug-cc-pVQZ & $1.743$ & $2.319$ & $4.138$ \\[0.1cm]
|
||||
& AVDZ & $1.745$ & $2.323$ & $4.140$ \\
|
||||
& AVTZ & $1.742$ & $2.318$ & $4.138$ \\
|
||||
& AVQZ & $1.743$ & $2.319$ & $4.138$ \\[0.1cm]
|
||||
SF-TD-$\omega$B97X-V & 6-31+G(d) & $1.810$ & $2.377$ & $4.220$ \\
|
||||
& aug-cc-pVDZ & $1.800$ & $2.356$ & $4.217$ \\
|
||||
& aug-cc-pVTZ & $1.797$ & $2.351$ & $4.213$ \\[0.1cm]
|
||||
& AVDZ & $1.800$ & $2.356$ & $4.217$ \\
|
||||
& AVTZ & $1.797$ & $2.351$ & $4.213$ \\[0.1cm]
|
||||
SF-TD-M11 & 6-31+G(d) & $1.566$ & $2.687$ & $4.292$ \\
|
||||
& aug-cc-pVDZ & $1.546$ & $2.640$ & $4.267$ \\
|
||||
& aug-cc-pVTZ & $1.559$ & $2.651$ & $4.300$ \\
|
||||
& aug-cc-pVQZ & $1.557$ & $2.650$ & $4.299$ \\[0.1cm]
|
||||
& AVDZ & $1.546$ & $2.640$ & $4.267$ \\
|
||||
& AVTZ & $1.559$ & $2.651$ & $4.300$ \\
|
||||
& AVQZ & $1.557$ & $2.650$ & $4.299$ \\[0.1cm]
|
||||
SF-TD-LC-$\omega $PBE08 & 6-31+G(d) & $1.917$ & $2.445$ & $4.353$ \\
|
||||
& aug-cc-pVDZ & $1.897$ & $2.415$ & $4.346$ \\
|
||||
& aug-cc-pVTZ & $1.897$ & $2.415$ & $4.348$ \\
|
||||
& aug-cc-pVQZ & $1.897$ & $2.415$ & $4.348$ \\[0.1cm]
|
||||
SF-TD-B2PLYP & 6-31+G(d) & $1.538$ & $2.827$ & $4.462$ \\
|
||||
& aug-cc-pVDZ & $1.531$ & $2.788$ & $4.455$ \\
|
||||
& aug-cc-pVTZ & $1.525$ & $2.776$ & $4.448$ \\[0.1cm]
|
||||
SF-TD-B2GPPLYP & 6-31+G(d) & \\
|
||||
& aug-cc-pVDZ & \\
|
||||
& aug-cc-pVTZ & \\
|
||||
& AVDZ & $1.897$ & $2.415$ & $4.346$ \\
|
||||
& AVTZ & $1.897$ & $2.415$ & $4.348$ \\
|
||||
& AVQZ & $1.897$ & $2.415$ & $4.348$ \\[0.1cm]
|
||||
|
||||
\end{tabular}
|
||||
\end{ruledtabular}
|
||||
\end{table}
|
||||
@ -344,24 +357,24 @@ SF-TD-B2PLYP & 6-31+G(d) & $1.538$ & $2.827$ & $4.462$ \\
|
||||
Method & Basis & $1\,{}^3B_{1g}$ & $1\,{}^1B_{1g}$ & $2\,{}^1A_{g}$ \\
|
||||
\hline
|
||||
SF-CIS & 6-31+G(d) & $1.514$ & $3.854$ & $5.379$ \\
|
||||
& aug-cc-pVDZ & $1.487$ & $3.721$ & $5.348$ \\
|
||||
& aug-cc-pVTZ & $1.472$ & $3.701$ & $5.342$ \\
|
||||
& aug-cc-pVQZ & $1.471$ & $3.702$ & $5.342$ \\[0.1cm]
|
||||
& AVDZ & $1.487$ & $3.721$ & $5.348$ \\
|
||||
& AVTZ & $1.472$ & $3.701$ & $5.342$ \\
|
||||
& AVQZ & $1.471$ & $3.702$ & $5.342$ \\[0.1cm]
|
||||
SF-ADC(2)-s & 6-31+G(d) & $1.577$ & $3.303$ & $4.196$ \\
|
||||
& aug-cc-pVDZ & $1.513$ & $3.116$ & $4.114$ \\
|
||||
& aug-cc-pVTZ & $1.531$ & $3.099$ & $4.131$ \\
|
||||
& aug-cc-pVQZ & $1.544$ & $3.101$ & $4.140$ \\[0.1cm]
|
||||
& AVDZ & $1.513$ & $3.116$ & $4.114$ \\
|
||||
& AVTZ & $1.531$ & $3.099$ & $4.131$ \\
|
||||
& AVQZ & $1.544$ & $3.101$ & $4.140$ \\[0.1cm]
|
||||
SF-ADC(2)-x & 6-31+G(d) & $1.557$ & $3.232$ & $3.728$ \\
|
||||
& aug-cc-pVDZ & $1.524$ & $3.039$ & $3.681$ \\
|
||||
& aug-cc-pVTZ & $1.539$ & $3.031$ & $3.703$ \\[0.1cm]
|
||||
& AVDZ & $1.524$ & $3.039$ & $3.681$ \\
|
||||
& AVTZ & $1.539$ & $3.031$ & $3.703$ \\[0.1cm]
|
||||
SF-ADC(3) & 6-31+G(d) & $1.435$ & $3.352$ & $4.242$ \\
|
||||
& aug-cc-pVDZ & $1.422$ & $3.180$ & $4.208$ \\
|
||||
& aug-cc-pVTZ & $1.419$ & $3.162$ & $4.224$ \\[0.1cm]
|
||||
& AVDZ & $1.422$ & $3.180$ & $4.208$ \\
|
||||
& AVTZ & $1.419$ & $3.162$ & $4.224$ \\[0.1cm]
|
||||
SF-EOM-CCSD & 6-31+G(d) & $1.663$ & $3.515$ & $4.275$ \\
|
||||
& aug-cc-pVDZ & $1.611$ & $3.315$ & $3.856$ \\
|
||||
& aug-cc-pVTZ & $1.609$ & $3.293$ & $4.245$ \\[0.1cm]
|
||||
& AVDZ & $1.611$ & $3.315$ & $3.856$ \\
|
||||
& AVTZ & $1.609$ & $3.293$ & $4.245$ \\[0.1cm]
|
||||
SF-EOM-CC(2,3) & 6-31+G(d) & $1.490$ & $3.333$ & $4.061$ \\
|
||||
& aug-cc-pVDZ & $1.464$ & $3.156$ & $4.027$ \\
|
||||
& AVDZ & $1.464$ & $3.156$ & $4.027$ \\
|
||||
\end{tabular}
|
||||
\end{ruledtabular}
|
||||
\end{table}
|
||||
@ -381,63 +394,63 @@ SF-EOM-CC(2,3) & 6-31+G(d) & $1.490$ & $3.333$ & $4.061$ \\
|
||||
Method & Basis & $1\,{}^3B_{1g}$ & $1\,{}^1B_{1g}$ & $2\,{}^1A_{g}$ \\
|
||||
\hline
|
||||
CC3 &6-31+G(d)& $1.42$ & $3.341$ & $4.658$ \\
|
||||
& aug-cc-pVDZ & $1.396$ & $3.158$ & $4.711$ \\
|
||||
& aug-cc-pVTZ & $1.402$ & $3.119$ & $4.777$ \\
|
||||
& aug-cc-pVQZ & $1.409$ & $3.113$ & $4.774$ \\[0.1cm]
|
||||
& AVDZ & $1.396$ & $3.158$ & $4.711$ \\
|
||||
& AVTZ & $1.402$ & $3.119$ & $4.777$ \\
|
||||
& AVQZ & $1.409$ & $3.113$ & $4.774$ \\[0.1cm]
|
||||
CCSDT &6-31+G(d)& $1.442$ & $3.357$ & $4.311$ \\
|
||||
& aug-cc-pVDZ & $1.411$ & $3.175$ & $4.327$ \\
|
||||
& aug-cc-pVTZ & $1.411$ & $3.139$ & $4.429$ \\[0.1cm]
|
||||
& AVDZ & $1.411$ & $3.175$ & $4.327$ \\
|
||||
& AVTZ & $1.411$ & $3.139$ & $4.429$ \\[0.1cm]
|
||||
CC4 &6-31+G(d)& & $3.343$ & $4.067$ \\
|
||||
& aug-cc-pVDZ & & $3.164$ & $4.041$ \\[0.1cm]
|
||||
& AVDZ & & $3.164$ & $4.041$ \\[0.1cm]
|
||||
CCSDTQ &6-31+G(d)& & $3.34$ & $4.073$ \\[0.1cm]
|
||||
SA2-CASSCF(4,4) &6-31+G(d)& $1.662$ & $4.657$ & $4.439$ \\
|
||||
& aug-cc-pVDZ & $1.672$ & $4.563$ & $4.448$ \\
|
||||
& aug-cc-pVTZ & $1.67$ & $4.546$ & $4.441$ \\
|
||||
& aug-cc-pVQZ & $1.671$ & $4.549$ & $4.44$ \\[0.1cm]
|
||||
& AVDZ & $1.672$ & $4.563$ & $4.448$ \\
|
||||
& AVTZ & $1.67$ & $4.546$ & $4.441$ \\
|
||||
& AVQZ & $1.671$ & $4.549$ & $4.44$ \\[0.1cm]
|
||||
CASPT2(4,4) &6-31+G(d)& $1.44$ & $3.162$ & $4.115$ \\
|
||||
& aug-cc-pVDZ & $1.414$ & $2.971$ & $4.068$ \\
|
||||
& aug-cc-pVTZ & $1.412$ & $2.923$ & $4.072$ \\
|
||||
& aug-cc-pVQZ & $1.417$ & $2.911$ & $4.081$ \\[0.1cm]
|
||||
& AVDZ & $1.414$ & $2.971$ & $4.068$ \\
|
||||
& AVTZ & $1.412$ & $2.923$ & $4.072$ \\
|
||||
& AVQZ & $1.417$ & $2.911$ & $4.081$ \\[0.1cm]
|
||||
XMS-CASPT2(4,4) &6-31+G(d)& & & $4.151$ \\
|
||||
& aug-cc-pVDZ & & & $4.105$ \\
|
||||
& aug-cc-pVTZ & & & $4.114$ \\
|
||||
& aug-cc-pVQZ & & & $4.125$ \\[0.1cm]
|
||||
& AVDZ & & & $4.105$ \\
|
||||
& AVTZ & & & $4.114$ \\
|
||||
& AVQZ & & & $4.125$ \\[0.1cm]
|
||||
SC-NEVPT2(4,4) &6-31+G(d)& $1.407$ & $2.707$ & $4.145$ \\
|
||||
& aug-cc-pVDZ & $1.381$ & $2.479$ & $4.109$ \\
|
||||
& aug-cc-pVTZ & $1.379$ & $2.422$ & $4.108$ \\
|
||||
& aug-cc-pVQZ & $1.384$ & $2.408$ & $4.125$ \\[0.1cm]
|
||||
& AVDZ & $1.381$ & $2.479$ & $4.109$ \\
|
||||
& AVTZ & $1.379$ & $2.422$ & $4.108$ \\
|
||||
& AVQZ & $1.384$ & $2.408$ & $4.125$ \\[0.1cm]
|
||||
PC-NEVPT2(4,4) &6-31+G(d)& $1.409$ & $2.652$ & $4.12$ \\
|
||||
& aug-cc-pVDZ & $1.384$ & $2.424$ & $4.084$ \\
|
||||
& aug-cc-pVTZ & $1.382$ & $2.368$ & $4.083$ \\
|
||||
& aug-cc-pVQZ & $1.387$ & $2.353$ & $4.091$ \\[0.1cm]
|
||||
& AVDZ & $1.384$ & $2.424$ & $4.084$ \\
|
||||
& AVTZ & $1.382$ & $2.368$ & $4.083$ \\
|
||||
& AVQZ & $1.387$ & $2.353$ & $4.091$ \\[0.1cm]
|
||||
MRCI(4,4) &6-31+G(d)& $1.564$ & $3.802$ & $4.265$ \\
|
||||
& aug-cc-pVDZ & $1.558$ & $3.67$ & $4.254$ \\
|
||||
& aug-cc-pVTZ & $1.568$ & $3.678$ & $4.27$ \\
|
||||
& aug-cc-pVQZ & $1.574$ & $3.681$ & $4.28$ \\[0.1cm]
|
||||
& AVDZ & $1.558$ & $3.67$ & $4.254$ \\
|
||||
& AVTZ & $1.568$ & $3.678$ & $4.27$ \\
|
||||
& AVQZ & $1.574$ & $3.681$ & $4.28$ \\[0.1cm]
|
||||
SA2-CASSCF(12,12) &6-31+G(d)& $1.675$ & $3.924$ & $4.22$ \\
|
||||
& aug-cc-pVDZ & $1.685$ & $3.856$ & $4.221$ \\
|
||||
& aug-cc-pVTZ & $1.686$ & $3.844$ & $4.217$ \\
|
||||
& aug-cc-pVQZ & $1.687$ & $3.846$ & $4.216$ \\[0.1cm]
|
||||
& AVDZ & $1.685$ & $3.856$ & $4.221$ \\
|
||||
& AVTZ & $1.686$ & $3.844$ & $4.217$ \\
|
||||
& AVQZ & $1.687$ & $3.846$ & $4.216$ \\[0.1cm]
|
||||
CASPT2(12,12) &6-31+G(d)& $1.508$ & $3.407$ & $4.099$ \\
|
||||
& aug-cc-pVDZ & $1.489$ & $3.256$ & $4.044$ \\
|
||||
& aug-cc-pVTZ & $1.48$ & $3.183$ & $4.043$ \\
|
||||
& aug-cc-pVQZ & $1.482$ & $3.163$ & $4.047$ \\[0.1cm]
|
||||
& AVDZ & $1.489$ & $3.256$ & $4.044$ \\
|
||||
& AVTZ & $1.48$ & $3.183$ & $4.043$ \\
|
||||
& AVQZ & $1.482$ & $3.163$ & $4.047$ \\[0.1cm]
|
||||
XMS-CASPT2(12,12) &6-31+G(d)& & & $4.111$ \\
|
||||
& aug-cc-pVDZ & & & $4.056$ \\
|
||||
& aug-cc-pVTZ & & & $4.059$ \\
|
||||
& aug-cc-pVQZ & & & $4.065$ \\[0.1cm]
|
||||
& AVDZ & & & $4.056$ \\
|
||||
& AVTZ & & & $4.059$ \\
|
||||
& AVQZ & & & $4.065$ \\[0.1cm]
|
||||
SC-NEVPT2(12,12) &6-31+G(d)& $1.522$ & $3.409$ & $4.13$ \\
|
||||
& aug-cc-pVDZ & $1.511$ & $3.266$ & $4.093$ \\
|
||||
& aug-cc-pVTZ & $1.501$ & $3.188$ & $4.086$ \\
|
||||
& aug-cc-pVQZ & $1.503$ & $3.167$ & $4.088$ \\[0.1cm]
|
||||
& AVDZ & $1.511$ & $3.266$ & $4.093$ \\
|
||||
& AVTZ & $1.501$ & $3.188$ & $4.086$ \\
|
||||
& AVQZ & $1.503$ & $3.167$ & $4.088$ \\[0.1cm]
|
||||
PC-NEVPT2(12,12) &6-31+G(d)& $1.487$ & $3.296$ & $4.103$ \\
|
||||
& aug-cc-pVDZ & $1.472$ & $3.141$ & $4.064$ \\
|
||||
& aug-cc-pVTZ & $1.462$ & $3.063$ & $4.056$ \\
|
||||
& aug-cc-pVQZ & $1.464$ & $3.043$ & $4.059$ \\[0.1cm]
|
||||
& AVDZ & $1.472$ & $3.141$ & $4.064$ \\
|
||||
& AVTZ & $1.462$ & $3.063$ & $4.056$ \\
|
||||
& AVQZ & $1.464$ & $3.043$ & $4.059$ \\[0.1cm]
|
||||
MRCI(12,12) &6-31+G(d)& & & $4.125$ \\[0.1cm]
|
||||
CIPSI &6-31+G(d)& $1.486\pm 0.005$ & $3.348\pm 0.024$ & $4.084\pm 0.012$ \\
|
||||
& aug-cc-pVDZ & $1.458\pm 0.009$ & $3.187\pm 0.035$ & $4.04\pm 0.04$ \\
|
||||
& aug-cc-pVTZ & $1.461\pm 0.030$ & $3.142\pm 0.035$ & $4.03\pm 0.09$ \\
|
||||
& AVDZ & $1.458\pm 0.009$ & $3.187\pm 0.035$ & $4.04\pm 0.04$ \\
|
||||
& AVTZ & $1.461\pm 0.030$ & $3.142\pm 0.035$ & $4.03\pm 0.09$ \\
|
||||
\end{tabular}
|
||||
\end{ruledtabular}
|
||||
\end{table}
|
||||
@ -457,57 +470,57 @@ CIPSI &6-31+G(d)& $1.486\pm 0.005$ & $3.348\pm 0.024$ & $4.084\pm 0.012$ \\
|
||||
Method & Basis & $1\,{}^3A_{2g}$ & $2\,{}^1A_{1g}$ & $1\,{}^1B_{2g}$ \\
|
||||
\hline
|
||||
CCSD & 6-31+G(d) & $0.148$ & $1.788$ & \\
|
||||
& aug-cc-pVDZ & $0.1$ & $1.65$ & \\
|
||||
& aug-cc-pVTZ & $0.085$ & $1.6$ & \\
|
||||
& aug-cc-pVQZ & $0.084$ & $1.588$ & \\[0.1cm]
|
||||
& AVDZ & $0.1$ & $1.65$ & \\
|
||||
& AVTZ & $0.085$ & $1.6$ & \\
|
||||
& AVQZ & $0.084$ & $1.588$ & \\[0.1cm]
|
||||
CC3 & 6-31+G(d) & & $1.809$ & $2.836$ \\
|
||||
& aug-cc-pVDZ & & $1.695$ & $2.646$ \\
|
||||
& aug-cc-pVTZ & & $1.662$ & $2.72$ \\[0.1cm]
|
||||
& AVDZ & & $1.695$ & $2.646$ \\
|
||||
& AVTZ & & $1.662$ & $2.72$ \\[0.1cm]
|
||||
CCSDT & 6-31+G(d) & $0.21$ & $1.751$ & $2.565$ \\
|
||||
& aug-cc-pVDZ & $0.165$ & $1.659$ & $2.45$ \\
|
||||
& aug-cc-pVTZ & $0.149$ & $1.631$ & $2.537$ \\[0.1cm]
|
||||
& AVDZ & $0.165$ & $1.659$ & $2.45$ \\
|
||||
& AVTZ & $0.149$ & $1.631$ & $2.537$ \\[0.1cm]
|
||||
CC4 & 6-31+G(d) & & $1.604$ & $2.121$ \\
|
||||
& aug-cc-pVDZ & & $1.539$ & $1.934$ \\[0.1cm]
|
||||
& AVDZ & & $1.539$ & $1.934$ \\[0.1cm]
|
||||
CCSDTQ & 6-31+G(d) & $0.205$ & $1.593$ & $2.134$ \\[0.1cm]
|
||||
SA2-CASSCF(4,4) & 6-31+G(d) & $0.447$ & $2.257$ & $3.549$ \\
|
||||
& aug-cc-pVDZ & $0.438$ & $2.24$ & $3.443$ \\
|
||||
& aug-cc-pVTZ & $0.434$ & $2.234$ & $3.424$ \\
|
||||
& aug-cc-pVQZ & $0.435$ & $2.235$ & $3.427$ \\[0.1cm]
|
||||
& AVDZ & $0.438$ & $2.24$ & $3.443$ \\
|
||||
& AVTZ & $0.434$ & $2.234$ & $3.424$ \\
|
||||
& AVQZ & $0.435$ & $2.235$ & $3.427$ \\[0.1cm]
|
||||
CASPT2(4,4) & 6-31+G(d) & $0.176$ & $1.588$ & $1.899$ \\
|
||||
& aug-cc-pVDZ & $0.137$ & $1.54$ & $1.708$ \\
|
||||
& aug-cc-pVTZ & $0.128$ & $1.506$ & $1.635$ \\
|
||||
& aug-cc-pVQZ & $0.128$ & $1.498$ & $1.612$ \\[0.1cm]
|
||||
& AVDZ & $0.137$ & $1.54$ & $1.708$ \\
|
||||
& AVTZ & $0.128$ & $1.506$ & $1.635$ \\
|
||||
& AVQZ & $0.128$ & $1.498$ & $1.612$ \\[0.1cm]
|
||||
SC-NEVPT2(4,4) & 6-31+G(d) & $0.083$ & $1.52$ & $1.38$ \\
|
||||
& aug-cc-pVDZ & $0.037$ & $1.465$ & $1.14$ \\
|
||||
& aug-cc-pVTZ & $0.024$ & $1.428$ & $1.055$ \\
|
||||
& aug-cc-pVQZ & $0.024$ & $1.42$ & $1.03$ \\[0.1cm]
|
||||
& AVDZ & $0.037$ & $1.465$ & $1.14$ \\
|
||||
& AVTZ & $0.024$ & $1.428$ & $1.055$ \\
|
||||
& AVQZ & $0.024$ & $1.42$ & $1.03$ \\[0.1cm]
|
||||
PC-NEVPT2(4,4) & 6-31+G(d) & $0.085$ & $1.496$ & $1.329$ \\
|
||||
& aug-cc-pVDZ & $0.039$ & $1.44$ & $1.088$ \\
|
||||
& aug-cc-pVTZ & $0.026$ & $1.403$ & $1.003$ \\
|
||||
& aug-cc-pVQZ & $0.026$ & $1.395$ & $0.977$ \\[0.1cm]
|
||||
& AVDZ & $0.039$ & $1.44$ & $1.088$ \\
|
||||
& AVTZ & $0.026$ & $1.403$ & $1.003$ \\
|
||||
& AVQZ & $0.026$ & $1.395$ & $0.977$ \\[0.1cm]
|
||||
MRCI(4,4) & 6-31+G(d) & $0.297$ & $1.861$ & $2.571$ \\
|
||||
& aug-cc-pVDZ & $0.273$ & $1.823$ & $2.419$ \\
|
||||
& aug-cc-pVTZ & $0.271$ & $1.824$ & $2.415$ \\
|
||||
& aug-cc-pVQZ & $0.273$ & $1.825$ & $2.413$ \\[0.1cm]
|
||||
& AVDZ & $0.273$ & $1.823$ & $2.419$ \\
|
||||
& AVTZ & $0.271$ & $1.824$ & $2.415$ \\
|
||||
& AVQZ & $0.273$ & $1.825$ & $2.413$ \\[0.1cm]
|
||||
SA2-CASSCF(12,12) & 6-31+G(d) & $0.386$ & $1.974$ & $2.736$ \\
|
||||
& aug-cc-pVDZ & $0.374$ & $1.947$ & $2.649$ \\
|
||||
& aug-cc-pVTZ & $0.37$ & $1.943$ & $2.634$ \\
|
||||
& aug-cc-pVQZ & $0.371$ & $1.945$ & $2.637$ \\[0.1cm]
|
||||
& AVDZ & $0.374$ & $1.947$ & $2.649$ \\
|
||||
& AVTZ & $0.37$ & $1.943$ & $2.634$ \\
|
||||
& AVQZ & $0.371$ & $1.945$ & $2.637$ \\[0.1cm]
|
||||
CASPT2(12,12) & 6-31+G(d) & $0.235$ & $1.635$ & $2.17$ \\
|
||||
& aug-cc-pVDZ & $0.203$ & $1.588$ & $2.015$ \\
|
||||
& aug-cc-pVTZ & $0.183$ & $1.538$ & $1.926$ \\
|
||||
& aug-cc-pVQZ & $0.179$ & $1.522$ & $1.898$ \\[0.1cm]
|
||||
& AVDZ & $0.203$ & $1.588$ & $2.015$ \\
|
||||
& AVTZ & $0.183$ & $1.538$ & $1.926$ \\
|
||||
& AVQZ & $0.179$ & $1.522$ & $1.898$ \\[0.1cm]
|
||||
SC-NEVPT2(12,12) & 6-31+G(d) & $0.218$ & $1.644$ & $2.143$ \\
|
||||
& aug-cc-pVDZ & $0.189$ & $1.6$ & $1.991$ \\
|
||||
& aug-cc-pVTZ & $0.165$ & $1.546$ & $1.892$ \\
|
||||
& aug-cc-pVQZ & $0.16$ & $1.529$ & $1.862$ \\[0.1cm]
|
||||
& AVDZ & $0.189$ & $1.6$ & $1.991$ \\
|
||||
& AVTZ & $0.165$ & $1.546$ & $1.892$ \\
|
||||
& AVQZ & $0.16$ & $1.529$ & $1.862$ \\[0.1cm]
|
||||
PC-NEVPT2(12,12) & 6-31+G(d) & $0.189$ & $1.579$ & $2.02$ \\
|
||||
& aug-cc-pVDZ & $0.156$ & $1.53$ & $1.854$ \\
|
||||
& aug-cc-pVTZ & $0.131$ & $1.476$ & $1.756$ \\
|
||||
& aug-cc-pVQZ & $0.126$ & $1.46$ & $1.727$ \\[0.1cm]
|
||||
& AVDZ & $0.156$ & $1.53$ & $1.854$ \\
|
||||
& AVTZ & $0.131$ & $1.476$ & $1.756$ \\
|
||||
& AVQZ & $0.126$ & $1.46$ & $1.727$ \\[0.1cm]
|
||||
CIPSI & 6-31+G(d) & $0.2010\pm 0.0030$ & $1.602\pm 0.007$ & $2.13\pm 0.04$ \\
|
||||
& aug-cc-pVDZ & $0.1570\pm 0.0030$ & $1.587\pm 0.005$ & $2.102\pm 0.027$ \\
|
||||
& aug-cc-pVTZ & $0.169\pm 0.029$ & $1.63\pm 0.05$ & \\
|
||||
& AVDZ & $0.1570\pm 0.0030$ & $1.587\pm 0.005$ & $2.102\pm 0.027$ \\
|
||||
& AVTZ & $0.169\pm 0.029$ & $1.63\pm 0.05$ & \\
|
||||
\end{tabular}
|
||||
\end{ruledtabular}
|
||||
|
||||
@ -528,50 +541,50 @@ CIPSI & 6-31+G(d) & $0.2010\pm 0.0030$ & $1.602\pm 0.007$ & $2.13\pm 0.04$ \\
|
||||
Method & Basis & $1\,{}^3A_{2g}$ & $2\,{}^1A_{1g}$ & $1\,{}^1B_{2g}$ \\
|
||||
\hline
|
||||
SF-CIS & 6-31+G(d) & $0.355$ & $2.742$ & $3.101$ \\
|
||||
& aug-cc-pVDZ & $0.318$ & $2.593$ & $3.052$ \\
|
||||
& aug-cc-pVTZ & $0.305$ & $2.576$ & $3.053$ \\
|
||||
& aug-cc-pVQZ & $0.306$ & $2.577$ & $3.056$ \\[0.1cm]
|
||||
& AVDZ & $0.318$ & $2.593$ & $3.052$ \\
|
||||
& AVTZ & $0.305$ & $2.576$ & $3.053$ \\
|
||||
& AVQZ & $0.306$ & $2.577$ & $3.056$ \\[0.1cm]
|
||||
SF-TD-B3LYP & 6-31+G(d) & $-0.016$ & $0.487$ & $0.542$ \\
|
||||
& aug-cc-pVDZ & $-0.019$ & $0.477$ & $0.536$ \\
|
||||
& aug-cc-pVTZ & $-0.020$ & $0.472$ & $0.533$ \\
|
||||
& aug-cc-pVQZ & $-0.020$ & $0.473$ & $0.533$ \\[0.1cm]
|
||||
& AVDZ & $-0.019$ & $0.477$ & $0.536$ \\
|
||||
& AVTZ & $-0.020$ & $0.472$ & $0.533$ \\
|
||||
& AVQZ & $-0.020$ & $0.473$ & $0.533$ \\[0.1cm]
|
||||
SF-TD-PBE0 & 6-31+G(d) & $-0.012$ & $0.618$ & $0.689$ \\
|
||||
& aug-cc-pVDZ & $-0.016$ & $0.602$ & $0.680$ \\
|
||||
& aug-cc-pVTZ & $-0.019$ & $0.597$ & $0.677$ \\
|
||||
& aug-cc-pVQZ & $-0.018$ & $0.597$ & $0.677$ \\[0.1cm]
|
||||
& AVDZ & $-0.016$ & $0.602$ & $0.680$ \\
|
||||
& AVTZ & $-0.019$ & $0.597$ & $0.677$ \\
|
||||
& AVQZ & $-0.018$ & $0.597$ & $0.677$ \\[0.1cm]
|
||||
SF-TD-BH\&HLYP& 6-31+G(d) & $0.064$ & $1.305$ & $1.458$ \\
|
||||
& aug-cc-pVDZ & $0.051$ & $1.260$ & $1.437$ \\
|
||||
& aug-cc-pVTZ & $0.045$ & $1.249$ & $1.431$ \\
|
||||
& aug-cc-pVQZ & $0.046$ & $1.250$ & $1.432$ \\[0.1cm]
|
||||
& AVDZ & $0.051$ & $1.260$ & $1.437$ \\
|
||||
& AVTZ & $0.045$ & $1.249$ & $1.431$ \\
|
||||
& AVQZ & $0.046$ & $1.250$ & $1.432$ \\[0.1cm]
|
||||
SF-TD-M06-2X & 6-31+G(d) & $0.102$ & $1.476$ & $1.640$ \\
|
||||
& aug-cc-pVDZ & $0.086$ & $1.419$ & $1.611$ \\
|
||||
& aug-cc-pVTZ & $0.078$ & $1.403$ & $1.602$ \\
|
||||
& aug-cc-pVQZ & $0.079$ & $1.408$ & $1.607$ \\[0.1cm]
|
||||
& AVDZ & $0.086$ & $1.419$ & $1.611$ \\
|
||||
& AVTZ & $0.078$ & $1.403$ & $1.602$ \\
|
||||
& AVQZ & $0.079$ & $1.408$ & $1.607$ \\[0.1cm]
|
||||
SF-TD-CAM-B3LYP & 6-31+G(d) & $0.021$ & $0.603$ & $0.672$ \\
|
||||
& aug-cc-pVDZ & $0.012$ & $0.585$ & $0.666$ \\
|
||||
& aug-cc-pVTZ & $0.010$ & $0.580$ & $0.664$ \\
|
||||
& aug-cc-pVQZ & $0.010$ & $0.580$ & $0.664$ \\[0.1cm]
|
||||
& AVDZ & $0.012$ & $0.585$ & $0.666$ \\
|
||||
& AVTZ & $0.010$ & $0.580$ & $0.664$ \\
|
||||
& AVQZ & $0.010$ & $0.580$ & $0.664$ \\[0.1cm]
|
||||
SF-TD-$\omega $B97X-V & 6-31+G(d) & $0.040$ & $0.600$ & $0.670$ \\
|
||||
& aug-cc-pVDZ & $0.029$ & $0.576$ & $0.664$ \\
|
||||
& aug-cc-pVTZ & $0.026$ & $0.572$ & $0.662$ \\
|
||||
& aug-cc-pVQZ & $0.026$ & $0.572$ & $0.662$ \\[0.1cm]
|
||||
& AVDZ & $0.029$ & $0.576$ & $0.664$ \\
|
||||
& AVTZ & $0.026$ & $0.572$ & $0.662$ \\
|
||||
& AVQZ & $0.026$ & $0.572$ & $0.662$ \\[0.1cm]
|
||||
SF-TD-M11 & 6-31+G(d) & $0.102$ & $1.236$ & $1.374$ \\
|
||||
& aug-cc-pVDZ & $0.087$ & $1.196$ & $1.362$ \\
|
||||
& aug-cc-pVTZ & $0.081$ & $1.188$ & $1.359$ \\
|
||||
& aug-cc-pVQZ & $0.080$ & $1.185$ & $1.357$ \\[0.1cm]
|
||||
& AVDZ & $0.087$ & $1.196$ & $1.362$ \\
|
||||
& AVTZ & $0.081$ & $1.188$ & $1.359$ \\
|
||||
& AVQZ & $0.080$ & $1.185$ & $1.357$ \\[0.1cm]
|
||||
SF-TD-LC-$\omega $PBE08 & 6-31+G(d) & $0.078$ & $0.593$ & $0.663$ \\
|
||||
& aug-cc-pVDZ & $0.060$ & $0.563$ & $0.659$ \\
|
||||
& aug-cc-pVTZ & $0.058$ & $0.561$ & $0.658$ \\
|
||||
& aug-cc-pVQZ & $0.058$ & $0.561$ & $0.659$ \\[0.1cm]
|
||||
& AVDZ & $0.060$ & $0.563$ & $0.659$ \\
|
||||
& AVTZ & $0.058$ & $0.561$ & $0.658$ \\
|
||||
& AVQZ & $0.058$ & $0.561$ & $0.659$ \\[0.1cm]
|
||||
SF-ADC(2)-s & 6-31+G(d) & $0.345$ & $1.760$ & $2.096$ \\
|
||||
& aug-cc-pVDZ & $0.269$ & $1.656$ & $1.894$ \\
|
||||
& aug-cc-pVTZ & $0.256$ & $1.612$ & $1.844$ \\[0.1cm]
|
||||
& AVDZ & $0.269$ & $1.656$ & $1.894$ \\
|
||||
& AVTZ & $0.256$ & $1.612$ & $1.844$ \\[0.1cm]
|
||||
SF-ADC(2)-x & 6-31+G(d) & $0.264$ & $1.181$ & $1.972$ \\
|
||||
& aug-cc-pVDZ & $0.216$ & $1.107$ & $1.760$ \\
|
||||
& aug-cc-pVTZ & $0.212$ & $1.091$ & $1.731$ \\[0.1cm]
|
||||
& AVDZ & $0.216$ & $1.107$ & $1.760$ \\
|
||||
& AVTZ & $0.212$ & $1.091$ & $1.731$ \\[0.1cm]
|
||||
SF-ADC(3) & 6-31+G(d) & $0.123$ & $1.650$ & $2.078$ \\
|
||||
& aug-cc-pVDZ & $0.088$ & $1.571$ & $1.878$ \\
|
||||
& aug-cc-pVTZ & $0.079$ & $1.575$ & $1.853$ \\
|
||||
& AVDZ & $0.088$ & $1.571$ & $1.878$ \\
|
||||
& AVTZ & $0.079$ & $1.575$ & $1.853$ \\
|
||||
\end{tabular}
|
||||
\end{ruledtabular}
|
||||
|
||||
@ -586,17 +599,18 @@ SF-ADC(3) & 6-31+G(d) & $0.123$ & $1.650$ & $2.078$ \\
|
||||
|
||||
%================================================
|
||||
\subsection{Autoisomerization barrier}
|
||||
\label{sec:auto}
|
||||
The autoisomerization barrier for the CBD molecule is defined as the energy difference between the singlet ground state of the square ($D_{4h}$) structure and the singlet ground state of the rectangular ($D_{2h}$) geometry. Results for the calculation of the automerization barrier are shown in Tables \ref{tab:auto_standard} and \ref{tab:auto_spin_flip}. As said in \ref{sec:intro} the range for this barrier is quite large. So again, it is important to define our reference in this work in order to be able to compare our results. Table \ref{tab:auto_standard} gives standard methods results, we can observe a large difference for the autoisomerization barrier between the multi-configurational methods. Indeed, for the CASSCF(12,12) we have a difference of the order of 3 kcal.mol$^{-1}$ with CASPT2(12,12) and NEVPT2(12,12) for all the basis. However, the difference between CASPT2(12,12) and NEVPT2(12,12) is much smaller, of the order of 0.2 kcal.mol$^{-1}$ for all the basis.
|
||||
|
||||
%%% TABLE VI %%%
|
||||
\begin{squeezetable}
|
||||
\begin{table}
|
||||
\caption{Autoisomerization barrier for standard methods in kcal.mol$^{-1}$.}
|
||||
\caption{Autoisomerization barrier for standard methods in \kcalmol.}
|
||||
|
||||
\label{tab:auto_standard}
|
||||
\begin{ruledtabular}
|
||||
\begin{tabular}{llrrrr}
|
||||
Method & 6-31+G(d) & aug-cc-pVDZ& aug-cc-pVTZ & aug-cc-pVQZ\\
|
||||
Method & 6-31+G(d) & AVDZ& AVTZ & AVQZ\\
|
||||
\hline
|
||||
CASSCF(12,12) & $10.19$ & $10.75$ & $11.59$ & $11.62$ \\
|
||||
CASPT2(12,12) & $7.24$ & $7.53$ & $8.51$ & $8.71$ \\
|
||||
@ -619,7 +633,7 @@ CIPSI & $7.91\pm 0.21$ & $8.58\pm 0.14$ & & \\
|
||||
%%% TABLE VII %%%
|
||||
\begin{squeezetable}
|
||||
\begin{table}
|
||||
\caption{Autoisomerization barrier for spin-flip methods in kcal.mol$^{-1}$.}
|
||||
\caption{Autoisomerization barrier for spin-flip methods in \kcalmol.}
|
||||
|
||||
\label{tab:auto_spin_flip}
|
||||
\begin{ruledtabular}
|
||||
@ -635,8 +649,6 @@ SF-TD-M06-2X & $9.34$ & $9.68$ & $10.39$ & $10.40$ \\
|
||||
SF-TD-CAM-B3LYP & $18.21$ & $18.30$ & $18.98$ & $18.97$ \\
|
||||
SF-TD-$\omega$B97X-V & $18.46$ & $18.48$ & $19.14$ & $19.12$ \\
|
||||
SF-TD-M11 & $11.13$ & $10.38$ & $11.28$ & $11.19$ \\
|
||||
SF-TD-B2PLYP & $11.37$ & $11.54$ & $12.19$ & \\
|
||||
SF-TD-B2GPPLYP & $8.86$ & $9.05$ & $9.70$ & \\
|
||||
SF-ADC2-s & $6.69$ & $7.15$ & $8.64$ & $8.85$ \\
|
||||
SF-ADC2-x & $8.66$ & $9.15$ & $10.40$ & \\
|
||||
SF-ADC3 & $8.06$ & $8.76$ & $9.58$ & \\
|
||||
@ -651,7 +663,7 @@ SF-ADC3 & $8.06$ & $8.76$ & $9.58$ & \\
|
||||
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
\section{Conclusion}
|
||||
\label{sec:res}
|
||||
\label{sec:conclusion}
|
||||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
||||
|
||||
|
||||
|
||||
148
Notebooks/CBD.nb
148
Notebooks/CBD.nb
@ -10,10 +10,10 @@
|
||||
NotebookFileLineBreakTest
|
||||
NotebookFileLineBreakTest
|
||||
NotebookDataPosition[ 158, 7]
|
||||
NotebookDataLength[ 385478, 9480]
|
||||
NotebookOptionsPosition[ 379867, 9372]
|
||||
NotebookOutlinePosition[ 380260, 9388]
|
||||
CellTagsIndexPosition[ 380217, 9385]
|
||||
NotebookDataLength[ 386938, 9514]
|
||||
NotebookOptionsPosition[ 381326, 9406]
|
||||
NotebookOutlinePosition[ 381719, 9422]
|
||||
CellTagsIndexPosition[ 381676, 9419]
|
||||
WindowFrame->Normal*)
|
||||
|
||||
(* Beginning of Notebook Content *)
|
||||
@ -37,7 +37,8 @@ Cell[BoxData[{
|
||||
mhchem}\>\"", ",", "\"\<\\\\definecolor{darkgreen}{RGB}{0, 180, 0}\>\""}],
|
||||
"}"}]}]}], "]"}], ";"}]}], "Input",
|
||||
InitializationCell->True,
|
||||
CellLabel->"In[57]:=",ExpressionUUID->"73548dbc-1321-4b7c-9133-54658820615c"]
|
||||
CellLabel->
|
||||
"In[222]:=",ExpressionUUID->"73548dbc-1321-4b7c-9133-54658820615c"]
|
||||
}, Closed]],
|
||||
|
||||
Cell[CellGroupData[{
|
||||
@ -9145,6 +9146,27 @@ Cell[BoxData[
|
||||
RowBox[{"-", "154.35746882"}], "+", "154.37689835"}], ")"}], "*",
|
||||
"627.509180"}], ",", "\"\<\>\""}], "}"}], ",", "\[IndentingNewLine]",
|
||||
"\[IndentingNewLine]",
|
||||
RowBox[{"(*",
|
||||
RowBox[{"SF", "-", "TD", "-", "B2GPPLYP"}], "*)"}],
|
||||
"\[IndentingNewLine]",
|
||||
RowBox[{"{",
|
||||
RowBox[{"\"\<SF-TD-B2GPPLYP\>\"", ",",
|
||||
RowBox[{
|
||||
RowBox[{"(",
|
||||
RowBox[{
|
||||
RowBox[{"-", "154.23096101"}], "+", "154.24508025"}], ")"}], "*",
|
||||
"627.509180"}], ",",
|
||||
RowBox[{
|
||||
RowBox[{"(",
|
||||
RowBox[{
|
||||
RowBox[{"-", "154.24355480"}], "+", "154.25798320"}], ")"}], "*",
|
||||
"627.509180"}], ",",
|
||||
RowBox[{
|
||||
RowBox[{"(",
|
||||
RowBox[{
|
||||
RowBox[{"-", "154.27910591"}], "+", "154.29456080"}], ")"}], "*",
|
||||
"627.509180"}], ",", "\"\<\>\""}], "}"}], ",", "\[IndentingNewLine]",
|
||||
"\[IndentingNewLine]",
|
||||
RowBox[{"(*",
|
||||
RowBox[{"SF", "-", "ADC2", "-", "s"}], "*)"}], "\[IndentingNewLine]",
|
||||
RowBox[{"{",
|
||||
@ -9244,8 +9266,13 @@ Cell[BoxData[
|
||||
3.8296465750794897`*^9}, {3.8296469007826138`*^9, 3.829646931299312*^9}, {
|
||||
3.829647054357782*^9, 3.829647083143505*^9}, {3.829990027448785*^9,
|
||||
3.8299900565963984`*^9}, {3.8299901013780727`*^9,
|
||||
3.8299901026680927`*^9}},ExpressionUUID->"9cdf6e6d-a218-450a-95f8-\
|
||||
e3b0048775a3"],
|
||||
3.8299901026680927`*^9}, {3.8302274662543163`*^9, 3.830227534528288*^9}, {
|
||||
3.8302275899525843`*^9, 3.83022759067994*^9}, {3.830227637511037*^9,
|
||||
3.830227639797657*^9}, {3.830227796192348*^9, 3.830227812945609*^9}, {
|
||||
3.8302278836137953`*^9, 3.830227932282262*^9}, 3.830228057477911*^9,
|
||||
3.8302294159571657`*^9},
|
||||
CellLabel->
|
||||
"In[232]:=",ExpressionUUID->"9cdf6e6d-a218-450a-95f8-e3b0048775a3"],
|
||||
|
||||
Cell[CellGroupData[{
|
||||
|
||||
@ -9268,7 +9295,8 @@ Cell[BoxData[
|
||||
3.827829559473291*^9}, {3.827830123216207*^9, 3.827830160115258*^9}, {
|
||||
3.8278301906043577`*^9, 3.8278303147342863`*^9}, {3.827831100880259*^9,
|
||||
3.8278311034111967`*^9}, {3.828427008356379*^9, 3.828427014104094*^9}},
|
||||
CellLabel->"In[66]:=",ExpressionUUID->"224e51d4-fb80-4d27-8453-1ee6ceb061d4"],
|
||||
CellLabel->
|
||||
"In[233]:=",ExpressionUUID->"224e51d4-fb80-4d27-8453-1ee6ceb061d4"],
|
||||
|
||||
Cell[BoxData[
|
||||
TagBox[
|
||||
@ -9321,7 +9349,11 @@ Cell[BoxData[
|
||||
"11.2805631013730369716`4.052330779232006",
|
||||
"11.1895930955612428193`4.048814293895153"},
|
||||
{"\<\"SF-TD-B2PLYP\"\>", "11.3660675022429717984`4.055610230959512",
|
||||
"11.5435898492662101233`4.062340887613776", "\<\"\"\>", "\<\"\"\>"},
|
||||
"11.5435898492662101233`4.062340887613776",
|
||||
"12.1922084380831883266`4.0860823787605804", "\<\"\"\>"},
|
||||
{"\<\"SF-TD-B2GPPLYP\"\>", "8.859952714631754489`3.94743140407336",
|
||||
"9.0539534527208527237`3.956838257413658",
|
||||
"9.6980853508907571836`3.986686001934359", "\<\"\"\>"},
|
||||
{"\<\"SF-ADC2-s\"\>", "6.6935807565410261333`3.8256585074208167",
|
||||
"7.1483462608356180823`3.8542055812605818",
|
||||
"8.6424470638195618477`3.9366367281089816",
|
||||
@ -9360,10 +9392,12 @@ Cell[BoxData[
|
||||
3.827892204100655*^9, 3.827892577618187*^9, {3.827892825440534*^9,
|
||||
3.8278928432874327`*^9}, 3.827893206303258*^9, 3.8278935359926977`*^9, {
|
||||
3.8284270036463346`*^9, 3.8284270145700912`*^9}, 3.829639409858081*^9, {
|
||||
3.829646924755624*^9, 3.829646936556965*^9}, 3.829647085850748*^9},
|
||||
3.829646924755624*^9, 3.829646936556965*^9}, 3.829647085850748*^9,
|
||||
3.8302274558301*^9, 3.830227643490749*^9, 3.8302278162360373`*^9,
|
||||
3.830227935178979*^9, 3.830228060839126*^9, 3.830229419650457*^9},
|
||||
CellLabel->
|
||||
"Out[66]//TableForm=",ExpressionUUID->"65811671-3e8e-40fe-871b-\
|
||||
5c8799cb80df"]
|
||||
"Out[233]//TableForm=",ExpressionUUID->"c93ef115-5cf8-4b73-b414-\
|
||||
5465f153ebdb"]
|
||||
}, Open ]]
|
||||
}, Open ]]
|
||||
}, Open ]]
|
||||
@ -9389,94 +9423,94 @@ CellTagsIndex->{}
|
||||
Notebook[{
|
||||
Cell[CellGroupData[{
|
||||
Cell[580, 22, 157, 3, 98, "Title",ExpressionUUID->"f9fac6ad-b4b0-4fc4-b662-ff26b4fcb383"],
|
||||
Cell[740, 27, 533, 12, 89, "Input",ExpressionUUID->"73548dbc-1321-4b7c-9133-54658820615c",
|
||||
Cell[740, 27, 537, 13, 89, "Input",ExpressionUUID->"73548dbc-1321-4b7c-9133-54658820615c",
|
||||
InitializationCell->True]
|
||||
}, Closed]],
|
||||
Cell[CellGroupData[{
|
||||
Cell[1310, 44, 208, 4, 72, "Title",ExpressionUUID->"d98052b5-b8bf-49c7-976f-70f5883fce25"],
|
||||
Cell[1314, 45, 208, 4, 72, "Title",ExpressionUUID->"d98052b5-b8bf-49c7-976f-70f5883fce25"],
|
||||
Cell[CellGroupData[{
|
||||
Cell[1543, 52, 302, 8, 68, "Section",ExpressionUUID->"0d5f8085-b423-4bab-b579-a185275ba81f"],
|
||||
Cell[1547, 53, 302, 8, 68, "Section",ExpressionUUID->"0d5f8085-b423-4bab-b579-a185275ba81f"],
|
||||
Cell[CellGroupData[{
|
||||
Cell[1870, 64, 156, 3, 54, "Subsection",ExpressionUUID->"010a9483-a730-4fa6-9b2a-13873895483d"],
|
||||
Cell[2029, 69, 422, 11, 33, "Text",ExpressionUUID->"203dfa49-b391-422b-90c4-2b6a46c41acd"],
|
||||
Cell[2454, 82, 618, 12, 196, "Text",ExpressionUUID->"cd19f4a5-5b55-4f6b-ab28-a67e3893faaf"]
|
||||
Cell[1874, 65, 156, 3, 54, "Subsection",ExpressionUUID->"010a9483-a730-4fa6-9b2a-13873895483d"],
|
||||
Cell[2033, 70, 422, 11, 33, "Text",ExpressionUUID->"203dfa49-b391-422b-90c4-2b6a46c41acd"],
|
||||
Cell[2458, 83, 618, 12, 196, "Text",ExpressionUUID->"cd19f4a5-5b55-4f6b-ab28-a67e3893faaf"]
|
||||
}, Closed]],
|
||||
Cell[CellGroupData[{
|
||||
Cell[3109, 99, 152, 3, 38, "Subsection",ExpressionUUID->"146f0b45-c767-44b4-adf3-b56aaf7be25a"],
|
||||
Cell[3113, 100, 152, 3, 38, "Subsection",ExpressionUUID->"146f0b45-c767-44b4-adf3-b56aaf7be25a"],
|
||||
Cell[CellGroupData[{
|
||||
Cell[3286, 106, 167, 3, 45, "Subsubsection",ExpressionUUID->"85eac05c-0c0d-4570-9562-016c59c123c0"],
|
||||
Cell[3456, 111, 447, 12, 30, "Input",ExpressionUUID->"5e934c70-0bb3-4c28-a094-c1ba13c9aa1b"],
|
||||
Cell[3906, 125, 16123, 362, 1627, "Input",ExpressionUUID->"283f461f-fa73-44be-9b2e-3fccbe85020e"]
|
||||
Cell[3290, 107, 167, 3, 45, "Subsubsection",ExpressionUUID->"85eac05c-0c0d-4570-9562-016c59c123c0"],
|
||||
Cell[3460, 112, 447, 12, 30, "Input",ExpressionUUID->"5e934c70-0bb3-4c28-a094-c1ba13c9aa1b"],
|
||||
Cell[3910, 126, 16123, 362, 1627, "Input",ExpressionUUID->"283f461f-fa73-44be-9b2e-3fccbe85020e"]
|
||||
}, Open ]],
|
||||
Cell[CellGroupData[{
|
||||
Cell[20066, 492, 168, 3, 45, "Subsubsection",ExpressionUUID->"496a2a44-bf08-47c1-928f-5e9378cdcbe6"],
|
||||
Cell[20237, 497, 34689, 895, 1606, "Input",ExpressionUUID->"1f47c0ab-65f0-4e3d-b3bf-78a937e27006"]
|
||||
Cell[20070, 493, 168, 3, 45, "Subsubsection",ExpressionUUID->"496a2a44-bf08-47c1-928f-5e9378cdcbe6"],
|
||||
Cell[20241, 498, 34689, 895, 1606, "Input",ExpressionUUID->"1f47c0ab-65f0-4e3d-b3bf-78a937e27006"]
|
||||
}, Open ]]
|
||||
}, Open ]],
|
||||
Cell[CellGroupData[{
|
||||
Cell[54975, 1398, 271, 4, 54, "Subsection",ExpressionUUID->"e9b6b4b3-0323-418a-aaa5-63000ac657a3"],
|
||||
Cell[54979, 1399, 271, 4, 54, "Subsection",ExpressionUUID->"e9b6b4b3-0323-418a-aaa5-63000ac657a3"],
|
||||
Cell[CellGroupData[{
|
||||
Cell[55271, 1406, 167, 3, 45, "Subsubsection",ExpressionUUID->"c2a68803-eb09-4af4-95e0-ded4fd26bba2"],
|
||||
Cell[55275, 1407, 167, 3, 45, "Subsubsection",ExpressionUUID->"c2a68803-eb09-4af4-95e0-ded4fd26bba2"],
|
||||
Cell[CellGroupData[{
|
||||
Cell[55463, 1413, 3016, 47, 56, "Input",ExpressionUUID->"517fe565-dbd4-4f72-9bce-30280a372548"],
|
||||
Cell[58482, 1462, 7766, 161, 1056, "Output",ExpressionUUID->"2be1012c-883f-42eb-b9e3-ba443f8e5e1d"]
|
||||
Cell[55467, 1414, 3016, 47, 56, "Input",ExpressionUUID->"517fe565-dbd4-4f72-9bce-30280a372548"],
|
||||
Cell[58486, 1463, 7766, 161, 1056, "Output",ExpressionUUID->"2be1012c-883f-42eb-b9e3-ba443f8e5e1d"]
|
||||
}, Open ]]
|
||||
}, Closed]],
|
||||
Cell[CellGroupData[{
|
||||
Cell[66297, 1629, 168, 3, 37, "Subsubsection",ExpressionUUID->"e347b647-085c-4a7a-aacb-16cb902d6b69"],
|
||||
Cell[66301, 1630, 168, 3, 37, "Subsubsection",ExpressionUUID->"e347b647-085c-4a7a-aacb-16cb902d6b69"],
|
||||
Cell[CellGroupData[{
|
||||
Cell[66490, 1636, 1559, 31, 77, "Input",ExpressionUUID->"9326eaa0-341b-4e7e-b996-418c884cb254"],
|
||||
Cell[68052, 1669, 109780, 2764, 1850, "Output",ExpressionUUID->"97c1e6af-851b-4074-9392-0292b07355a3"]
|
||||
Cell[66494, 1637, 1559, 31, 77, "Input",ExpressionUUID->"9326eaa0-341b-4e7e-b996-418c884cb254"],
|
||||
Cell[68056, 1670, 109780, 2764, 1850, "Output",ExpressionUUID->"97c1e6af-851b-4074-9392-0292b07355a3"]
|
||||
}, Open ]]
|
||||
}, Open ]]
|
||||
}, Open ]],
|
||||
Cell[CellGroupData[{
|
||||
Cell[177893, 4440, 152, 3, 54, "Subsection",ExpressionUUID->"a0eea47c-0908-49f4-87f8-8a7d458525f2"],
|
||||
Cell[178048, 4445, 2808, 53, 115, "Input",ExpressionUUID->"c94fe3f6-bb2b-425b-9436-cceb389f4f54"],
|
||||
Cell[180859, 4500, 8151, 180, 283, "Input",ExpressionUUID->"ce9a37c9-f290-4ff9-8e20-fb2ac069191f"],
|
||||
Cell[189013, 4682, 129, 3, 30, "Input",ExpressionUUID->"afe9eaa6-d491-4b7a-be7d-51d4faa79596"]
|
||||
Cell[177897, 4441, 152, 3, 54, "Subsection",ExpressionUUID->"a0eea47c-0908-49f4-87f8-8a7d458525f2"],
|
||||
Cell[178052, 4446, 2808, 53, 115, "Input",ExpressionUUID->"c94fe3f6-bb2b-425b-9436-cceb389f4f54"],
|
||||
Cell[180863, 4501, 8151, 180, 283, "Input",ExpressionUUID->"ce9a37c9-f290-4ff9-8e20-fb2ac069191f"],
|
||||
Cell[189017, 4683, 129, 3, 30, "Input",ExpressionUUID->"afe9eaa6-d491-4b7a-be7d-51d4faa79596"]
|
||||
}, Closed]]
|
||||
}, Open ]],
|
||||
Cell[CellGroupData[{
|
||||
Cell[189191, 4691, 351, 9, 68, "Section",ExpressionUUID->"031beb2c-3612-4398-bbed-82451a91c66a"],
|
||||
Cell[189195, 4692, 351, 9, 68, "Section",ExpressionUUID->"031beb2c-3612-4398-bbed-82451a91c66a"],
|
||||
Cell[CellGroupData[{
|
||||
Cell[189567, 4704, 156, 3, 54, "Subsection",ExpressionUUID->"e797c147-2261-4ad9-bf30-a876fe039a42"],
|
||||
Cell[189726, 4709, 475, 11, 33, "Text",ExpressionUUID->"adf832dd-e1af-42ec-b58c-cbbbe9c3192b"],
|
||||
Cell[190204, 4722, 480, 12, 196, "Text",ExpressionUUID->"c1f07e0a-8d48-41b9-a295-c059851cd566"]
|
||||
Cell[189571, 4705, 156, 3, 54, "Subsection",ExpressionUUID->"e797c147-2261-4ad9-bf30-a876fe039a42"],
|
||||
Cell[189730, 4710, 475, 11, 33, "Text",ExpressionUUID->"adf832dd-e1af-42ec-b58c-cbbbe9c3192b"],
|
||||
Cell[190208, 4723, 480, 12, 196, "Text",ExpressionUUID->"c1f07e0a-8d48-41b9-a295-c059851cd566"]
|
||||
}, Closed]],
|
||||
Cell[CellGroupData[{
|
||||
Cell[190721, 4739, 152, 3, 38, "Subsection",ExpressionUUID->"0b0b3775-588f-4d6e-8521-a48a2543bc88"],
|
||||
Cell[190876, 4744, 15753, 350, 1522, "Input",ExpressionUUID->"0fe9f196-cce7-4233-9666-614981e2b593"],
|
||||
Cell[206632, 5096, 31172, 807, 1732, "Input",ExpressionUUID->"86ef0a50-daa3-492e-bff4-0c102d15e915"]
|
||||
Cell[190725, 4740, 152, 3, 38, "Subsection",ExpressionUUID->"0b0b3775-588f-4d6e-8521-a48a2543bc88"],
|
||||
Cell[190880, 4745, 15753, 350, 1522, "Input",ExpressionUUID->"0fe9f196-cce7-4233-9666-614981e2b593"],
|
||||
Cell[206636, 5097, 31172, 807, 1732, "Input",ExpressionUUID->"86ef0a50-daa3-492e-bff4-0c102d15e915"]
|
||||
}, Open ]],
|
||||
Cell[CellGroupData[{
|
||||
Cell[237841, 5908, 271, 4, 54, "Subsection",ExpressionUUID->"29d9a581-07ec-4f13-a8a0-7566f2d2f79d"],
|
||||
Cell[237845, 5909, 271, 4, 54, "Subsection",ExpressionUUID->"29d9a581-07ec-4f13-a8a0-7566f2d2f79d"],
|
||||
Cell[CellGroupData[{
|
||||
Cell[238137, 5916, 1122, 23, 56, "Input",ExpressionUUID->"e8d3b26f-3e58-4a8e-a660-5c2cf37c7202"],
|
||||
Cell[239262, 5941, 94254, 2369, 1596, "Output",ExpressionUUID->"6f94d4c2-5a26-49ac-aa70-234050744f9f"]
|
||||
Cell[238141, 5917, 1122, 23, 56, "Input",ExpressionUUID->"e8d3b26f-3e58-4a8e-a660-5c2cf37c7202"],
|
||||
Cell[239266, 5942, 94254, 2369, 1596, "Output",ExpressionUUID->"6f94d4c2-5a26-49ac-aa70-234050744f9f"]
|
||||
}, Open ]],
|
||||
Cell[CellGroupData[{
|
||||
Cell[333553, 8315, 3156, 50, 56, "Input",ExpressionUUID->"a64c1c75-5181-46a3-916b-6c0d8d53cd2a"],
|
||||
Cell[336712, 8367, 13306, 291, 966, "Output",ExpressionUUID->"e2124d33-2d89-4e1c-8d3a-3faa59b02399"]
|
||||
Cell[333557, 8316, 3156, 50, 56, "Input",ExpressionUUID->"a64c1c75-5181-46a3-916b-6c0d8d53cd2a"],
|
||||
Cell[336716, 8368, 13306, 291, 966, "Output",ExpressionUUID->"e2124d33-2d89-4e1c-8d3a-3faa59b02399"]
|
||||
}, Open ]],
|
||||
Cell[350033, 8661, 147, 3, 52, "Input",ExpressionUUID->"91058d2c-e268-4c8b-8ff4-d01b8a02c12e"]
|
||||
Cell[350037, 8662, 147, 3, 52, "Input",ExpressionUUID->"91058d2c-e268-4c8b-8ff4-d01b8a02c12e"]
|
||||
}, Open ]]
|
||||
}, Open ]],
|
||||
Cell[CellGroupData[{
|
||||
Cell[350229, 8670, 164, 3, 67, "Section",ExpressionUUID->"5e574607-3b56-4320-867c-484c37c5330b"],
|
||||
Cell[350233, 8671, 164, 3, 67, "Section",ExpressionUUID->"5e574607-3b56-4320-867c-484c37c5330b"],
|
||||
Cell[CellGroupData[{
|
||||
Cell[350418, 8677, 152, 3, 54, "Subsection",ExpressionUUID->"1c76f8e7-0977-47de-b4f9-95febd6ffe98"],
|
||||
Cell[350573, 8682, 2921, 65, 241, "Input",ExpressionUUID->"1606e71a-fdd3-4bff-84a4-4cf2d64faef5"],
|
||||
Cell[350422, 8678, 152, 3, 54, "Subsection",ExpressionUUID->"1c76f8e7-0977-47de-b4f9-95febd6ffe98"],
|
||||
Cell[350577, 8683, 2921, 65, 241, "Input",ExpressionUUID->"1606e71a-fdd3-4bff-84a4-4cf2d64faef5"],
|
||||
Cell[CellGroupData[{
|
||||
Cell[353519, 8751, 3030, 47, 52, "Input",ExpressionUUID->"4b815da6-1ca4-49fb-9f5e-99820f6623b2"],
|
||||
Cell[356552, 8800, 3424, 79, 200, "Output",ExpressionUUID->"37fab670-976c-41e6-a633-c9de5012fcd2"]
|
||||
Cell[353523, 8752, 3030, 47, 52, "Input",ExpressionUUID->"4b815da6-1ca4-49fb-9f5e-99820f6623b2"],
|
||||
Cell[356556, 8801, 3424, 79, 200, "Output",ExpressionUUID->"37fab670-976c-41e6-a633-c9de5012fcd2"]
|
||||
}, Open ]],
|
||||
Cell[CellGroupData[{
|
||||
Cell[360013, 8884, 168, 3, 45, "Subsubsection",ExpressionUUID->"39bfcb3c-82e9-40f1-a609-e1fd21238104"],
|
||||
Cell[360184, 8889, 14019, 358, 1144, "Input",ExpressionUUID->"9cdf6e6d-a218-450a-95f8-e3b0048775a3"],
|
||||
Cell[360017, 8885, 168, 3, 45, "Subsubsection",ExpressionUUID->"39bfcb3c-82e9-40f1-a609-e1fd21238104"],
|
||||
Cell[360188, 8890, 15103, 384, 1228, "Input",ExpressionUUID->"9cdf6e6d-a218-450a-95f8-e3b0048775a3"],
|
||||
Cell[CellGroupData[{
|
||||
Cell[374228, 9251, 968, 19, 52, "Input",ExpressionUUID->"224e51d4-fb80-4d27-8453-1ee6ceb061d4"],
|
||||
Cell[375199, 9272, 4604, 93, 266, "Output",ExpressionUUID->"65811671-3e8e-40fe-871b-5c8799cb80df"]
|
||||
Cell[375316, 9278, 972, 20, 52, "Input",ExpressionUUID->"224e51d4-fb80-4d27-8453-1ee6ceb061d4"],
|
||||
Cell[376291, 9300, 4971, 99, 330, "Output",ExpressionUUID->"c93ef115-5cf8-4b73-b414-5465f153ebdb"]
|
||||
}, Open ]]
|
||||
}, Open ]]
|
||||
}, Open ]]
|
||||
|
||||
Loading…
Reference in New Issue
Block a user