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

Holger Fehske

Holger Fehske appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elcond-mat.mes-hallphysics.plasm-phDistributed, Parallel, and Cluster ComputingPerformancecond-mat.quant-gasMathematical Softwarecond-mat.othercond-mat.dis-nncond-mat.mtrl-sciNumerical Analysisphysics.comp-phphysics.opticsquant-phComputational Engineering, Finance, and Sciencemath.NAphysics.chem-ph

Catalog footprint

What is connected

62works
17topics
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

62 published item(s)

preprint2022arXiv

Nonequilibrium dynamics in pumped Mott insulators

We use time-evolution techniques for (infinite) matrix-product-states to calculate, directly in the thermodynamic limit, the time-dependent photoemission spectra and dynamic structure factors of the half-filled Hubbard chain after pulse irradiation. These quantities exhibit clear signatures of the photoinduced phase transition from insulator to metal that occurs because of the formation of so-called $η$ pairs. In addition, the spin dynamic structure factor loses spectral weight in the whole momentum space, reflecting the suppression of antiferromagnetic correlations due to the buildup of $η$-pairing states. The numerical method demonstrated in this work can be readily applied to other one-dimensional models driven out of equilibrium by optical pumping.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Photoinduced metallization of excitonic insulators

Utilizing the time-dependent density-matrix renormalization group technique, we numerically prove photoinduced pairing states in the extended Falicov-Kimball model (EFKM) at half filling, both with and without internal SU(2) symmetry. In the time-dependent photoemission spectra an extra band appears above the Fermi energy after pulse irradiation, indicating an insulator-to-metal transition. Even in the absence of the SU(2) structure, the pair correlations are enhanced during the pump, and afterwards decrease over time. This implies the possible metallization of Ta$_2$NiSe$_5$, a strong candidate for an excitonic insulator material, for which the EFKM is considered to be the minimal theoretical model. Simulating the photoemission with optimized pulse parameters, we demonstrate a photoinduced quantum phase transition, in accord with recent findings in time- and angle-resolved photoemission spectroscopy experiments on Ta$_2$NiSe$_5$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Controlling the direction of topological transport in a non-Hermitian time-reversal symmetric Floquet ladder

We propose a one-dimensional Floquet ladder that possesses two distinct topological transport channels with opposite directionality. The transport channels occur due to a $\mathbb Z_2$ non-Hermitian Floquet topological phase that is protected by time-reversal symmetry. The signatures of this phase are two pairs of Kramers degenerate Floquet quasienergy bands that are separated by an imaginary gap. We discuss how the Floquet ladder can be implemented in a photonic waveguide lattice and show that the direction of transport in the resulting waveguide structure can be externally controlled by focusing two light beams into adjacent waveguides. The relative phase between the two light beams selects which of the two transport channels is predominantly populated, while the angles of incidence of the two light beams determine which of the transport channels is suppressed by non-Hermitian losses. We identify the optimal lattice parameters for the external control of transport and demonstrate the robustness of this mechanism against disorder.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Scrutinizing the Debye plasma model: Rydberg excitons unravel the properties of low-density plasmas in semiconductors

For low-density plasmas, the classical limit described by the Debye-Hückel theory is still considered as an appropriate description even though a clear experimental proof of this paradigm is lacking due to the problems in determining the plasma-induced shift of single-particle energies in atomic systems. We show that Rydberg excitons in states with a high principal quantum number are highly sensitive probes for their surrounding making it possible to unravel accurately the basic properties of low-density non-degenerate electron-hole plasmas. To this end, we accurately measure the parameters of Rydberg excitons such as energies and linewidths in absorption spectra of bulk cuprous oxide crystals in which a tailored electron-hole plasma has been generated optically. Since from the absorption spectra exciton energies, as well as the shift of the single-particle energies given by the band edge, can be directly derived, the measurements allow us to determine the plasma density and temperature independently, which has been a notoriously hard problem in semiconductor physics. Our analysis shows unambiguously that the impact of the plasma cannot be described by the classical Debye model, but requires a quantum many-body theory, not only for the semiconductor plasma investigated here, but in general. Furthermore, it reveals a new exciton scattering mechanism with coupled plasmon-phonon modes becoming important even at very low plasma densities.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cutting off the non-Hermitian boundary from an anomalous Floquet topological insulator

In two-dimensional anomalous Floquet insulators, non-Hermitian boundary state engineering can be used to completely separate chiral boundary states from bulk bands in the quasienergy spectrum. The topological properties of such spectrally separated boundary states are no longer restricted by the strict bulk-boundary correspondence of Hermitian systems. We show that this additional topological freedom enables one to faithfully transfer the topological properties of a boundary attached to a Floquet insulator to a non-Hermitian Floquet chain obtained by physically cutting off the boundary from the bulk. We implement this scenario for a simple model of an anomalous Floquet insulator with Hermitian and non-Hermitian boundaries, and discuss the relevance of our construction for the experimental realization of non-Hermitian topological phases that connect dimensions one and two.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electronic Properties of $α-\mathcal{T}_3$ Quantum Dots in Magnetic Fields

We address the electronic properties of quantum dots in the two-dimensional $α-\mathcal{T}_3$ lattice when subjected to a perpendicular magnetic field. Implementing an infinite mass boundary condition, we first solve the eigenvalue problem for an isolated quantum dot in the low-energy, long-wavelength approximation where the system is described by an effective Dirac-like Hamiltonian that interpolates between the graphene (pseudospin 1/2) and Dice (pseudospin 1) limits. Results are compared to a full numerical (finite-mass) tight-binding lattice calculation. In a second step we analyse charge transport through a contacted $α-\mathcal{T}_3$ quantum dot in a magnetic field by calculating the local density of states and the conductance within the kernel polynomial and Landauer-Büttiker approaches. Thereby the influence of a disordered environment is discussed as well.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Photoinduced $η$-pairing at finite temperatures

We numerically prove photoinduced $η$-pairing in a half-filled fermionic Hubbard chain at both zero and finite temperature. The result, obtained by combining the matrix-product-state based infinite time-evolving block decimation technique and the purification method, applies to the thermodynamic limit. Exciting the Mott insulator by a laser electric field docked on via the Peierls phase, we track the time-evolution of the correlated many-body system and determine the optimal parameter set for which the nonlocal part of the $η$-pair correlation function becomes dominant during the laser pump at zero and low temperatures. These correlations vanish at higher temperatures and long times after pulse irradiation. In the high laser frequency strong Coulomb coupling regime we observe a remnant enhancement of the Brillouin-zone boundary pair-correlation function also at high temperatures, if the Hubbard interaction is about a multiple of the laser frequency, which can be attributed to an enhanced double occupancy in the virtual Floquet state.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Photoinduced $η$-pairing in One-dimensional Mott Insulators

