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

Timo A. Lähde

Timo A. Lähde appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-latnucl-thhep-phcond-mat.str-elnucl-excond-mat.quant-gascond-mat.stat-mechhep-thastro-ph.COastro-ph.SRcond-mat.mes-hallphysics.comp-ph

Catalog footprint

What is connected

31works
12topics
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

31 published item(s)

preprint2021arXiv

Alpha-alpha scattering in the Multiverse

We investigate the phase shifts of low-energy alpha-alpha scattering under variations of the fundamental parameters of the Standard Model, namely the light quark mass, the electromagnetic fine-structure constant as well as the QCD theta-angle. As a first step, we recalculate alpha-alpha scattering in our Universe utilizing various improvements in the adiabatic projection method, which leads to an improved, parameter-free prediction of the S- and D-wave phase shifts for laboratory energies below 10 MeV. We find that positive shifts in the pion mass have a small effect on the S-wave phase shift, whereas lowering the pion mass adds some repulsion in the two-alpha system. The effect on the D-wave phase shift turns out to be more pronounced as signaled by the D-wave resonance parameters. Variations of the fine-structure constant have almost no effect on the low-energy alpha-alpha phase shifts. We further show that up-to-and-including next-to-leading order in the chiral expansion, variations of these phase shifts with respect to the QCD theta-angle can be expressed in terms of the theta-dependent pion mass.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Avoiding Ergodicity Problems in Lattice Discretizations of the Hubbard Model

The Hubbard model arises naturally when electron-electron interactions are added to the tight-binding descriptions of many condensed matter systems. For instance, the two-dimensional Hubbard model on the honeycomb lattice is central to the ab initio description of the electronic structure of carbon nanomaterials, such as graphene. Such low-dimensional Hubbard models are advantageously studied with Markov chain Monte Carlo methods, such as Hybrid Monte Carlo (HMC). HMC is the standard algorithm of the lattice gauge theory community, as it is well suited to theories of dynamical fermions. As HMC performs continuous, global updates of the lattice degrees of freedom, it provides superior scaling with system size relative to local updating methods. A potential drawback of HMC is its susceptibility to ergodicity problems due to so-called exceptional configurations, for which the fermion operator cannot be inverted. Recently, ergodicity problems were found in some formulations of HMC simulations of the Hubbard model. Here, we address this issue directly and clarify under what conditions ergodicity is maintained or violated in HMC simulations of the Hubbard model. We study different lattice formulations of the fermion operator and provide explicit, representative calculations for small systems, often comparing to exact results. We show that a fermion operator can be found which is both computationally convenient and free of ergodicity problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Semimetal-Mott Insulator Quantum Phase Transition of the Hubbard Model on the Honeycomb Lattice

We take advantage of recent improvements in the grand canonical Hybrid Monte Carlo algorithm, to perform a precision study of the single-particle gap in the hexagonal Hubbard model, with on-site electron-electron interactions. After carefully controlled analyses of the Trotter error, the thermodynamic limit, and finite-size scaling with inverse temperature, we find a critical coupling of $U_c/κ=3.834(14)$ and the critical exponent $zν=1.185(43)$. Under the assumption that this corresponds to the expected anti-ferromagnetic Mott transition, we are also able to provide a preliminary estimate $β=1.095(37)$ for the critical exponent of the order parameter. We consider our findings in view of the $SU(2)$ Gross-Neveu, or chiral Heisenberg, universality class. We also discuss the computational scaling of the Hybrid Monte Carlo algorithm, and possible extensions of our work to carbon nanotubes, fullerenes, and topological insulators.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An update on fine-tunings in the triple-alpha process

