Localized helium excitations in 4He_N-benzene clusters

Patrick Huang,K. Birgitta Whaley

DOI: https://doi.org/10.1103/PhysRevB.67.155419

2003-01-28

Abstract:We compute ground and excited state properties of small helium clusters 4He_N containing a single benzene impurity molecule. Ground-state structures and energies are obtained for N=1,2,3,14 from importance-sampled, rigid-body diffusion Monte Carlo (DMC). Excited state energies due to helium vibrational motion near the molecule surface are evaluated using the projection operator, imaginary time spectral evolution (POITSE) method. We find excitation energies of up to ~23 K above the ground state. These states all possess vibrational character of helium atoms in a highly anisotropic potential due to the aromatic molecule, and can be categorized in terms of localized and collective vibrational modes. These results appear to provide precursors for a transition from localized to collective helium excitations at molecular nanosubstrates of increasing size. We discuss the implications of these results for analysis of anomalous spectral features in recent spectroscopic studies of large aromatic molecules in helium clusters.

Chemical Physics,Atomic and Molecular Clusters

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to understand the ground - state and excited - state properties of small - sized helium clusters (4HeN) containing a single benzene molecule, especially the vibration modes of helium atoms near the benzene molecule surface. Specifically, the authors aim to:

1. **Calculate the ground - state structure and energy**: For small helium clusters with N = 1, 2, 3, 14, calculate their ground - state structures and energies by the importance - sampling rigid - body diffusion Monte Carlo (DMC) method.

2. **Evaluate the excited - state energy**: Use the projection - operator imaginary - time spectral evolution (POITSE) method to evaluate the excited - state energy caused by the vibration motion of helium atoms near the benzene molecule surface. The energy of these excited states can be up to about 23 K.

3. **Analyze the excited - state properties**: Classify these excited states according to their spatial properties into local vibration modes and collective vibration modes, and explore their potential mechanisms of interaction with the aromatic molecule surface.

4. **Explain experimental phenomena**: The research results provide a theoretical basis for explaining the abnormal spectral features of large aromatic molecules in helium clusters observed in recent spectroscopic studies. In particular, these results reveal the transition from local to collective helium excited states, and this transition may be of great significance for understanding the helium adsorption behavior in larger molecular systems as the size of the molecular nano - substrate increases.

### Formula summary

- **Hamiltonian**:
\[
\hat{H}=\hat{T}^{(\text{imp})}+\hat{T}^{(\text{He})}+\hat{V}
\]
where:
- \(\hat{T}^{(\text{imp})}\) is the kinetic energy of the impurity rigid body
- \(\hat{T}^{(\text{He})}\) is the kinetic energy of helium atoms
- \(\hat{V}\) is the total potential energy

- **Kinetic energy of the impurity rigid body**:
\[
\hat{T}^{(\text{imp})}=-D_I\nabla_I^2 - d_a\left(\frac{\partial^2}{\partial\phi_x^2}+\frac{\partial^2}{\partial\phi_y^2}\right)-d_c\frac{\partial^2}{\partial\phi_z^2}
\]
where:
- \(D_I = \frac{\hbar^2}{2m_I}\)
- \(d_a=\frac{\hbar^2}{2I_a}\)
- \(d_c=\frac{\hbar^2}{2I_c}\)

- **Kinetic energy of helium atoms**:
\[
\hat{T}^{(\text{He})}=-D_4\sum_{j = 1}^N\nabla_j^2,\quad D_4=\frac{\hbar^2}{2m_4}
\]

- **Model potential energy**:
\[
\hat{V}=\sum_{j = 1}^N V_{\text{He - I}}(r_j)+\sum_{i < j}V_{\text{He - He}}(r_{ij})
\]

- **Helium - hydrogen interaction**:
\[
V_H(r_{\beta j})=4\epsilon_H\left[\left(\frac{\sigma_H}{r_{\beta j}}\right)^{12}-\left(\frac{\sigma_H}{r_{\beta j}}\right)^6\right]
\]

- **Helium - carbon interaction**:
\[
V_C(r_{\alpha j})=4\epsilon'_C\left[\left(\frac{\sigma'_C}{r_{\alpha j}}\right)^{14}-\left(\frac{\sigma'_C}{r_{\alpha j}}\right)^6\right]