We study the zero-temperature phase diagram of the 2D spin-imbalanced Fermi gas with short-ranged attractive interactions using the recently developed neural network variational Monte Carlo method with the AGPs FermiNet Ansatz. The Fulde-Ferrell-Larkin-Ovchinnikov phase is observed in the weakly interacting BCS limit and a polarised superfluid is seen in the strongly interacting BEC limit. When the interactions are strong, the minority-spin momentum density is reduced almost to zero in the momentum-space region occupied by the unpaired majority-spin electrons. When the interactions are very strong, phase separation occurs, with regions containing bosonic pairs and unpaired regions occupied by the remaining majority-spin particles. In addition, we observe translational symmetry breaking at intermediate interaction strengths, where the system forms an exotic crystal of Cooper pairs in a Fermi fluid of unpaired majority-spin particles. We provide a possible explanation for the formation of the crystalline phase, explain the origins of the k-space momentum-density hole when the pairs are tightly bound, and discuss how our approach opens new directions for future work.
Deep-learning-based variational Monte Carlo has emerged as a highly accurate method for solving the many-electron Schrödinger equation. Despite favorable scaling with the number of electrons, $${\mathcal{O}}({{n}_{{\rm{el}}}}^{4})$$ O ( n el 4 ) , the practical value of deep-learning-based variational Monte Carlo is limited by the high cost of optimizing the neural network weights for every system studied. Recent research has proposed optimizing a single neural network across multiple systems, reducing the cost per system. Here we extend this approach to solids, which require numerous calculations across different geometries, boundary conditions and supercell sizes. We demonstrate that optimization of a single ansatz across these variations significantly reduces optimization steps. Furthermore, we successfully transfer a network trained on 2 × 2 × 2 supercells of LiH, to 3 × 3 × 3 supercells, reducing the number of optimization steps required to simulate the large system by a factor of 50 compared with previous work.
This paper investigates interaction-induced symmetry breaking in circular quantum dots. We start by explaining what is known about symmetry breaking in quantum dots, pointing out that the anisotropic "static Wigner molecule" ground states frequently observed in simulations are created by interference effects that occur even in the non-interacting limit. They have nothing in common with the interaction-driven crystallization of the uniform electron gas described by Wigner. This leads us to define the term Wigner molecule more carefully via a finite analog of the spontaneous symmetry breaking that arises in the homogeneous electron gas when the interactions are strong. According to this definition, the charge density patterns characteristic of true interaction-induced Wigner molecules can only be seen if a small symmetry-breaking perturbation is applied to a strongly interacting quantum dot. A simple argument based on separation of variables into center-of-mass and internal coordinates shows that the strength of the perturbation required to produce a finite effect on the density tends to zero in the limit as the strength of the interaction tends to infinity. We confirm computationally that interaction-induced Wigner molecules satisfying this definition exist. The neural-network variational Monte Carlo method used in our simulations proves more accurate than the coupled cluster and diffusion Monte Carlo methods employed in previous benchmark calculations of quantum dots at small to intermediate interaction strengths. For high interaction strengths, our neural-network variational Monte Carlo energies agree very well with existing fixed-node diffusion Monte Carlo benchmarks, proving ∼0.01% better for small values of the total spin projection Sz but ∼0.01% worse for fully spin-polarized systems.
Understanding superfluidity remains a major goal of condensed matter physics. Here we tackle this challenge utilizing the recently developed Fermionic neural network (FermiNet) wave function Ansatz for variational Monte Carlo calculations. We study the unitary Fermi gas, a system with strong, short-range, two-body interactions known to possess a superfluid ground state but difficult to describe quantitatively. We demonstrate key limitations of the FermiNet Ansatz in studying the unitary Fermi gas and propose a simple modification that outperforms the original FermiNet significantly, giving highly accurate results. We prove mathematically that the new Ansatz, which only differs from the original Ansatz by the method of antisymmetrization, is a strict generalization of the original FermiNet architecture, despite the use of fewer parameters. Our approach shares several advantages with the FermiNet: the use of a neural network removes the need for an underlying basis set; and the flexibility of the network yields extremely accurate results within a variational quantum Monte Carlo framework that provides access to unbiased estimates of arbitrary ground-state expectation values. We discuss how the method can be extended to study other superfluids.
We present an algorithm to estimate the excited states of a quantum system by variational Monte Carlo, which has no free parameters and requires no orthogonalization of the states, instead transforming the problem into that of finding the ground state of an expanded system. Arbitrary observables can be calculated, including off-diagonal expectations, such as the transition dipole moment. The method works particularly well with neural network ansätze, and by combining this method with the FermiNet and Psiformer ansätze, we can accurately recover excitation energies and oscillator strengths on a range of molecules. We achieve accurate vertical excitation energies on benzene-scale molecules, including challenging double excitations. Beyond the examples presented in this work, we expect that this technique will be of interest for atomic, nuclear, and condensed matter physics.
Deep neural networks have been very successful as highly accurate wave function Ansätze for variational Monte Carlo calculations of molecular ground states. We present an extension of one such Ansatz, FermiNet, to calculations of the ground states of periodic Hamiltonians, and study the homogeneous electron gas. FermiNet calculations of the ground-state energies of small electron gas systems are in excellent agreement with previous initiator full configuration interaction quantum Monte Carlo and diffusion Monte Carlo calculations. We investigate the spin-polarized homogeneous electron gas and demonstrate that the same neural network architecture is capable of accurately representing both the delocalized Fermi liquid state and the localized Wigner crystal state. The network converges on the translationally invariant ground state at high density and spontaneously breaks the symmetry to produce the crystalline ground state at low density, despite being given no a priori knowledge that a phase transition exists.
Grid cells are part of a widespread network which supports navigation and spatial memory. Stable grid patterns appear late in development, in concert with extracellular matrix aggregates termed perineuronal nets (PNNs) that condense around inhibitory neurons. It has been suggested that PNNs stabilize synaptic connections and long-term memories, but their role in the grid cell network remains elusive. We show that removal of PNNs leads to lower inhibitory spiking activity, and reduces grid cells' ability to create stable representations of a novel environment. Furthermore, in animals with disrupted PNNs, exposure to a novel arena corrupted the spatiotemporal relationships within grid cell modules, and the stored representations of a familiar arena. Finally, we show that PNN removal in entorhinal cortex distorted spatial representations in downstream hippocampal neurons. Together this work suggests that PNNs provide a key stabilizing element for the grid cell network.