Search

Quantum chemical calculations of the ground-state properties of positron-molecule complexes are challenging. The main difficulty lies in employing an appropriate basis set for representing the coalescence between electrons and a positron. Here, we tackle this problem with the recently developed Fermionic neural network (FermiNet) wavefunction, which does not depend on a basis set. We find that FermiNet produces highly accurate, in some cases state-of-the-art, ground-state energies across a range of atoms and small molecules with a wide variety of qualitatively distinct positron binding characteristics. We calculate the binding energy of the challenging non-polar benzene molecule, finding good agreement with the experimental value, and obtain annihilation rates which compare favourably with those obtained with explicitly correlated Gaussian wavefunctions. Our results demonstrate a generic advantage of neural network wavefunction-based methods and broaden their applicability to systems beyond the standard molecular Hamiltonian., Comment: 14 pages, 3 figures

Deep neural networks have been extremely successful as highly accurate wave function ans\"atze 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 is given no \emph{a priori} knowledge that a phase transition exists, but converges on the translationally invariant ground state at high density and spontaneously breaks the symmetry to produce the crystalline ground state at low density., Comment: 12 pages, 3 figures

Density matrix quantum Monte Carlo (DMQMC) is used to sample exact-on-average $N$-body density matrices for uniform electron gas systems of up to 10$^{124}$ matrix elements via a stochastic solution of the Bloch equation. The results of these calculations resolve a current debate over the accuracy of the data used to parametrize finite-temperature density functionals. Exchange-correlation energies calculated using the real-space restricted path-integral formalism and the $k$-space configuration path-integral formalism disagree by up to $\sim$$10$\% at certain reduced temperatures $T/T_F \le 0.5$ and densities $r_s \le 1$. Our calculations confirm the accuracy of the configuration path-integral Monte Carlo results available at high density and bridge the gap to lower densities, providing trustworthy data in the regime typical of planetary interiors and solids subject to laser irradiation. We demonstrate that DMQMC can calculate free energies directly and present exact free energies for $T/T_F \ge 1$ and $r_s \le 2$., Comment: Accepted version: added free energy data and restructured text. Now includes supplementary material

Within Full Configuration Interaction Quantum Monte Carlo, we investigate how the statistical error behaves as a function of the parameters which control the stochastic sampling. We define the inefficiency as a measure of the statistical error per particle sampling the space and per timestep and show there is a sizeable parameter regime where this is minimised. We find that this inefficiency increases sublinearly with Hilbert space size and can be reduced by localising the canonical Hartree--Fock molecular orbitals, suggesting that the choice of basis impacts the method beyond that of the sign problem., Comment: 7 pages, 4 figures, Supplementary material available as 'Ancillary Files'

Hubbard-like Hamiltonians are widely used to describe on-site Coulomb interactions in magnetic and strongly-correlated solids, but there is much confusion in the literature about the form these Hamiltonians should take for shells of p and d orbitals. This paper derives the most general s, p and d orbital Hubbard-like Hamiltonians consistent with the relevant symmetries, and presents them in ways convenient for practical calculations. We use the full configuration interaction method to study p and d orbital dimers and compare results obtained using the correct Hamiltonian and the collinear and vector Stoner Hamiltonians. The Stoner Hamiltonians can fail to describe properly the nature of the ground state, the time evolution of excited states, and the electronic heat capacity., Comment: Updated the paper to make some clarifications and include colour figures

The recently developed density matrix quantum Monte Carlo (DMQMC) algorithm stochastically samples the N -body thermal density matrix and hence provides access to exact properties of many-particle quantum systems at arbitrary temperatures. We demonstrate that moving to the interaction picture provides substantial benefits when applying DMQMC to interacting fermions. In this first study, we focus on a system of much recent interest: the uniform electron gas in the warm dense regime. The basis set incompleteness error at finite temperature is investigated and extrapolated via a simple Monte Carlo sampling procedure. Finally, we provide benchmark calculations for a four-electron system, comparing our results to previous work where possible., Comment: Minor text modifications and added deviation subplot to Fig. 9 following referee's review comments. Include link to supplementary material in references

We expand upon the recent semi-stochastic adaptation to full configuration interaction quantum Monte Carlo (FCIQMC). We present an alternate method for generating the deterministic space without a priori knowledge of the wave function and present stochastic efficiencies for a variety of both molecular and lattice systems. The algorithmic details of an efficient semi-stochastic implementation are presented, with particular consideration given to the effect that the adaptation has on parallel performance in FCIQMC. We further demonstrate the benefit for calculation of reduced density matrices in FCIQMC through replica sampling, where the semi-stochastic adaptation seems to have even larger efficiency gains. We then combine these ideas to produce explicitly correlated corrected FCIQMC energies for the Beryllium dimer, for which stochastic errors on the order of wavenumber accuracy are achievable., Comment: 10 pages, 7 figures

