Let $\Delta$ be the Laplacian on a smooth domain $\Omega\subset \mathbb{R}^2$ with Dirichlet boundary conditions. I am interested in whether the implication
\begin{align}
\Omega \text{ is asymmetric } \Rightarrow \Delta \text{ only has single eigenvalues}
\end{align}
holds, where by asymmetric I mean that the only planar rotation $\varphi$ such that $\varphi(\Omega) = \Omega$ is the identity.
It is well known that the set of Dirichlet eigenvalues $\{\lambda_i\}_{i\in \mathbb{N}}$ forms a sequence
\begin{align}
0<\lambda_1< \lambda_2 \leq \lambda_3 \leq \dots \nearrow \infty
\end{align}
with the first eigenvalue being single.
Examples of symmetric domains with multiple eigenvalues (see [1]) include a square, an equilateral triangle or a circular disk.
An example of a symmetric domain with only simple eigenvalues would be a rectangle with side lengths $l_1$, $l_2$ such that $\left(\frac{l_1}{l_2}\right)^2$ is irrational [1, 3.1].
I have carried out some numerical experiments with asymmetric shapes:
https://i.sstatic.net/uwroB.png
https://i.sstatic.net/6VJk6.png
https://i.sstatic.net/v1suq.png
but up to the 200th eigenvalue there are only single ones.
Does anyone know whether the above claim is true or false, or has ideas for counterexample shapes to check (i.e., asymmetric shapes with multiple eigenvalues)?
[1] Grebenkov, Denis S., and B-T. Nguyen. "Geometrical structure of Laplacian eigenfunctions." siam REVIEW 55.4 (2013): 601-667.