Employing the density-matrix renormalization group technique in the matrix-product-state representation, we investigate the photoexcited superconducting correlations induced by the $η$-pairing mechanism in the half-filled Hubbard chain. We estimate the characteristic pair correlation function and verify the accuracy of our numerical results by comparison with exact-diagonalization data for small systems. The optimal parameter set of the pump that most enhances the $η$-pair correlations, is calculated in the strong-coupling regime. For such a pump, we explore the possibility of quasi-long-range order.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Real and imaginary edge states in stacked Floquet honeycomb lattices

We present a non-Hermitian Floquet model with topological edge states in real and imaginary band gaps. The model utilizes two stacked honeycomb lattices which can be related via four different types of non-Hermitian time-reversal symmetry. Implementing the correct time-reversal symmetry provides us with either two counterpropagating edge states in a real gap, or a single edge state in an imaginary gap. The counterpropagating edge states allow for either helical or chiral transport along the lattice perimeter. In stark contrast, we find that the edge state in the imaginary gap does not propagate. Instead, it remains spatially localized while its amplitude continuously increases. Our model is well-suited for realizing these edge states in photonic waveguide lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Resonant scattering of Dice quasiparticles on oscillating quantum dots

We consider a Dice model with Dirac cones intersected by a topologically flat band at the charge neutrality point and analyze the inelastic scattering of massless pseudospin-1 particles on a circular, gate-defined, oscillating barrier. Focusing on the resonant scattering regime at small energy of the incident wave, we calculate the reflection and transmission coefficients and derive explicit expressions for the time-dependent particle probability, current density and scattering efficiency within (Floquet) Dirac-Weyl theory, both in the near-field and the far-field. We discuss the importance of sideband scattering and Fano resonances in the quantum limit. When resonance conditions are fulfilled, the particle is temporarily trapped in vortices located close to edge of the quantum dot before it gets resubmitted with strong angular dependence. Interestingly even periodically alternating forward and backward radiation may occur. We also demonstrate the revival of resonant scattering related to specific fusiform boundary trapping profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Topological origin of quantized transport in non-Hermitian Floquet chains

We show that non-Hermiticity enables topological phases with unidirectional transport in one-dimensional Floquet chains. The topological signatures of these phases are non-contractible loops in the spectrum of the Floquet propagator that are separated by an imaginary gap. Such loops occur exclusively in non-Hermitian Floquet systems. We define the corresponding topological invariant as the winding number of the Floquet propagator relative to the imaginary gap. To relate topology to transport, we introduce the concept of regularized dynamics of non-Hermitian chains. We establish that, under the conditions of regularized dynamics, transport is quantized in so far as the charge transferred over one period equals the topological winding number. We illustrate these theoretical findings with the example of a Floquet chain that features a topological phase transition and acts as a charge pump in the non-trivial topological phase. We finally discuss whether these findings justify the notion that non-Hermitian Floquet chains support topological transport.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Understanding HPC Benchmark Performance on Intel Broadwell and Cascade Lake Processors

Hardware platforms in high performance computing are constantly getting more complex to handle even when considering multicore CPUs alone. Numerous features and configuration options in the hardware and the software environment that are relevant for performance are not even known to most application users or developers. Microbenchmarks, i.e., simple codes that fathom a particular aspect of the hardware, can help to shed light on such issues, but only if they are well understood and if the results can be reconciled with known facts or performance models. The insight gained from microbenchmarks may then be applied to real applications for performance analysis or optimization. In this paper we investigate two modern Intel x86 server CPU architectures in depth: Broadwell EP and Cascade Lake SP. We highlight relevant hardware configuration settings that can have a decisive impact on code performance and show how to properly measure on-chip and off-chip data transfer bandwidths. The new victim L3 cache of Cascade Lake and its advanced replacement policy receive due attention. Finally we use DGEMM, sparse matrix-vector multiplication, and the HPCG benchmark to make a connection to relevant application scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Chebyshev Filter Diagonalization on Modern Manycore Processors and GPGPUs

Chebyshev filter diagonalization is well established in quantum chemistry and quantum physics to compute bulks of eigenvalues of large sparse matrices. Choosing a block vector implementation, we investigate optimization opportunities on the new class of high-performance compute devices featuring both high-bandwidth and low-bandwidth memory. We focus on the transparent access to the full address space supported by both architectures under consideration: Intel Xeon Phi "Knights Landing" and Nvidia "Pascal." We propose two optimizations: (1) Subspace blocking is applied for improved performance and data access efficiency. We also show that it allows transparently handling problems much larger than the high-bandwidth memory without significant performance penalties. (2) Pipelining of communication and computation phases of successive subspaces is implemented to hide communication costs without extra memory traffic. As an application scenario we use filter diagonalization studies on topological insulator materials. Performance numbers on up to 512 nodes of the OakForest-PACS and Piz Daint supercomputers are presented, achieving beyond 100 Tflop/s for computing 100 inner eigenvalues of sparse matrices of dimension one billion.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Fermionic time-reversal symmetry in a photonic topological insulator

Much of the recent enthusiasm directed towards topological insulators as a new state of matter is motivated by their hallmark feature of protected chiral edge states. In fermionic systems, Kramers degeneracy gives rise to these entities in the presence of time-reversal symmetry (TRS). In contrast, bosonic systems obeying TRS are generally assumed to be fundamentally precluded from supporting edge states. In this work, we dispel this perception and experimentally demonstrate counter-propagating chiral states at the edge of a time-reversal-symmetric photonic waveguide structure. The pivotal step in our approach is encoding the effective spin of the propagating states as a degree of freedom of the underlying waveguide lattice, such that our photonic topological insulator is characterised by a $\mathbb{Z}_2$-type invariant. Our findings allow for fermionic properties to be harnessed in bosonic systems, thereby opening new avenues for topological physics in photonics as well as acoustics, mechanics and even matter waves.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Electron confinement in graphene with gate-defined quantum dots

