Source author record

J. E. Gubernatis

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

cond-mat.str-el cond-mat.stat-mech cond-mat cond-mat.mtrl-sci physics.comp-ph cond-mat.supr-con physics.chem-ph

Catalog footprint

Research overview

13works

7topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2015arXiv

Triplet $p+ip$ paring correlation in doped Kane-Mele-Hubbard model: A quantum Monte Carlo study

By using the constrained-phase quantum Monte Carlo method, we performed a systematic study of the pairing correlations in the ground state of the doped Kane-Mele-Hubbard model on a honeycomb lattice. We find that pairing correlations with $d+id$ symmetry dominate close to half filling, but pairing correlations with $p+ip$ symmetry dominate as hole doping moves the system below three-quarters filling. We correlate these behaviors of the pairing correlations with the topology of the Fermi surfaces of the non-interacting problem. We also find that the effective pairing correlation is enhanced greatly as the interaction increases, and these superconducting correlations are robust against varying the spin-orbit coupling strength. Our numerical results suggest a possible way to realize spin triplet superconductivity in doped honeycomb-like materials or ultracold atoms in optical traps.

preprint2013arXiv

Magnetic Impurity Affected by Spin-Orbit Coupling: Behavior near a Topological Phase Transition

We investigate the effect of spin-orbit coupling on the behavior of magnetic impurity at the edge of a zigzag graphene ribbon by means of quantum Monte Carlo simulations. A peculiar interplay of Kane-Mele type spin-orbit and impurity-host coupling is found to affect local properties such as the impurity magnetic moment and spectral densities. The special helical nature of the topological insulator on the edge is found to affect nonlocal quantities, such as the two-particle and spin-spin correlation functions linking electrons on the impurity with those in the conduction band.

preprint2012arXiv

Behavior of a Magnetic Impurity in Graphene in the Presence of a Vacancy

With quantum Monte Carlo methods, we investigate the consequences of placing a magnetic adatom adjacent to a vacancy in a graphene sheet. We find that instead of the adatom properties depending on the energy of the adatom orbital, as in a single impurity problem, they develop a dependence on the energy of the split localized state associated with the single vacancy problem. Shifting the chemical potential through this experimentally more accessible energy scale reveals novel behavior in the spectral density, magnetic susceptibility, and the correlations of the adatom spin and charge with those of the conduction electrons. In general, the behavior of the adatom in the presence of a vacancy differs significantly from its behavior in the absence of a vacancy.

preprint2011arXiv

An Infinite Swapping Approach to the Rare-Event Sampling Problem

We describe a new approach to the rare-event Monte Carlo sampling problem. This technique utilizes a symmetrization strategy to create probability distributions that are more highly connected and thus more easily sampled than their original, potentially sparse counterparts. After discussing the formal outline of the approach and devising techniques for its practical implementation, we illustrate the utility of the technique with a series of numerical applications to Lennard-Jones clusters of varying complexity and rare-event character.

preprint2011arXiv

Magnetic Impurities in Graphene

We used a quantum Monte Carlo method to study the magnetic impurity adatoms on graphene. We found that by tuning the chemical potential we could switch the values of the impurity's local magnet moment between relatively large and small values. Our computations of the impurity's spectral density found its behavior to differ significantly from that of an impurity in a normal metal and our computations of the charge-charge and spin-spin correlations between the impurity and the conduction band electrons found them to be strongly suppressed. In general our results are consistent with those from poor man's scaling and numerical renormalization group methods.

preprint2010arXiv

Band structure of SnTe studied by Photoemission Spectroscopy

We present an angle-resolved photoemission spectroscopy study of the electronic structure of SnTe, and compare the experimental results to ab initio band structure calculations as well as a simplified tight-binding model of the p-bands. Our study reveals the conjectured complex Fermi surface structure near the L-points showing topological changes in the bands from disconnected pockets, to open tubes, and then to cuboids as the binding energy increases, resolving lingering issues about the electronic structure. The chemical potential at the crystal surface is found to be 0.5eV below the gap, corresponding to a carrier density of p =1.14x10^{21} cm^{-3} or 7.2x10^{-2} holes per unit cell. At a temperature below the cubic-rhombohedral structural transition a small shift in spectral energy of the valance band is found, in agreement with model predictions.

preprint2010arXiv

Comparative Monte Carlo Efficiency by Monte Carlo Analysis

