BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Joaquín E. Drut

Joaquín E. Drut appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-latnucl-thcond-mat.quant-gascond-mat.stat-mechcond-mat.str-elquant-phhep-phcond-mat.mes-hallcond-mat.othercond-mat.supr-conhep-thphysics.chem-phphysics.comp-ph

Catalog footprint

What is connected

32works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

32 published item(s)

preprint2021arXiv

Pairing patterns in polarized unitary Fermi gases above the superfluid transition

We non-perturbatively study pairing in the high-temperature regime of polarized unitary two-component Fermi gases by extracting the pair-momentum distribution and shot-noise correlations. Whereas the pair-momentum distribution allows us to analyze the propagation of pairs composed of one spin-up and one spin-down fermion, shot-noise correlations provide us with a tomographic insight into pairing correlations around the Fermi surfaces associated with the two species. Assuming that the dominant pairing patterns right above the superfluid transition also govern the formation of condensates in the low-temperature regime, our analysis suggests that the superfluid ground state is homogeneous and of the Bardeen-Cooper-Schrieffer-type over a wide range of polarizations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Emergence of a pseudogap in the BCS-BEC crossover

Strongly correlated Fermi systems with pairing interactions become superfluid below a critical temperature $T_c$. The extent to which such pairing correlations alter the behavior of the liquid at temperatures $T > T_c$ is a subtle issue that remains an area of debate, in particular regarding the appearance of the so-called pseudogap in the BCS-BEC crossover of unpolarized spin-$1/2$ nonrelativistic matter. To shed light on this, we extract several quantities of crucial importance at and around the unitary limit, namely: the odd-even staggering of the total energy, the spin susceptibility, the pairing correlation function, the condensate fraction, and the critical temperature $T_c$, using a non-perturbative, constrained-ensemble quantum Monte Carlo algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pairing patterns in one-dimensional spin- and mass-imbalanced Fermi gases

We study spin- and mass-imbalanced mixtures of spin-$\tfrac{1}{2}$ fermions interacting via an attractive contact potential in one spatial dimension. Specifically, we address the influence of unequal particle masses on the pair formation by means of the complex Langevin method. By computing the pair-correlation function and the associated pair-momentum distribution we find that inhomogeneous pairing is present for all studied spin polarizations and mass imbalances. To further characterize the pairing behavior, we analyze the density-density correlations in momentum space, the so-called shot noise, which is experimentally accessible through time-of-flight imaging. At finite spin polarization, the latter is known to show distinct maxima at momentum configurations associated with the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) instability. Besides those maxima, we find that additional features emerge in the noise correlations when mass imbalance is increased, revealing the stability of FFLO-type correlations against mass imbalance and furnishing an experimentally accessible signature to probe this type of pairing.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Thermodynamics of spin-orbit coupled bosons in two dimensions from complex Langevin

We investigate the thermal properties of interacting spin-orbit coupled bosons with contact interactions in two spatial dimensions. To that end, we implement the complex Langevin method, motivated by the appearance of a sign problem, on a square lattice with periodic boundary conditions. We calculate the density equation of state non-perturbatively in a range of spin-orbit couplings and chemical potentials. Our results show that mean-field solutions tend to underestimate the average density, especially for stronger values of the spin-orbit coupling. Additionally, the finite nature of the simulation volume induces the formation of pseudo-condensates. These have been observed to be destroyed by the spin-orbit interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Entanglement, noise, and the cumulant expansion

We put forward a simpler and improved variation of a recently proposed method to overcome the signal-to-noise problem found in Monte Carlo calculations of the entanglement entropy of interacting fermions. The present method takes advantage of the approximate lognormal distributions that characterize the signal-to-noise properties of other approaches. In addition, we show that a simple rewriting of the formalism allows circumvention of the inversion of the restricted one-body density matrix in the calculation of the $n$-th Rényi entanglement entropy for $n>2$. We test our technique by implementing it in combination with the hybrid Monte Carlo algorithm and calculating the $n=2,3,4, \dots, 10$ Rényi entropies of the 1D attractive Hubbard model. We use that data to extrapolate to the von Neumann ($n=1$) and $n\to\infty$ cases.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The entanglement spectrum and Rényi entropies of non-relativistic conformal fermions