We theoretically analyse the possibility to electrostatically confine electrons in circular quantum dot arrays, impressed on contacted graphene nanoribbons by top gates. Utilising exact numerical techniques, we compute the scattering efficiency of a single dot and demonstrate that for small-sized scatterers the cross-sections are dominated by quantum effects, where resonant scattering leads to a series of quasi-bound dot states. Calculating the conductance and the local density of states for quantum dot superlattices we show that the resonant carrier transport through such graphene-based nanostructures can be easily tuned by varying the gate voltage.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

GHOST: Building blocks for high performance sparse linear algebra on heterogeneous systems

While many of the architectural details of future exascale-class high performance computer systems are still a matter of intense research, there appears to be a general consensus that they will be strongly heterogeneous, featuring "standard" as well as "accelerated" resources. Today, such resources are available as multicore processors, graphics processing units (GPUs), and other accelerators such as the Intel Xeon Phi. Any software infrastructure that claims usefulness for such environments must be able to meet their inherent challenges: massive multi-level parallelism, topology, asynchronicity, and abstraction. The "General, Hybrid, and Optimized Sparse Toolkit" (GHOST) is a collection of building blocks that targets algorithms dealing with sparse matrix representations on current and future large-scale systems. It implements the "MPI+X" paradigm, has a pure C interface, and provides hybrid-parallel numerical kernels, intelligent resource management, and truly heterogeneous parallelism for multicore CPUs, Nvidia GPUs, and the Intel Xeon Phi. We describe the details of its design with respect to the challenges posed by modern heterogeneous supercomputers and recent algorithmic developments. Implementation details which are indispensable for achieving high efficiency are pointed out and their necessity is justified by performance measurements or predictions based on performance models. The library code and several applications are available as open source. We also provide instructions on how to make use of GHOST in existing software packages, together with a case study which demonstrates the applicability and performance of GHOST as a component within a larger software stack.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ground-state and spectral signatures of cavity exciton-polariton condensates

We propose a projector-based renormalization framework to study exciton-polariton Bose-Einstein condensation in a microcavity matter-light system. Treating Coulomb interaction and electron-hole/photon coupling effects on an equal footing we analyze the ground-state properties of the exciton polariton model according to the detuning and the excitation density. We demonstrate that the condensate by its nature shows a crossover from an excitonic insulator (of Bose-Einstein respectively BCS type) to a polariton and finally photonic condensed state as the excitation density increases at large detuning. If the detuning is weak polariton or photonic phases dominate. While in both cases a notable renormalization of the quasiparticle band structure occurs that strongly affects the coherent part of the excitonic luminescence, the incoherent wavevector-resolved luminescence spectrum develops a flat bottom only for small detuning.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

High-performance implementation of Chebyshev filter diagonalization for interior eigenvalue computations

We study Chebyshev filter diagonalization as a tool for the computation of many interior eigenvalues of very large sparse symmetric matrices. In this technique the subspace projection onto the target space of wanted eigenvectors is approximated with filter polynomials obtained from Chebyshev expansions of window functions. After the discussion of the conceptual foundations of Chebyshev filter diagonalization we analyze the impact of the choice of the damping kernel, search space size, and filter polynomial degree on the computational accuracy and effort, before we describe the necessary steps towards a parallel high-performance implementation. Because Chebyshev filter diagonalization avoids the need for matrix inversion it can deal with matrices and problem sizes that are presently not accessible with rational function methods based on direct or iterative linear solvers. To demonstrate the potential of Chebyshev filter diagonalization for large-scale problems of this kind we include as an example the computation of the $10^2$ innermost eigenpairs of a topological insulator matrix with dimension $10^9$ derived from quantum physics applications.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Topological insulators in random potentials

We investigate the effects of magnetic and nonmagnetic impurities on the two-dimensional surface states of three-dimensional topological insulators (TIs). Modeling weak and strong TIs using a generic four-band Hamiltonian, which allows for a breaking of inversion and time-reversal symmetries and takes into account random local potentials as well as the Zeeman and orbital effects of external magnetic fields, we compute the local density of states, the single-particle spectral function, and the conductance for a (contacted) slab geometry by numerically exact techniques based on kernel polynomial expansion and Green's function approaches. We show that bulk disorder refills the suface-state Dirac gap induced by a homogeneous magnetic field with states, whereas orbital (Peierls-phase) disorder perserves the gap feature. The former effect is more pronounced in weak TIs than in strong TIs. At moderate randomness, disorder-induced conducting channels appear in the surface layer, promoting diffusive metallicity. Random Zeeman fields rapidly destroy any conducting surface states. Imprinting quantum dots on a TI's surface, we demonstrate that carrier transport can be easily tuned by varying the gate voltage, even to the point where quasi-bound dot states may appear.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Unconventional superconductivity and interaction induced Fermi surface reconstruction in the two-dimensional Edwards model

We study the competition between unconventional superconducting pairing and charge density wave (CDW) formation for the two-dimensional Edwards Hamiltonian at half filling, a very general two-dimensional transport model in which fermionic charge carriers couple to a correlated background medium. Using the projective renormalization method we find that a strong renormalization of the original fermionic band causes a new hole-like Fermi surface to emerge near the center of the Brillouin zone, before it eventually gives rise to the formation of a charge density wave. On the new, disconnected parts of the Fermi surface superconductivity is induced with a sign-changing order parameter. We discuss these findings in the light of recent experiments on iron-based oxypnictide superconductors.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Absorption of an electron by a dielectric wall

We introduce a method for calculating the probability with which a low-energy electron hitting the wall of a plasma gets stuck in it and apply the method to a dielectric wall with positive electron affinity smaller than the bandgap using MgO as an example. In accordance with electron beam scattering data we obtain energy-dependent sticking probabilities significantly less than unity and question thereby for electrons the perfect absorber assumption unisono used in plasma modeling.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Comparative density-matrix renormalization group study of symmetry-protected topological phases in spin-1 chain and Bose-Hubbard models

