Multiphase magnetism in Yb2Ti2O7.

We use neutron scattering to show that ferromagnetism and antiferromagnetism coexist in the low T state of the pyrochlore quantum magnet [Formula: see text] While magnetic Bragg peaks evidence long-range static ferromagnetic order, inelastic scattering shows that short-range correlated antiferromagnetism is also present. Small-angle neutron scattering provides direct evidence for mesoscale magnetic structure that we associate with metastable antiferromagnetism. Classical Monte Carlo simulations based on exchange interactions inferred from [Formula: see text]-oriented high-field spin wave measurements confirm that antiferromagnetism is metastable within the otherwise ferromagnetic ground state. The apparent lack of coherent spin wave excitations and strong sensitivity to quenched disorder characterizing [Formula: see text] is a consequence of this multiphase magnetism.

Spin-lattice Coupling and the Emergence of the Trimerized Phase in the S=1 Kagome Antiferromagnet Na_{2}Ti_{3}Cl_{8}.

Spin-1 antiferromagnets are abundant in nature, but few theories exist to understand their properties and behavior when geometric frustration is present. Here we study the S=1 kagome compound Na_{2}Ti_{3}Cl_{8} using a combination of density functional theory, exact diagonalization, and density matrix renormalization group approaches to achieve a first principles supported explanation of its exotic magnetic phases. We find that the effective magnetic Hamiltonian includes essential non-Heisenberg terms that do not stem from spin-orbit coupling, and both trimerized and spin-nematic magnetic phases are relevant. The experimentally observed structural transition to a breathing kagome phase is driven by spin-lattice coupling, which favors the trimerized magnetic phase against the quadrupolar one. We thus show that lattice effects can be necessary to understand the magnetism in frustrated magnetic compounds and surmise that Na_{2}Ti_{3}Cl_{8} is a compound that cannot be understood from only electronic or only lattice Hamiltonians, very much like VO_{2}.

Tunable Magnon Interactions in a Ferromagnetic Spin-1 Chain.

NiNb_{2}O_{6} is an almost ideal realization of a 1D spin-1 ferromagnetic Heisenberg chain compound with weak unidirectional anisotropy. Using time-domain THz spectroscopy, we measure the low-energy electrodynamic response of NiNb_{2}O_{6} as a function of temperature and external magnetic field. At low temperatures, we find a magnonlike spin excitation, which corresponds to the lowest energy excitation at q∼0. At higher temperatures, we unexpectedly observe a temperature-dependent renormalization of the spin-excitation energy, which has a strong dependence on field direction. Using theoretical arguments, exact diagonalizations, and finite temperature dynamical Lanczos calculations, we construct a picture of magnon-magnon interactions that naturally explains the observed renormalization. We show how magnetic field strength and direction may be used to directly tune the sign of the magnon-magnon interaction. This unique scenario is a consequence of the spin-1 nature and has no analog in the more widely studied spin-1/2 systems.

Dynamical Structure Factor of the Three-Dimensional Quantum Spin Liquid Candidate NaCaNi_{2}F_{7}.

We study the dynamical structure factor of the spin-1 pyrochlore material NaCaNi_{2}F_{7}, which is well described by a weakly perturbed nearest-neighbour Heisenberg Hamiltonian, Our three approaches-molecular dynamics simulations, stochastic dynamical theory, and linear spin wave theory-reproduce remarkably well the momentum dependence of the experimental inelastic neutron scattering intensity as well as its energy dependence with the exception of the lowest energies. We discuss two surprising aspects and their implications for quantum spin liquids in general: the complete lack of sharp quasiparticle excitations in momentum space and the success of the linear spin wave theory in a regime where it would be expected to fail for several reasons.

Macroscopically Degenerate Exactly Solvable Point in the Spin-1/2 Kagome Quantum Antiferromagnet.

Frustrated quantum magnets are a central theme in condensed matter physics due to the richness of their phase diagrams. They support a panoply of phases including various ordered states and topological phases. Yet, this problem has defied a solution for a long time due to the lack of controlled approximations which make it difficult to distinguish between competing phases. Here we report the discovery of a special quantum macroscopically degenerate point in the XXZ model on the spin-1/2 kagome quantum antiferromagnet for the ratio of Ising to antiferromagnetic transverse coupling J_{z}/J=-1/2. This point is proximate to many competing phases explaining the source of the complexity of the phase diagram. We identify five phases near this point including both spin-liquid and broken-symmetry phases and give evidence that the kagome Heisenberg antiferromagnet is close to a transition between two phases.

Density-matrix based determination of low-energy model Hamiltonians from ab initio wavefunctions.

We propose a way of obtaining effective low energy Hubbard-like model Hamiltonians from ab initio quantum Monte Carlo calculations for molecular and extended systems. The Hamiltonian parameters are fit to best match the ab initio two-body density matrices and energies of the ground and excited states, and thus we refer to the method as ab initio density matrix based downfolding. For benzene (a finite system), we find good agreement with experimentally available energy gaps without using any experimental inputs. For graphene, a two dimensional solid (extended system) with periodic boundary conditions, we find the effective on-site Hubbard U(∗)/t to be 1.3 ± 0.2, comparable to a recent estimate based on the constrained random phase approximation. For molecules, such parameterizations enable calculation of excited states that are usually not accessible within ground state approaches. For solids, the effective Hamiltonian enables large-scale calculations using techniques designed for lattice models.

Emergent spin excitations in a Bethe lattice at percolation.