We propose a modified power method for computing the subdominant eigenvalue $λ_2$ of a matrix or continuous operator. Here we focus on defining simple Monte Carlo methods for its application. The methods presented use random walkers of mixed signs to represent the subdominant eigenfuction. Accordingly, the methods must cancel these signs properly in order to sample this eigenfunction faithfully. We present a simple procedure to solve this sign problem and then test our Monte Carlo methods by computing the $λ_2$ of various Markov chain transition matrices. We first computed ${λ_2}$ for several one and two dimensional Ising models, which have a discrete phase space, and compared the relative efficiencies of the Metropolis and heat-bath algorithms as a function of temperature and applied magnetic field. Next, we computed $λ_2$ for a model of an interacting gas trapped by a harmonic potential, which has a mutidimensional continuous phase space, and studied the efficiency of the Metropolis algorithm as a function of temperature and the maximum allowable step size $Δ$ . Based on the $λ_2$ criterion, we found for the Ising models that small lattices appear to give an adequate picture of comparative efficiency and that the heat-bath algorithm is more efficient than the Metropolis algorithm only at low temperatures where both algorithms are inefficient. For the harmonic trap problem, we found that the traditional rule-of-thumb of adjusting $Δ$ so the Metropolis acceptance rate is around 50% range is often sub-optimal. In general, as a function of temperature or $Δ$ , $λ_2$ for this model displayed trends defining optimal efficiency that the acceptance ratio does not. The cases studied also suggested that Monte Carlo simulations for a continuum model are likely more efficient than those for a discretized version of the model.

preprint2010arXiv

Tin telluride: a weakly co-elastic metal

We report resonant ultrasound spectroscopy (RUS), dilatometry/magnetostriction, magnetotransport, magnetization, specific heat, and $^{119}$ Sn Mössbauer spectroscopy measurements on SnTe and Sn $_{0.995}$ Cr $_{0.005}$ Te. Hall measurements at $T=77$ K indicate that our Bridgman-grown single crystals have a $p$ -type carrier concentration of $3.4 \times 10^{19}$ cm $^{-3}$ and that our Cr-doped crystals have an $n$ -type concentration of $5.8 \times 10^{22}$ cm $^{-3}$ . Although our SnTe crystals are diamagnetic over the temperature range $2\, \text{K} \leq T \leq 1100\, \text{K}$ , the Cr-doped crystals are room temperature ferromagnets with a Curie temperature of 294 K. For each sample type, three-terminal capacitive dilatometry measurements detect a subtle 0.5 micron distortion at $T_c \approx 85$ K. Whereas our RUS measurements on SnTe show elastic hardening near the structural transition, pointing to co-elastic behavior, similar measurements on Sn $_{0.995}$ Cr $_{0.005}$ Te show a pronounced softening, pointing to ferroelastic behavior. Effective Debye temperature, $θ_D$ , values of SnTe obtained from $^{119}$ Sn Mössbauer studies show a hardening of phonons in the range 60--115K ( $θ_D$ = 162K) as compared with the 100--300K range ( $θ_D$ = 150K). In addition, a precursor softening extending over approximately 100 K anticipates this collapse at the critical temperature, and quantitative analysis over three decades of its reduced modulus finds $ΔC_{44}/C_{44}=A|(T-T_0)/T_0|^{-κ}$ with $κ= 0.50 \pm 0.02$ , a value indicating a three-dimensional softening of phonon branches at a temperature $T_0 \sim 75$ K, considerably below $T_c$ . We suggest that the differences in these two types of elastic behaviors lie in the absence of elastic domain wall motion in the one case and their nucleation in the other.

preprint2008arXiv

Monte Carlo Determination of Multiple Extremal Eigenpairs

We present a Monte Carlo algorithm that allows the simultaneous determination of a few extremal eigenpairs of a very large matrix without the need to compute the inner product of two vectors or store all the components of any one vector. The new algorithm, a Monte Carlo implementation of a deterministic one we recently benchmarked, is an extension of the power method. In the implementation presented, we used a basic Monte Carlo splitting and termination method called the comb, incorporated the weight cancellation method of Arnow {\it et al.}, and exploited a new sampling method, the sewing method, that does a large state space sampling as a succession of small state space samplings. We illustrate the effectiveness of the algorithm by its determination of the two largest eigenvalues of the transfer matrices for variously-sized two-dimensional, zero field Ising models. While very likely useful for other transfer matrix problems, the algorithm is however quite general and should find application to a larger variety of problems requiring a few dominant eigenvalues of a matrix.

preprint1998arXiv

Effects of doping on spin correlations in the periodic Anderson model

We studied the effects of hole doping on spin correlations in the periodic Anderson model, mainly at the full and three-quarters-full lower bands cases. In the full lower band case, strong anti-ferromagnetic correlations develop when the on-site repulsive interaction strength $U$ becomes comparable to the quasi-particle band width. In the three-quarters full case, a novel kind of spin correlation develops that is consistent with the resonance between a $(π,0)$ and a $(0,π)$ spin-density wave. In this state the spins on different sublattices appear uncorrelated. Hole doping away from the completely full case rapidly destroys the long-range anti-ferromagnetic correlations, in a manner reminiscent of the destruction of anti-ferromagnetism in the Hubbard model. In contrast to the Hubbard model, the doping does not shift the peak in the magnetic structure factor from the $(π,π)$ position. At dopings intermediate to the full and three-quarters full cases, only weak spin correlations exist.

