Source author record

Topi Siro

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

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

10 published item(s)

preprint2016arXiv

Exact diagonalization of quantum lattice models on coprocessors

We implement the Lanczos algorithm on an Intel Xeon Phi coprocessor and compare its performance to a multi-core Intel Xeon CPU and an NVIDIA graphics processor. The Xeon and the Xeon Phi are parallelized with OpenMP and the graphics processor is programmed with CUDA. The performance is evaluated by measuring the execution time of a single step in the Lanczos algorithm. We study two quantum lattice models with different particle numbers, and conclude that for small systems, the multi-core CPU is the fastest platform, while for large systems, the graphics processor is the clear winner, reaching speedups of up to 7.6 compared to the CPU. The Xeon Phi outperforms the CPU with sufficiently large particle number, reaching a speedup of 2.5.

preprint2016arXiv

Topological phase transitions in the repulsively interacting Haldane-Hubbard model

Using dynamical mean-field theory and exact diagonalization we study the phase diagram of the repulsive Haldane-Hubbard model, varying the interaction strength and the sublattice potential difference. In addition to the quantum Hall phase with Chern number $C=2$ and the band insulator with $C=0$ present already in the noninteracting model, the system also exhibits a $C=0$ Mott insulating phase, and a $C=1$ quantum Hall phase. We explain the latter phase by a spontaneous symmetry breaking where one of the spin-components is in the Hall state and the other in the band insulating state.

preprint2014arXiv

Impurities and Landau level mixing in a fractional quantum Hall state in a flatband lattice model

We study the toplogical checkerboard lattice model around the $ν=\frac{1}{3}$ fractional quantum Hall phase using numerical exact diagonalization without Landau level projections. We add local perturbations, modified hoppings and on-site potentials, and observe phase transitions from the fractional quantum Hall phase to metallic and insulating phases when the strength and number of impurities is increased. In addition to evaluating the energy spectrum, we identify the phase diagrams by computing the topological Chern number of the many-body ground state manifold, and we show how the ground states lose their correlations due to the impurities by evaluating the spectrum of the one-body reduced density matrix. Our results show that the phase transition from the fractional quantum Hall phase to the metallic phase occurs for both impurity hoppings and potentials. Strong impurity hoppings cause a further transition into the insulating state, regardless of the sign of the hopping, when their density is high enough. In contrast, the same happens only for attractive potentials. Furthermore, the mixing to the higher band in a two-band model, generally denoted as Landau level mixing, is measured concluding that the lowest Landau level projection works well even with remarkably strong interactions and in the presence of impurities.

preprint2013arXiv

Accelerated molecular dynamics force evaluation on graphics processing units for thermal conductivity calculations

In this paper, we develop a highly efficient molecular dynamics code fully implemented on graphics processing units for thermal conductivity calculations using the Green-Kubo formula. We compare two different schemes for force evaluation, a previously used thread-scheme where a single thread is used for one particle and each thread calculates the total force for the corresponding particle, and a new block-scheme where a whole block is used for one particle and each thread in the block calculates one or several pair forces between the particle associated with the given block and its neighbor particle(s) associated with the given thread. For both schemes, two different classical potentials, namely, the Lennard-Jones potential and the rigid-ion potential are implemented. While the thread-scheme performs a little better for relatively large systems, the block-scheme performs much better for relatively small systems. The relative performance of the block-scheme over the thread-scheme also increases with the increasing of the cutoff radius. We validate the implementation by calculating lattice thermal conductivities of solid argon and lead telluride.

preprint2013arXiv

Computational Physics on Graphics Processing Units

The use of graphics processing units for scientific computations is an emerging strategy that can significantly speed up various different algorithms. In this review, we discuss advances made in the field of computational physics, focusing on classical molecular dynamics, and on quantum simulations for electronic structure calculations using the density functional theory, wave function techniques, and quantum field theory.

preprint2013arXiv

Efficient linear-scaling quantum transport calculations on graphics processing units and applications on electron transport in graphene

We implement, optimize, and validate the linear-scaling Kubo-Greenwood quantum transport simulation on graphics processing units by examining resonant scattering in graphene. We consider two practical representations of the Kubo-Greenwood formula: a Green-Kubo formula based on the velocity auto-correlation and an Einstein formula based on the mean square displacement. The code is fully implemented on graphics processing units with a speedup factor of up to 16 (using double-precision) relative to our CPU implementation. We compare the kernel polynomial method and the Fourier transform method for the approximation of the Dirac delta function and conclude that the former is more efficient. In the ballistic regime, the Einstein formula can produce the correct quantized conductance of one-dimensional graphene nanoribbons except for an overshoot near the band edges. In the diffusive regime, the Green-Kubo and the Einstein formalisms are demonstrated to be equivalent. A comparison of the length-dependence of the conductance in the localization regime obtained by the Einstein formula with that obtained by the non-equilibrium Green's function method reveals the challenges in defining the length in the Kubo-Greenwood formalism at the strongly localized regime.

preprint2013arXiv

Enhancement models of momentum densities of annihilating electron-positron pairs: the many-body picture of natural geminals

The correlated motion of a positron surrounded by electrons is a fundamental many-body problem. We approach this by modeling the momentum density of annihilating electron-positron pairs using the framework of reduced density matrices, natural orbitals and natural geminals (electron-positron pair wave functions) of the quantum theory of many-particle systems. We find that an expression based on the natural geminals provides an exact, unique and compact expression for the momentum density. The natural geminals can be used to define and determine enhancement factors for enhancement models going beyond the independent-particle model for a better understanding of results of positron annihilation experiments.

preprint2013arXiv

Lagrange mesh and exact diagonalization for numerical study of semiconductor quantum dot systems with application in singlet-triplet qubits

We present a highly flexible computational scheme for studying correlated electrons confined by an arbitrary external potential in two-dimensional semiconductor quantum dots. The method starts by a Lagrange mesh calculation for the single-particle states, followed by the calculation of the Coulomb interaction matrix elements between these, and combining both in the exact diagonalization of the many-body Hamiltonian. We apply the method in simulation of double quantum dot singlet-triplet qubits. We simulate the full quantum control and dynamics of one singlet-triplet qubit. We also use our method to provide an exact diagonalization based first-principles model for studying two singlet-triplet qubits and their capacitative coupling via the long-distance Coulomb interaction.

preprint2013arXiv

Non-adiabatic charge state transitions in singlet-triplet qubits

In double quantum dot singlet-triplet qubits, the exchange interaction is used in both quantum gate operation and the measurement of the state of the qubit. The exchange can be controlled electronically by applying gate voltage pulses. We simulate the exchange induced charge state transitions in one and two singlet-triplet qubit systems using the exact diagonalization method. We find that fast detuning pulses may result in leakage between different singlet charge states. The leakage could cause measurement errors and hinder quantum gate operation for example in the case of the two-qubit Coulomb gate.

preprint2012arXiv

Exact diagonalization of the Hubbard model on graphics processing units

We solve the Hubbard model with the exact diagonalization method on a graphics processing unit (GPU). We benchmark our GPU program against a sequential CPU code by using the Lanczos algorithm to solve the ground state energy in two cases: a one-dimensional ring and a two-dimensional square lattice. In the one-dimensional case, we obtain speedups of over 100 and 60 in single and double precision arithmetic, respectively. In the two-dimensional case, the corresponding speedups are over 110 and 70.