Search

We explore different ways of incorporating accurate trial wave functions into free projection auxiliary-field quantum Monte Carlo (fp-AFQMC). States employed include coupled-cluster singles and doubles, multi-Slater, and symmetry-projected mean-field wave functions. We adapt a recently proposed fast multi-Slater local energy evaluation algorithm for fp-AFQMC, making the use of long expansions from selected configuration interaction methods feasible. We demonstrate how these wave functions serve to mitigate the sign problem and accelerate convergence in quantum chemical problems, allowing the application of fp-AFQMC to systems of substantial sizes. Our calculations on the widely studied model Cu2O22+ system show that many previously reported isomerization energies differ substantially from the near-exact fp-AFQMC value.

A conjecture is made that the weight space for 4D, $\cal N$-extended supersymmetrical representations is embedded within the permutahedra associated with permutation groups ${\mathbb{S}}{}_{d}$. Adinkras and Coxeter Groups associated with minimal representations of 4D, $\cal N$ = 1 supersymmetry provide evidence supporting this conjecture. It is shown the appearance of the mathematics of 4D, $\cal N$ = 1 minimal off-shell supersymmetry representations is equivalent to solving a four color problem on the truncated octahedron. This observation suggest an entirely new way to approach the off-shell SUSY auxiliary field problem based on IT algorithms probing the properties of ${\mathbb{S}}{}_{d}$., LaTeX twice, 39 pages, 14 figures, 28 tables (v2 added one table and associated content, v3 note added in proof comment)

We present a new methodology for polymer self-consistent field theory (SCFT) that has spectral accuracy in the contour dimension while retaining linear scaling of computational effort with system size. In contrast, traditional linear-scaling algorithms only have polynomial order accuracy. The improved accuracy allows for faster simulations and lower memory costs compared to traditional algorithms. The new spectral methods are enabled by converting from an auxiliary field representation to a recently developed “polymer coherent states” framework.

The spin susceptibility in pure neutron matter is computed from auxiliary field diffusion Monte Carlo calculations over a wide range of densities. The calculations are performed for different spin asymmetries, while using twist-averaged boundary conditions to reduce finite-size effects. The employed nuclear interactions include both the phenomenological Argonne AV8$^\prime$+UIX potential and local interactions that are derived from chiral effective field theory up to next-to-next-to-leading order., 15 pages, 10 figures

We discuss the construction of a nuclear Energy Density Functional (EDF) from ab initio calculations, and we advocate the need of a methodical approach that is free from ad hoc assumptions. The equations of state (EoS) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO$_{\rm sat}$ and the phenomenological AV4$^\prime$+UIX$_{c}$ Hamiltonians as inputs in the Self-consistent Green's Function (SCGF) and Auxiliary Field Diffusion Monte Carlo (AFDMC) methods, respectively. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio-based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO$_{\rm sat}$-based EDF produces encouraging results, whereas the AV4$^\prime$+UIX$_{c}$-based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our work paves the way for a practical and systematic way to merge ab initio nuclear theory and DFT, while at the same time it sheds light on some of the critical aspects of this procedure., Comment: 14 pages, 6 figures

Resolving the heavy fermion band in the conduction electron momentum resolved spectral function of the Kondo lattice model is challenging since, in the weak coupling limit, its spectral weight is exponentially small. In this article we consider a composite fermion operator, consisting of a conduction electron dressed by spin fluctuations that shares the same quantum numbers as the electron operator. Using approximation free auxiliary field quantum Monte Carlo simulations we show that for the SU(2) spin-symmetric model on the square lattice at half filling, the composite fermion acts as a magnifying glass for the heavy fermion band. In comparison to the conduction electron residue that scales as $e^{-W/J_k}$ with $W$ the bandwidth and $J_k$ the Kondo coupling, the residue of the composite fermion tracks $J_k$. This result holds down to $J_k/W = 0.05$, and confirms the point of view that magnetic ordering, present below $J_k/W = 0.18$, does not destroy the heavy quasiparticle. We furthermore investigate the spectral function of the composite fermion in the ground state and at finite temperatures, for SU($N$) generalizations of the Kondo lattice model, as well as for ferromagnetic Kondo couplings, and compare our results to analytical calculations in the limit of high temperatures, large-$N$, large-$S$, and large $J_k$. Based on these calculations, we conjecture that the composite fermion operator provides a unique tool to study the destruction of the heavy fermion quasiparticle in Kondo breakdown transitions. The relation of our results to scanning tunneling spectroscopy and photoemission experiments is discussed., 16 pages, 14 figures, typos corrected