preprint1997arXiv

Quantum Monte Carlo Study of Hole Binding and Pairing Correlations in the Three-Band Hubbard Model

We simulated the 3-band Hubbard model using the Constrained Path Monte Carlo (CPMC) method in search for a possible superconducting ground state. The CPMC is a ground state method which is free of the exponential scaling of computing time with system size. We calculated the binding energy of a pair of holes for systems up to $6 \times 4$ unit cells. We also studied the pairing correlation functions versus distance for both the d-wave and extended s-wave channels in systems up to $6 \times 6$ . We found that holes bind for a wide range of parameters and that the binding increased as the system size is increased. However, the pairing correlation functions decay quickly with distance. For the extended s channel, we found that as the Coulomb interaction $U_d$ on the Cu sites is increased, the long-range part of the correlation functions is suppressed and fluctuates around zero. For the $d_{x^2 - y^2}$ channel, we found that the correlations decay rapidly with distance towards a small positive value. However, this value becomes smaller as the interaction $U_d$ or the system size is increased.

preprint1995arXiv

A Constrained Path Quantum Monte Carlo Method for Fermion Ground States

We propose a new quantum Monte Carlo algorithm to compute fermion ground-state properties. The ground state is projected from an initial wavefunction by a branching random walk in an over-complete basis space of Slater determinants. By constraining the determinants according to a trial wavefunction $|Ψ_T \rangle$ , we remove the exponential decay of signal-to-noise ratio characteristic of the sign problem. The method is variational and is exact if $|Ψ_T\rangle$ is exact. We report results on the two-dimensional Hubbard model up to size $16\times 16$ , for various electron fillings and interaction strengths.

preprint1995arXiv

Generalization of the Fortuin-Kasteleyn transformation and its application to quantum spin simulations,

We generalize the Fortuin-Kasteleyn (FK) cluster representation of the partition function of the Ising model to represent the partition function of quantum spin models with an arbitrary spin magnitude in arbitrary dimensions. This generalized representation enables us to develop a new cluster algorithm for the simulation of quantum spin systems by the worldline Monte Carlo method. Because the Swendsen-Wang algorithm is based on the FK representation, the new cluster algorithm naturally includes it as a special case. As well as the general description of the new representation, we present an illustration of our new algorithm for some special interesting cases: the Ising model, the antiferromagnetic Heisenberg model with $S=1$ , and a general Heisenberg model. The new algorithm is applicable to models with any range of the exchange interaction, any lattice geometry, and any dimensions.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

J. E. Gubernatis

Research overview

Is this your work?

13 publications

Triplet $p+ip$ paring correlation in doped Kane-Mele-Hubbard model: A quantum Monte Carlo study

Magnetic Impurity Affected by Spin-Orbit Coupling: Behavior near a Topological Phase Transition

Behavior of a Magnetic Impurity in Graphene in the Presence of a Vacancy

An Infinite Swapping Approach to the Rare-Event Sampling Problem

Magnetic Impurities in Graphene

Band structure of SnTe studied by Photoemission Spectroscopy

Comparative Monte Carlo Efficiency by Monte Carlo Analysis

Tin telluride: a weakly co-elastic metal

Monte Carlo Determination of Multiple Extremal Eigenpairs

Effects of doping on spin correlations in the periodic Anderson model

Quantum Monte Carlo Study of Hole Binding and Pairing Correlations in the Three-Band Hubbard Model

A Constrained Path Quantum Monte Carlo Method for Fermion Ground States

Generalization of the Fortuin-Kasteleyn transformation and its application to quantum spin simulations,

See the researcher in context

Building this map preview

J. E. Gubernatis

Research overview

Is this your work?

13 publications

Triplet p+ipp+ipp+ip paring correlation in doped Kane-Mele-Hubbard model: A quantum Monte Carlo study

Magnetic Impurity Affected by Spin-Orbit Coupling: Behavior near a Topological Phase Transition

Behavior of a Magnetic Impurity in Graphene in the Presence of a Vacancy

An Infinite Swapping Approach to the Rare-Event Sampling Problem

Magnetic Impurities in Graphene

Band structure of SnTe studied by Photoemission Spectroscopy

Comparative Monte Carlo Efficiency by Monte Carlo Analysis

Tin telluride: a weakly co-elastic metal

Monte Carlo Determination of Multiple Extremal Eigenpairs

Effects of doping on spin correlations in the periodic Anderson model

Quantum Monte Carlo Study of Hole Binding and Pairing Correlations in the Three-Band Hubbard Model

A Constrained Path Quantum Monte Carlo Method for Fermion Ground States

Generalization of the Fortuin-Kasteleyn transformation and its application to quantum spin simulations,

See the researcher in context

Building this map preview

Triplet $p+ip$ paring correlation in doped Kane-Mele-Hubbard model: A quantum Monte Carlo study