Effects of substrate corrugation during helium adsorption on graphene in the grand canonical ensemble

Gage Erwin,Adrian Del Maestro

2024-04-29

Abstract:Adsorption of 4He on atomically flat substrates such as graphene provides a route towards the engineering of low dimensional quantum phases including superfluids and strongly interacting insulators. In this study, we explore the effects of graphene corrugation on the helium adsorption process via quantum Monte Carlo simulations in the grand canonical ensemble. We utilize an empirical adsorption potential based on the superposition of individual helium-carbon interactions and systematically control corrugation, from a smooth membrane to the fully rough potential, via the implementation of a cutoff in reciprocal space. The results highlight the importance of using a fully corrugated potential to understand the plethora of commensurate and incommensurate phases in the first adsorbed layer. Surprisingly, some residual effects of corrugation are still present before and during the promotion and onset of a second adsorbed layer, where a smooth adsorption can lead to enhanced particle fluctuations in a helium-interaction dominated regime.

Mesoscale and Nanoscale Physics

What problem does this paper attempt to address?

### What problems does this paper attempt to solve?

This paper aims to study how the wrinkling effect of the substrate affects the phase - transition behavior in low - dimensional systems when helium (4He) is adsorbed on a graphene substrate. Specifically, the authors used the worm - algorithm quantum Monte Carlo method to simulate the adsorption process of helium on graphene substrates with different degrees of wrinkling in the grand - canonical ensemble. In this way, they explored the following problems:

1. **The influence of the wrinkling effect on the first - and second - layer helium adsorption**: The authors constructed different helium - graphene potential models from smooth to fully wrinkled to test the influence of the wrinkling effect on helium adsorption. The study found that ignoring the substrate wrinkling may lead to finite - size effects and behave differently during the first - and second - layer adsorption processes.

2. **Verify previous research results**: Previous experimental and theoretical studies have shown that helium adsorption on a graphene substrate can exhibit a rich phase diagram, including the possibility of coexistence of two - dimensional superfluidity and density - wave order. However, it is still controversial whether these phenomena can be realized in simulations. In this paper, through quantum Monte Carlo simulations in the grand - canonical ensemble, the existence conditions of these phenomena are verified.

3. **Determine the appropriate potential model**: The authors compared helium - graphene potential models at different approximation levels, especially the differences between the smooth potential and the fully wrinkled potential. The research results show that in order to accurately simulate helium adsorption in low - dimensional systems, a fully wrinkled potential model must be used.

4. **The influence of the finite - size effect**: By simulating different system sizes, the authors analyzed the influence of the finite - size effect on the adsorption process. The results show that in certain chemical potential ranges, the finite - size effect is more significant, especially at high filling fractions.

### Main conclusions

- Ignoring the substrate wrinkling will lead to the finite - size effect of the first - and second - layer helium adsorption.
- The fully wrinkled helium - graphene potential model is the key to simulating helium adsorption in low - dimensional systems.
- The wrinkling effect also has a significant influence on the second - layer helium adsorption, especially at high chemical potentials.
- The influence of the finite - size effect on certain specific filling fractions is more obvious, especially at high filling fractions.

### Formula summary

The main formulas involved in the paper include:

- **Adsorption potential**:
\[
V_{\text{He - s}}(\mathbf{r}_i)=\frac{\epsilon\sigma^2}{4\pi A}\left[\left(\frac{2}{5}\left(\frac{\sigma}{z}\right)^{10}-\left(\frac{\sigma}{z}\right)^4\right)+\sum_{g\neq0}q\sum_{k = 1}e^{i\mathbf{g}\cdot[\mathbf{m}_k+\mathbf{r}_i]}\left(\frac{1}{60}\left(\frac{g\sigma^2}{2z}\right)^5K_5(gz)-\left(\frac{g\sigma^2}{2z}\right)^2K_2(gz)\right)\right]
\]

- **N - body Hamiltonian**:
\[
H =-\frac{\hbar^2}{2m}\sum_{i = 1}^N\nabla_i^2+\sum_{i = 1}^N V_{\text{He - s}}(\mathbf{r}_i)+\sum_{i<j}V_{\text{He - He}}(\mathbf{r}_i-\mathbf{r}_j)
\]

- **Average number of particles**:
\[
n=\frac{\langle\hat{N}\rangle}{N_s}
\]

- **Energy per particle**:
\[
\frac{E}{N}=\frac{\langle\hat{H}\rangle}{N}
\]

- **Compressibility**:
\[
\kappa=\fra