Source author record

L. Covaci

L. Covaci 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

22works
4topics
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

22 published item(s)

preprint2020arXiv

Wave-packet scattering at a normal-superconductor interface in two-dimensional materials: a generalized theoretical approach

A wave-packet time evolution method, based on the split-operator technique, is developed to investigate the scattering of quasi-particles at a normal-superconductor interface of arbitrary profile and shape. As a practical application, we consider a system where low energy electrons can be described as Dirac particles, which is the case for most two-dimensional materials, such as graphene and transition metal dichalcogenides. However the method is easily adapted for other cases such as electrons in few layer black phosphorus, or any Schrödinger quasi-particles within the effective mass approximation in semiconductors. We employ the method to revisit Andreev reflection in graphene, where specular and retro reflection cases are observed for electrons scattered by a step-like superconducting region. The effect of opening a zero-gap channel across the superconducting region on the electron and hole scattering is also addressed, as an example of the versatility of the technique proposed here.

preprint2019arXiv

Double moiré with a twist: super-moiré in encapsulated graphene

A periodic spatial modulation, as created by a moiré pattern, has been extensively studied with the view to engineer and tune the properties of graphene. Graphene encapsulated by hexagonal boron nitride (hBN) when slightly misaligned with the top and bottom hBN layers experiences two interfering moiré patterns, resulting in a so-called super-moiré (SM). This leads to a lattice and electronic spectrum reconstruction. A geometrical construction of the non-relaxed SM patterns allows us to indicate qualitatively the induced changes in the electronic properties and to locate the SM features in the density of states and in the conductivity. To emphasize the effect of lattice relaxation, we report band gaps at all Dirac-like points in the hole doped part of the reconstructed spectrum, which are expected to be enhanced when including interaction effects. Our result is able to distinguish effects due to lattice relaxation and due to the interfering SM and provides a clear picture on the origin of recently experimentally observed effects in such trilayer heterostuctures.

preprint2015arXiv

Vortex detection and quantum transport in mesoscopic graphene Josephson junction arrays

We investigate mesoscopic Josephson junction arrays created by patterning superconducting disks on monolayer graphene, concentrating on the high-$T/T_c$ regime of these devices and the phenomena which contribute to the superconducting glass state in diffusive arrays. We observe features in the magnetoconductance at rational fractions of flux quanta per array unit cell, which we attribute to the formation of flux-quantized vortices. The applied fields at which the features occur are well described by Ginzburg-Landau simulations that take into account the number of unit cells in the array. We find that the mean conductance and universal conductance fluctuations are both enhanced below the critical temperature and field of the superconductor, with greater enhancement away from the graphene Dirac point.

preprint2014arXiv

Disordered graphene Josephson junctions

A tight-binding approach based on the Chebyshev-Bogoliubov-de Gennes method is used to describe disordered single-layer graphene Josephson junctions. Scattering by vacancies, ripples or charged impurities is included. We compute the Josephson current and investigate the nature of multiple Andreev reflections, which induce bound states appearing as peaks in the density of states for energies below the superconducting gap. In the presence of single atom vacancies, we observe a strong suppression of the supercurrent that is a consequence of strong inter-valley scattering. Although lattice deformations should not induce inter-valley scattering, we find that the supercurrent is still suppressed, which is due to the presence of pseudo-magnetic barriers. For charged impurities, we consider two cases depending on whether the average doping is zero, i.e. existence of electron-hole puddles, or finite. In both cases, short range impurities strongly affect the supercurrent, similar to the vacancies scenario.

preprint2014arXiv

Electronic structure of graphene hexagonal flake subjected to triaxial stress

The electronic properties of a triaxially strained hexagonal graphene flake with either armchair or zig-zag edges are investigated using molecular dynamics simulations and tight-binding calculations. We found that: i) the pseudo-magnetic field in the strained graphene flakes is not uniform neither in the center nor at the edge of zig-zag terminated flakes, ii) the pseudo-magnetic field is almost zero in the center of armchair terminated flakes but increases dramatically near the edges, iii) the pseudo-magnetic field increases linearly with strain, for strains lower than 15$%$ while growing non-linearly beyond this threshold, iv) the local density of states in the center of the zig-zag hexagon exhibits pseudo-Landau levels with broken sub-lattice symmetry in the zero'th pseudo-Landau level, and in addition there is a shift in the Dirac cone due to strain induced scalar potentials. This study provides a realistic model of the electronic properties of inhomogeneously strained graphene where the relaxation of the atomic positions is correctly included together with strain induced modifications of the hopping terms up to next-nearest neighbors.

preprint2014arXiv

Non-linear response to electric field in extended Hubbard models