The triple-alpha process, whereby evolved stars create carbon and oxygen, is believed to be fine-tuned to a high degree. Such fine-tuning is suggested by the unusually strong temperature dependence of the triple-alpha reaction rate at stellar temperatures. This sensitivity is due to the resonant character of the triple-alpha process, which proceeds through the so-called "Hoyle state" of $^{12}$C with spin-parity $0^+$. The question of fine-tuning can be studied within the {\it ab initio} framework of nuclear lattice effective field theory, which makes it possible to relate {\it ad hoc} changes in the energy of the Hoyle state to changes in the fundamental parameters of the nuclear Hamiltonian, which are the light quark mass $m_q$ and the electromagnetic fine-structure constant. Here, we update the effective field theory calculation of the sensitivity of the triple-alpha process to small changes in the fundamental parameters. In particular, we consider recent high-precision lattice QCD calculations of the nucleon axial coupling $g_A$, as well as new and more comprehensive results from stellar simulations of the production of carbon and oxygen. While the updated stellar simulations allow for much larger {\it ad hoc} shifts in the Hoyle state energy than previously thought, recent lattice QCD results for the nucleon S-wave singlet and triplet scattering lengths now disfavor the scenario of no fine-tuning in the light quark mass $m_q$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nuclear binding near a quantum phase transition

How do protons and neutrons bind to form nuclei? This is the central question of ab initio nuclear structure theory. While the answer may seem as simple as the fact that nuclear forces are attractive, the full story is more complex and interesting. In this work we present numerical evidence from ab initio lattice simulations showing that nature is near a quantum phase transition, a zero-temperature transition driven by quantum fluctuations. Using lattice effective field theory, we perform Monte Carlo simulations for systems with up to twenty nucleons. For even and equal numbers of protons and neutrons, we discover a first-order transition at zero temperature from a Bose-condensed gas of alpha particles (4He nuclei) to a nuclear liquid. Whether one has an alpha-particle gas or nuclear liquid is determined by the strength of the alpha-alpha interactions, and we show that the alpha-alpha interactions depend on the strength and locality of the nucleon-nucleon interactions. This insight should be useful in improving calculations of nuclear structure and important astrophysical reactions involving alpha capture on nuclei. Our findings also provide a tool to probe the structure of alpha cluster states such as the Hoyle state responsible for the production of carbon in red giant stars and point to a connection between nuclear states and the universal physics of bosons at large scattering length.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Precise determination of lattice phase shifts and mixing angles

We introduce a general and accurate method for determining lattice phase shifts and mixing angles, which is applicable to arbitrary, non-cubic lattices. Our method combines angular momentum projection, spherical wall boundaries and an adjustable auxiliary potential. This allows us to construct radial lattice wave functions and to determine phase shifts at arbitrary energies. For coupled partial waves, we use a complex-valued auxiliary potential that breaks time-reversal invariance. We benchmark our method using a system of two spin-1/2 particles interacting through a finite-range potential with a strong tensor component. We are able to extract phase shifts and mixing angles for all angular momenta and energies, with precision greater than that of extant methods. We discuss a wide range of applications from nuclear lattice simulations to optical lattice experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantum Monte Carlo Calculations for Carbon Nanotubes

We show how lattice Quantum Monte Carlo can be applied to the electronic properties of carbon nanotubes in the presence of strong electron-electron correlations. We employ the path-integral formalism and use methods developed within the lattice QCD community for our numerical work. Our lattice Hamiltonian is closely related to the hexagonal Hubbard model augmented by a long-range electron-electron interaction. We apply our method to the single-quasiparticle spectrum of the (3,3) armchair nanotube configuration, and consider the effects of strong electron-electron correlations. Our approach is equally applicable to other nanotubes, as well as to other carbon nanostructures. We benchmark our Monte Carlo calculations against the two- and four-site Hubbard models, where a direct numerical solution is feasible.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ab initio alpha-alpha scattering

Processes involving alpha particles and alpha-like nuclei comprise a major part of stellar nucleosynthesis and hypothesized mechanisms for thermonuclear supernovae. In an effort towards understanding alpha processes from first principles, we describe in this letter the first ab initio calculation of alpha-alpha scattering. We use lattice effective field theory to describe the low-energy interactions of nucleons and apply a technique called the adiabatic projection method to reduce the eight-body system to an effective two-cluster system. We find good agreement between lattice results and experimental phase shifts for S-wave and D-wave scattering. The computational scaling with particle number suggests that alpha processes involving heavier nuclei are also within reach in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Breaking and restoration of rotational symmetry for irreducible tensor operators on the lattice