We characterize non-perturbatively the Rényi entropies of degree n=2,3,4, and 5 of three-dimensional, strongly coupled many-fermion systems in the scale-invariant regime of short interaction range and large scattering length, i.e. in the unitary limit. We carry out our calculations using lattice methods devised recently by us. Our results show the effect of strong pairing correlations on the entanglement entropy, which modify the sub-leading behavior for large subsystem sizes (as characterized by the dimensionless parameter x=kF L_A, where kF is the Fermi momentum and L_A the linear subsystem size), but leave the leading order unchanged relative to the non-interacting case. Moreover, we find that the onset of the sub-leading asymptotic regime is at surprisingly small x=2-4. We provide further insight into the entanglement properties of this system by analyzing the spectrum of the entanglement Hamiltonian of the two-body problem from weak to strong coupling. The low-lying entanglement spectrum displays clear features as the strength of the coupling is varied, such as eigenvalue crossing, a sharp change in the Schmidt gap, and scale invariance at unitarity. Beyond the low-lying component, the spectrum appears as a quasi-continuum distribution, for which we present a statistical characterization; we find, in particular, that the mean shifts to infinity as the coupling is turned off, which indicates that that part of the spectrum represents non-perturbative contributions to the entanglement Hamiltonian. In contrast, the low-lying entanglement spectrum evolves to finite values in the noninteracting limit. The scale invariance of the unitary regime guarantees that our results are universal features intrinsic to 3D quantum mechanics and represent a well-defined prediction for ultracold atom experiments, which were recently shown to have direct access to the entanglement entropy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A hybrid Monte Carlo approach to the entanglement entropy of interacting fermions

The Monte Carlo calculation of Rényi entanglement entropies $S^{}_n$ of interacting fermions suffers from a well-known signal-to-noise problem, even for a large number of situations in which the infamous sign problem is absent. A few methods have been proposed to overcome this issue, such as ensemble switching and the use of auxiliary partition-function ratios. Here, we present an approach that builds on the recently proposed free-fermion decomposition method; it incorporates entanglement in the probability measure in a natural way; it takes advantage of the hybrid Monte Carlo algorithm (an essential tool in lattice quantum chromodynamics and other gauge theories with dynamical fermions); and it does not suffer from noise problems. This method displays no sign problem for the same cases as other approaches and is therefore useful for a wide variety of systems. As a proof of principle, we calculate $S_2^{}$ for the one-dimensional, half-filled Hubbard model and compare with results from exact diagonalization and the free-fermion decomposition method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Convexity of the entanglement energy of SU($2N$)-symmetric fermions with attractive interactions

The positivity of the probability measure of attractively interacting systems of $2N$-component fermions enables the derivation of an exact convexity property for the ground-state energy of such systems. Using analogous arguments, applied to path-integral expressions for the $n$-th Rényi entanglement entropy $S^{}_n$ derived recently, we prove non-perturbative analytic relations for the entanglement energies $\mathcal E^{}_n$ of those systems defined via $\mathcal{E}^{}_n \equiv \frac{n-1}{n}T S^{}_n + F$ where $β= 1/T$ is the extent of the imaginary time direction and $-βF = \ln \mathcal Z$ where $\mathcal Z$ is the partition sum appropriate to the temperature. These relations are valid for all sub-system sizes, particle numbers and dimensions, and in arbitrary external trapping potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hard-Wall and Non-Uniform Lattice Monte Carlo Approaches to One-Dimensional Fermi Gases in a Harmonic Trap

We present in detail two variants of the lattice Monte Carlo method aimed at tackling systems in external trapping potentials: a uniform-lattice approach with hard-wall boundary conditions, and a non-uniform Gauss-Hermite lattice approach. Using those two methods, we compute the ground-state energy and spatial density profile for systems of N=4 - 8 harmonically trapped fermions in one dimension. From the favorable comparison of both energies and density profiles (particularly in regions of low density), we conclude that the trapping potential is properly resolved by the hard-wall basis. Our work paves the way to higher dimensions and finite temperature analyses, as calculations with the hard-wall basis can be accelerated via fast Fourier transforms, the cost of unaccelerated methods is otherwise prohibitive due to the unfavorable scaling with system size.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase structure of mass- and spin-imbalanced unitary Fermi gases