We reexamine the one-dimensional spin-1 $XXZ$ model with on-site uniaxial single-ion anisotropy as to the appearance and characterization of the symmetry-protected topological Haldane phase. By means of large-scale density-matrix renormalization group (DMRG) calculations the central charge can be determined numerically via the von Neumann entropy, from which the ground-sate phase diagram of the model can be derived with high precision. The nontrivial gapped Haldane phase shows up in between the trivial gapped even Haldane and N{é}el phases, appearing at large single-ion and spin--exchange interaction anisotropies, respectively. We furthermore carve out a characteristic degeneracy of the lowest entanglement level in the topological Haldane phase, which is determined using a conventional finite-system DMRG technique with both periodic and open boundary conditions. Defining the spin and neutral gaps in analogy to the single-particle and neutral gaps in the intimately connected extended Bose-Hubbard model, we show that the excitation gaps in the spin model qualitatively behave just as for the bosonic system. We finally compute the dynamical spin structure factor in the three different gapped phases and find significant differences in the intensity maximum which might be used to distinguish these phases experimentally.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Competition between excitonic charge and spin density waves: Influence of electron-phonon and Hund's rule couplings

We analyze the stability of excitonic ground states in the two-band Hubbard model with additional electron-phonon and Hund's rule couplings using a combination of mean-field and variational cluster approaches. We show that both the interband Coulomb interaction and the electron-phonon interaction will cooperatively stabilize a charge density wave (CDW) state which typifies an "excitonic" CDW if predominantly triggered by the effective interorbital electron-hole attraction or a "phononic" CDW if mostly caused by the coupling to the lattice degrees of freedom. By contrast, the Hund's rule coupling promotes an excitonic spin density wave. We determine the transition between excitonic charge and spin density waves and comment on a fixation of the phase of the excitonic order parameter that would prevent the formation of a superfluid condensate of excitons. The implications for exciton condensation in several material classes with strongly correlated electrons are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Criticality at the Haldane-insulator charge-density-wave quantum phase transition

Exploiting the entanglement concept within a matrix-product-state based infinite density-matrix renormalization group approach, we show that the spin-density-wave and bond-order-wave ground states of the one-dimensional half-filled extended Hubbard model give way to a symmetry-protected topological Haldane state in case an additional alternating ferromagnetic spin interaction is added. In the Haldane insulator the lowest entanglement level features a characteristic twofold degeneracy. Increasing the ratio between nearest-neighbor and local Coulomb interaction $V/U$, the enhancement of the entanglement entropy, the variation of the charge, spin and neutral gaps, and the dynamical spin/density response signal a quantum phase transition to a charge-ordered state. Below a critical point, which belongs to the universality class of the tricritical Ising model with central charge 7/10, the model is critical with $c=1/2$ along the transition line. Above this point, the transition between the Haldane insulator and charge-density-wave phases becomes first order.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Entangled light from driven dissipative microcavities

We study the generation of entangled light in planar semiconductor microcavities. The focus is on a particular pump configuration where the dissipative internal polariton dynamics leads to the emission of entangled light in a W-state. Our study is based on the nonlinear equations of motion for the polariton operators derived within the dynamics-controlled truncation formalism. They include the losses through the cavity mirrors, the interaction with lattice vibrations, and the external laser driving in a Langevin approach. We find that the generated entanglement is robust against decoherence under realistic experimental conditions. The results show that pair correlations in solid-state devices can be used to stabilize the nonlocal properties of the emitted radiation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Performance Engineering of the Kernel Polynomial Method on Large-Scale CPU-GPU Systems

The Kernel Polynomial Method (KPM) is a well-established scheme in quantum physics and quantum chemistry to determine the eigenvalue density and spectral properties of large sparse matrices. In this work we demonstrate the high optimization potential and feasibility of peta-scale heterogeneous CPU-GPU implementations of the KPM. At the node level we show that it is possible to decouple the sparse matrix problem posed by KPM from main memory bandwidth both on CPU and GPU. To alleviate the effects of scattered data access we combine loosely coupled outer iterations with tightly coupled block sparse matrix multiple vector operations, which enables pure data streaming. All optimizations are guided by a performance analysis and modelling process that indicates how the computational bottlenecks change with each optimization step. Finally we use the optimized node-level KPM with a hybrid-parallel framework to perform large scale heterogeneous electronic structure calculations for novel topological materials on a petascale-class Cray XC30 system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A unified sparse matrix data format for efficient general sparse matrix-vector multiply on modern processors with wide SIMD units

Sparse matrix-vector multiplication (spMVM) is the most time-consuming kernel in many numerical algorithms and has been studied extensively on all modern processor and accelerator architectures. However, the optimal sparse matrix data storage format is highly hardware-specific, which could become an obstacle when using heterogeneous systems. Also, it is as yet unclear how the wide single instruction multiple data (SIMD) units in current multi- and many-core processors should be used most efficiently if there is no structure in the sparsity pattern of the matrix. We suggest SELL-C-sigma, a variant of Sliced ELLPACK, as a SIMD-friendly data format which combines long-standing ideas from General Purpose Graphics Processing Units (GPGPUs) and vector computer programming. We discuss the advantages of SELL-C-sigma compared to established formats like Compressed Row Storage (CRS) and ELLPACK and show its suitability on a variety of hardware platforms (Intel Sandy Bridge, Intel Xeon Phi and Nvidia Tesla K20) for a wide range of test matrices from different application areas. Using appropriate performance models we develop deep insight into the data transfer properties of the SELL-C-sigma spMVM kernel. SELL-C-sigma comes with two tuning parameters whose performance impact across the range of test matrices is studied and for which reasonable choices are proposed. This leads to a hardware-independent ("catch-all") sparse matrix format, which achieves very high efficiency for all test matrices across all hardware platforms.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Infrared light extinction by charged dielectric core-coat particles

We study the effect of surplus electrons on the infrared extinction of dielectric particles with a core-coat structure and propose to use it for an optical measurement of the particle charge in a dusty plasma. The particles consist of an inner core with negative and an outer coat with positive electron affinity. Both the core and the coat give rise to strong transverse optical phonon resonances, leading to anomalous light scattering in the infrared. Due to the radial profile of the electron affinity electrons accumulate in the coat region making the infrared extinction of this type of particles very charge-sensitive, in particular, the extinction due to a resonance arising solely due to the core-coat structure. The maximum of this resonance is in the far-infrared and responds to particle charges realizable in ordinary dusty laboratory plasmas.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Linear response within the projection-based renormalization method: Many-body corrections beyond the random phase approximation