We study the spin-1/2 quantum Heisenberg antiferromagnet on a Bethe lattice diluted to the percolation threshold. Dilution creates areas of even or odd sublattice imbalance resulting in "dangling spins" [L. Wang and A. W. Sandvik, Phys. Rev. Lett. 97, 117204 (2006); Phys. Rev. B 81, 054417 (2010)]. These collectively act as "emergent" spin-1/2 degrees of freedom and are responsible for the creation of a set of low-lying "quasidegenerate states." Using density matrix renormalization group calculations, we detect the presence and location of these emergent spins. We find an effective Hamiltonian of these emergent spins, with Heisenberg interactions that decay exponentially with the distance between them.

Semistochastic projector Monte Carlo method.

We introduce a semistochastic implementation of the power method to compute, for very large matrices, the dominant eigenvalue and expectation values involving the corresponding eigenvector. The method is semistochastic in that the matrix multiplication is partially implemented numerically exactly and partially stochastically with respect to expectation values only. Compared to a fully stochastic method, the semistochastic approach significantly reduces the computational time required to obtain the eigenvalue to a specified statistical uncertainty. This is demonstrated by the application of the semistochastic quantum Monte Carlo method to systems with a sign problem: the fermion Hubbard model and the carbon dimer.

Investigation of the monopole magneto-chemical potential in spin ices using capacitive torque magnetometry.

The single-ion anisotropy and magnetic interactions in spin-ice systems give rise to unusual non-collinear spin textures, such as Pauling states and magnetic monopoles. The effective spin correlation strength (Jeff) determines the relative energies of the different spin-ice states. With this work, we display the capability of capacitive torque magnetometry in characterizing the magneto-chemical potential associated with monopole formation. We build a magnetic phase diagram of Ho2Ti2O7, and show that the magneto-chemical potential depends on the spin sublattice (α or ß), i.e., the Pauling state, involved in the transition. Monte Carlo simulations using the dipolar-spin-ice Hamiltonian support our findings of a sublattice-dependent magneto-chemical potential, but the model underestimates the Jeff for the ß-sublattice. Additional simulations, including next-nearest neighbor interactions (J2), show that long-range exchange terms in the Hamiltonian are needed to describe the measurements. This demonstrates that torque magnetometry provides a sensitive test for Jeff and the spin-spin interactions that contribute to it.

Continuum of quantum fluctuations in a three-dimensional S = 1 Heisenberg magnet.

Conventional crystalline magnets are characterized by symmetry breaking and normal modes of excitation called magnons, with quantized angular momentum ℏ. Neutron scattering correspondingly features extra magnetic Bragg diffraction at low temperatures and dispersive inelastic scattering associated with single magnon creation and annihilation. Exceptions are anticipated in so-called quantum spin liquids, as exemplified by the one-dimensional spin-1/2 chain, which has no magnetic order and where magnons accordingly fractionalize into spinons with angular momentum ℏ/2. This is spectacularly revealed by a continuum of inelastic neutron scattering associated with two-spinon processes. Here, we report evidence for these key features of a quantum spin liquid in the three-dimensional antiferromagnet NaCaNi2F7. We show that despite the complication of random Na1+-Ca2+ charge disorder, NaCaNi2F7 is an almost ideal realization of the spin-1 antiferromagnetic Heisenberg model on a pyrochlore lattice. Magnetic Bragg diffraction is absent and 90% of the neutron spectral weight forms a continuum of magnetic scattering with low-energy pinch points, indicating NaCaNi2F7 is in a Coulomb-like phase. Our results demonstrate that disorder can act to freeze only the lowest-energy magnetic degrees of freedom; at higher energies, a magnetic excitation continuum characteristic of fractionalized excitations persists.

Efficient Heat-Bath Sampling in Fock Space.

We introduce an algorithm for sampling many-body quantum states in Fock space. The algorithm efficiently samples states with probability approximately proportional to an arbitrary function of the second-quantized Hamiltonian matrix element connecting the sampled state to the current state. We apply the new sampling algorithm to the recently developed semistochastic full configuration interaction quantum Monte Carlo (S-FCIQMC) method, a semistochastic implementation of the power method for projecting out the ground state energy in a basis of Slater determinants. Our new sampling method requires modest additional computational time and memory compared to uniform sampling but results in newly spawned weights that are approximately of the same magnitude, thereby greatly improving the efficiency of projection. A comparison in efficiency between our sampling algorithm and uniform sampling is performed on the all-electron nitrogen dimer at equilibrium in Dunning's cc-pVXZ basis sets with X ∈ {D, T, Q, 5}, demonstrating a large gain in efficiency that increases with basis set size. In addition, a comparison in efficiency is performed on three all-electron first-row dimers, B2, N2, and F2, in a cc-pVQZ basis, demonstrating that the gain in efficiency compared to uniform sampling also increases dramatically with the number of electrons.

Charge density waves in disordered media circumventing the Imry-Ma argument.

Two powerful theoretical predictions, Anderson localization and the Imry-Ma argument, impose significant restrictions on the phases of matter that can exist in the presence of even the smallest amount of disorder in one-dimensional systems. These predictions forbid electrically conducting states and ordered states respectively. It was thus remarkable that a mechanism to circumvent Anderson localization relying on the presence of correlated disorder was found, that is also realized in certain biomolecular systems. In a similar manner, we show that the Imry-Ma argument can be circumvented, resulting in the formation of stable ordered states with discrete broken symmetries in disordered one dimensional systems. We then investigate other mechanisms by which disorder can destroy an ordered state.