The electric-field response of a one-dimensional ring of interacting fermions, where the interactions are described by the extended Hubbard model, is investigated. By using an accurate real-time propagation scheme based on the Chebyshev expansion of the evolution operator, we uncover various non-linear regimes for a range of interaction parameters that allows modeling of metallic and insulating (either charge density wave or spin density wave insulators) rings. The metallic regime appears at the phase boundary between the two insulating phases and provides the opportunity to describe either weakly or strongly correlated metals. We find that the {\it fidelity susceptibility} of the ground state as a function of magnetic flux piercing the ring provides a very good measure of the short-time response. Even completely different interacting regimes behave in a similar manner at short time-scales as long as the fidelity susceptibility is the same. Depending on the strength of the electric field we find various types of responses: persistent currents in the insulating regime, dissipative regime or damped Bloch-like oscillations with varying frequencies or even irregular in nature. Furthermore, we also consider the dimerization of the ring and describe the response of a correlated band insulator. In this case the distribution of the energy levels is more clustered and the Bloch-like oscillations become even more irregular.

preprint2014arXiv

The Tomasch effect in nanoscale superconductors

The Tomasch effect (TE) is due to quasiparticle interference (QPI) as induced by a nonuniform superconducting order parameter, which results in oscillations in the density of states (DOS) at energies above the superconducting gap. Quantum confinement in nanoscale superconductors leads to an inhomogenerous distribution of the Cooper-pair condensate, which we found triggers the manifestation of a new TE. We investigate the electronic structure of the nanoscale superconductors by solving the Bogoliubov-de Gennes (BdG) equations self-consistently and obtain the TE determined by two types of processes, involving two or three particle QPIs. Both types of QPIs result in additional BCS-like Bogoliubov-quasiparticles and BCS-like energy gaps leading to oscillations in the DOS and modulated wave patterns in the local density of states. These effects are strongly related to the symmetries of the system. A reduced $4 \times 4$ inter-subbands BdG Hamiltonian is established in order to describe analytically the TE of two-QPIs. Our study is relevant to nanoscale superconductors, either nanowires or thin films, Bose-Einsten condensates and confined systems such as two-dimensional electron gas interface superconductivity.

preprint2013arXiv

Field effect on surface states in a doped Mott-Insulator thin film

Surface effects of a doped thin film made of a strongly correlated material are investigated both in the absence and presence of a perpendicular electric field. We use an inhomogeneous Gutzwiller approximation for a single band Hubbard model in order to describe correlation effects. For low doping, the bulk value of the quasiparticle weight is recovered exponentially deep into the slab, but with increasing doping, additional Friedel oscillations appear near the surface. We show that the inverse correlation length has a power-law dependence on the doping level. In the presence of an electrical field, considerable changes in the quasiparticle weight can be realized throughout the system. We observe a large difference (as large as five orders of magnitude) in the quasiparticle weight near the opposite sides of the slab. This effect can be significant in switching devices that use the surface states for transport.

preprint2013arXiv

Superconducting current and proximity effect in ABA and ABC multilayer graphene Josephson junctions

Using a numerical tight-binding approach based on the Chebyshev-Bogoliubov-de Gennes method we describe Josephson junctions made of multilayer graphene contacted by top superconducting gates. Both Bernal (ABA) and rhombohedral (ABC) stacking are considered and we find that the type of stacking has a strong effect on the proximity effect and the supercurrent flow. For both cases the pair amplitude shows a polarization between dimer and non-dimer atoms, being more pronounced for rhombohedral stacking. Even though the proximity effect in non-dimer sites is enhanced when compared to single layer graphene, we find that the supercurrent is suppressed. The spatial distribution of the supercurrent shows that for Bernal stacking the current flows only in the top-most layers while for rhombohedral stacking the current flows throughout the whole structure.

preprint2013arXiv

Tight-binding description of intrinsic superconducting correlations in multilayer graphene

Using highly efficient GPU-based simulations of the tight-binding Bogoliubov-de Gennes equations we solve self-consistently for the pair correlation in rhombohedral (ABC) and Bernal (ABA) multilayer graphene by considering a finite intrinsic s-wave pairing potential. We find that the two different stacking configurations have opposite bulk/surface behavior for the order parameter. Surface superconductivity is robust for ABC stacked multilayer graphene even at very low pairing potentials for which the bulk order parameter vanishes, in agreement with a recent analytical approach. In contrast, for Bernal stacked multilayer graphene, we find that the order parameter is always suppressed at the surface and that there exists a critical value for the pairing potential below which no superconducting order is achieved. We considered different doping scenarios and find that homogeneous doping strongly suppresses surface superconductivity while non-homogeneous field-induced doping has a much weaker effect on the superconducting order parameter. For multilayer structures with hybrid stacking (ABC and ABA) we find that when the thickness of each region is small (few layers), high-temperature surface superconductivity survives throughout the bulk due to the proximity effect between ABC/ABA interfaces where the order parameter is enhanced.