The explicit evaluation of linear response coefficients for interacting many-particle systems still poses a considerable challenge to theoreticians. In this work we use a novel many-particle renormalization technique, the so-called projector-based renormalization method, to show how such coefficients can systematically be evaluated. To demonstrate the prospects and power of our approach we consider the dynamical wave-vector dependent spin susceptibility of the two-dimensional Hubbard model and also determine the subsequent magnetic phase diagram close to half-filling. We show that the superior treatment of (Coulomb) correlation and fluctuation effects within the projector-based renormalization method significantly improves the standard random phase approximation results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spectral and Entanglement Properties of the Bosonic Haldane Insulator

We discuss the existence of a nontrivial topological phase in one-dimensional interacting systems described by the extended Bose-Hubbard model with a mean filling of one boson per site. Performing large-scale density-matrix renormalization group calculations we show that the presence of nearest-neighbor repulsion enriches the ground-state phase diagram of the paradigmatic Bose-Hubbard model by stabilizing a novel gapped insulating state, the so-called Haldane insulator, which, embedded into superfluid, Mott insulator, and density wave phases, is protected by the lattice inversion symmetry. The quantum phase transitions between the different insulating phases were determined from the central charge via the von Neumann entropy. The Haldane phase reveals a characteristic fourfold degeneracy of the entanglement spectrum. We finally demonstrate that the intensity maximum of the dynamical charge structure factor, accessible by Bragg spectroscopy, features the gapped dispersion known from the spin-1 Heisenberg chain.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Charge-density-wave formation in the Edwards fermion-boson model at one-third band filling

We examine the ground-state properties of the one-dimensional Edwards spinless fermion transport model by means of large-scale density-matrix renormalization-group calculations. Determining the single-particle gap and the Tomonaga-Luttinger liquid parameter ($K_ρ$) at zero temperature, we prove the existence of a metal-to-insulator quantum phase transition at one-third band filling. The insulator---established by strong correlation in the background medium---typifies a charge density wave (CDW) that is commensurate with the band filling. $K_ρ=2/9$ is very small at the quantum critical point, and becomes $K_ρ^{\rm CDW}=1/9$ in the infinitesimally doped three-period CDW, as predicted by the bosonization approach.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Surface electrons at plasma walls

In this chapter we introduce a microscopic modelling of the surplus electrons on the plasma wall which complements the classical description of the plasma sheath. First we introduce a model for the electron surface layer to study the quasistationary electron distribution and the potential at an unbiased plasma wall. Then we calculate sticking coefficients and desorption times for electron trapping in the image states. Finally we study how surplus electrons affect light scattering and how charge signatures offer the possibility of a novel charge measurement for dust grains.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Characterization of Mott-insulating and superfluid phases in the one-dimensional Bose--Hubbard model

We use strong-coupling perturbation theory, the variational cluster approach (VCA), and the dynamical density-matrix renormalization group (DDMRG) method to investigate static and dynamical properties of the one-dimensional Bose--Hubbard model in both the Mott-insulating and superfluid phases. From the von Neumann entanglement entropy we determine the central charge and the transition points for the first two Mott lobes. Our DMRG results for the ground-state energy, momentum distribution function, boson correlation function decay, Mott gap, and single particle-spectral function are reproduced very well by the strong-coupling expansion to fifth order, and by VCA with clusters up to 12 sites as long as the ratio between the hopping amplitude and on-site repulsion, t/U, is smaller than 0.15 and 0.25, respectively. In addition, in the superfluid phase VCA captures well the ground-state energy and the sound velocity of the linear phonon modes. This comparison provides an authoritative estimate for the range of applicability of these methods. In strong-coupling theory for the Mott phase, the dynamical structure factor is obtained from the solution of an effective single-particle problem with an attractive potential. The resulting resonances show up as double-peak structure close to the Brillouin zone boundary. These high-energy features also appear in the superfluid phase which is characterized by a pronounced phonon mode at small momenta and energies, as predicted by Bogoliubov and field theory. In one dimension, there are no traces of an amplitude mode in the dynamical single-particle and two-particle correlation functions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Charge carrier localisation in disordered graphene nanoribbons

We study the electronic properties of actual-size graphene nanoribbons subjected to substitutional disorder particularly with regard to the experimentally observed metal-insulator transition. Calculating the local, mean and typical density of states, as well as the time-evolution of the particle density we comment on a possible disorder-induced localisation of charge carriers at and close to the Dirac point within a percolation transition scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Condensation of Excitons in Cu2O at Ultracold Temperatures: Experiment and Theory

We present experiments on the luminescence of excitons confined in a potential trap at milli-Kelvin bath temperatures under cw-excitation. They reveal several distinct features like a kink in the dependence of the total integrated luminescence intensity on excitation laser power and a bimodal distribution of the spatially resolved luminescence. Furthermore, we discuss the present state of the theoretical description of Bose-Einstein condensation of excitons with respect to signatures of a condensate in the luminescence. The comparison of the experimental data with theoretical results with respect to the spatially resolved as well as the integrated luminescence intensity shows the necessity of taking into account a Bose-Einstein condensed excitonic phase in order to understand the behaviour of the trapped excitons.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Coulomb interaction effects in graphene bilayers: electron-hole pairing and plasmaron formation

We report a theoretical study of the many-body effects of electron-electron interaction on the ground-state and spectral properties of double-layer graphene. Using a projector-based renormalization method we show that if a finite voltage difference is applied between the graphene layers electron-hole pairs can be formed and---at very low temperatures---an excitonic instability might emerge in a double-layer graphene structure. The single-particle spectral function near the Fermi surface exhibits a prominent quasiparticle peak, different from neutral (undoped) graphene bilayers. Away from the Fermi surface, we find that the charge carriers strongly interact with plasmons, thereby giving rise to a broad plasmaron peak in the angle-resolved photoemission spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effect of Electron-Phonon Interaction Range for a Half-Filled Band in One Dimension

We demonstrate that fermion-boson models with nonlocal interactions can be simulated at finite band filling with the continuous-time quantum Monte Carlo method. We apply this method to explore the influence of the electron-phonon interaction range for a half-filled band in one dimension, covering the full range from the Holstein to the Fröhlich regime. The phase diagram contains metallic, Peierls, and phase-separated regions, which we characterize in terms of static and dynamical correlation functions. In particular, our results reveal a suppression of $2k_F$ charge correlations with increasing interaction range, allowing for a power-law decay comparable to the pairing correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Electron surface layer at the interface of a plasma and a dielectric wall