We study the phase diagram of mass- and spin-imbalanced unitary Fermi gases, in search for the emergence of spatially inhomogeneous phases. To account for fluctuation effects beyond the mean-field approximation, we employ renormalization group techniques. We thus obtain estimates for critical values of the temperature, mass and spin imbalance, above which the system is in the normal phase. In the unpolarized, equal-mass limit, our result for the critical temperature is in accordance with state-of-the-art Monte Carlo calculations. In addition, we estimate the location of regions in the phase diagram where inhomogeneous phases are likely to exist. We show that an intriguing relation exists between the general structure of the many-body phase diagram and the binding energies of the underlying two-body bound-state problem, which further supports our findings. Our results suggest that inhomogeneous condensates form for mass ratios of the spin-down and spin-up fermions greater than three. The extent of the inhomogeneous phase in parameter space increases with increasing mass imbalance.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phases of spin- and mass-imbalanced ultracold Fermi gases in harmonic traps

We analyze the phase structure of mass- and spin-imbalanced unitary Fermi gases in harmonic traps. To this end, we employ Density Functional Theory in the local density approximation. Depending on the values of the control parameters measuring mass and spin imbalance, we observe that three regions exist in the trap, namely: a superfluid region at the center, surrounded by a mixed region of resonantly interacting spin-up and spin-down fermions, and finally a fully polarized phase surrounding the previous two regions. We also find regimes in the phase diagram where the existence of a superfluid region at the center of the trap is not energetically favored. We point out the limitations of our approach at the present stage, and call for more detailed (ab initio) studies of the equation of state of uniform, mass-imbalanced unitary Fermi gases.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Velocity renormalization in graphene from lattice Monte Carlo

We compute the Fermi velocity of the Dirac quasiparticles in clean graphene at the charge neutrality point for strong Coulomb coupling alpha_g. We perform a Lattice Monte Carlo calculation within the low-energy Dirac theory, which includes an instantaneous, long-range Coulomb interaction. We find a renormalized Fermi velocity v_FR > v_F, where v_F = c/300. Our results are consistent with a momentum-independent v_FR which increases approximately linearly with alpha_g, although a logarithmic running with momentum cannot be excluded at present. At the predicted critical coupling alpha_gc for the semimetal-insulator transition due to excitonic pair formation, we find v_FR/v_F = 3.3, which we discuss in light of experimental findings for v_FR/v_F at the charge neutrality point in ultra-clean suspended graphene.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Zero-temperature equation of state of mass-imbalanced resonant Fermi gases

We calculate the zero-temperature equation of state of mass-imbalanced resonant Fermi gases in an ab initio fashion, by implementing the recent proposal of imaginary-valued mass difference to bypass the sign problem in lattice Monte Carlo calculations. The fully non-perturbative results thus obtained are analytically continued to real mass imbalance to yield the physical equation of state, providing predictions for upcoming experiments with mass-imbalanced atomic Fermi gases. In addition, we present an exact relation for the rate of change of the equation of state at small mass imbalances, showing that it is fully determined by the energy of the mass-balanced system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cooper pairing above the critical temperature in a unitary Fermi gas

We present an ab initio determination of spin responses of the unitary Fermi gas, based on finite temperature quantum Monte Carlo calculations and the Kubo linear-response formalism. We determine the temperature dependence of the spin susceptibility and the spin conductivity. We show that both quantities exhibit suppression above the critical temperature of the superfluid to normal phase transition due to presence of the Cooper pairs. The spin diffusion transport coefficient does not display the existence of a minimum in the vicinity of the critical temperature and it drops to very low values D_s approx 0.8hbar/m in the superfluid phase. All these spin observables show a smooth and monotonic behavior with temperature when crossing the critical temperature T_c, until the Fermi liquid regime is attained at the temperature T*, where the pseudogap regime disappears.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fermi gases with imaginary mass imbalance and the sign problem in Monte Carlo calculations

Fermi gases in strongly coupled regimes, such as the unitary limit, are inherently challenging for many-body methods. Although much progress has been made with purely analytic methods, quantitative results require ab initio numerical approaches, such as Monte Carlo (MC) calculations. However, mass-imbalanced and spin-imbalanced gases are not accessible to MC calculations due to the infamous sign problem. It was recently pointed out that the sign problem, for finite spin imbalance, can be circumvented by resorting to imaginary polarizations and analytic continuation. Large parts of the phase diagram spanned by temperature and polarization then become accessible to MC calculations. We propose to apply a similar strategy to the mass-imbalanced case, which opens up the possibility to study the associated phase diagram with MC calculations. In particular, our analysis suggests that a detection of a (tri-)critical point in this phase diagram is possible. We also discuss calculations in the zero-temperature limit with our approach.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Imaginary polarization as a way to surmount the sign problem in ab initio calculations of spin-imbalanced Fermi gases