preprint2013arXiv

Vortex states in nanoscale superconducting squares: the influence of quantum confinement

Bogoliubov-de Gennes theory is used to investigate the effect of the size of a superconducting square on the vortex states in the quantum confinement regime. When the superconducting coherence length is comparable to the Fermi wavelength, the shape resonances of the superconducting order parameter have strong influence on the vortex configuration. Several unconventional vortex states, including asymmetric ones, giant multi-vortex combinations, and states comprising giant antivortex, were found as ground states and their stability was found to be very sensitive on the value of $k_Fξ_0$, the size of the sample $W$, and the magnetic flux $Φ$. By increasing the temperature and/or enlarging the size of the sample, quantum confinement is suppressed and the conventional mesoscopic vortex states as predicted by the Ginzburg-Laudau (GL) theory are recovered. However, contrary to the GL results we found that the states containing symmetry-induced vortex-antivortex pairs are stable over the whole temperature range. It turns out that the inhomogeneous order parameter induced by quantum confinement favors vortex-antivortex molecules, as well as giant vortices with a rich structure in the vortex core - unattainable in the GL domain.

preprint2012arXiv

Impact of Dresselhaus vs. Rashba spin-orbit coupling on the Holstein polaron

We utilize an exact variational numerical procedure to calculate the ground state properties of a polaron in the presence of Rashba and linear Dresselhaus spin-orbit coupling. We find that when the linear Dresselhaus spin-orbit coupling approaches the Rashba spin-orbit coupling, the Van-Hove singularity in the density of states will be shifted away from the bottom of the band and finally disappear when the two spin-orbit couplings are tuned to be equal. The effective mass will be suppressed; the trend will become more significant for low phonon frequency. The presence of two dominant spin-orbit couplings will make it possible to tune the effective mass with more varied observables.

preprint2012arXiv

Nano-engineered non-uniform strain in graphene

Recent experiments showed that non-uniform strain can be produced by depositing graphene over pillars. We employed atomistic calculations to study the non-uniform strain and the induced pseudo-magnetic field up to 5000 Tesla in graphene on top of nano-pillars. By decreasing the distance between the nano-pillars a complex distribution for the pseudo-magnetic field can be generated. Furthermore, we performed tight-binding calculations of the local density of states (LDOS) by using the relaxed graphene configuration obtained from the atomistic calculations. We find that the quasiparticle LDOS are strongly modified near the pillars, both at low energies showing sub-lattice polarization, and at high energies showing shifts of the van Hove singularity. Our study shows that changing the specific pattern of the nano-pillars allows us to create a desired shape of the pseudo-magnetic field profile while the LDOS maps provide an input for experimental verifications by scanning tunneling microscopy.

preprint2012arXiv

Partially unzipped carbon nanotubes as magnetic field sensors

The conductance, $G(E)$, through graphene nanoribbons (GNR) connected to a partially unzipped carbon nanotube (CNT) is studied in the presence of an external magnetic field applied parallel to the long axis of the tube by means of non-equilibrium Green's function technique. We consider (z)igzag and (a)rmchair CNTs that are partially unzipped to form aGNR/zCNT/aGNR or zGNR/aCNT/zGNR junctions. We find that the inclusion of a longitudinal magnetic field affects the electronic states only in the CNT region, leading to the suppression of the conductance at low energies. Unlike previous studies, for the zGNR/aCNT/zGNR junction in zero field, we find a sharp dip in the conductance as the energy approaches the Dirac point and we attribute this non-trivial behavior to the peculiar band dispersion of the constituent subsystems. We demonstrate that both types of junctions can be used as magnetic field sensors.

preprint2012arXiv

Tight-binding study of bilayer graphene Josephson junctions

Using highly efficient simulations of the tight-binding Bogoliubov-de Gennes model we solved self-consistently for the pair correlation and the Josephson current in a Superconducting-Bilayer graphene-Superconducting Josephson junction. Different doping levels for the non-superconducting link are considered in the short and long junction regime. Self-consistent results for the pair correlation and superconducting current resemble those reported previously for single layer graphene except in the Dirac point where remarkable differences in the proximity effect are found as well as a suppression of the superconducting current in long junction regime. Inversion symmetry is broken by considering a potential difference between the layers and we found that the supercurrent can be switched if junction length is larger than the Fermi length.

preprint2012arXiv

Unconventional vortex states in nanoscale superconductors due to shape-induced resonances in the inhomogeneous Cooper-pair condensate