We study the potential and the charge distribution across the interface of a plasma and a dielectric wall. For this purpose, the charge bound to the wall is modelled as a quasi-stationary electron surface layer which satisfies Poisson's equation and minimizes the grand canonical potential of the wall-thermalized excess electrons constituting the wall charge. Based on an effective model for a graded interface taking into account the image potential and the offset of the conduction band to the potential just outside the dielectric, we specifically calculate the potential and the electron distribution for magnesium oxide, silicon dioxide and sapphire surfaces in contact with a helium discharge. Depending on the electron affinity of the surface, we find two vastly different behaviors. For negative electron affinity, electrons do not penetrate into the wall and an external surface charge is formed in the image potential, while for positive electron affinity, electrons penetrate into the wall and a space charge layer develops in the interior of the dielectric. We also investigate how the electron surface layer merges with the bulk of the dielectric.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pseudo-particle approach for charge-transferring molecule-surface collisions

Based on a semi-empirical generalized Anderson-Newns model we construct a pseudo-particle description for electron emission due to de-excitation of metastable molecules at surfaces. The pseudo-particle approach allows us to treat resonant charge-transfer and Auger processes on an equal footing, as it is necessary when both channels are open. This is for instance the case when a metastable nitrogen molecule hits a diamond surface. Using non-equilibrium Green functions and physically motivated approximations to the self-energies of the Dyson equations we derive a system of rate equations for the probabilities with which the metastable molecule, the molecular ground state, and the negative ion can be found in the course of the scattering event. From the rate equations we also obtain the spectrum of the emitted electron and the secondary electron emission coefficient. Our numerical results indicate the resonant tunneling process undermining the source of the Auger channel which therefore contributes only a few percent to the secondary electron emission.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Resonant charge transfer at dielectric surfaces

We report on the theoretical description of secondary electron emission due to resonant charge transfer occurring during the collision of metastable nitrogen molecules with dielectric surfaces. The emission is described as a two step process consisting of electron capture to form an intermediate shape resonance and subsequent electron emission by decay of this ion, either due to its natural life time or its interaction with the surface. The electron capture is modeled using the Keldysh Green's function technique and the negative ion decay is described by a combination of the Keldysh technique and a rate equation approach. We find the resonant capture of electrons to be very efficient and the natural decay to be clearly dominating over the surface-induced decay. Secondary electron emission coefficients are calculated for aluminum oxide, magnesium oxide, silicon oxide, and diamond at several kinetic energies of the projectile. With the exception of magnesium oxide the coefficients turn out to be of the order of 0.1 over the whole range of kinetic energies. This rather large value is a direct consequence of the shape resonance acting as a relay state for electron emission.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sparse matrix-vector multiplication on GPGPU clusters: A new storage format and a scalable implementation

Sparse matrix-vector multiplication (spMVM) is the dominant operation in many sparse solvers. We investigate performance properties of spMVM with matrices of various sparsity patterns on the nVidia "Fermi" class of GPGPUs. A new "padded jagged diagonals storage" (pJDS) format is proposed which may substantially reduce the memory overhead intrinsic to the widespread ELLPACK-R scheme. In our test scenarios the pJDS format cuts the overall spMVM memory footprint on the GPGPU by up to 70%, and achieves 95% to 130% of the ELLPACK-R performance. Using a suitable performance model we identify performance bottlenecks on the node level that invalidate some types of matrix structures for efficient multi-GPGPU parallelization. For appropriate sparsity patterns we extend previous work on distributed-memory parallel spMVM to demonstrate a scalable hybrid MPI-GPGPU code, achieving efficient overlap of communication and computation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The fate of topological-insulator surface states under strong disorder

Three-dimensional topological insulators feature Dirac-like surface states which are topologically protected against the influence of weak quenched disorder. Here we investigate the effect of surface disorder beyond the weak-disorder limit using large-scale numerical simulations. We find two qualitatively distinct regimes: Moderate disorder destroys the Dirac cone and induces diffusive metallic behavior at the surface. Even more remarkably, for strong surface disorder a Dirac cone reappears, as new weakly disordered "surface" states emerge in the sample beneath the disordered surface layer, which can be understood in terms of an interface between a topological and an Anderson insulator. Together, this demonstrates the drastic effect of disorder on topological surface states, which cannot be captured within effective two-dimensional models for the surface states alone.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Aharonov-Bohm interference for a hole in a two-dimensional Ising antiferromagnet in a transverse magnetic field

We show that a proper consideration of the contribution of Trugman loops leads to a fairly low effective mass for a hole moving in a square lattice Ising antiferromagnet, if the bare hopping and the exchange energy scales are comparable. This contradicts the general view that because of the absence of spin fluctuations, this effective mass must be extremely large. Moreover, in the presence of a transverse magnetic field, we show that the effective hopping integrals acquire an unusual dependence on the magnetic field, through Aharonov-Bohm interference, in addition to significant retardation effects. The effect of the Aharonov-Bohm interference on the cyclotron frequency (for small magnetic fields) and the Hofstadter butterfly (for large magnetic fields) is analyzed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Auger de-excitation of metastable molecules at metallic surfaces

We study secondary electron emission from metallic surfaces due to Auger de-excitation of diatomic metastable molecules. Our approach is based on an effective model for the two active electrons involved in the process -- a molecular electron described by a linear combination of atomic orbitals when it is bound and a two-center Coulomb wave when it is not and a metal electron described by the eigenfunctions of a step potential -- and employs Keldysh Green's functions. Solving the Dyson equation for the retarded Green's function by exponential resummation we are able to treat time-nonlocal self-energies and to avoid the wide-band approximation.Results are presented for the de-excitation of \NitrogenDominantMetastableState\ on aluminum and tungsten and discussed in view of previous experimental and theoretical investigations. We find quantitative agreement with experimental data for tungsten indicating that the effective model captures the physics of the process quite well. For aluminum we predict secondary electron emission due to Auger de-excitation to be one to two orders of magnitude smaller than the one found for resonant charge-transfer and subsequent auto-detachment.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamic charge correlations near the Peierls transition