The HANDE quantum Monte Carlo project offers accessible stochastic algorithms for general use for scientists in the field of quantum chemistry. HANDE is an ambitious and general high-performance code developed by a geographically-dispersed team with a variety of backgrounds in computational science. In the course of preparing a public, open-source release, we have taken this opportunity to step back and look at what we have done and what we hope to do in the future. We pay particular attention to development processes, the approach taken to train students joining the project, and how a flat hierarchical structure aids communication, Comment: 6 pages. Submission to WSSSPE2

We show that Full Configuration Interaction Quantum Monte Carlo (FCIQMC) is a Markov Chain in its present form. We construct the Markov matrix of FCIQMC for a two determinant system and hence compute the stationary distribution. These solutions are used to quantify the dependence of the population dynamics on the parameters defining the Markov chain. Despite the simplicity of a system with only two determinants, it still reveals a population control bias inherent to the FCIQMC algorithm. We investigate the effect of simulation parameters on the population control bias for the neon atom and suggest simulation setups to in general minimise the bias. We show a reweighting scheme to remove the bias caused by population control commonly used in Diffusion Monte Carlo [J. Chem. Phys. 99, 2865 (1993)] is effective and recommend its use as a post processing step., Comment: Supplementary material available as 'Ancillary Files'

We present a quantum Monte Carlo method capable of sampling the full density matrix of a many-particle system at finite temperature. This allows arbitrary reduced density matrix elements and expectation values of complicated non-local observables to be evaluated easily. The method resembles full configuration interaction quantum Monte Carlo but works in the space of many-particle operators instead of the space of many-particle wave functions. One simulation provides the density matrix at all temperatures simultaneously, from $T=\infty$ to $T=0$, allowing the temperature dependence of expectation values to be studied. The direct sampling of the density matrix also allows the calculation of some previously inaccessible entanglement measures. We explain the theory underlying the method, describe the algorithm, and introduce an importance-sampling procedure to improve the stochastic efficiency. To demonstrate the potential of our approach, the energy and staggered magnetization of the isotropic antiferromagnetic Heisenberg model on small lattices, the concurrence of one-dimensional spin rings, and the Renyi $S_2$ entanglement entropy of various sublattices of the $6 \times 6$ Heisenberg model are calculated. The nature of the sign problem in the method is also investigated., Comment: 13 pages, 6 figures

The recently proposed full configuration interaction quantum Monte Carlo method allows access to essentially exact ground-state energies of systems of interacting fermions substantially larger than previously tractable without knowledge of the nodal structure of the ground-state wave function. We investigate the nature of the sign problem in this method and how its severity depends on the system studied. We explain how cancelation of the positive and negative particles sampling the wave function ensures convergence to a stochastic representation of the many-fermion ground state and accounts for the characteristic population dynamics observed in simulations., Comment: 11 pages. 6 figures

The worksphere, a user controlled computer workstation enclosure, was expanded in scope to an engineering workstation suitable for use on the Space Station as a crewmember desk in orbit. The concept was also explored as a module control station capable of enclosing enough equipment to control the station from each module. The concept has commercial potential for the Space Station and surface workstation applications. The central triangular beam interior configuration was expanded and refined to seven different beam configurations. These included triangular on center, triangular off center, square, hexagonal small, hexagonal medium, hexagonal large and the H beam. Each was explored with some considerations as to the utilities and a suggested evaluation factor methodology was presented. Scale models of each concept were made. The models were helpful in researching the seven beam configurations and determining the negative residual (unused) volume of each configuration. A flexible hardware evaluation factor concept is proposed which could be helpful in evaluating interior space volumes from a human factors point of view. A magnetic version with all the graphics is available from the author or the technical monitor.

A study was conducted to define the architectural elements and issues of the Space Station. The objective of the study was to identify those questions which require further research and suggest ways in which the research can be undertaken. The study examined five primary topics, asked salient questions and described the merits of alternative solutions.

Lessons learned from the Alaskan North Slope construction camps and the Western Regional Undersea Laboratory are analyzed with respect to possible improvements for space station interior space utilization and living areas. The human factors engineering aspects have a direct bearing on the condition of crew and occupants.

Application of the connected-element radio interferometer of the National Radio Astronomy Observation in Green Bank, West Virginia to the determination of improved source coordinates, astronomical constants, and variations in Earth rotation parameters is discussed. It is concluded that because of the brevity and discontinuity of the data so far no reliable conclusions regarding the accuracy of the data can be drawn.