Vortex matter in mesoscopic superconductors is known to be strongly affected by the geometry of the sample. Here we show that in nanoscale superconductors with coherence length comparable to the Fermi wavelength the shape resonances of the order parameter results in an additional contribution to the quantum topological confinement - leading to unconventional vortex configurations. Our Bogoliubov-de Gennes calculations in a square geometry reveal a plethora of asymmetric, giant multi-vortex, and vortex-antivortex structures, stable over a wide range of parameters and which are very different from those predicted by the Ginzburg-Landau theory. These unconventional states are relevant for high-Tc nanograins, confined Bose-Einstein condensates, and graphene flakes with proximity-induced superconductivity.

preprint2012arXiv

Wavepacket scattering on graphene edges in the presence of a (pseudo) magnetic field

The scattering of a Gaussian wavepacket in armchair and zigzag graphene edges is theoretically investigated by numerically solving the time dependent Schrödinger equation for the tight-binding model Hamiltonian. Our theory allows to investigate scattering in reciprocal space, and depending on the type of graphene edge we observe scattering within the same valley, or between different valleys. In the presence of an external magnetic field, the well know skipping orbits are observed. However, our results demonstrate that in the case of a pseudo-magnetic field, induced by non-uniform strain, the scattering by an armchair edge results in a non-propagating edge state.

preprint2011arXiv

Electrical field induced shift of the Mott Metal-Insulator transition in thin films

The ground state properties of a paramagnetic Mott insulator are investigated in the presence of an external electrical field using the inhomogeneous Gutzwiller approximation for a single band Hubbard model in a slab geometry. The metal insulator transition is shifted towards higher Hubbard repulsions by applying an electric field perpendicular to the slab. The spatial distribution of site dependent quasiparticle weight shows that the quasiparticle weight is maximum in few layers beneath the surface. Moreover only at higher Hubbard repulsion, larger than the bulk critical U, the electric field will be totally screened only for centeral cites. Our results show that by presence of an electric field perpendicular to a thin film made of a strongly correlated material, states near the surface will remain metallic while the bulk becomes insulating after some critical U. In contrast, in the absence of the electric field the surface becomes insulating before the bulk.

preprint2011arXiv

Impact of spin-orbit coupling on the Holstein polaron

We utilize an exact variational numerical procedure to calculate the ground state properties of a polaron in the presence of a Rashba-like spin orbit interaction. Our results corroborate with previous work performed with the Momentum Average approximation and with weak coupling perturbation theory. We find that spin orbit coupling increases the effective mass in the regime with weak electron phonon coupling, and decreases the effective mass in the intermediate and strong electron phonon coupling regime. Analytical strong coupling perturbation theory results confirm our numerical results in the small polaron regime. A large amount of spin orbit coupling can lead to a significant lowering of the polaron effective mass.

preprint2011arXiv

Rectification by imprinted phase in a Josephson junction

A Josephson phase shift can be induced in a Josephson junction by a strategically nearby pinned Abrikosov vortex (AV). For an asymmetric distribution of imprinted phase along the junction (controlled by the position of the AV) such a simple system is capable of rectification of $ac$ current in a broad and tunable frequency range. The resulting rectified voltage is a consequence of the directed motion of a Josephson antivortex which forms a pair with the AV when at local equilibrium. The proposed realization of the ratchet potential by imprinted phase is more efficient than the asymmetric geometry of the junction itself, it is easily realizable experimentally, and provides rectification even in the absence of applied magnetic field.

preprint2011arXiv

Wave packet dynamics and valley filter in strained graphene

The time evolution of a wavepacket in strained graphene is studied within the tight-binding model and continuum model. The effect of an external magnetic field, as well as a strain-induced pseudo-magnetic field, on the wave packet trajectories and zitterbewegung are analyzed. Combining the effects of strain with those of an external magnetic field produces an effective magnetic field which is large in one of the Dirac cones, but can be practically zero in the other. We construct an efficient valley filter, where for a propagating incoming wave packet consisting of momenta around the K and K' Dirac points, the outgoing wave packet exhibits momenta in only one of these Dirac points, while the components of the packet that belong to the other Dirac point are reflected due to the Lorentz force. We also found that the zitterbewegung is permanent in time in the presence of either external or strain-induced magnetic fields, but when both the external and strain-induced magnetic fields are present, the zitterbewegung is transient in one of the Dirac cones, whereas in the other cone the wave packet exhibits permanent spatial oscillations.

preprint2010arXiv

Chebyshev-BdG: an efficient numerical approach to inhomogeneous superconductivity

We propose a highly efficient numerical method to describe inhomogeneous superconductivity by using the kernel polynomial method in order to calculate the Green's functions of a superconductor. Broken translational invariance of any type (impurities, surfaces or magnetic fields) can be easily incorporated. We show that limitations due to system size can be easily circumvented and therefore this method opens the way for the study of scenarios and/or geometries that were unaccessible before. The proposed method is highly efficient and amenable to large scale parallel computation. Although we only use it in the context of superconductivity, it is applicable to other inhomogeneous mean-field theories.