標題: | 小型過渡金屬簇之精確電子結構計算 Accurate electronic structure calculations on small transition metal clusters |
作者: | 凱富 Cristopher Camacho Leandro 魏恆理 Henryk Witek 應用化學系碩博士班 |
關鍵字: | Perturbation theory;Intruder states;Transition metals;CASSCF;Perturbation theory;Intruder states;Transition metals;CASSCF |
公開日期: | 2010 |
摘要: | The potential energy curve of the ground state of Mn$_2$ has been
studied using a systematic sequence of complete active spaces.
Deficiencies of the routinely used active space, built from atomic $4s$
and $3d$ orbitals, has been identified and discussed. It is shown that
an additional $\sigma_g$ orbital, originating from the atomic virtual
$4p_z$ orbitals, is essential for a proper description of static
correlation in the $^1\Sigma_{g}^{+}$ state of Mn$_2$. The calculated
spectroscopic parameters of the $^1\Sigma_{g}^{+}$ state agree well with
available experimental data. The calculated equilibrium bond lengths
are located between 3.24 and~3.50~{\AA}, the harmonic vibrational
frequencies, between 44 and~72~cm$^{-1}$, and the dissociation energies,
between 0.05 and~0.09~eV.
A detailed analysis of a severe intruder state problem in the multistate
multireference perturbation theory (MS-MRPT) calculations on the ground
state of manganese dimer is presented. An enormous number of detected
intruder states ($>$5000) do not permit finding even an approximate
shape of the $X^1\Sigma_{g}^{+}$ potential energy curve. The intruder
states are explicitly demonstrated to originate from quasidegeneracies
in the zeroth-order Hamiltonian spectrum. The electronic configurations
responsible for appearance of the quasidegeneracies are identified as
single and double excitations from the active orbitals to the external
orbitals. It is shown that the quasidegeneracy problem can be
completely eliminated using shift techniques despite of its severity.
The resultant curves are smooth and continuous. Unfortunately, strong
dependence of the spectroscopic parameters of the $X^1\Sigma_{g}^{+}$
state on the shift parameter is observed. This finding rises serious
controversies regarding validity of employing shift techniques for
solving the intruder state problem in multistate multireference
perturbation theory.
Prediction of a false ground state with popular variants of
multireference perturbation theory (CASPT2 and MRMP) is reported. The
failure occurs for a remarkably simple chemical system: the Sc$_2$
molecule. Reasons for the failure are discussed and appropriate
remedies are suggested. The presented finding has far-reaching
consequences for all the chemical community giving a serious warning on
the applicability of multireference perturbation theory in the presence
of intruder states.
A systematic investigation of low-lying states of Sc$_2$ using
multireference perturbation theory (NEVPT2 and NEVPT3) indicates that
the ground state of this system is $^5\Sigma_u^-$ with
$r_e=2.611$~{\AA}, $\omega_e=241.8$~cm$^{-1}$, and $D_e=1.78$~eV. This
state is closely followed by other low-lying states of Sc$_2$:
$^3\Sigma_u^-$, $^5\Delta_u$, $^3\Pi_g$, $^1\Pi_g$, and $^1\Sigma_u^-$.
Our energy ordering of the $^5\Sigma_u^-$ and $^3\Sigma_u^-$ states
confirms the recent MRCI results of Kalemos \textit{et al.}
[\textit{J.Chem.Phys.} \textbf{132}, 024309 (2010)] and is at variance
with the earlier DMC predictions of Matxain \textit{et al.}
[\textit{J.Chem.Phys.} \textbf{128}, 194315 (2008)]. An excellent
agreement between the second- and third-order NEVPT results and between
the computed and experimental values of $\omega_e$ (241.8 vs.
238.9~cm$^{-1}$) for the $^5\Sigma_u^-$ state suggests high accuracy of
our predictions. The potential energy curve of the ground state of Mn$_2$ has been studied using a systematic sequence of complete active spaces. Deficiencies of the routinely used active space, built from atomic $4s$ and $3d$ orbitals, has been identified and discussed. It is shown that an additional $\sigma_g$ orbital, originating from the atomic virtual $4p_z$ orbitals, is essential for a proper description of static correlation in the $^1\Sigma_{g}^{+}$ state of Mn$_2$. The calculated spectroscopic parameters of the $^1\Sigma_{g}^{+}$ state agree well with available experimental data. The calculated equilibrium bond lengths are located between 3.24 and~3.50~{\AA}, the harmonic vibrational frequencies, between 44 and~72~cm$^{-1}$, and the dissociation energies, between 0.05 and~0.09~eV. A detailed analysis of a severe intruder state problem in the multistate multireference perturbation theory (MS-MRPT) calculations on the ground state of manganese dimer is presented. An enormous number of detected intruder states ($>$5000) do not permit finding even an approximate shape of the $X^1\Sigma_{g}^{+}$ potential energy curve. The intruder states are explicitly demonstrated to originate from quasidegeneracies in the zeroth-order Hamiltonian spectrum. The electronic configurations responsible for appearance of the quasidegeneracies are identified as single and double excitations from the active orbitals to the external orbitals. It is shown that the quasidegeneracy problem can be completely eliminated using shift techniques despite of its severity. The resultant curves are smooth and continuous. Unfortunately, strong dependence of the spectroscopic parameters of the $X^1\Sigma_{g}^{+}$ state on the shift parameter is observed. This finding rises serious controversies regarding validity of employing shift techniques for solving the intruder state problem in multistate multireference perturbation theory. Prediction of a false ground state with popular variants of multireference perturbation theory (CASPT2 and MRMP) is reported. The failure occurs for a remarkably simple chemical system: the Sc$_2$ molecule. Reasons for the failure are discussed and appropriate remedies are suggested. The presented finding has far-reaching consequences for all the chemical community giving a serious warning on the applicability of multireference perturbation theory in the presence of intruder states. A systematic investigation of low-lying states of Sc$_2$ using multireference perturbation theory (NEVPT2 and NEVPT3) indicates that the ground state of this system is $^5\Sigma_u^-$ with $r_e=2.611$~{\AA}, $\omega_e=241.8$~cm$^{-1}$, and $D_e=1.78$~eV. This state is closely followed by other low-lying states of Sc$_2$: $^3\Sigma_u^-$, $^5\Delta_u$, $^3\Pi_g$, $^1\Pi_g$, and $^1\Sigma_u^-$. Our energy ordering of the $^5\Sigma_u^-$ and $^3\Sigma_u^-$ states confirms the recent MRCI results of Kalemos \textit{et al.} [\textit{J.Chem.Phys.} \textbf{132}, 024309 (2010)] and is at variance with the earlier DMC predictions of Matxain \textit{et al.} [\textit{J.Chem.Phys.} \textbf{128}, 194315 (2008)]. An excellent agreement between the second- and third-order NEVPT results and between the computed and experimental values of $\omega_e$ (241.8 vs. 238.9~cm$^{-1}$) for the $^5\Sigma_u^-$ state suggests high accuracy of our predictions. |
URI: | http://140.113.39.130/cdrfb3/record/nctu/#GT079425816 http://hdl.handle.net/11536/40850 |
顯示於類別: | 畢業論文 |