Ab Initio Study of the Beryllium Isotopes ^{7}Be to ^{12}Be

We present a systematic ab initio study of the low-lying states in beryllium isotopes from ^{7}Be to ^{12}Be using nuclear lattice effective field theory with the N^{3}LO interaction. Our calculations achieve good agreement with experimental data for energies, radii, and electromagnetic properties. We introduce a novel, model-independent method to quantify nuclear shapes, uncovering a distinct pattern in the interplay between positive and negative parity states across the isotopic chain. By combining Monte Carlo sampling of the many-body density operator with a novel nucleon-grouping algorithm, the prominent two-center cluster structures, the emergence of one-neutron halo, complex nuclear molecular dynamics such as π orbital and σ orbital, emerge naturally.