From ultracold atoms to quantum chromodynamics, reliable ab initio studies of strongly interacting fermions require numerical methods, typically in some form of quantum Monte Carlo calculation. Unfortunately, (non)relativistic systems at finite density (spin polarization) generally have a sign problem, such that those ab initio calculations are impractical. It is well-known, however, that in the relativistic case imaginary chemical potentials solve this problem, assuming the data can be analytically continued to the real axis. Is this feasible for nonrelativistic systems? Are the interesting features of the phase diagram accessible in this manner? By introducing complex chemical potentials, for real total particle number and imaginary polarization, the sign problem is avoided in the nonrelativistic case. To give a first answer to the above questions, we perform a mean-field study of the finite-temperature phase diagram of spin-1/2 fermions with imaginary polarization.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Inhomogeneous phases in one-dimensional mass- and spin-imbalanced Fermi gases

We compute the phase diagram of strongly interacting fermions in one dimension at finite temperature, with mass and spin imbalance. By including the possibility of the existence of a spatially inhomogeneous ground state, we find regions where spatially varying superfluid phases are favored over homogeneous phases. We obtain estimates for critical values of the temperature, mass and spin imbalance, above which these phases disappear. Finally, we show that an intriguing relation exists between the general structure of the phase diagram and the binding energies of the underlying two-body bound-state problem.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lattice methods for strongly interacting many-body systems

Lattice field theory methods, usually associated with non-perturbative studies of quantum chromodynamics, are becoming increasingly common in the calculation of ground-state and thermal properties of strongly interacting non-relativistic few- and many-body systems, blurring the interfaces between condensed matter, atomic and low-energy nuclear physics. While some of these techniques have been in use in the area of condensed matter physics for a long time, others, such as hybrid Monte Carlo and improved effective actions, have only recently found their way across areas. With this topical review, we aim to provide a modest overview and a status update on a few notable recent developments. For the sake of brevity we focus on zero-temperature, non-relativistic problems. After a short introduction, we lay out some general considerations and proceed to discuss sampling algorithms, observables, and systematic effects. We show selected results on ground- and excited-state properties of fermions in the limit of unitarity. The appendix contains details on group theory on the lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Precise determination of the structure factor and contact in a unitary Fermi gas

We present a high-precision determination of the universal contact parameter in a strongly interacting Fermi gas. In a trapped gas at unitarity we find the contact to be $3.06 \pm 0.08$ at a temperature of 0.08 of the Fermi temperature in a harmonic trap. The contact governs the high-momentum (short-range) properties of these systems and this low temperature measurement provides a new benchmark for the zero temperature homogeneous contact. The experimental measurement utilises Bragg spectroscopy to obtain the dynamic and static structure factors of ultracold Fermi gases at high momentum in the unitarity and molecular Bose-Einstein condensate (BEC) regimes. We have also performed quantum Monte Carlo calculations of the static properties, extending from the weakly coupled Bardeen-Cooper-Schrieffer (BCS) regime to the strongly coupled BEC case, which show agreement with experiment at the level of a few percent.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Improved lattice operators for non-relativistic fermions

In this work I apply a recently proposed improvement procedure, originally conceived to reduce finite lattice spacing effects in transfer matrices for dilute Fermi systems, to tuning operators for the calculation of observables. I construct, in particular, highly improved representations for the energy and the contact, as a first step in an improvement program for finite-temperature calculations. I illustrate the effects of improvement on those quantities with a ground-state lattice calculation at unitarity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Shear Viscosity of a Unitary Fermi Gas

We present the first ab initio determination of the shear viscosity eta of the Unitary Fermi Gas, based on finite temperature quantum Monte Carlo calculations and the Kubo linear-response formalism. We determine the temperature dependence of the shear viscosity to entropy density ratio eta/s. The minimum of eta/s appears to be located above the critical temperature for the superfluid-to-normal phase transition with the most probable value being eta/s approx 0.2 hbar/kB, which almost saturates the Kovtun-Son-Starinets universal value hbar/(4 pi kB).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Equation of State of the Unitary Fermi Gas: An Update on Lattice Calculations