We consider fermions on a honeycomb lattice supplemented by a spin invariant interaction that dynamically generates a quantum spin Hall insulator. This lattice model provides an instance of Gross-Neveu Heisenberg criticality, as realized for example by the Hubbard model on the honeycomb lattice. Using auxiliary field quantum Monte Carlo simulations we show that we can compute with unprecedented precision susceptibilities of the order parameter. In O(N) Gross-Neveu transitions, the anomalous dimension of the bosonic mode grows as a function of N such that in the large-N limit it is of particular importance to consider susceptibilities rather than equal time correlations so as to minimize contributions from the background. For the N=3 case, we obtain $1/\nu=1.11(4)$, $\eta_{\phi}=0.80(9)$, and $\eta_{\psi}=0.29(2)$ for respectively the correlation length exponent, bosonic and fermionic anomalous dimensions., Comment: 11 pages,13 figures

We analyze the interaction of antineutrinos with water particle in electrochemical detectors. We postulate that some biological products generate a scalar auxiliary field B0 which breaks low-energy antineutrinos into boson vb and fermion v f particles. Low-energy anti neutrinos are suggested to interact with water molecules and produce output currents. We examine the output current of neutrino interactions in the electrochemical apparatus with chemical-reaction equations and half-cell model under postulated influence of weak interaction. The environmental neutrino is analysed. The output currents are treated to be generated by hydrogen ion and oxygen with the half-cell model with inclusion of weak interaction effect on hydrooxide ion recombination.

The bond dissociation energies of a set of 44 3 d transition metal-containing diatomics are computed with phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing a correlated sampling technique. We investigate molecules with H, N, O, F, Cl, and S ligands, including those in the 3dMLBE20 database first compiled by Truhlar and co-workers with calculated and experimental values that have since been revised by various groups. In order to make a direct comparison of the accuracy of our ph-AFQMC calculations with previously published results from 10 DFT functionals, CCSD(T), and icMR-CCSD(T), we establish an objective selection protocol which utilizes the most recent experimental results except for a few cases with well-specified discrepancies. With the remaining set of 41 molecules, we find that ph-AFQMC gives robust agreement with experiment superior to that of all other methods, with a mean absolute error (MAE) of 1.4(4) kcal/mol and maximum error of 3(3) kcal/mol (parentheses account for reported experimental uncertainties and the statistical errors of our ph-AFQMC calculations). In comparison, CCSD(T) and B97, the best performing DFT functional considered here, have MAEs of 2.8 and 3.7 kcal/mol, respectively, and maximum errors in excess of 17 kcal/mol (for the CoS diatomic). While a larger and more diverse data set would be required to demonstrate that ph-AFQMC is truly a benchmark method for transition metal systems, our results indicate that the method has tremendous potential, exhibiting unprecedented consistency and accuracy compared to other approximate quantum chemical approaches.

We report on an auxiliary field approach for solving nonlinear eigenvalue problems occurring in nano-optical systems with material dispersion. The material dispersion can be described by a rational function for the frequency-dependent permittivity, e.g., by a Drude-Lorentz model or a rational function fit to measured material data. The approach is applied to compute plasmonic resonances of a metallic grating.

The aim of this work is to discuss some aspects of the reduction of order formalism in the context of the Fadeev–Jackiw symplectic formalism, both at the classical and the quantum level. We start by reviewing the symplectic analysis in a regular theory (a higher derivative massless scalar theory), both using the Ostrogradsky prescription and also by reducing the order of the Lagrangian with an auxiliary field, showing the equivalence of these two approaches. The interpretation of the degrees of freedom is discussed in some detail. Finally, we perform the similar analysis in a singular higher derivative gauge theory (the Podolsky electrodynamics), in the reduced order formalism: we claim that this approach have the advantage of clearly separating the symplectic structure of the model into a Maxwell and a Proca (ghost) sector, thus complementing the understanding of the degrees of freedom of the theory and simplifying calculations involving matrices.