We study the breaking of rotational symmetry on the lattice for irreducible tensor operators and practical methods for suppressing this breaking. We illustrate the features of the general problem using an $α$ cluster model for $^{8}$Be. We focus on the lowest states with non-zero angular momentum and examine the matrix elements of multipole moment operators. We show that the physical reduced matrix element is well reproduced by averaging over all possible orientations of the quantum state, and this is expressed as a sum of matrix elements weighted by the corresponding Clebsch-Gordan coefficients. For our $α$ cluster model we find that the effects of rotational symmetry breaking can be largely eliminated for lattice spacings of $a\leq 1.7$ fm, and we expect similar improvement for actual lattice Monte Carlo calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nuclear Lattice Simulations using Symmetry-Sign Extrapolation

Projection Monte Carlo calculations of lattice Chiral Effective Field Theory suffer from sign oscillations to a varying degree dependent on the number of protons and neutrons. Hence, such studies have hitherto been concentrated on nuclei with equal numbers of protons and neutrons, and especially on the alpha nuclei where the sign oscillations are smallest. Here, we introduce the "symmetry-sign extrapolation" method, which allows us to use the approximate Wigner SU(4) symmetry of the nuclear interaction to systematically extend the Projection Monte Carlo calculations to nuclear systems where the sign problem is severe. We benchmark this method by calculating the ground-state energies of the $^{12}$C, $^6$He and $^6$Be nuclei, and discuss its potential for studies of neutron-rich halo nuclei and asymmetric nuclear matter.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ab initio calculation of the spectrum and structure of $^{16}$O

We present ab initio lattice calculations of the low-energy even-parity states of $^{16}$O using chiral nuclear effective field theory. We find good agreement with the empirical energy spectrum, and with the electromagnetic properties and transition rates. For the ground state, we find that the nucleons are arranged in a tetrahedral configuration of alpha clusters. For the first excited spin-0 state, we find that the predominant structure is a square configuration of alpha clusters, with rotational excitations that include the first spin-2 state.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Breaking and restoration of rotational symmetry on the lattice for bound state multiplets

We explore the breaking of rotational symmetry on the lattice for bound state energies and practical methods for suppressing this breaking. We demonstrate the general problems associated with lattice discretization errors and finite-volume errors using an $α$ cluster model for $^8$Be and $^{12}$C. We consider the two and three $α$-particle systems and focus on the lowest states with non-zero angular momentum which split into multiplets corresponding to different irreducible representations of the cubic group. We examine the dependence of such splittings on the lattice spacing and box size. We find that lattice spacing errors are closely related to the commensurability of the lattice with the intrinsic length scales of the system. We also show that rotational symmetry breaking effects can be significantly reduced by using improved lattice actions, and that the physical energy levels are accurately reproduced by the weighted average of a given spin multiplets.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Lattice Effective Field Theory for Medium-Mass Nuclei

We extend Nuclear Lattice Effective Field Theory (NLEFT) to medium-mass nuclei, and present results for the ground states of alpha nuclei from $^4$He to $^{28}$Si, calculated up to next-to-next-to-leading order (NNLO) in the EFT expansion. This computational advance is made possible by extrapolations of lattice data using multiple initial and final states. For our soft two-nucleon interaction, we find that the overall contribution from multi-nucleon forces must change sign from attractive to repulsive with increasing nucleon number. This effect is not produced by three-nucleon forces at NNLO, but it can be approximated by an effective four-nucleon interaction. We discuss the convergence of the EFT expansion and the broad significance of our findings for future ab initio calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Hoyle state in Nuclear Lattice EFT