The thermodynamic properties of the unitary Fermi gas (UFG) have recently been measured to unprecedented accuracy at the MIT. In particular, these measurements provide an improved understanding of the regime below T/eF ~ 0.20, where a transition into a superfluid phase occurs. In light of this development, we present an overview of state-of-the-art auxiliary field quantum Monte Carlo (AFQMC) results for the UFG at finite temperature and compare them with the MIT data for the energy, chemical potential, and density. These AFQMC results have been obtained using methods based on the hybrid Monte Carlo (HMC) algorithm, which was first introduced within the context of lattice QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Exact-exchange density functional theory for neutron drops

We compute the ground-state properties of finite systems of neutrons in an external harmonic trap, interacting via the Minnesota potential, using the "exact-exchange" form of orbital-dependent density functional theory. We compare our results with Hartree-Fock calculations and find very close agreement. Within the context of the interaction studied, we conclude that this simple orbital-dependent functional brings conventional nuclear density functional theory to the level of Hartree-Fock in an ab initio fashion. Our work is a first step towards higher-order ab initio nuclear functionals based on realistic nucleon-nucleon interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Momentum Distribution and Contact of the Unitary Fermi gas

We calculate the momentum distribution n(k) of the Unitary Fermi Gas using Quantum Monte Carlo calculations at finite temperature T/ε_F as well as in the ground state. At large momenta k/k_F, we find that n(k) falls off as C/k^4, in agreement with the Tan relations. From the asymptotics of n(k), we determine the contact C as a function of T/ε_F and present a comparison with theory. At low T/ε_F, we find that C increases with temperature, and we tentatively identify a maximum around T/ε_F \simeq 0.4. Our calculations are performed on lattices of spatial extent up to N_x = 14 with a particle number per unit volume of \simeq 0.03 - 0.07.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Strongly Coupled Graphene on the Lattice

The two-dimensional carbon allotrope graphene has recently attracted a lot of attention from researchers in the disciplines of Lattice Field Theory, Lattice QCD and Monte Carlo calculations. This interest has been prompted by several remarkable properties of the conduction electrons in graphene. For instance, the conical band structure of graphene at low energies is strongly reminiscent of relativistic Dirac fermions. Also, due the low Fermi velocity of v_F = c/300, where c is the speed of light in vacuum, the physics of the conduction electrons in graphene is qualitatively similar to Quantum Electrodynamics in a strongly coupled regime. In turn, this opens up the prospect of the experimental realization of gapped, strongly correlated states in the electronic phase diagram of graphene. Here, we review the experimental and theoretical motivations for Lattice Field Theory studies of graphene, and describe the directions that such research is likely to progress in during the next few years. We also give a brief overview of the two main lattice theories of graphene, the hexagonal Hubbard theory and the low-energy Dirac theory. Finally, we describe the prospect of extracting response functions, such as the electric conductivity, using Lattice Field Theory calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The unitary Fermi gas at finite temperature: momentum distribution and contact

The Unitary Fermi Gas (UFG) is one of the most strongly interacting systems known to date, as it saturates the unitarity bound on the quantum mechanical scattering cross section. The UFG corresponds to a two-component Fermi gas in the limit of short interaction range and large scattering length, and is currently realized in ultracold-atom experiments via Feshbach resonances. While easy to define, the UFG poses a challenging quantum many-body problem, as it lacks any characteristic scale other than the density. As a consequence, accurate quantitative predictions of the thermodynamic properties of the UFG require Monte Carlo calculations. However, significant progress has also been made with purely analytical methods. Notably, in 2005 Tan derived a set of exact thermodynamic relations in which a universal quantity known as the "contact" C plays a crucial role. Recently, C has also been found to determine the prefactor of the high- frequency power-law decay of correlators as well as the right-hand-sides of shear- and bulk viscosity sum rules. The contact is therefore a central piece of information on the UFG in equilibrium as well as away from equilibrium. In this talk we describe some of the known aspects of Fermi gases at and around unitarity, show our latest Monte Carlo results for the contact at finite temperature, and summarize the open questions in the field, some of which we are starting to answer using large-scale Monte Carlo calculations by adapting methods from Lattice QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

First-order chiral transition in the compact lattice theory of graphene and the case for improved actions

A comparison of the compact and non-compact lattice versions of the low-energy theory of graphene is presented. The compact theory is found to exhibit a chiral phase transition which appears to be of first order, at a critical coupling of beta_c = 0.42 +/- 0.01. We confirm that the non-compact theory exhibits a second-order transition at beta_c = 0.072 +/- 0.003, and determine the effects of UV-divergent tadpole contributions in both cases. Upon tadpole improvement of the non-compact theory we find beta_c^{TI} = 0.163 +/- 0.002, which strengthens the case for a semimetal-insulator transition in graphene at strong Coulomb coupling. Finally, we highlight the need for systematic studies using improved lattice actions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Graphene: from materials science to particle physics