The quantum phase transition between a repulsive Luttinger liquid and an insulating Peierls state is studied in the framework of the one-dimensional spinless Holstein model. We focus on the adiabatic regime but include the full quantum dynamics of the phonons. Using continuous-time quantum Monte Carlo simulations, we track in particular the dynamic charge structure factor and the single-particle spectrum across the transition. With increasing electron-phonon coupling, the dynamic charge structure factor reveals the emergence of a charge gap, and a clear signature of phonon softening at the zone boundary. The single-particle spectral function evolves continuously across the transition. Hybridization of the charge and phonon modes of the Luttinger liquid description leads to two modes, one of which corresponds to the coherent polaron band. This band acquires a gap upon entering the Peierls phase, whereas the other mode constitutes the incoherent, high-energy spectrum with backfolded shadow bands. Coherent polaronic motion is a direct consequence of quantum lattice fluctuations. In the strong-coupling regime, the spectrum is described by the static, mean-field limit. Importantly, whereas finite electron density in general leads to screening of polaron effects, the latter reappear at half filling due to charge ordering and lattice dimerization.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Excitonic resonances in the 2D extended Falicov-Kimball model

Using the projector-based renormalization method we investigate the formation of the excitonic insulator phase in the two-dimensional (2D) spinless Falicov-Kimball model with dispersive $f$ electrons and address the existence of excitonic bound states at high temperatures on the semiconductor side of the semimetal-semiconductor transition. To this end we calculate the imaginary part of the dynamical electron-hole pair susceptibility and analyze the wave-vector and energy dependence of excitonic resonances emerging in the band gap. We thereby confirm the existence of the exciton insulator and its exciton environment within a generic two-band lattice model with local Coulomb attraction.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hybrid-parallel sparse matrix-vector multiplication with explicit communication overlap on current multicore-based systems

We evaluate optimized parallel sparse matrix-vector operations for several representative application areas on widespread multicore-based cluster configurations. First the single-socket baseline performance is analyzed and modeled with respect to basic architectural properties of standard multicore chips. Beyond the single node, the performance of parallel sparse matrix-vector operations is often limited by communication overhead. Starting from the observation that nonblocking MPI is not able to hide communication cost using standard MPI implementations, we demonstrate that explicit overlap of communication and computation can be achieved by using a dedicated communication thread, which may run on a virtual core. Moreover we identify performance benefits of hybrid MPI/OpenMP programming due to improved load balancing even without explicit communication overlap. We compare performance results for pure MPI, the widely used "vector-like" hybrid programming strategies, and explicit overlap on a modern multicore-based cluster and a Cray XE6 system.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Metal-to-insulator transition and electron-hole puddle formation in disordered graphene nanoribbons

The experimentally observed metal-to-insulator transition in hydrogenated graphene is numerically confirmed for actual sized graphene samples and realistic impurity concentrations. The eigenstates of our tight-binding model with substitutional disorder corroborate the formation of electron-hole-puddles with characteristic length scales comparable to the ones found in experiments. The puddles cause charge inhomogeneities and tend to suppress Anderson localization. Even though, monitoring the charge carrier quantum dynamics and performing a finite-size scaling of the local density of states distribution, we find strong evidence for the existence of localized states in graphene nanoribbons with short-range but also correlated long-range disorder.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Plasma walls beyond the perfect absorber approximation for electrons

Plasma walls accumulate electrons more efficiently than ions leading to wall potentials which are negative with respect to the plasma potential. Theoretically, walls are usually treated as perfect absorber for electrons and ions implying perfect sticking of the particles to the wall and infinitely long desorption times for particles stuck to the wall. For electrons we question the perfect absorber model and calculate, specifically for a planar dielectric wall, the electron sticking coefficient $s_e$ and the electron desorption time $τ_e$. For the uncharged wall we find $s_e\ll 1$ and $τ_e\approx 10^{-4}s$. Thus, in the early stage of the build-up of the wall potential, when the wall is essentially uncharged, the wall is not a perfect absorber for electrons. For the charged wall we find $τ_e^{-1}\approx 0$. Thus, $τ_e$ approaches the perfect absorber value. But $s_e$ is still only of the order of $10^{-1}$. Calculating $s_e$ as a function of the wall potential and combining this expression with the quasi-stationary balance equations for the electron and ion surface densities we find the selfconsistent wall potential, including surface effects, to be 30% of the perfect absorber value.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Distribution of the local density of states as a criterion for Anderson localization: Numerically exact results for various lattices in two and three dimensions

Numerical approaches to Anderson localization face the problem of having to treat large localization lengths while being restricted to finite system sizes. We show that by finite-size scaling of the probability distribution of the local density of states (LDOS) this long-standing problem can be overcome. To this end we reexamine the approach, propose numerical refinements, and apply it to study the dependence of the distribution of the LDOS on the dimensionality and coordination number of the lattice. Particular attention is given to the graphene lattice. We show that the system-size dependence of the LDOS distribution is indeed an unambiguous sign of Anderson localization, irrespective of the dimension and lattice structure. The numerically exact LDOS data obtained by us agree with a log-normal distribution over up to ten orders of magnitude and thereby fulfill a nontrivial symmetry relation previously derived for the non-linear $σ$-model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Momentum average approximation for models with boson-modulated hopping: Role of closed loops in the dynamical generation of a finite quasiparticle mass

We generalize the momentum average approximation to study the properties of single polarons in models with boson affected hopping, where the fermion-boson scattering depends explicitly on both the fermion's and the boson's momentum. As a specific example, we investigate the Edwards fermion-boson model in both one and two dimensions. In one dimension, this allows us to compare our results with exact diagonalization results, to validate the accuracy of our approximation. The generalization to two-dimensional lattices allows us to calculate the polaron's quasiparticle weight and dispersion throughout the Brillouin zone and to demonstrate the importance of Trugman loops in generating a finite effective mass even when the free fermion has an infinite mass.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Parallel sparse matrix-vector multiplication as a test case for hybrid MPI+OpenMP programming

We evaluate optimized parallel sparse matrix-vector operations for two representative application areas on widespread multicore-based cluster configurations. First the single-socket baseline performance is analyzed and modeled with respect to basic architectural properties of standard multicore chips. Going beyond the single node, parallel sparse matrix-vector operations often suffer from an unfavorable communication to computation ratio. Starting from the observation that nonblocking MPI is not able to hide communication cost using standard MPI implementations, we demonstrate that explicit overlap of communication and computation can be achieved by using a dedicated communication thread, which may run on a virtual core. We compare our approach to pure MPI and the widely used "vector-like" hybrid programming strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase separation of multicomponent excitonic Bose-Einstein condensates