We review the calculation of the Hoyle state of $^{12}$C in Nuclear Lattice Effective Field Theory (NLEFT) and its anthropic implications for the nucleosynthesis of $^{12}$C and $^{16}$O in red giant stars. We also review the extension of NLEFT to the regime of medium-mass nuclei, with emphasis on the determination of the ground-state energies of the alpha nuclei $^{16}$O, $^{20}$Ne, $^{24}$Mg and $^{28}$Si by means of Euclidean time projection. Finally, we review recent NLEFT results for the spectrum, electromagnetic properties, and alpha-cluster structure of $^{16}$O.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Uncertainties of Euclidean Time Extrapolation in Lattice Effective Field Theory

Extrapolations in Euclidean time form a central part of Nuclear Lattice Effective Field Theory (NLEFT) calculations using the Projection Monte Carlo method, as the sign problem in many cases prevents simulations at large Euclidean time. We review the next-to-next-to-leading order NLEFT results for the alpha nuclei up to $^{28}$Si, with emphasis on the Euclidean time extrapolations, their expected accuracy and potential pitfalls. We also discuss possible avenues for improving the reliability of Euclidean time extrapolations in NLEFT.

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
preprint2013arXiv

Dependence of the triple-alpha process on the fundamental constants of nature

We present an ab initio calculation of the quark mass dependence of the ground state energies of ^4He, ^8Be and ^{12}C, and of the energy of the Hoyle state in ^{12}C. These investigations are performed within the framework of lattice chiral Effective Field Theory. We address the sensitivity of the production rate of carbon and oxygen in red giant stars to the fundamental constants of nature by considering the impact of variations in the light quark masses and the electromagnetic fine-structure constant on the reaction rate of the triple-alpha process. As carbon and oxygen are essential to life as we know it, we also discuss the implications of our findings for an anthropic view of the Universe. We find strong evidence that the physics of the triple-alpha process is driven by alpha clustering, and that shifts in the fundamental parameters at the \simeq 2 - 3 % level are unlikely to be detrimental to the development of life. Tolerance against much larger changes cannot be ruled out at present, given the relatively limited knowledge of the quark mass dependence of the two-nucleon S-wave scattering parameters. Lattice QCD is expected to provide refined estimates of the scattering parameters in the future.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Two-loop Sunset Integrals at Finite Volume

We show how to compute the two-loop sunset integrals at finite volume, for non-degenerate masses and non-zero momentum. We present results for all integrals that appear in the Chiral Perturbation Therory ($χ$PT) calculation of the pseudoscalar meson masses and decay constants at NNLO, including the case of Partially Quenched $χ$PT. We also provide numerical implementations of the finite-volume sunset integrals, and review the results for one-loop integrals at finite volume.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Viability of carbon-based life as a function of the light quark mass

The Hoyle state plays a crucial role in the helium burning of stars that have reached the red giant stage. The close proximity of this state to the triple-alpha threshold is needed for the production of carbon, oxygen, and other elements necessary for life. We investigate whether this life-essential condition is robust or delicately fine-tuned by measuring its dependence on the fundamental constants of nature, specifically the light quark mass and the strength of the electromagnetic interaction. We show that there exist strong correlations between the alpha-particle binding energy and the various energies relevant to the triple-alpha process. We derive limits on the variation of these fundamental parameters from the requirement that sufficient amounts of carbon and oxygen be generated in stars. We also discuss the implications of these results for an anthropic view of the universe.

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

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
preprint2007arXiv

Initial-State Coulomb Interaction in the dd -> alpha pi^0 Reaction

The effects of initial-state Coulomb interactions in the charge-symmetry-breaking reaction dd -> alpha pi^0 are investigated within a previously published formalism. This is a leading order effect in which the Coulomb interaction between the two initial state protons leads to the breakup of the two deuterons into a continuum state that is well connected to the final alpha pi^0 state by the strong emission of a pion. As a first step, we use a simplified set of d and alpha wave functions and a plane-wave approximation for the initial dd state. This Coulomb mechanism, by itself, yields cross sections that are much larger than the experimental ones, and which are comparable in size to the contributions from other mechanisms. Inclusion of this mechanism is therefore necessary in a realistic calculation.

0 citations0 reviewsNo code linkedOpen access