Comment on “Anab initiocluster study of the structure of the Si(001) surface” [J. Chem. Phys.112, 2994 (2000)]

On the structure of Si(100) surface: importance of higher order correlations for buckled dimer

We revisit a dangling theoretical question of whether the surface reconstruction of the Si(100) surface would energetically favor the symmetric or buckled dimers on the intrinsic potential energy surfaces at 0 K. This seemingly simple question is still unanswered definitively since all existing density functional based calculations predict the dimers to be buckled, while most wavefunction based correlated treatments prefer the symmetric configurations. Here, we use the doubly hybrid density functional (DHDF) geometry optimizations, in particular, XYGJ-OS, complete active space self-consistent field theory, multi-reference perturbation theory, multi-reference configuration interaction (MRCI), MRCI with the Davidson correction (MRCI + Q), multi-reference average quadratic CC (MRAQCC), and multi-reference average coupled pair functional (MRACPF) methods to address this question. The symmetric dimers are still shown to be lower in energy than the buckled dimers when using the CASPT2 method on the DHDF optimized geometries, consistent with the previous results using B3LYP geometries [Y. Jung, Y. Shao, M. S. Gordon, D. J. Doren, and M. Head-Gordon, J. Chem. Phys. 119, 10917 (2003)]. Interestingly, however, the MRCI + Q, MRAQCC, and MRACPF results (which give a more refined description of electron correlation effects) suggest that the buckled dimer is marginally more stable than its symmetric counterpart. The present study underlines the significance of having an accurate description of the electron-electron correlation as well as proper multi-reference wave functions when exploring the extremely delicate potential energy surfaces of the reconstructed Si(100) surface.

The structure of the Si9H12 cluster: A coupled cluster and multi-reference perturbation theory study

Full geometry optimizations using both singles and doubles coupled cluster theory with perturbative triple excitations, CCSD(T), and second order multi-reference perturbation theory, MRMP2, have been employed to predict the structure of Si9H12, a cluster commonly used in calculations to represent the Si(100) surface. Both levels of theory predict the structure of this cluster to be symmetric (not buckled), and no evidence for a buckled (asymmetric) structure is found at either level of theory.

Electronic Excited States of Si(100) and Organic Molecules Adsorbed on Si(100)

The electronically excited states of the Si(100) surface and acetylene, benzene, and 9,10-phenanthrenequinone adsorbed on Si(100) are studied with time-dependent density functional theory. The computational cost of these calculations can be reduced through truncation of the single excitation space. This allows larger cluster models of the surface in conjunction with large adsorbates to be studied. On clean Si(100), the low-lying excitations correspond to transitions between the pi orbitals of the silicon-silicon dimers. These excitations are predicted to occur in the range 0.4-2 eV. When organic molecules are adsorbed on the surface, surface --> molecule, molecule --> surface, and electronic excitations localized within the adsorbate are also observed at higher energies. For acetylene and benzene, the remaining pipi* excitations are found to lie at lower energies than in the corresponding gas-phase species. Even though the aromaticity of 9,10-phenanthrenequinone is retained, significant shifts in the pipi* excitations of the aromatic rings are predicted. This is in part due to structural changes that occur upon adsorption.

Density-functional study of the cycloaddition of acrylonitrile on the Si(100) surface

Using a density functional approach, we have explored the cycloaddition of acrylonitrile on the Si(100) surface. The buckling of the surface dimers characteristic for the (2x1) reconstructed surface is shown to favor structures with a dipolar moment such as the resonant form of acrylonitrile with cumulative double bonds. The bond of acrylonitrile via a single C atom is a possible intermediate leading to the nitrile structure of the adsorbed molecule.

Role of electronic correlation in the Si(100) reconstruction: a quantum Monte Carlo study

Recent low-temperature scanning tunneling experiments have questioned the generally accepted picture of buckled silicon dimers as the ground state reconstruction of the Si(100) surface, undermining the ability of density functional theory to accurately describe electronic correlations at surfaces. We present quantum Monte Carlo calculations on large cluster models of the surface, and conclude that buckling remains energetically favorable even when the present-day best treatment of electronic correlation is employed. The implications for experimental interpretation are discussed.