For the observation of Bose-Einstein condensation, excitons in cuprous oxide are regarded as promising candidates due to their large binding energy and long lifetime. High particle densities may be achieved by entrapment in a stress induced potential. We consider a multi-component gas of interacting para- and orthoexcitons in cuprous oxide confined in a three-dimensional potential trap. Based on the Hartree-Fock-Bogoliubov theory, we calculate density profiles as well as decay luminescence spectra which exhibit signatures of the separation of the Bose-condensed phases.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Polarons and slow quantum phonons

We describe the formation and properties of Holstein polarons in the entire parameter regime. Our presentation focuses on the polaron mass and radius, which we obtain with an improved numerical technique. It is based on the combination of variational exact diagonalization with an improved construction of phonon states, providing results even for the strong coupling adiabatic regime. In particular we can describe the formation of large and heavy adiabatic polarons. A comparison of the polaron mass for the one and three dimensional situation explains how the different properties in the static oscillator limit determine the behavior in the adiabatic regime. The transport properties of large and small polarons are characterized by the f-sum rule and the optical conductivity. Our calculations are approximation-free and have negligible numerical error. This allows us to give a conclusive and impartial description of polaron formation. We finally discuss the implications of our results for situations beyond the Holstein model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spectral signatures of the BCS-BEC crossover in the excitonic insulator phase of the extended Falicov-Kimball model

We explore the spontaneous formation of an excitonic insulator state at the semimetal-semiconductor transition of mixed-valence materials in the framework of the spinless Falicov-Kimball model with direct $f$-$f$ electron hopping. Adapting the projector-based renormalization method, we obtain a set of renormalization differential equations for the extended Falicov-Kimball model parameters and finally derive analytical expressions for the order parameter, as well as for the renormalized $c$- and $f$-electron dispersions, momentum distributions, and wave-vector resolved single-particle spectral functions. Our numerical results proved the valence transition picture, related to the appearance of the excitonic insulator phase, in the case of overlapping $c$ and $f$ bands. Thereby the photoemission spectra show significant differences between the weak-to-intermediate and intermediate-to-strong Coulomb attraction regimes, indicating a BCS-BEC transition of the excitonic condensate.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stability of edge states and edge magnetism in graphene nanoribbons

We critically discuss the stability of edge states and edge magnetism in zigzag edge graphene nanoribbons (ZGNRs). We point out that magnetic edge states might not exist in real systems, and show that there are at least three very natural mechanisms - edge reconstruction, edge passivation, and edge closure - which dramatically reduce the effect of edge states in ZGNRs or even totally eliminate them. Even if systems with magnetic edge states could be made, the intrinsic magnetism would not be stable at room temperature. Charge doping and the presence of edge defects further destabilize the intrinsic magnetism of such systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Anderson disorder in graphene nanoribbons: A local distribution approach

Disorder effects strongly influence the transport properties of graphene based nanodevices even to the point of Anderson localization. Focusing on the local density of states and its distribution function, we analyze the localization properties of actual size graphene nanoribbons. In particular we determine the time evolution and localization length of the single particle wave function in dependence on the ribbon extension and edge geometry, as well as on the disorder type and strength.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Numerical approaches to time evolution of complex quantum systems

We examine several numerical techniques for the calculation of the dynamics of quantum systems. In particular, we single out an iterative method which is based on expanding the time evolution operator into a finite series of Chebyshev polynomials. The Chebyshev approach benefits from two advantages over the standard time-integration Crank-Nicholson scheme: speedup and efficiency. Potential competitors are semiclassical methods such as the Wigner-Moyal or quantum tomographic approaches. We outline the basic concepts of these techniques and benchmark their performance against the Chebyshev approach by monitoring the time evolution of a Gaussian wave packet in restricted one-dimensional (1D) geometries. Thereby the focus is on tunnelling processes and the motion in anharmonic potentials. Finally we apply the prominent Chebyshev technique to two highly non-trivial problems of current interest: (i) the injection of a particle in a disordered 2D graphene nanoribbon and (ii) the spatiotemporal evolution of polaron states in finite quantum systems. Here, depending on the disorder/electron-phonon coupling strength and the device dimensions, we observe transmission or localisation of the matter wave.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Performance limitations for sparse matrix-vector multiplications on current multicore environments

The increasing importance of multicore processors calls for a reevaluation of established numerical algorithms in view of their ability to profit from this new hardware concept. In order to optimize the existent algorithms, a detailed knowledge of the different performance-limiting factors is mandatory. In this contribution we investigate sparse matrix-vector multiplication, which is the dominant operation in many sparse eigenvalue solvers. Two conceptually different storage schemes and computational kernels have been conceived in the past to target cache-based and vector architectures, respectively. Starting from a series of microbenchmarks we apply the gained insight on optimized sparse MVM implementations, whose serial and OpenMP-parallel performance we review on state-of-the-art multicore systems.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Comparative numerical study of Anderson localization in disordered electron systems

Taking into account that a proper description of disordered systems should focus on distribution functions, the authors develop a powerful numerical scheme for the determination of the probability distribution of the local density of states (LDOS), which is based on a Chebyshev expansion with kernel polynomial refinement and allows the study of large finite clusters (up to $100^3$). For the three-dimensional Anderson model it is demonstrated that the distribution of the LDOS shows a significant change at the disorder induced delocalisation-localisation transition. Consequently, the so-called typical density of states, defined as the geometric mean of the LDOS, emerges as a natural order parameter. The calculation of the phase diagram of the Anderson model proves the efficiency and reliability of the proposed approach in comparison to other localisation criteria, which rely, e.g., on the decay of the wavefunction or the inverse participation number.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Anderson localization of polaron states

Using the vanishing of the typical polaron tunneling rate as an indicator of the breakdown of itinerancy, we study the localization of polaron states in a generic model for a disordered polaronic material. We find that extremely small disorder causes an Anderson localization of small polaron states. However, the ratio between the critical disorder strength needed to localize all states in the polaron band and the renormalized bandwidth is not necessarily smaller than for a bare electron.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Density-Matrix Algorithm for Phonon Hilbert Space Reduction in the Numerical Diagonalization of Quantum Many-Body Systems

Combining density-matrix and Lanczos algorithms we propose a new optimized phonon approach for finite-cluster diagonalizations of interacting electron-phonon systems. To illustrate the efficiency and reliability of our method, we investigate the problem of bipolaron band formation in the extended Holstein Hubbard model.

0 citations0 reviewsNo code linkedOpen access