Asymptotic equivalence of the shell-model and local-density descriptions of Coulombic systems confined by radially symmetric potentials in two and three dimensions

Solution to the Thomson Problem for Clifford Tori with an Application to Wigner Crystals

In its original version, the Thomson problem consists of the search for the minimum-energy configuration of a set of point-like electrons that are confined to the surface of a two-dimensional sphere (S 2 ) that repel each other according to Coulomb's law, in which the distance is the Euclidean distance in the embedding space of the sphere, i.e., R 3 . In this work, we consider the analogous problem where the electrons are confined to an n-dimensional flat Clifford torus T n with n = 1, 2, 3. Since the torus T n can be embedded in the complex manifold C n , we define the distance in the Coulomb law as the Euclidean distance in C n , in analogy to what is done for the Thomson problem on the sphere. The Thomson problem on a Clifford torus is of interest because supercells with the topology of a Clifford torus can be used to describe periodic systems such as Wigner crystals. In this work, we numerically solve the Thomson problem on a square Clifford torus. To illustrate the usefulness of our approach, we apply it to Wigner crystals. We demonstrate that the equilibrium configurations we obtain for large numbers of electrons are consistent with the predicted structures of Wigner crystals. Finally, in the one-dimensional case, we analytically obtain the energy spectrum and the phonon dispersion law.

Chemical Langevin equation: A path-integral view of Gillespie's derivation

In 2000, Gillespie rehabilitated the chemical Langevin equation (CLE) by describing two conditions that must be satisfied for it to yield a valid approximation of the chemical master equation (CME). In this work, we construct an original path-integral description of the CME and show how applying Gillespie's two conditions to it directly leads to a path-integral equivalent to the CLE. We compare this approach to the path-integral equivalent of a large system size derivation and show that they are qualitatively different. In particular, both approaches involve converting many sums into many integrals, and the difference between the two methods is essentially the difference between using the Euler-Maclaurin formula and using Riemann sums. Our results shed light on how path integrals can be used to conceptualize coarse-graining biochemical systems and are readily generalizable.

Self-organization of charged particles in circular geometry

The basic principles of self-organization of one-component charged particles, confined in disk and circular parabolic potentials, are proposed. A system of equations is derived, which allows us to determine equilibrium configurations for an arbitrary, but finite, number of charged particles that are distributed over several rings. Our approach reduces significantly the computational effort in minimizing the energy of equilibrium configurations and demonstrates a remarkable agreement with the values provided by molecular dynamics calculations. With the increase of particle number n>180 we find a steady formation of a centered hexagonal lattice that smoothly transforms to valence circular rings in the ground-state configurations for both potentials.

Harmonium atoms at weak confinements: The formation of the Wigner molecules

The formation of the Wigner molecules in three-dimensional assemblies of equicharged particles upon weakening of harmonic confinements with spherical symmetry is investigated using five electronic states of the four-electron harmonium atom as a representative example. Numerical results of accurate explicitly-correlated calculations, carried out for a wide range of confinement strengths ω, reveal nonmonotonic convergence of the total energies and one-electron densities to those predicted by the formalism asymptotically exact at the ω→0 limit. The convergence rates are strongly state- and property-dependent, which implies the nonexistence of a single measure quantifying the extent of the formation of the Wigner molecule in a particular state of the harmonium atom. In light of these observations, it is concluded that, like in two dimensions, the emergence of the Wigner molecules in Coulombic systems confined by spherically symmetric harmonic potentials is a complex and gradual process that takes place over a range of confinement strengths spanning several orders of magnitude.

Rovibrational states of Wigner molecules in spherically symmetric confining potentials

The strong-localization limit of three-dimensional Wigner molecules, in which repulsively interacting particles are confined by a weak spherically symmetric potential, is investigated. An explicit prescription for computation of rovibrational wavefunctions and energies that are asymptotically exact at this limit is presented. The prescription is valid for systems with arbitrary angularly-independent interparticle and confining potentials, including those involving Coulombic and screened (i.e., Yukawa/Debye) interactions. The necessary derivations are greatly simplified by explicit constructions of the Eckart frame and the parity-adapted primitive wavefunctions. The performance of the new formalism is illustrated with the three- and four-electron harmonium atoms at their strong-correlation limits. In particular, the involvement of vibrational modes with the E symmetry is readily pinpointed as the origin of the "anomalous" weak-confinement behavior of the (1)S+ state of the four-electron species that is absent in its (1)D+ companion of the strong-confinement regime.