Since its discovery in 2004, graphene, a two-dimensional hexagonal carbon allotrope, has generated great interest and spurred research activity from materials science to particle physics and vice versa. In particular, graphene has been found to exhibit outstanding electronic and mechanical properties, as well as an unusual low-energy spectrum of Dirac quasiparticles giving rise to a fractional quantum Hall effect when freely suspended and immersed in a magnetic field. One of the most intriguing puzzles of graphene involves the low-temperature conductivity at zero density, a central issue in the design of graphene-based nanoelectronic components. While suspended graphene experiments have shown a trend reminiscent of semiconductors, with rising resistivity at low temperatures, most theories predict a constant or even decreasing resistivity. However, lattice field theory calculations have revealed that suspended graphene is at or near the critical coupling for excitonic gap formation due to strong Coulomb interactions, which suggests a simple and straightforward explanation for the experimental data. In this contribution we review the current status of the field with emphasis on the issue of gap formation, and outline recent progress and future points of contact between condensed matter physics and Lattice QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Preconditioning the non-relativistic many-fermion problem

Preconditioning is at the core of modern many-fermion Monte Carlo algorithms, such as Hybrid Monte Carlo, where the repeated solution of a linear problem involving an ill-conditioned matrix is needed. We report on a performance comparison of three preconditioning strategies, namely Chebyshev polynomials, strong-coupling approximation and weak-coupling expansion. We use conjugate gradient (CG) on the normal equations as well as stabilized biconjugate gradient (BiCGStab) as solvers and focus on the fermion matrix of the unitary Fermi gas. Our results indicate that BiCGStab is by far the most efficient strategy, both in terms of the number of iterations and matrix-vector operations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Signatures of a gap in the conductivity of graphene

The electrical conductivity of suspended graphene has recently been measured for the first time, and found to behave as σ~ \sqrt{|n|} as expected for Dirac quasiparticles at large carrier density. The charge inhomogeneity is strongly reduced in suspended samples, which revealed an unexpected insulating trend in σ(T). Above a transitional density n*, the temperature dependence was found to revert to metallic. We show that these features of the DC conductivity are consistent with a simple model of gapped Dirac quasiparticles, with specific signatures of a gap in the vicinity of the charge-neutrality point. Our findings are reminiscent of the conductivity profile in semiconducting materials, exhibiting a thermal activation for T \geq \tilde T and a weakly T-dependent background for T \leq \tilde T, where \tilde T is given by the saturation density \tilde n associated with the residual charge inhomogeneity. We discuss possible origins of a bandgap in graphene as well as alternative scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Toward large-scale Hybrid Monte Carlo simulations of the Hubbard model on graphics processing units

The performance of the Hybrid Monte Carlo algorithm is determined by the speed of sparse matrix-vector multiplication within the context of preconditioned conjugate gradient iteration. We study these operations as implemented for the fermion matrix of the Hubbard model in d+1 space-time dimensions, and report a performance comparison between a 2.66 GHz Intel Xeon E5430 CPU and an NVIDIA Tesla C1060 GPU using double-precision arithmetic. We find speedup factors ranging between 30-350 for d = 1, and in excess of 40 for d = 3. We argue that such speedups are of considerable impact for large-scale simulational studies of quantum many-body systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Critical exponents of the semimetal-insulator transition in graphene: A Monte Carlo study

The low-energy theory of graphene exhibits spontaneous chiral symmetry breaking due to pairing of quasiparticles and holes, corresponding to a semimetal-insulator transition at strong Coulomb coupling. We report a Lattice Monte Carlo study of the critical exponents of this transition as a function of the number of Dirac flavors $N_f^{}$, finding $δ= 1.25 \pm 0.05$ for $N_f^{} = 0$, $δ= 2.26 \pm 0.06$ for $N_f^{} = 2$ and $δ= 2.62 \pm 0.11$ for $N_f^{} = 4$, with $γ\simeq 1$ throughout. We compare our results with recent analytical work for graphene and closely related systems, and discuss scenarios for the fate of the chiral transition at finite temperature and carrier density, an issue of relevance for upcoming experiments with suspended graphene samples.

0 citations0 reviewsNo code linkedOpen access