Source author record

F. M. Peeters

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

ResearcherUnclaimed source record

cond-mat.mes-hall cond-mat.supr-con cond-mat.mtrl-sci cond-mat.soft cond-mat.str-el cond-mat.stat-mech quant-ph cond-mat.quant-gas Biological Physics cond-mat.other hep-ph nlin.PS physics.class-ph physics.comp-ph physics.optics

Catalog footprint

What is connected

230works

15topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Alternating superconducting and charge density wave monolayers within bulk 6R-TaS2

Van der Waals (vdW) heterostructures continue to attract intense interest as a route of designing materials with novel properties that cannot be found in naturally occurring materials. Unfortunately, this approach is currently limited to only a few layers that can be stacked on top of each other. Here we report a bulk material consisting of superconducting monolayers interlayered with monolayers displaying charge density waves (CDW). This bulk vdW heterostructure is created by phase transition of 1T-TaS2 to 6R at 800 °C in an inert atmosphere. Electron microscopy analysis directly shows the presence of alternating 1T and 1H monolayers within the resulting bulk 6R phase. Its superconducting transition (Tc) is found at 2.6 K, exceeding the Tc of the bulk 2H phase of TaS2. The superconducting temperature can be further increased to 3.6 K by exfoliating 6R-TaS2 and then restacking its layers. Using first-principles calculations, we argue that the coexistence of superconductivity and CDW within 6R-TaS2 stems from amalgamation of the properties of adjacent 1H and 1T monolayers, where the former dominates the superconducting state and the latter the CDW behavior.

preprint2022arXiv

Anisotropic and tunable optical conductivity of a two-dimensional semi-Dirac system in the presence of elliptically polarized radiation

We investigate the effect of ellipticity ratio of the polarized radiation field on optoelectronic properties of a two-dimensional (2D) semi-Dirac (SD) system. The optical conductivity is calculated within the energy balance equation approach derived from the semiclassical Boltzmann equation. We find that there exists the anisotropic optical absorption induced via both the intra- and interband electronic transition channels in the perpendicular $xx$ and $yy$ directions. Furthermore, we examine the effects of the ellipticity ratio, the temperature, the carrier density, and the band-gap parameter on the optical conductivity of the 2D SD system placed in transverse and vertical directions, respectively. It is shown that the ellipticity ratio, temperature, carrier density, and band-gap parameter can play the important roles in tuning the strength, peak position, and shape of the optical conductivity spectrum. The results obtained from this study indicate that the 2D SD system can be a promising anisotropic and tunable optical and optoelectronic material for applications in innovative 2D optical and optoelectronic devices, which are active in the infrared and terahertz bandwidths.

preprint2022arXiv

Electronic Mach-Zehnder interference in a bipolar hybrid monolayer-bilayer graphene junction

Graphene matter in a strong magnetic field, realizing one-dimensional quantum Hall channels, provides a unique platform for studying electron interference. Here, using the Landauer-Büttiker formalism along with the tight-binding model, we investigate the quantum Hall (QH) effects in unipolar and bipolar monolayer-bilayer graphene (MLG-BLG) junctions. We find that a Hall bar made of an armchair MLG-BLG junction in the bipolar regime results in valley-polarized edge-channel interferences and can operate a fully tunable Mach-Zehnder (MZ) interferometer device. Investigation of the bar-width and magnetic-field dependence of the conductance oscillations shows that the MZ interference in such structures can be drastically affected by the type of (zigzag) edge termination of the second layer in the BLG region [composed of vertical dimer or non-dimer atoms]. Our findings reveal that both interfaces exhibit a double set of Aharonov-Bohm interferences, with the one between two oppositely valley-polarized edge channels dominating and causing a large-amplitude conductance oscillation ranging from 0 to $ 2 e^2 / h$. We explain and analyze our findings by analytically solving the Dirac-Weyl equation for a gated semi-infinite MLG-BLG junction.

preprint2022arXiv

Isolated and hybrid bilayer graphene rings

Using the continuum model, we investigate the electronic properties of two types of bilayer graphene (BLG) quantum ring (QR) geometries: (i) an isolated BLG QR and (ii) a monolayer graphene (MLG) with a QR put on top of an infinite graphene sheet (hybrid BLG QR). Solving the Dirac-Weyl equation in the presence of a perpendicular magnetic field and applying the infinite-mass boundary condition at the ring boundaries, we obtain analytical results for the energy levels and corresponding wave spinors for both structures. In the case of isolated BLG QR, we observe a sizeable and magnetically tunable band gap which agrees with the tight-binding transport simulations. Our analytical results also show the intervalley symmetry $ E^K_e (m) = -E^{K'}_h(m) $ between the electron (e) and hole (h) states ($ m $ being the angular momentum quantum number) for the energy spectrum of the isolated BLG QR. The presence of interface boundary in a hybrid BLG QR modifies drastically the energy levels as compared to that of an isolated BLG QR. Its energy levels are tunable from MLG dot, to isolated BLG QR, and to MLG Landau energy levels as magnetic field is varied. Our predictions can be verified experimentally using different techniques such as by magnetotransport measurements.

preprint2022arXiv

Strong anisotropic optical properties of 8-Pmmn borophene: a many-body perturbation study

Using first-principle many-body perturbation theory, we investigate the optical properties of 8- borophene at two levels of approximations; the GW method considering only the electron-electron interaction and the GW in combination with the Bethe-Salpeter equation including electron-hole coupling. The band structure exhibits anisotropic Dirac cones with semimetallic character. The optical absorption spectra are obtained for different light polarizations and we predict strong optical absorbance anisotropy. The absorption peaks undergo a global redshift when the electron-hole interaction is taken into account due to the formation of bound excitons which have an anisotropic excitonic wave function.

preprint2021arXiv

Reversible ratchet effects in a narrow superconducting ring

We study the ratchet effect in a narrow pinning-free superconductive ring based on time-dependent Ginzburg-Landau (TDGL) equations. Voltage responses to external dc an ac currents at various magnetic fields are studied. Due to asymmetric barriers for flux penetration and flux exit in the ring-shaped superconductor, the critical current above which the flux-flow state is reached, as well as the critical current for the transition to the normal state, are different for the two directions of applied current. These effects cooperatively cause ratchet signal reversal at high magnetic fields, which has not been reported to date in a pinning-free system. The ratchet signal found here is larger than those induced by asymmetric pinning potentials. Our results also demonstrate the feasibility of using mesoscopic superconductors to employ superconducting diode effect in versatile superconducting devices.

preprint2020arXiv

Hematite at its thinnest limit

Motivated by the recent synthesis of two-dimensional $α$-Fe$_2$O$_3$ [Balan $et$ $al.$ Nat. Nanotech. 13, 602 (2018)], we analyze the structural, vibrational, electronic and magnetic properties of single- and few-layer $α$-Fe$_2$O$_3$ compared to bulk, by $ab-initio$ and Monte-Carlo simulations. We reveal how monolayer $α$-Fe$_2$O$_3$ (hematene) can be distinguished from the few-layer structures, and how they all differ from bulk through observable Raman spectra. The optical spectra exhibit gradual shift of the prominent peak to higher energy, as well as additional features at lower energy when $α$-Fe$_2$O$_3$ is thinned down to a monolayer. Both optical and electronic properties have strong spin asymmetry, meaning that lower-energy optical and electronic activities are allowed for the single-spin state. Finally, our considerations of magnetic properties reveal that 2D hematite has anti-ferromagnetic ground state for all thicknesses, but the critical temperature for Morin transition increases with decreasing sample thickness. On all accounts, the link to available experimental data is made, and further measurements are prompted.

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.

preprint2019arXiv

Three-dimensional electron-hole superfluidity in a superlattice close to room temperature

Although there is strong theoretical and experimental evidence for electron-hole superfluidity in separated sheets of electrons and holes at low $T$, extending superfluidity to high $T$ is limited by strong 2D fluctuations and Kosterlitz-Thouless effects. We show this limitation can be overcome using a superlattice of alternating electron- and hole-doped semiconductor monolayers. The superfluid transition in a 3D superlattice is not topological, and for strong electron-hole pair coupling, the transition temperature $T_c$ can be at room temperature. As a quantitative illustration, we show $T_c$ can reach $270$ K for a superfluid in a realistic superlattice of transition metal dichalcogenide monolayers.

preprint2019arXiv

Tight-Binding Studio: A Technical Software Package to Find the Parameters of Tight-Binding Hamiltonian

We present the Tight-Binding Studio (TBStudio) software package for calculating tight-binding Hamiltonian from a set of Bloch energy bands obtained from first principle theories such as density functional theory, Hartree-Fock calculations or Semi-empirical band structure theory. This will be helpful for scientists who are interested in studying the electronic properties of structures using Green's function theory in tight-binding approximation. TBStudio is a cross-platform application written in C++ with a graphical user interface design that is user-friendly and easy to work with. This software is powered by Linear Algebra Package C interface library for solving the eigenvalue problems and the standard high-performance OpenGL graphic library for real-time plotting. TBStudio and its examples together with the tutorials are available for download from tight-binding.com.

preprint2016arXiv

All-strain based valley filter in graphene nanoribbons using snake states

A pseudo-magnetic field kink can be realized along a graphene nanoribbon using strain engineering. Electron transport along this kink is governed by snake states that are characterized by a single propagation direction. Those pseudo-magnetic fields point towards opposite directions in the K and K' valleys, leading to valley polarized snake states. In a graphene nanoribbon with armchair edges this effect results in a valley filter that is based only on strain engineering. We discuss how to maximize this valley filtering by adjusting the parameters that define the stress distribution along the graphene ribbon.

preprint2016arXiv

Anomalous Dynamical Behavior of Freestanding Graphene Membranes

We report subnanometer, high-bandwidth measurements of the out-of-plane (vertical) motion of atoms in freestanding graphene using scanning tunneling microscopy. By tracking the vertical position over a long time period, a 1000-fold increase in the ability to measure space-time dynamics of atomically thin membranes is achieved over the current state-of-the-art imaging technologies. We observe that the vertical motion of a graphene membrane exhibits rare long-scale excursions characterized by both anomalous mean-squared displacements and Cauchy-Lorentz power law jump distributions.

preprint2016arXiv

Bilayer SnS$_{2}$: Easy-tunable Stacking Sequence by Charging and Loading Pressure

Employing density functional theory-based methods, we investigate monolayer and bilayer structures of hexagonal SnS$_{2}$, which is recently synthesized monolayer metal dichalcogenide. Comparison of 1H and 1T phases of monolayer SnS$_{2}$ confirms the ground state to be the 1T phase. In its bilayer structure we examine different stacking configurations of the two layers. It is found that the interlayer coupling in bilayer SnS$_{2}$ is weaker than that of typical transition-metal dichalcogenides (TMDs) so that alternative stacking orders have similar structural parameters and they are separated with low energy barriers. Possible signature of the stacking order in SnS$_{2}$ bilayer has been sought in the calculated absorbance and reflectivity spectra. We also study the effects of the external electric field, charging, and loading pressure on the characteristic properties of bilayer SnS$_{2}$. It is found that (i) the electric field increases the coupling between the layers at its prefered stacking order, so the barrier height increases, (ii) the bang gap value can be tuned by the external E-field and under sufficient E-field, the bilayer SnS$_{2}$ can become semi-metal, (iii) the most favorable stacking order can be switched by charging and (iv) a loading pressure exceeding 3 GPa changes the stacking order. E-field tunable bandgap and easy-tunable stacking sequence of SnS$_{2}$ layers make this 2D crystal structure a good candidate for field effect transistor and nanoscale lubricant applications.

preprint2016arXiv

Commensurability Effects in Viscosity of Nanoconfined Water

The rate of water flow through hydrophobic nanocapillaries is greatly enhanced as compared to that expected from macroscopic hydrodynamics. This phenomenon is usually described in terms of a relatively large slip length, which is in turn defined by such microscopic properties as the friction between water and capillary surfaces, and the viscosity of water. We show that the viscosity of water and, therefore, its flow rate are profoundly affected by the layered structure of confined water if the capillary size becomes less than 2 nm. To this end we study the structure and dynamics of water confined between two parallel graphene layers using equilibrium molecular dynamics simulations. We find that the shear viscosity is not only greatly enhanced for subnanometer capillaries, but also exhibits large oscillations that originate from commensurability between the capillary size and the size of water molecules. Such oscillating behavior of viscosity and, consequently, the slip length should be taken into account in designing and studying graphene-based and similar membranes for desalination and filtration.

preprint2016arXiv

Comment on "Creating in-plane pseudomagnetic fields in excess of 1000 T by misoriented stacking in a graphene bilayer"

In a recent paper [Phys. Rev. B 89, 125418 (2014)], the authors argue that it is possible to map the electronic properties of twisted bilayer graphene to those of bilayer graphene in an in-plane magnetic field. However, their description of the low-energy dynamics of twisted bilayer graphene is restricted to the extended zone scheme and therefore neglects the effects of the superperiodic structure. If the energy spectrum is studied in the supercell Brillouin zone, we find that the comparison with an in-plane magnetic field fails because (i) the energy spectra of the two situations exhibit different symmetries and (ii) the low-energy spectra are very different.

preprint2016arXiv

Computing Optical Properties of Ultra-thin Crystals

An overview is given of recent advances in experimental and theoretical understanding of optical properties of ultra-thin crystal structures (graphene, phosphorene, silicene, MoS2, MoSe2 , WS2 , WSe2 , h-AlN, h-BN, fluorographene, graphane). Ultra-thin crystals are atomically-thick layered crystals that have unique properties which differ from their 3D counterpart. Because of the difficulties in the synthesis of few-atom-thick crystal structures, which are thought to be the main building blocks of future nanotechnology, reliable theoretical predictions of their electronic, vibrational and optical properties are of great importance. Recent studies revealed the reliable predictive power of existing theoretical approaches based on density functional theory (DFT).

preprint2016arXiv

Electric field induced structural changes of water confined between two graphene layers

An external electric field changes the physical properties of polar-liquids due to the reorientation of their permanent dipoles. For example it should affect significantly the physical properties of water confined in a nanochannel. The latter effect is profoundly enhanced, if the field is applied along the nanochannel. Using molecular dynamics simulations, we predict that an in-plane electric field applied parallel to the channel polarizes water molecules which are confined between two graphene layers, resulting in distinct-ferroelectricity and electrical hysteresis. We found that electric fields alter the in-plane order of the hydrogen bonds: reversing the electric field does not restore the system to the non-polar initial state, instead a residual dipole moment remains in the system. Our study provides insights into the ferroelectric state of water when confined in nanochannels and shows how this can be tuned by an electric field.

preprint2016arXiv

Large gap electron-hole superfluidity and shape resonances in coupled graphene nanoribbons

We predict enhanced electron-hole superfluidity in two coupled electron-hole armchair-edge terminated graphene nanoribbons separated by a thin insulating barrier. In contrast to graphene monolayers, the multiple subbands of the nanoribbons are parabolic at low energy with a gap between the conduction and valence bands, and with lifted valley degeneracy. These properties make screening of the electron-hole interaction much weaker than for coupled electron-hole monolayers, thus boosting the pairing strength and enhancing the superfluid properties. The pairing strength is further boosted by the quasi-one-dimensional quantum confinement of the carriers, as well as by the large density of states near the bottom of each subband. The latter magnifies the superfluid shape resonances caused by the quantum confinement. Several superfluid partial condensates are present for finite-width nanoribbons with multiple subbands. We find that superfluidity is predominately in the strongly-coupled BEC and BCS-BEC crossover regimes, with large superfluid gaps up to 100 meV and beyond. When the gaps exceed the subband spacing, there is significant mixing of the subbands, a rounding of the shape resonances, and a resulting reduction in the one-dimensional nature of the system.

preprint2016arXiv

Mg(OH)2 -WS2 Heterobilayer: Electric Field Tunable Bandgap Crossover

Magnesium hydroxide (Mg(OH)2) has a layered brucite-like structure in its bulk form and was recently isolated as a new member of 2D monolayer materials. We investigated the electronic and optical properties of monolayer crystals of Mg(OH)2 and WS2 and their possible heterobilayer structure by means of first principles calculations. It was found that both monolayers of Mg(OH)2 and WS 2 are direct-gap semiconductors and these two monolayers form a typical van der Waals heterostructure with a weak interlayer interaction and a type-II band alignment with a staggered gap that spatially seperates electrons and holes. We also showed that an out-of-plane electric field induces a transition from a staggered to a straddling type heterojunction. Moreover, by solving the Bethe-Salpeter equation on top of single shot G0 W0 calculations, we show that the oscillator strength of the intralayer excitons of the heterostructure is an order of magnitude larger than the oscillator strength of the interlayer excitons. Because of the staggered interfacial gap and the field- tunable energy band structure, the Mg(OH)2 -WS2 heterobilayer can become an important candidate for various optoelectronic device applications in nanoscale.

preprint2016arXiv

Optical Properties of GaS-Ca(OH)$_2$ bilayer heterostructure

Finding novel atomically-thin heterostructures and understanding their characteristic properties are critical for developing better nanoscale optoelectronic devices. In this study, we investigate the electronic and optical properties of GaS-Ca(OH)$_2$ heterostructure using first-principle calculations. The band gap of the GaS-Ca(OH)$_2$ heterostructure is significantly reduced when compared with those of the isolated constituent layers. Our calculations show that the GaS-Ca(OH)$_2$ heterostructure is a type-II heterojunction which can be used to separate photoinduced charge carriers where electrons are localized in GaS and holes in the Ca(OH)$_2$ layer. This leads to spatially indirect excitons which are important for solar energy and optoelectronic applications due to their long lifetime. By solving the Bethe-Salpeter equation on top of single shot GW calculation (G$_0$W$_0$) the dielectric function and optical oscillator strength of the constituent monolayers and the heterostructure are obtained. The oscillator strength of the optical transition for GaS monolayer is an order of magnitude larger than Ca(OH)$_2$ monolayer. We also found that the calculated optical spectra of different stacking types of the heterostructure show dissimilarities, although their electronic structures are rather similar. This prediction can be used to determine the stacking type of ultra-thin heterostructures.

preprint2016arXiv

Piezoelectricity in asymmetrically strained bilayer graphene

We study the electronic properties of commensurate faulted bilayer graphene by diagonalizing the one-particle Hamiltonian of the bilayer system in a complete basis of Bloch states of the individual graphene layers. Our novel approach is very general and can be easily extended to any commensurate graphene-based heterostructure. Here, we consider three cases: i) twisted bilayer graphene, ii) bilayer graphene where triaxial stress is applied to one layer, and iii) bilayer graphene where uniaxial stress is applied to one layer. We show that the resulting superstructures can be divided into distinct classes, depending on the twist angle or the magnitude of the induced strain. The different classes are distinguished from each other by the interlayer coupling mechanism, resulting in fundamentally different low-energy physics. For the cases of triaxial and uniaxial stress, the individual graphene layers tend to decouple and we find significant charge transfer between the layers. In addition, this piezoelectric effect can be tuned by applying a perpendicular electric field. Finally, we show how our approach can be generalized to multilayer systems.

preprint2016arXiv

Quantum Transport Characteristics of Lateral pn-Junction of Single Layer TiS3

Using density functional theory and nonequilibrium Greens functions-based methods we investigated the electronic and transport properties of monolayer TiS3 pn-junction. We constructed a lateral pn-junction in monolayer TiS3 by using Li and F adatoms. An applied bias voltage caused significant variability in the electronic and transport properties of the TiS3 pn-junction. In addition, spin dependent current-voltage characteristics of the constructed TiS3 pn-junction were analyzed. Important device characteristics were found such as negative differential resistance and rectifying diode behaviors for spin-polarized currents in the TiS3 pn-junction. These prominent conduction properties of TiS3 pn-junction offer remarkable opportunities for the design of nanoelectronic devices based on a recently synthesized single-layered material.

preprint2016arXiv

Quantum transport in graphene Hall bars: Effects of vacancy disorder

Using the tight-binding model, we investigate the influence of vacancy disorder on electrical transport in graphene Hall bars in the presence of quantizing magnetic fields. Disorder, induced by a random distribution of monovacancies, breaks the graphene sublattice symmetry and creates states localized on the vacancies. These states are observable in the bend resistance, as well as in the total DOS. Their energy is proportional to the square root of the magnetic field, while their localization length is proportional to the cyclotron radius. At the energies of these localized states, the electron current flows around the monovacancies and, as we show, it can follow unexpected paths depending on the particular arrangement of vacancies. We study how these localized states change with the vacancy concentration, and what are the effects of including the next nearest neighbor hopping term. Our results are also compared with the situation when double vacancies are present in the system. Double vacancies also induce localized states, but their energy and magnetic field dependences are different. Their localization energy scales linearly with the magnetic field, and their localization length appears not to depend on the field strength.

preprint2016arXiv

Spin transitions driven by electric dipole spin resonance in fluorinated single- and bilayer-graphene quantum dots

Spin transitions driven by a periodically varying electric potential in dilute fluorinated graphene quantum dots are investigated. Flakes of monolayer graphene are considered as well as electrostatic electron traps induced in bilayer graphene. The stationary states are obtained within the tight-binding approach and are used to the basis of eigenstates to describe the system dynamics. The dilute fluorination of the top layer lifts the valley degeneracy of the confined states and attenuates the orbital magnetic dipole moments due to current circulation within the flake. Moreover, the spin-orbit coupling introduced by the surface deformation of the top layer induced by the adatoms allows spin flips to be driven by the AC electric field. For the bilayer quantum dots the spin flip times is substantially shorter than the experimental spin relaxation. Dynamical effects including many-photon and multilevel transitions are also discussed.

preprint2016arXiv

Strain controlled valley filtering in multi-terminal graphene structures

Valley-polarized currents can be generated by local straining of multi-terminal graphene devices. The pseudo-magnetic field created by the deformation allows electrons from only one valley to transmit and a current of electrons from a single valley is generated at the opposite side of the locally strained region. We show that valley filtering is most effective with bumps of a certain height and width. Despite the fact that the highest contribution to the polarized current comes from electrons from the lowest sub-band, contributions of other sub-bands are not negligible and can significantly enhance the output current.

preprint2016arXiv

Strain-induced topological phase transition in phosphorene and phosphorene nanoribbons

Using the tight-binding (TB) approximation with inclusion of the spin-orbit interaction, we predict a topological phase transition in the electronic band structure of phosphorene in the presence of axial strains. We derive a low-energy TB Hamiltonian that includes the spin-orbit interaction for bulk phosphorene. Applying a compressive biaxial in-plane strain and perpendicular tensile strain in ranges where the structure is still stable leads to a topological phase transition. We also examine the influence of strain on zigzag phosphorene nanoribbons (zPNRs) and the formation of the corresponding protected edge states when the system is in the topological phase. For zPNRs up to a width of 100 nm the energy gap is at least three orders of magnitude larger than the thermal energy at room temperature.

preprint2016arXiv

Strained graphene Hall bar

The effects of strain, induced by a Gaussian bump, on the magnetic field dependent transport properties of a graphene Hall bar are investigated. The numerical simulations are performed using both classical and quantum mechanical transport theory and we found that both approaches exhibit similar characteristic features. The effects of the Gaussian bump are manifested by a decrease of the bend resistance, $R_B$, around zero-magnetic field and the occurrence of side-peaks in $R_B$. These features are explained as a consequence of bump-assisted scattering of electrons towards different terminals of the Hall bar. Using these features we are able to give an estimate of the size of the bump. Additional oscillations in $R_B$ are found in the quantum description that are due to the population/depopulation of Landau levels. The bump has a minor influence on the Hall resistance even for very high values of the pseudo-magnetic field. When the bump is placed outside the center of the Hall bar valley polarized electrons can be collected in the leads.

preprint2016arXiv

The split-operator technique for the study of spinorial wavepacket dynamics

The split-operator technique for wave packet propagation in quantum systems is expanded here to the case of propagating wave functions describing Schrödinger particles, namely, charge carriers in semiconductor nanostructures within the effective mass approximation, in the presence of Zeeman effect, as well as of Rashba and Dresselhaus spin-orbit interactions. We also demonstrate that simple modifications to the expanded technique allow us to calculate the time evolution of wave packets describing Dirac particles, which are relevant for the study of transport properties in graphene.

preprint2016arXiv

Theoretical investigation of electron-hole complexes in anisotropic two-dimensional materials

Trions and biexcitons in anisotropic two-dimensional materials are investigated within an effective mass theory. Explicit results are obtained for phosphorene and arsenene, materials that share features such as a direct quasi-particle gap and anisotropic conduction and valence bands. Trions are predicted to have remarkably high binding energies and an elongated electron-hole structure with a preference for alignment along the armchair direction, where the effective masses are lower. We find that biexciton binding energies are also notably large, especially for monolayer phosphorene, where they are found to be twice as large as those for typical monolayer transition metal dichalcogenides.

preprint2016arXiv

Transport properties of bilayer graphene in a strong in-plane magnetic field

A strong in-plane magnetic field drastically alters the low-energy spectrum of bilayer graphene by separating the parabolic energy dispersion into two linear Dirac cones. The effect of this dramatic change on the transport properties strongly depends on the orientation of the in-plane magnetic field with respect to the propagation direction of the charge carriers and the angle at which they impinge on the electrostatic potentials. For magnetic fields oriented parallel to the potential boundaries an additional propagating mode that results from the splitting into Dirac cones enhances the transmission probability for charge carriers tunneling through the potentials and increases the corresponding conductance. Our results show that the chiral suppression of transmission at normal incidence is turned into a chiral enhancement when the magnetic field increases, thus indicating a transition from a bilayer to a monolayer-like system at normal incidence. Further, we find that the typical transmission resonances stemming from confinement in a potential barrier are shifted to higher energy and are eventually transformed into anti-resonances with increasing magnetic field. For magnetic fields oriented perpendicular to the potential boundaries we find a very pronounced transition from a bilayer system to two separated monolayer-like systems with Klein tunneling emerging at certain incident angles symmetric around 0, which also leaves a signature in the conductance. For both orientations of the magnetic field, the transmission probability is still correctly described by pseudospin conservation. Finally, to motivate the large in-plane magnetic field, we show that its energy spectrum can be mimicked by specific lattice deformations such as a relative shift of one of the layers. With this equivalence we introduce the notion of an in-plane pseudo-magnetic field.

preprint2016arXiv

Valley filtering using electrostatic potentials in bilayer graphene

Propagation of an electron wave packet through a quantum point contact (QPC) defined by electrostatic gates in bilayer graphene is investigated. The gates provide a bias between the layers, in order to produce an energy gap. If the gates on both sides of the contact produce the same bias, steps in the electron transmission probability are observed, as in the usual QPC. However, if the bias is inverted on one of the sides of the QPC, only electrons belonging to one of the Dirac valleys are allowed to pass, which provides a very efficient valley filtering.

preprint2015arXiv

AA-stacked bilayer square ice between graphene layers?

Water confined between two layers with separation of a few Angstrom forms layered two- dimensional ice structure. Using large scale molecular dynamics simulations with the adoptable ReaxFF interatomic potential we found that flat monolayer ice with a rhombic-square structure nucleates between graphene layers which is non-polar and non-ferroelectric. Two layers of water are found to crystallize into a square lattice close to the experimental found AA-stacking [G. Algara- Siller et al. Nature 519, 443445 (2015)]. Each layer has a net dipole moment which are in opposite direction. Bilayer ice is also non-polar and non-ferroelectric. For three layer ice we found that each layer has a crystal structure similar to monolayer ice.

preprint2015arXiv

Ag and Au Atoms Intercalated in Bilayer Heterostructures of Transition Metal Dichalcogenides and Graphene

The diffusive motion of metal nanoparticles Au and Ag on monolayer and between bilayer heterostructures of transition metal dichalcogenides and graphene are investigated in the framework of density functional theory. We found that the minimum energy barriers for diffusion and the possibility of cluster formation depend strongly on both the type of nanoparticle and the type of monolayers and bilayers. Moreover, the tendency to form clusters of Ag and Au can be tuned by creating various bilayers. Tunability of the diffusion characteristics of adatoms in van der Waals heterostructures holds promise for controllable growth of nanostructures.

preprint2015arXiv

Aharonov-Bohm oscillations in a mesoscopic ring with asymmetric arm-dependent injection

Electron transport through mesoscopic, one-dimensional rings with asymmetric injection into the arms of the ring is studied, in the presence of a Aharonov-Bohm flux, by means of an appropriate $\boldsymbol{S}$ matrix. This matrix is expressed in terms of two parameters one of which ($λ$) accounts phenomenologically for this asymmetric injection into the arms of the ring. In addition, the effect of a scatterer placed in one arm of the ring is considered. Explicit expressions are obtained for the transmission as a function of the incident electron energy, the magnetic field, the asymmetry parameter $λ$, and the strength of the scatterer. Results of the literature for symmetric rings are described by $λ=1$ and readily recovered. We relate our results to rings of finite width.

preprint2015arXiv

An efficient finite-difference scheme for computation of electron states in free-standing and core-shell quantum wires

The electron states in axially symmetric quantum wires are computed by means of the effective-mass Schroedinger equation, which is written in cylindrical coordinates phi, rho, and z. We show that a direct discretization of the Schroedinger equation by central finite differences leads to a non-symmetric Hamiltonian matrix. Because diagonalization of such matrices is more complex it is advantageous to transform it in a symmetric form. This can be done by the Liouville-like transformation proposed by Rizea et al. (Comp. Phys. Comm. 179 (2008) 466-478), which replaces the wave function psi(rho) with the function F(rho)=psi(rho)sqrt(rho) and transforms the Hamiltonian accordingly. Even though a symmetric Hamiltonian matrix is produced by this procedure, the computed wave functions are found to be inaccurate near the origin, and the accuracy of the energy levels is not very high. In order to improve on this, we devised a finite-difference scheme which discretizes the Schroedinger equation in the first step, and then applies the Liouville-like transformation to the difference equation. Such a procedure gives a symmetric Hamiltonian matrix, resulting in an accuracy comparable to the one obtained with the finite element method. The superior efficiency of the new finite-difference (FDM) scheme is demonstrated for a few rho-dependent one-dimensional potentials which are usually employed to model the electron states in free-standing and core-shell quantum wires. The new scheme is compared with the other FDM schemes for solving the effective-mass Schroedinger equation, and is found to deliver energy levels with much smaller numerical error for all the analyzed potentials. Moreover, the PT symmetry is invoked to explain similarities and differences between the considered FDM schemes.

preprint2015arXiv

Anisotropic exciton Stark shift in black phosphorus

We calculate the excitonic spectrum of few-layer black phosphorus by direct diagonalization of the effective mass Hamiltonian in the presence of an applied in-plane electric field. The strong attractive interaction between electrons and holes in this system allows one to investigate the Stark effect up to very high ionizing fields, including also the excited states. Our results show that the band anisotropy in black phosphorus becomes evident in the direction dependent field induced polarizability of the exciton.

preprint2015arXiv

Chain-like transitions in Wigner crystals: Sequential or non-sequential?

The structural transitions of the ground state of a system of repulsively interacting particles confined in a quasi-one-dimensional channel, and the effect of the interparticle interaction as well as the functional form of the confinement potential on those transitions are investigated. Although the non-sequential ordering of transitions (non-SOT), i.e. $1$-$2$-$4$-$3$-$4$-$5$-$6$-$...$ sequence of chain configurations with increasing density, is widely robust as predicted in a number of theoretical studies, the sequential ordering of transitions (SOT), i.e.~$1$-$2$-$3$-$4$-$5$-$6$-$...$ chain, is found as the ground state for long-ranged interparticle interaction and hard-wall-like confinement potentials. We found an energy barrier between every two different phases around its transition point, which plays an important role in the preference of the system to follow either a SOT or a non-SOT. However, that preferential transition requires also the stability of the phases during the transition. Additionally, we analyze the effect of a small structural disorder on the transition between the two phases around its transition point. Our results show that a small deformation of the triangular structure, change dramatically the picture of the transition between two phases, removing in a considerable region the non-SOT in the system. This feature could explain the fact that the non-SOT is, up to now, not observed in experimental systems, and suggests a more advanced experimental set-up to detect the non-SOT.

preprint2015arXiv

Distinct magnetic signatures of fractional vortex configurations in multiband superconductors

Vortices carrying fractions of a flux quantum are predicted to exist in multiband superconductors, where vortex core can split between multiple band-specific components of the superconducting condensate. Using the two-component Ginzburg-Landau model, we examine such vortex configurations in a two-band superconducting slab in parallel magnetic field. The fractional vortices appear due to the band-selective vortex penetration caused by different thresholds for vortex entry within each band-condensate, and stabilize near the edges of the sample. We show that the resulting fractional vortex configurations leave distinct fingerprints in the static measurements of the magnetization, as well as in ac dynamic measurements of the magnetic susceptibility, both of which can be readily used for the detection of these fascinating vortex states in several existing multiband superconductors.

preprint2015arXiv

Electron polarization function and plasmons in metallic armchair graphene nanoribbons

We calculate the polarization function of Dirac fermions in metallic armchair graphene nanoribbons for an arbitrary temperature and doping. We find that at finite temperatures due to the phase space redistribution among inter-band and intra-band electronic transitions in the conduction and valence bands, the full polarization function becomes independent of the temperature and the position of the chemical potential. As a result, for a given width of nanoribbons there exists a single plasmon mode, with the energy dispersion determined by the graphene's fine structure constant. In Coulomb-coupled nanoribbons, this plasmon splits into the basic in-phase and out-of-phase plasmon modes, with the splitting energy determined additionally by the inter-ribbon spacing.

preprint2015arXiv

Electronic and Magnetic Properties of 1T-TiSe2 Nanoribbons

Motivated by the recent synthesis of single layer TiSe2 , we used state-of-the-art density functional theory calculations, to investigate the structural and electronic properties of zigzag and armchair- edged nanoribbons of this material. Our analysis reveals that, differing from ribbons of other ultra-thin materials such as graphene, TiSe2 nanoribbons have some distinctive properties. The electronic band gap of the nanoribbons decreases exponentially with the width and vanishes for ribbons wider than 20 Angstroms. For ultranarrow zigzag-edged nanoribbons we find odd-even oscillations in the band gap width, although their band structures show similar features. Moreover, our detailed magnetic-ground-state analysis reveals that zigzag and armchair edged ribbons have nonmagnetic ground states. Passivating the dangling bonds with hydrogen at the edges of the structures influences the band dispersion. Our results shed light on the characteristic properties of T phase nanoribbons of similar crystal structures.

preprint2015arXiv

Giant paramagnetic Meissner effect in multiband superconductors

Superconductors, ideally diamagnetic when in the Meissner state, can also exhibit paramagnetic behavior due to trapped magnetic flux. In the absence of pinning such paramagnetic response is weak, and ceases with increasing sample thickness. Here we show that in multiband superconductors paramagnetic response can be observed even in slab geometries, and can be far larger than any previous estimate - even multiply larger than the diamagnetic Meissner response for the same applied magnetic field. We link the appearance of this giant paramagnetic response to the broad crossover between conventional Type-I and Type-II superconductors, where Abrikosov vortices interact non-monotonically and multibody effects become important, causing unique flux configurations and their locking in the presence of surfaces.

preprint2015arXiv

GraphAne: From Synthesis to Applications

Atomically thin crystals have recently been the focus of attention in particular after the synthesis of graphene, a monolayer hexagonal crystal structure of carbon. In this novel material class the chemically derived graphenes have attracted tremendous interest. It was shown that although bulk graphite is a chemically inert material, the surface of single layer graphene is rather reactive against individual atoms. So far, synthesis of several graphene derivatives have been reported such as hydrogenated graphene "graphane" (CH), fluorographene (CF) and chlorographene (CCl). Moreover, the stability of bromine and iodine covered graphene were predicted using computational tools. Among these derivatives, easy synthesis, insulating electronic behavior and reversibly tunable crystal structure of graphane make this material special for future ultra-thin device applications. This overview, surveys structural, electronic, magnetic, vibrational and mechanical properties of graphane. We also present a detailed overview of research efforts devoted to the computational modeling of graphane and its derivatives. Furthermore recent progress in synthesis techniques and possible applications of graphane are reviewed as well.

preprint2015arXiv

Graphene Ripples as a Realization of a Two-Dimensional Ising Model: A Scanning Tunneling Microscope Study

Ripples in pristine freestanding graphene naturally orient themselves in an array that is alternately curved-up and curved-down; maintaining an average height of zero. Using scanning tunneling microscopy (STM) to apply a local force, the graphene sheet will reversibly rise and fall in height until the height reaches 60-70 percent of its maximum at which point a sudden, permanent jump occurs. We successfully model the ripples as a spin-half Ising magnetic system, where the height of the graphene is the spin. The permanent jump in height, controlled by the tunneling current, is found to be equivalent to an antiferromagnetic-to-ferromagnetic phase transition. The thermal load underneath the STM tip alters the local tension and is identified as the responsible mechanism for the phase transition. Four universal critical exponents are measured from our STM data, and the model provides insight into the statistical role of graphenes unusual negative thermal expansion coefficient.

preprint2015arXiv

Hamiltonian of a many-electron system with single-electron and electron-pair states in a two-dimensional periodic potential

Based on the metastable electron-pair energy band in a two-dimensional (2D) periodic potential obtained previously by Hai and Castelano [J. Phys.: Condens. Matter 26, 115502 (2014)], we present in this work a Hamiltonian of many electrons consisting of single electrons and electron pairs in the 2D system. The electron-pair states are metastable of energies higher than those of the single-electron states at low electron density. We assume two different scenarios for the single-electron band. When it is considered as the lowest conduction band of a crystal, we compare the obtained Hamiltonian with the phenomenological model Hamiltonian of a boson-fermion mixture proposed by Friedberg and Lee [Phys. Rev. B 40, 6745 (1989)]. Single-electron-electron-pair and electron-pair-electron-pair interaction terms appear in our Hamiltonian and the interaction potentials can be determined from the electron-electron Coulomb interactions. When we consider the single-electron band as the highest valence band of a crystal, we show that holes in this valence band are important for stabilization of the electron-pair states in the system.

preprint2015arXiv

Hexagonal AlN: Dimensional-Crossover-Driven Bandgap Transition

Motivated by a recent experiment that reported the successful synthesis of hexagonal (h) AlN [Tsipas et al. Appl. Phys. Lett. 103, 251605 (2013)] we investigate structural, electronic and vibrational properties of bulk, bilayer and monolayer structures of h-AlN by using first-principles calculations. We show that the hexagonal phase of the bulk h-AlN is a stable direct-bandgap semiconductor. Calculated phonon spectrum displays a rigid-layer shear mode at 274 cm-1 and an Eg mode at 703 cm-1 which are observable by Raman measurements. In addition, single layer h-AlN is an indirect-bandgap semiconductor with a nonmagnetic ground state. For the bilayer structure, AA' type stacking is found to be the most favorable one and interlayer interaction is strong. While N-layered h-AlN is an indirect bandgap semiconductor for N=1-10, we predict that thicker structures (N>10) have a direct-bandgap at the Gamma-point. The number-of-layer-dependent bandgap transitions in h-AlN is interesting in that it is significantly different from the indirect-to- direct crossover obtained in the transition metal dichalcogenides.

preprint2015arXiv

High temperature conductance of a two-dimensional superlattice controlled by spin-orbit interaction

Rashba-type spin-orbit interaction (SOI) controlled band structure of a two-dimensional superlattice allows for the modulation of the conductance of finite size devices by changing the strength of the SOI. We consider rectangular arrays and find that the temperature dependence of the conductance disappears for high temperatures, but the strength of the SOI still affects the conductance at these temperatures. The modulation effect can be seen even in the presence of strong dephasing, which can be important in case of practical applications.

preprint2015arXiv

Magneto-optical transport properties of monolayer phosphorene

The electronic properties of monolayer phosphorene are exotic due to its puckered structure and large intrinsic direct band gap. We derive and discuss its band structure in the presence of a perpendicular magnetic field. Further, we evaluate the magneto-optical Hall and longitudinal optical conductivities, as functions of temperature, magnetic field, and Fermi energy, and show that they are strongly influenced by the magnetic field. The imaginary part of the former and the real part of the latter exhibit regular {\it interband} oscillations as functions of the frequency $ω$ in the range $\hslashω\sim 1.5-2$ eV. Strong {\it intraband} responses in the latter and week ones in the former occur at much lower frequencies. The magneto-optical response can be tuned in the microwave-to-terahertz and visible frequency ranges in contrast with a conventional two-dimensional electron gas or graphene in which the response is limited to the terahertz regime. This ability to isolate carriers in an anisotropic structure may make phosphorene a promising candidate for new optical devices.

preprint2015arXiv

Oscillating spin-orbit interaction in two-dimensional superlattices: sharp transmission resonances and time-dependent spin polarized currents

We consider ballistic transport through a lateral, two-dimensional superlattice with experimentally realizable, sinusoidally oscillating Rashba-type spin-orbit interaction. The periodic structure of the rectangular lattice produces a spin-dependent miniband structure for static SOI. Using Floquet theory, transmission peaks are shown to appear in the mini-bandgaps as a consequence of the additional, time-dependent SOI. A detailed analysis shows that this effect is due to the generation of harmonics of the driving frequency, via which e.g., resonances that cannot be excited in the case of static SOI become available. Additionally, the transmitted current shows space and time-dependent partial spin-polarization, in other words, polarization waves propagate through the superlattice.

preprint2015arXiv

Pentagonal Monolayer Crystals of Carbon, Boron Nitride, and Silver Azide

In this study we present a theoretical investigation of structural, electronic and mechanical properties of pentagonal monolayers of carbon (p-graphene), boron nitride (p-B$_{2}$N$_{4}$ and p-B$_{4}$N$_{2}$) and silver azide (p-AgN$_{3}$) by performing state-of-the-art first principles calculations. Our total energy calculations suggest feasible formation of monolayer crystal structures composed entirely of pentagons. In addition, electronic band dispersion calculations indicate that while p-graphene and p-AgN$_{3}$ are semiconductors with indirect bandgaps, p-BN structures display metallic behavior. We also investigate the mechanical properties (in-plane stiffness and the Poisson's ratio) of four different pentagonal structures under uniaxial strain. p-graphene is found to have the highest stiffness value and the corresponding Poisson's ratio is found to be negative. Similarly, p-B$_{2}$N$_{4}$ and p-B$_{4}$N$_{2}$ have negative Poisson's ratio values. On the other hand, the p-AgN$_{3}$ has a large and positive Poisson's ratio. In dynamical stability tests based on calculated phonon spectra of these pentagonal monolayers, we find that only p-graphene and p-B$_{2}$N$_{4}$ are stable, but p-AgN$_{3}$ and p-B$_{4}$N$_{2}$ are vulnerable against vibrational excitations.

preprint2015arXiv

Spin- and valley-dependent transport through arrays of ferromagnetic silicene junctions

We study ballistic transport of Dirac fermions in silicene through arrays of barriers, of width $d$, in the presence of an exchange field $M$ and a tunable potential of height $U$ or depth $-U$. The spin- and valley-resolved conductances as functions of $U$ or $M$, exhibit resonances away from the Dirac point (DP) and close to it a pronounced dip that becomes a gap when a critical electric field $E_z$ is applied. This gap widens by increasing the number of barriers and can be used to realize electric field-controlled switching of the current. The spin $p_s$ and valley $p_v$ polarizations of the current near the DP increase with $E_z$ or $M$ and can reach 100\% for certain of their values. These field ranges widen significantly by increasing the number of barriers. Also, $p_s$ and $p_v$ oscillate nearly periodically with the separation between barriers or wells and can be inverted by reversing $M$.

preprint2015arXiv

Spin-orbit interaction controlled properties of two-dimensional superlattices: Spintronic crystals

The band structure of two-dimensional artificial superlattices in the presence of (Rashba-type) spin-orbit interaction (SOI) is presented. The position and shape of the energy bands in these spintronic crystals depend on the geometry as well as the strength of the SOI, which can be tuned by external gate voltages. For finite mesoscopic arrays we show that their conductance properties can be understood from these spin-dependent band diagrams.

preprint2015arXiv

Stable Half-Metallic Monolayers of FeCl$_2$

The structural, electronic and magnetic properties of single layers of Iron Dichloride (FeCl$_{2}$) were calculated using first principles calculations. We found that the 1T phase of the single layer FeCl$_{2}$ is 0.17 eV/unit cell more favorable than its 1H phase. The structural stability is confirmed by phonon calculations. We found that 1T-FeCl$_{2}$ possess three Raman-active (130, 179 and 237 cm$^{-1}$) and one Infrared-active (279 cm$^{-1}$) phonon branches. The electronic band dispersion of the 1T-FeCl$_{2}$ is calculated using both GGA-PBE and DFT-HSE06 functionals. Both functionals reveal that the 1T-FeCl$_{2}$ has a half-metallic ground state with a Curie temperature of 17 K.

preprint2015arXiv

Theory of anharmonic phonons in 2D crystals

Anharmonic effects in an atomic monolayer thin crystal with honeycomb lattice structure are investigated by analytical and numerical lattice dynamical methods. Starting from a semi-empirical model for anharmonic couplings of third and fourth order, we study the in-plane and out-of-plane (flexural) mode components of the generalized wave vector dependent Grüneisen parameters, the thermal tension and the thermal expansion coefficients as function of temperature and crystal size. From the resonances of the displacement-displacement correlation functions we obtain the renormalization and decay rate of in-plane and flexural phonons as function of temperature, wave vector and crystal size in the classical and in the quantum regime. Quantitative results are presented for graphene. There we find that the transition temperature from negative to positive thermal expansion is lowered with smaller system size and by renormalization of the flexural mode frequency. Extensive numerical analysis throughout the Brillouin zone explores various decay and scattering channels. The relative importance of Normal and Umklapp processes is investigated. The work is complementary to crystalline membrane theory and computational studies of anharmonic effects in two-dimensional crystals.

preprint2015arXiv

Vacancy Formation and Oxidation Characteristics of Single Layer TiS3

The structural, electronic, and magnetic properties of pristine, defective, and oxidized monolayer TiS3 are investigated using first-principles calculations in the framework of density functional theory. We found that a single layer of TiS3 is a direct band gap semiconductor, and the bonding nature of the crystal is fundamentally different from other transition metal chalcogenides. The negatively charged surfaces of single layer TiS3 makes this crystal a promising material for lubrication applications. The formation energies of possible vacancies, i.e. S, Ti, TiS, and double S, are investigated via total energy optimization calculations. We found that the formation of a single S vacancy was the most likely one among the considered vacancy types. While a single S vacancy results in a nonmagnetic, semiconducting character with an enhanced band gap, other vacancy types induce metallic behavior with spin polarization of 0.3-0.8 μB. The reactivity of pristine and defective TiS3 crystals against oxidation was investigated using conjugate gradient calculations where we considered the interaction with atomic O, O2, and O3. While O2 has the lowest binding energy with 0.05-0.07 eV, O3 forms strong bonds stable even at moderate temperatures. The strong interaction (3.9-4.0 eV) between atomic O and TiS3 results in dissociative adsorption of some O-containing molecules. In addition, the presence of S-vacancies enhances the reactivity of the surface with atomic O, whereas it had a negative effect on the reactivity with O2 and O3 molecules.

preprint2015arXiv

Veselago lensing in graphene with a p-n junction: classical versus quantum effects

The feasibility of Veselago lensing in graphene with a p-n junction is investigated numerically for realistic injection leads. Two different set-ups with two narrow leads are considered with absorbing or reflecting side edges. This allows us to separately determine the influence of scattering on electron focusing for the edges and the p-n interface. Both semiclassical and tight-binding simulations show a distinctive peak in the transmission probability that is attributed to the Veselago lensing effect. We investigate the robustness of this peak on the width of the injector, the position of the p-n interface and different gate potential profiles. Furthermore, the influence of scattering by both short- and long-range impurities is considered.

preprint2015arXiv

Vortical versus skyrmionic states in mesoscopic \emph{p}-wave superconductors

We investigate the superconducting states that arise as a consequence of mesoscopic confinement and a multi-component order parameter in the Ginzburg-Landau model for $p\,$-wave superconductivity. Conventional vortices, but also half-quantum vortices and skyrmions are found as the applied magnetic field and the anisotropy parameters of the Fermi surface are varied. The solutions are well differentiated by a topological charge that for skyrmions is given by the Hopf invariant and for vortices by the circulation of the superconducting velocity. We revealed several unique states combining vortices and skyrmions, their possible reconfiguration with varied magnetic field, as well as the novel temporal and field-induced transitions between vortical and skyrmionic states.

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

Electron energy and temperature relaxation in graphene on a piezoelectric substrate

We study the energy and temperature relaxation of electrons in graphene on a piezoelectric substrate. Scattering from the combined potential of extrinsic piezoelectric surface acoustical (PA) phonons of the substrate and intrinsic deformation acoustical (DA) phonons of graphene is considered in a (non-)degenerate gas of Dirac fermions. It is shown that in the regime of low energies or temperatures the PA phonons dominate the relaxation and change qualitatively its character. This prediction is relevant for quantum metrology and electronic applications using graphene devices and suggests a new experimental setup for probing electron-phonon coupling in graphene.

preprint2014arXiv

Electronic properties of graphene nano-flakes: Energy gap, permanent dipole, termination effect and Raman spectroscopy

The electronic properties of graphene nano-flakes (GNFs) with different edge passivation is investigated by using density functional theory. Passivation with F and H atoms are considered: C$_{N_c}$ X$_{N_x}$ (X=F or H). We studied GNFs with $10<N_c<56$ and limit ourselves to the lowest energy configurations. We found that: i) the energy difference $Δ$ between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) decreases with $N_c$, ii) topological defects (pentagon and heptagon) break the symmetry of the GNFs and enhance the electric polarization, iii) the mutual interaction of bilayer GNFs can be understood by dipole-dipole interaction which were found sensitive to the relative orientation of the GNFs, iv) the permanent dipoles depend on the edge terminated atom, while the energy gap is independent of it, and v) the presence of heptagon and pentagon defects in the GNFs results in the largest difference between the energy of the spin-up and spin-down electrons which is larger for the H-passivated GNFs as compared to F-passivated GNFs. Our study shows clearly the effect of geometry, size, termination and bilayer on the electronic properties of small GNFs.This study reveals important features of graphene nano-flakes which can be detected using Raman spectroscopy.

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

Enhancement of electron-hole superfluidity in double few-layer graphene

We propose two coupled electron-hole sheets of few-layer graphene as a new nanostructure to observe superfluidity at enhanced densities and enhanced transition temperatures. For ABC stacked few-layer graphene we show that the strongly correlated electron-hole pairing regime is readily accessible experimentally using current technologies. We find for double trilayer and quadlayer graphene sheets spatially separated by a nano-thick hexagonal boron-nitride insulating barrier, that the transition temperature for electron-hole superfluidity can approach temperatures of 40 K.

preprint2014arXiv

Generic ordering of structural transitions in quasi-one-dimensional Wigner crystals

We investigate the dependence of the structural phase transitions in an infinite quasi-one-dimensional system of repulsively interacting particles on the profile of the confining channel. Three different functional expressions for the confinement potential related to real experimental systems are used that can be tuned continuously from a parabolic to a hard-wall potential in order to find a thorough understanding of the ordering of the chain-like structure transitions. We resolve the longstanding issue why the most theories predicted a 1-2-4-3-4 sequence of chain configurations with increasing density, while some experiments found the 1-2-3-4 sequence.

preprint2014arXiv

Graphene on boron-nitride: Moiré pattern in the van der Waals energy

The spatial dependence of the van der Waals (vdW) energy between graphene and hexagonal boron-nitride (h-BN) is investigated using atomistic simulations. The van der Waals energy between graphene and h-BN shows a hexagonal superlattice structure identical to the observed Moiré pattern in the local density of states (LDOS) which depends on the lattice mismatch and misorientation angle between graphene and h-BN. Our results provide atomistic features of the weak van der Waals interaction between graphene and BN which are in agreement with experiment and provide an analytical expression for the size of the spatial variation of the weak van der Waals interaction. We also found that the A-B-lattice symmetry of graphene is broken a long the armchair direction.

preprint2014arXiv

Graphene on hexagonal lattice substrate: Stress and Pseudo-magnetic field

Moir'e patterns in the pseudo-magnetic field and in the strain profile of graphene (GE) when put on top of a hexagonal lattice substrate are predicted from elasticity theory. %which are confirmed by atomistic simulations. The van der Waals (vdW) interaction between GE and the substrate induces out-of-plane deformations in graphene which results in a strain field, and consequently in a pseudo-magnetic field. When the misorientation angle is about 0.5 deg. a three-fold symmetric strain field is realized that results in a pseudo-magnetic field very similar to the one proposed by F. Guinea, M. I. Katsnelson, and A. K. Geim [Nat. Phys. 6, 30 (2010)]. Our results show that the periodicity and length of the pseudo-magnetic field can be tuned in GE by changing the misorientation angle and substrate adhesion parameters and a considerable energy gap (23 meV) can be obtained due to out-of-plane deformation of graphene which is in the range of recent experimental measurements (20-30 meV).

preprint2014arXiv

Interplay between snake and quantum edge states in a graphene Hall bar with a pn-junction

The magneto- and Hall resistance of a locally gated cross shaped graphene Hall bar is calculated. The edge of the top gate is placed diagonally across the center of the Hall cross. Four-probe resistance is calculated using the Landauer-Büttiker formalism, while the transmission coefficients are obtained using the non-equilibrium Green's function approach. The interplay between transport due to edge channels and snake states is investigated. When two edge channels are occupied we predict oscillations in the Hall and the bend resistance as function of the magnetic field which are a consequence of quantum interference between the occupied snake states.

preprint2014arXiv

Magnetic particles confined in a modulated channel: structural transitions tunable by tilting a magnetic field

The ground state of colloidal magnetic particles in a modulated channel are investigated as function of the tilt angle of an applied magnetic field. The particles are confined by a parabolic potential in the transversal direction while in the axial direction a periodic substrate potential is present. By using Monte Carlo (MC) simulations, we construct a phase diagram for the different crystal structures as a function of the magnetic field orientation, strength of the modulated potential and the commensurability factor of the system. Interestingly, we found first and second order phase transitions between different crystal structures, which can be manipulated by the orientation of the external magnetic field. A re-entrant behavior is found between two- and four-chain configurations, with continuous second order transitions. Novel configurations are found consisting of frozen in solitons. By changing the orientation and/or strength of the magnetic field and/or the strength and the spatial frequency of the periodic substrate potential, the system transits through different phases.

preprint2014arXiv

Membrane amplitude and triaxial stress in twisted bilayer graphene deciphered using first-principles directed elasticity theory and scanning tunneling microscopy

Twisted graphene layers produce a moiré pattern (MP) structure with a predetermined wavelength for given twist angle. However, predicting the membrane corrugation amplitude for any angle other than pure AB-stacked or AA-stacked graphene is impossible using first-principles density functional theory (DFT) due to the large supercell. Here, within elasticity theory we define the MP structure as the minimum energy configuration, thereby leaving the height amplitude as the only unknown parameter. The latter is determined from DFT calculations for AB and AA stacked bilayer graphene in order to eliminate all fitting parameters. Excellent agreement with scanning tunneling microscopy (STM) results across multiple substrates is reported as function of twist angle.

preprint2014arXiv

Multilayer graphene, Moiré patterns, grain boundaries and defects identified by scanning tunneling microscopy on the m-plane, non-polar surface of SiC

Epitaxial graphene is grown on a non-polar n+ 6H-SiC m-plane substrate and studied using atomic scale scanning tunneling microscopy. Multilayer graphene is found throughout the surface and exhibits rotational disorder. Moiré patterns of different spatial periodicities are found, and we found that as the wavelength increases, so does the amplitude of the modulations. This relationship reveals information about the interplay between the energy required to bend graphene and the interaction energy, i.e. van der Waals energy, with the graphene layer below. Our experiments are supported by theoretical calculations which predict that the membrane topographical amplitude scales with the Moiré pattern wavelength, L as L^-1 + αL^-2.

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

Oscillating spin-orbit interaction as a source of spin-polarized wave packets in two-terminal nanoscale devices

Ballistic transport through nanoscale devices with time-dependent Rashba-type spin-orbit interaction (SOI) can lead to spin-polarized wave packets that appear even for completely unpolarized input. The SOI that oscillates in a finite domain generates density and spin polarization fluctuations that leave the region as propagating waves. Particularly, spin polarization has space and time dependence even in regions without SOI. Our results are based on an analytic solution of the time-dependent Schrödinger equation. The relevant Floquet quasi-energies that are obtained appear in the energy spectrum of both the transmitted and reflected waves.

preprint2014arXiv

Plasmons and their interaction with electrons in trilayer graphene

The interaction between electrons and plasmons in trilayer graphene is investigated within the Overhauser approach resulting in the 'plasmaron' quasi-particle. This interaction is cast into a field theoretical problem, nd its effect on the energy spectrum is calculated using improved Wigner-Brillouin perturbation theory. The plasmaron spectrum is shifted with respect to the bare electron spectrum by $ΔE(\mathbf{k})\sim 50\div200\,{\rm meV}$ for ABC stacked trilayer graphene and for ABA trilayer graphene by $ΔE(\mathbf{k})\sim 30\div150\,{\rm meV}$ ($ΔE(\mathbf{k})\sim 1\div5\,{\rm meV}$) for the hyperbolic linear) part of the spectrum. The shift in general increases with the electron concentration $n_{e}$ and electron momentum. The dispersion of plasmarons is more pronounced in \textit{ABC} stacked than in ABA tacked trilayer graphene, because of the different energy band structure and their different plasmon dispersion.

preprint2014arXiv

Realization of free-standing silicene using bilayer graphene

The available synthesized silicene-like structures have been only realized on metallic substrates which are very different from the standalone buckled silicene, e.g. the Dirac cone of silicene is destroyed due to lattice distortion and the interaction with the substrate. Using graphene bilayer as a scaffold a novel route is proposed to synthesize silicene with electronic properties decoupled from the substrate. The buckled hexagonal arrangement of silicene between the graphene layers is found to be very similar to the theoretically predicted standalone buckled silicene which is only very weakly van der Waals coupled to the graphene layers with a graphite-like interlayer distance of 3.42\,Åand without only lattice distortion. We found that these stacked layers are stable well above room temperature.

preprint2014arXiv

Rippling, buckling and melting of single- and multi-layer MoS$_2$

Large-scale atomistic simulations using the reactive empirical bond order force field approach is implemented to investigate thermal and mechanical properties of single-layer (SL) and multi-layer (ML) molybdenum disulfide (MoS$_2$). The amplitude of the intrinsic ripples of SL-MoS$_2$ are found to be smaller than those exhibited by graphene (GE). Furthermore, because of the van der Waals interaction between layers, the out-of-plane thermal fluctuations of ML-MoS$_2$ decreases rapidly with increasing number of layers. This trend is confirmed by the buckling transition due to uniaxial stress which occurs for a significantly larger applied tension as compared to graphene. For SL-MoS$_2$, the melting temperature is estimated to be 3700~K which occurs through dimerization followed by the formation of small molecules consisting of 2 to 5 atoms. When different types of vacancies are inserted in the SL-MoS$_2$ it results in a decrease of both the melting temperature as well as the stiffness.

preprint2014arXiv

Spin and valley polarization of plasmons in silicene due to external fields

The electronic properties of the novel two dimensional (2D) material silicene are strongly influenced by the application of a perpendicular electric field $E_z$ and of an exchange field $M$ due to adatoms positioned on the surface or a ferromagnetic substrate. Within the random phase approximation, we investigate how electron-electron interactions are affected by these fields and present analytical and numerical results for the dispersion of plasmons, their lifetime, and their oscillator strength. We find that the combination of the fields $E_z$ and $M$ brings a spin and valley texture to the particle-hole excitation spectrum and allows the formation of spin- and valley-polarized plasmons. When the Fermi level lies in the gap of one spin in one valley, the intraband region of the corresponding spectrum disappears. For zero $E_z$ and finite $M$ the spin symmetry is broken and spin polarization is possible. The lifetime and oscillator strength of the plasmons are shown to depend strongly on the number of spin and valley type electrons that form the electron-hole pairs.

preprint2014arXiv

Stabilized silicene within bilayer graphene: A proposal based on molecular dynamics and density-functional tight-binding calculations

Free standing silicene is predicted to display comparable electronic properties as graphene. However, the yet synthesized silicene-like structures have been only realized on different substrates which turned out to exhibit versatile crystallographic structures that are very different from the theoretically predicted buckled phase of freestanding silicene. This calls for a different approach where silicene is stabilized using very weakly interacting surfaces. We propose here a novel route by using graphene bilayer as a scaffold. The confinement between the flat graphene layers results in a planar clustering of Si atoms with small buckling, which is energetically unfavorable in vacuum. Buckled hexagonal arrangement of Si atoms similar to free-standing silicene is observed for large clusters, which, in contrast to Si atoms on metallic surfaces, is only very weakly van der Waals coupled to the graphene layers. These clusters are found to be stable well above room temperature. Our findings, which are supported by density functional tight-binding calculations, show that intercalating bilayer graphene with Si is a favorable route to realize silicene.

preprint2014arXiv

Temporary cooling of quasiparticles and delay in voltage response of superconducting bridges after abrupt switching on the supercritical current

We revisit the problem of the dynamic response of a superconducting bridge after abruptly switching on the supercritical current $I>I_c$. In contrast to previous theoretical works we take into account spatial gradients and use both the local temperature approach and the kinetic equation for the distribution function of quasiparticles. In both models the finite delay time $t_d$ in the voltage response is connected with temporary cooling of quasiparticles due to the suppression of the superconducitng order parameter by current. We find that $t_d$ has different values and different temperature dependencies in the considered models. In turns out that the presence of even small inhomogeneities in the bridge or of bulk leads/contacts at the ends of the {\it homogenous} bridge favors a local suppression of the superconducting order parameter $|Δ|$ during the dynamic response. It results in a decrease of the delay time, in comparison with the spatially uniform model, due to the diffusion of nonequilibrium quasiparticles from the region with locally suppressed $|Δ|$. In case the current distribution is spatially nonuniform across the bridge the delay time is mainly connected with the time needed for the nucleation of the first vortex at the position where the current density is maximal (at $I\sim I_c$ and for not very wide films). We also find that a short alternating current pulse (sinusoid like) with zero time-average may result in a nonzero time-averaged voltage response where its sign depends on the phase of the ac current.

preprint2014arXiv

The role of atomic vacancies and boundary conditions on ballistic thermal transport in graphene nanoribbons

Quantum thermal transport in armchair and zig-zag graphene nanoribbons are investigated in the presence of single atomic vacancies and subject to different boundary conditions. We start with a full comparison of the phonon polarizations and energy dispersions as given by a fifth-nearest-neighbor force-constant model (5NNFCM) and by elasticity theory of continuum membranes (ETCM). For free-edges ribbons we discuss the behavior of an additional acoustic edge-localized flexural mode, known as fourth acoustic branch (4ZA), which has a small gap when it is obtained by the 5NNFCM. Then, we show that ribbons with supported-edges have a sample-size dependent energy gap in the phonon spectrum which is particularly large for in-plane modes. Irrespective to the calculation method and the boundary condition, the dependence of the energy gap for the low-energy optical phonon modes against the ribbon width W is found to be proportional to 1/W for in-plane, and 1/W$^2$ for out-of-plane phonon modes. Using the 5NNFCM, the ballistic thermal conductance and its contributions from every single phonon mode are then obtained by the non equilibrium Green's function technique. We found that, while edge and central localized single atomic vacancies do not affect the low-energy transmission function of in-plane phonon modes, they reduce considerably the contributions of the flexural modes. On the other hand, in-plane modes contributions are strongly dependent on the boundary conditions and at low temperatures can be highly reduced in supported-edges samples. These findings could open a route to engineer graphene based devices where it is possible to discriminate the relative contribution of polarized phonons and to tune the thermal transport on the nanoscale.

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.

preprint2014arXiv

Thermal mirror buckling in freestanding graphene locally controlled by scanning tunneling microscopy

Knowledge of and control over the curvature of ripples in freestanding graphene are desirable for fabricating and designing flexible electronic devices, and recent progress in these pursuits has been achieved using several advanced techniques such as scanning tunneling microscopy. The electrostatic forces induced through a bias voltage (or gate voltage) were used to manipulate the interaction of freestanding graphene with a tip (substrate). Such forces can cause large movements and sudden changes in curvature through mirror buckling. Here we explore an alternative mechanism, thermal load, to control the curvature of graphene. We demonstrate thermal mirror buckling of graphene by scanning tunneling microscopy and large-scale molecular dynamic simulations. The negative thermal expansion coefficient of graphene is an essential ingredient in explaining the observed effects. This new control mechanism represents a fundamental advance in understanding the influence of temperature gradients on the dynamics of freestanding graphene and future applications with electro-thermal-mechanical nanodevices.

preprint2014arXiv

Thermodynamic properties of the electron gas in multilayer graphene in the presence of a perpendicular magnetic field

The thermodynamic properties of the electron gas in multilayer graphene depend strongly on the number of layers and the type of stacking. Here we analyse how those properties change when we vary the number of layers for rhombohedral stacked multilayer graphene and compare our results with those from a conventional two dimensional electron gas. We show that the highly degenerate zero energy Landau level which is partly filled with electrons and partly with holes has a strong influence on the value of the different thermodynamic quantities.

preprint2014arXiv

Unusual ultralow frequency fluctuations in freestanding graphene

Intrinsic ripples in freestanding graphene have been exceedingly difficult to study. Individual ripple geometry was recently imaged using scanning tunneling microscopy, but these measurements are limited to static configurations. Thermally-activated flexural phonon modes should generate dynamic changes in curvature. Here we show how to track the vertical movement of a one-square-angstrom region of freestanding graphene using scanning tunneling microscopy, thereby allowing measurement of the out-of-plane time trajectory and fluctuations over long time periods. We also present a model from elasticity theory to explain the very-low-frequency oscillations. Unexpectedly, we sometimes detect a sudden colossal jump, which we interpret as due to mirror buckling. This innovative technique provides a much needed atomic-scale probe for the time-dependent behavior of intrinsic ripples. The discovery of this novel progenitor represents a fundamental advance in the use of scanning tunneling microscopy, which together with the application of a thermal load provides a low-frequency nano-resonator.

preprint2013arXiv

Anomalous Raman Spectra and Thickness Dependent Electronic properties of WSe2

Typical Raman spectra of transition metal dichalcogenides (TMDs) display two prominent peaks, E2g and A1g, that are well separated from each other. We find that these modes are degenerate in bulk WSe2 yielding one single Raman peak. As the dimensionality is lowered, the observed peak splits in two as a result of broken degeneracy. In contrast to our experimental findings, our phonon dispersion calculations reveal that these modes remain degenerate independent of the number of layers. Interestingly, for minuscule biaxial strain the degeneracy is preserved but once the crystal symmetry is broken by uniaxial strain, the degeneracy is lifted. Our calculated phonon dispersion for uniaxially strained WSe2 shows a perfect match to the measured Raman spectrum which suggests that uniaxial strain exists in WSe2 flakes possibly induced during the sample preparation and/or as a result of interaction between WSe2 and the substrate. Furthermore, we find that WSe2 undergoes an indirect to direct bandgap transition from bulk to monolayers which is ubiquitous for semiconducting TMDs. These results not only allow us to understand the vibrational properties of WSe2 but also provides detailed insight to their physical properties.

preprint2013arXiv

Bilayer graphene Hall bar with a pn-junction

We investigate the magnetic field dependence of the Hall and the bend resistances for a ballistic Hall bar structure containing a pn-junction sculptured from a bilayer of graphene. The electric response is obtained using the billiard model and we investigate the cases of bilayer graphene with and without a band gap. Two different conduction regimes are possible: $i$) both sides of the junction have the same carrier type, and $ii$) one side of the junction is n-type while the other one is p-type. The first case shows Hall plateau-like features in the Hall resistance that fade away as the band gap opens. The second case exhibits a bend resistance that is asymmetric in magnetic field as a consequence of snake states along the pn-interface, where the maximum is shifted away from zero magnetic field.

preprint2013arXiv

Boron Nitride Monolayer: A Strain-Tunable Nanosensor

The influence of triaxial in-plane strain on the electronic properties of a hexagonal boron-nitride sheet is investigated using density functional theory. Different from graphene, the triaxial strain localizes the molecular orbitals of the boron-nitride flake in its center depending on the direction of the applied strain. The proposed technique for localizing the molecular orbitals that are close to the Fermi level in the center of boron nitride flakes can be used to actualize engineered nanosensors, for instance, to selectively detect gas molecules. We show that the central part of the strained flake adsorbs polar molecules more strongly as compared with an unstrained sheet.

preprint2013arXiv

Braess paradox at the mesoscopic scale

We theoretically demonstrate that the transport inefficiency recently found experimentally for branched-out mesoscopic networks can also be observed in a quantum ring of finite width with an attached central horizontal branch. This is done by investigating the time evolution of an electron wave packet in such a system. Our numerical results show that the conductivity of the ring does not necessary improves if one adds an extra channel. This ensures that there exists a quantum analogue of the Braess Paradox, originating from quantum scattering and interference.

preprint2013arXiv

Dynamics of self-organized driven particles with competing range interaction

Non-equilibrium self-organized patterns formed by particles interacting through competing range interaction are driven over a substrate by an external force. We show that, with increasing driving force, the pre-existed static patterns evolve into dynamic patterns either via disordered phase or depinned patterns, or via the formation of non-equilibrium stripes. Strikingly, the stripes are formed either in the direction of the driving force or in the transverse direction, depending on the pinning strength. The revealed dynamical patterns are summarized in a dynamical phase diagram.

preprint2013arXiv

Effect of nonhomogenous dielectric background on the plasmon modes in graphene double-layer structures at finite temperatures

We have calculated the plasmon modes in graphene double layer structures at finite temperatures, taking into account the inhomogeneity of the dielectric background of the system. The effective dielectric function is obtained from the solution of the Poisson equation of three-layer dielectric medium with the graphene sheets located at the interfaces, separating the different materials. Due to the momentum dispersion of the effective dielectric function, the intra- and inter-layer bare Coulomb interactions in the graphene double layer system acquires an additional momentum dependence--an effect that is of the order of the inter-layer interaction itself. We show that the energies of the in-phase and out-of-phase plasmon modes are determined largely by different values of the spatially dependent effective dielectric function. The effect of the dielectric inhomogeneity increases with temperature and even at high temperatures the energy shift induced by the dielectric inhomogeneity and temperature itself remains larger than the broadening of the plasmon energy dispersions due to the Landau damping. The obtained new features of the plasmon dispersions can be observed in frictional drag measurements and in inelastic light scattering and electron energy-loss spectroscopies.

preprint2013arXiv

Electron-electron interactions in bilayer graphene quantum dots

A parabolic quantum dot (QD) as realized by biasing nanostructured gates on bilayer graphene is investigated in the presence of electron-electron interaction. The energy spectrum and the phase diagram reveal unexpected transitions as function of a magnetic field. For example, in contrast to semiconductor QDs, we find a novel valley transition rather than only the usual singlet-triplet transition in the ground state of the interacting system. The origin of these new features can be traced to the valley degree of freedom in bilayer graphene. These transitions have important consequences for cyclotron resonance experiments.

preprint2013arXiv

Electronic states in a graphene flake strained by a Gaussian bump

The effect of strain in graphene is usually modeled by a pseudo-magnetic vector potential which is, however, derived in the limit of small strain. In realistic cases deviations are expected in view of graphene's very high strain tolerance, which can be up to 25%. Here we investigate the pseudo-magnetic field generated by a Gaussian bump and we show that it exhibits significant differences with numerical tight-binding results. Furthermore, we calculate the electronic states in the strained region for a hexagon shaped flake with armchair edges. We find that the six-fold symmetry of the wave functions inside the Gaussian bump is directly related to the different effect of strain along the fundamental directions of graphene: zigzag and armchair. Low energy electrons are strongly confined in the armchair directions and are localized on the carbon atoms of a single sublattice.

preprint2013arXiv

Energy-Momentum dispersion relation of plasmarons in bilayer graphene

The relation between the energy and momentum of plasmarons in bilayer graphene is investigated within the Overhauser approach, where the electron-plasmon interaction is described as a field theoretical problem. We find that the Dirac-like spectrum is shifted by $ΔE(\mathbf{k})\sim 100\div150\,{\rm meV}$ depending on the electron concentration $n_{e}$ and electron momentum. The shift increases with electron concentration as the energy of plasmons becomes larger. The dispersion of plasmarons is more pronounced than in the case of single layer graphene, which is explained by the fact that the energy dispersion of electrons is quadratic and not linear. We expect that these predictions can be verified using angle-resolved photoemission spectroscopy (ARPES).

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

Four band tunneling in bilayer graphene

The conductance, the transmission and the reflection probabilities through rectangular potential barriers and pn-junctions are obtained for bilayer graphene taking into account the four bands of the energy spectrum. We have evaluated the importance of the skew hopping parameters γ3 and γ4 to these properties and show that for energies E>γ1/100 their effect is negligible. For high energies two modes of propagation exist and we investigate scattering between these modes. For perpendicular incidence both propagation modes are decoupled and scattering between them is forbidden. This extends the concept of pseudospin as defined within the two band approximation to a four band model and corresponds to the (anti)symmetry of the wavefunctions under in-plane mirroring. New transmission resonances are found that appear as sharp peaks in the conductance which are absent in the two band approximation. The application of an interlayer bias to the system: 1) breaks the pseudospin structure, 2) opens a bandgap that results in a distinct feature of suppressed transmission in the conductance, and 3) breaks the angular symmetry with respect to normal incidence in the transmission and reflection.

preprint2013arXiv

Ginzburg-Landau theory for multiband superconductors: microscopic derivation

A procedure to derive the Ginzburg-Landau (GL) theory from the multiband BCS Hamiltonian is developed in a general case with an arbitrary number of bands and arbitrary interaction matrix. It combines the standard Gor'kov truncation and a subsequent reconstruction in order to match accuracies of the obtained terms. This reconstruction recovers the phenomenological GL theory as obtained from the Landau model of phase transitions but offers explicit microscopic expressions for the relevant parameters. Detailed calculations are presented for a three-band system treated as a prototype multiband superconductor. It is demonstrated that the symmetry in the coupling matrix may lead to the chiral ground state with the phase frustration, typical for systems with broken time-reversal symmetry.

preprint2013arXiv

Graphene Hall bar with an asymmetric pn-junction

We investigated the magnetic field dependence of the Hall and the bend resistances in the ballistic regime for a single layer graphene Hall bar structure containing a pn-junction. When both regions are n-type the Hall resistance dominates and Hall type of plateaus are formed. These plateaus occur as a consequence of the restriction on the angle imposed by Snell's law allowing only electrons with a certain initial angles to transmit though the potential step. The size of the plateau and its position is determined by the position of the potential interface as well as the value of the applied potential. When the second region is p-type the bend resistance dominates which is asymmetric in field due to the presence of snake states. Changing the position of the pn-interface in the Hall bar strongly affects these states and therefore the bend resistance is also changed. Changing the applied potential we observe that the bend resistance exhibits a peak around the charge-neutrality point (CNP) which is independent of the position of the pn-interface, while the Hall resistance shows a sign reversal when the CNP is crossed, which is in very good agreement with a recent experiment [J. R. Williams et al., Phys. Rev. Lett. 107, 046602(2011)].

preprint2013arXiv

Helical liquid of snake states

We derive an exact solution to the problem of spin snake states induced in a nonhomogeneous magnetic field by a combined action of the Rashba spin-orbit and Zeeman fields. In an antisymmetric magnetic field the spin snake states are nonlocal composite particles, originating from spatially separated entangled spins. Adding an external homogeneous magnetic field breaks the spin-parity symmetry gapping out the spectral branches, which results in a regular beating pattern of the spin current. These new phenomena in a helical liquid of snake states are proposed for an experimental realization.

preprint2013arXiv

Klein paradox for a pn junction in multilayer graphene

Charge carriers in single and multilayered graphene systems behave as chiral particles due to the particular lattice symmetry of the crystal. We show that the interplay between the meta-material properties of graphene multilayers and the pseudospinorial properties of the charge carriers result in the occurrence of Klein and anti-Klein tunneling for rhombohedral stacked multilayers. We derive an algebraic formula predicting the angles at which these phenomena occur and support this with numerical calculations for systems up to four layers. We present a decomposition of an arbitrarily stacked multilayer into pseudospin doublets that have the same properties as rhombohedral systems with a lower number of layers.

preprint2013arXiv

Melting of graphene clusters

Density-functional tight-binding and classical molecular dynamics simulations are used to investigate the structural deformations and melting of planar carbon nano-clusters $C_{N}$ with N=2-55. The minimum energy configurations for different clusters are used as starting configuration for the study of the temperature effects on the bond breaking/rotation in carbon lines (N$<$6), carbon rings (5$<$N$<$19) and graphene nano-flakes. The larger the rings (graphene nano-flake) the higher the transition temperature (melting point) with ring-to-line (perfect-to-defective) transition structures. The melting point was obtained by using the bond energy, the Lindemann criteria, and the specific heat. We found that hydrogen-passivated graphene nano-flakes (C$_{N}$H$_M$) have a larger melting temperature with a much smaller dependence on its size. The edges in the graphene nano-flakes exhibit several different meta-stable configurations (isomers) during heating before melting occurs.

preprint2013arXiv

Melting of Partially Fluorinated Graphene: From Detachment of Fluorine Atoms to Large Defects and Random Coils

The melting of fluorographene is very unusual and depends strongly on the degree of fluorination. For temperatures below 1000 K, fully fluorinated graphene (FFG) is thermo-mechanically more stable than graphene but at T$_m\approx$2800 K FFG transits to random coils which is almost twice lower than the melting temperature of graphene, i.e. 5300 K. For fluorinated graphene (PFG) up to 30 % ripples causes detachment of individual F-atoms around 2000 K while for 40-60 % fluorination, large defects are formed beyond 1500 K and beyond 60% of fluorination F-atoms remain bonded to graphene until melting. The results agree with recent experiments on the dependence of the reversibility of the fluorination process on the percentage of fluorination.

preprint2013arXiv

Multiband tunneling in trilayer graphene

The electronic tunneling properties of the two stable forms of trilayer graphene (TLG), rhombohedral ABC and Bernal ABA, are examined for pn and pnp junctions as realized by using a single gate (SG) or a double gate (DG). For the rhombohedral form, due to the chirality of the electrons, the Klein paradox is found at normal incidence for SG devices while at high energy interband scattering between additional propagation modes can occur. The electrons in Bernal ABA TLG can have a monolayer- or bilayer-like character when incident on a SG device. Using a DG however both propagation modes will couple by breaking the mirror symmetry of the system which induces intermode scattering and resonances that depend on the width of the DG pnp junction. For ABC TLG the DG opens up a band gap which suppresses Klein tunneling. The DG induces also an unexpected asymmetry in the tunneling angle for single valley electrons.

preprint2013arXiv

Phonon Softening and Direct to Indirect Bandgap Crossover in Strained Single Layer MoSe2

Motivated by recent experimental observations of Tongay et al. [Tongay et al., Nano Letters, 12(11), 5576 (2012)] we show how the electronic properties and Raman characteristics of single layer MoSe2 are affected by elastic biaxial strain. We found that with increasing strain: (1) the E' and E" Raman peaks (E1g and E2g in bulk) exhibit significant red shifts (up to 30 cm-1), (2) the position of the A1' peak remains at 180 cm-1 (A1g in bulk) and does not change considerably with further strain, (3) the dispersion of low energy flexural phonons crosses over from quadratic to linear and (4) the electronic band structure undergoes a direct to indirect bandgap crossover under 3% biaxial tensile strain. Thus the application of strain appears to be a promising approach for a rapid and reversible tuning of the electronic, vibrational and optical properties of single layer MoSe2 and similar MX2 dichalcogenides.

preprint2013arXiv

Pseudo magnetic field in strained graphene: revisited

We revisit the theory of the pseudo magnetic field as induced by strain in graphene using the tight-binding approach. A systematic expansion of the hopping parameter and the deformation of the lattice vectors is presented from which we obtain an expression for the pseudo magnetic field for low energy electrons. We generalize and discuss previous results and propose a novel effective Hamiltonian. The contributions of the different terms to the pseudo magnetic field expression is investigated for a model triaxial strain profile and are compared with the full solution. Our work suggests that the previous proposed pseudo magnetic field expression is valid up to reasonably high strain 15 % and there is no K-dependent pseudo-magnetic field.

preprint2013arXiv

Realization of Artificial Ice Systems for Magnetic Vortices in a Superconducting MoGe Thin-film with Patterned Nanostructures

We report an anomalous matching effect in MoGe thin films containing pairs of circular holes arranged in such a way that four of those pairs meet at each vertex point of a square lattice. A remarkably pronounced fractional matching was observed in the magnetic field dependences of both the resistance and the critical current. At the half matching field the critical current can be even higher than that at zero field. This has never been observed before for vortices in superconductors with pinning arrays. Numerical simulations within the nonlinear Ginzburg-Landau theory reveal a square vortex ice configuration in the ground state at the half matching field and demonstrate similar characteristic features in the field dependence of the critical current, confirming the experimental realization of an artificial ice system for vortices for the first time.

preprint2013arXiv

Single-file and normal diffusion of magnetic colloids in modulated channels

Diffusive properties of interacting magnetic dipoles confined in a parabolic narrow channel and in the presence of a periodic modulated (corrugated) potential along the unconfined direction are studied using Brownian dynamics simulations. We compare our simulation results with the analytical result for the effective diffusion coefficient of a single-particle by Festa and d'Agliano [Physica A \textbf{90}, 229 (1978)] and show the importance of inter-particle interaction on the diffusion process. We present results for the diffusion of magnetic dipoles as a function of linear density, strength of the periodic modulation and commensurability factor.

preprint2013arXiv

Snake states in graphene quantum dots in the presence of a p-n junction

We investigate the magnetic interface states of graphene quantum dots that contain p-n junctions. Within a tight-binding approach, we consider rectangular quantum dots in the presence of a perpendicular magnetic field containing p-n, as well as p-n-p and n-p-n junctions. The results show the interplay between the edge states associated with the zigzag terminations of the sample and the snake states that arise at the p-n junction, due to the overlap between electron and hole states at the potential interface. Remarkable localized states are found at the crossing of the p-n junction with the zigzag edge having a dumb-bell shaped electron distribution. The results are presented as function of the junction parameters and the applied magnetic flux.

preprint2013arXiv

Soft vortex matter in a type-I/type-II superconducting bilayer

Magnetic flux patterns are known to strongly differ in the intermediate state of type-I and type-II superconductors. Using a type-I/type-II bilayer we demonstrate hybridization of these flux phases into a plethora of unique new ones. Owing to a complicated multi-body interaction between individual fluxoids, many different intriguing patterns are possible under applied magnetic field, such as few-vortex clusters, vortex chains, mazes or labyrinthal structures resembling the phenomena readily encountered in soft matter physics. However, in our system the patterns are tunable by sample parameters, magnetic field, current and temperature, which reveals transitions from short-range clustering to long-range ordered phases such as parallel chains, gels, glasses and crystalline vortex lattices, or phases where lamellar type-I flux domains in one layer serve as a bedding potential for type-II vortices in the other - configurations clearly beyond the soft-matter analogy.

preprint2013arXiv

Spectroscopy of snake states using a graphene Hall bar

An approach to observe snake states in a graphene Hall bar containing a pn-junction is proposed. The magnetic field dependence of the bend resistance in a ballistic graphene Hall bar structure containing a tilted pn-junction oscillates as a function of applied magnetic field. We show that each oscillation is due to a specific snake state that moves along the pn-interface. Furthermore depending on the value of the magnetic field and applied potential we can control the lead in which the electrons will end up and hence control the response of the system.

preprint2013arXiv

Spin-orbit interaction induced singularity of the charge density relaxation propagator

The charge density relaxation propagator of a two dimensional electron system, which is the slope of the imaginary part of the polarization function, exhibits singularities for bosonic momenta having the order of the spin-orbit momentum and depending on the momentum orientation. We have provided an intuitive understanding for this non-analytic behavior in terms of the inter chirality subband electronic transitions, induced by the combined action of Bychkov-Rashba (BR) and Dresselhaus (D) spin-orbit coupling. It is shown that the regular behavior of the relaxation propagator is recovered in the presence of only one BR or D spin-orbit field or for spin-orbit interaction with equal BR and D coupling strengths. This creates a new possibility to influence carrier relaxation properties by means of an applied electric field.

preprint2013arXiv

Spiral graphone and one sided fluorographene nano-ribbons

The instability of a free-standing one sided hydrogenated/fluorinated graphene nano-ribbon, i.e. graphone/fluorographene, is studied using ab-initio, semiempirical and large scale molecular dynamics simulations. Free standing semi-infinite arm-chair like hydrogenated/fluorinated graphene (AC-GO/AC-GF) and boat like hydrogenated/fluorinated graphene (B-GO/B-GF) (nano-ribbons which are periodic along the zig-zag direction) are unstable and spontaneously transform into spiral structures. We find that rolled, spiral B-GO and B-GF are energetically more favorable than spiral AC-GO and AC-GF which is opposite to the double sided flat hydrogenated/fluorinated graphene, i.e. graphane/fluorographene. We found that the packed, spiral structures exhibit unexpected localized HOMO-LUMO at the edges with increasing energy gap during rolling. These rolled hydrocarbon structures are stable beyond room temperature up to at least $T$=1000\,K.

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

The 30-band k.p modeling of electron and hole states in silicon quantum wells

We modeled the electron and hole states in Si/SiO2 quantum wells within a basis of standing waves using the 30-band k.p theory. The hard-wall confinement potential is assumed, and the influence of the peculiar band structure of bulk silicon on the quantum-well sub-bands is explored. Numerous spurious solutions in the conduction-band and valence-band energy spectra are found and are identified to be of two types: (1) spurious states which have large contributions of the bulk solutions with large wave-vectors (the high-k spurious solutions) and (2) states which originate mainly from the spurious valley outside the Brillouin zone (the extra-valley spurious solutions). An algorithm to remove all those nonphysical solutions from the electron and hole energy spectra is proposed. Furthermore, slow and oscillatory convergence of the hole energy levels with the number of basis functions is found and is explained by the peculiar band mixing and the confinement in the considered quantum well. We discovered that assuming the hard-wall potential leads to numerical instability of the hole states computation. Nonetheless, allowing the envelope functions to exponentially decay in a barrier of finite height is found to improve the accuracy of the computed hole states.

preprint2013arXiv

Thermal properties of fluorinated graphene

Large scale atomistic simulations using the reactive force field approach (ReaxFF) are implemented to investigate the thermomechanical properties of fluorinated graphene (FG). A new set of parameters for the reactive force field potential (ReaxFF) optimized to reproduce key quantum mechanical properties of relevant carbon-fluor cluster systems are presented. Molecular dynamics (MD) simulations are used to investigate the thermal rippling behavior of FG and its mechanical properties and compare them with graphene (GE), graphane (GA) and a sheet of BN. The mean square value of the height fluctuations $< h^2>$ and the height-height correlation function $H(q)$ for different system sizes and temperatures show that FG is an un-rippled system in contrast to the thermal rippling behavior of graphene (GE). The effective Young's modulus of a flake of fluorinated graphene is obtained to be 273 N/m and 250 N/m for a flake of FG under uniaxial strain along arm-chair and zig-zag direction, respectively.

preprint2013arXiv

Thermomechanical properties of a single hexagonal boron nitride sheet

Using atomistic simulations we investigate the thermodynamical properties of a single atomic layer of hexagonal boron nitride (h-BN). The thermal induced ripples, heat capacity, and thermal lattice expansion of large scale h-BN sheets are determined and compared to those found for graphene (GE) for temperatures up to 1000 K. By analyzing the mean square height fluctuations $< h^2>$ and the height-height correlation function $H(q)$ we found that the h-BN sheet is a less stiff material as compared to graphene. The bending rigidity of h-BN: i) is about 16% smaller than the one of GE at room temperature (300 K), and ii) increases with temperature as in GE. The difference in stiffness between h-BN and GE results in unequal responses to external uniaxial and shear stress and different buckling transitions. In contrast to a GE sheet, the buckling transition of a h-BN sheet depends strongly on the direction of the applied compression. The molar heat capacity, thermal expansion coefficient and the Gruneisen parameter are estimated to be 25.2 J\,mol$^{-1}$\,K$^{-1}$, 7.2$\times10^{-6}$K$^{-1}$ and 0.89, respectively.

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

Tunable diffusion of magnetic particles in a quasi-one-dimensional channel

The diffusion of a system of ferromagnetic dipoles confined in a quasi-one-dimensional parabolic trap is studied using Brownian dynamics simulations. We show that the dynamics of the system is tunable by an in-plane external homogeneous magnetic field. For a strong applied magnetic field, we find that the mobility of the system, the exponent of diffusion and the crossover time among different diffusion regimes can be tuned by the orientation of the magnetic field. For weak magnetic fields, the exponent of diffusion in the subdiffusive regime is independent of the orientation of the external field.

preprint2013arXiv

Vibrational properties of nanographene

The eigenmodes and the vibrational density of states of the ground state configuration of graphene clusters are calculated using atomistic simulations. The modified Brenner potential is used to describe the carbon-carbon interaction and carbon-hydrogen interaction in case of H-passivated edges. For a given configuration of the C-atoms the eigenvectors and eigenfrequencies of the normal modes are obtained after diagonalisation of the dynamical matrix whose elements are the second derivative of the potential energy. The compressional and shear properties are obtained from the divergence and rotation of the velocity field. For symmetric and defective clusters with pentagon arrangement on the edge, the highest frequency modes are shear modes. The specific heat of the clusters is also calculated within the harmonic approximation and the convergence to the result for bulk graphene is investigated.

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

Atypical BCS-BEC crossover induced by quantum-size effects

Quantum-size oscillations of the basic physical characteristics of a confined fermionic condensate are a well-known phenomenon. Its conventional understanding is based on the single-particle physics, whereby the oscillations follow the size-dependent changes in the single-particle density of states. Here we present a study of a cigar-shaped ultracold superfluid Fermi gas, which demonstrates an important many-body aspect of the quantum-size effects, overlooked previously. The many-body physics is revealed in the atypical crossover from the Bardeen-Cooper-Schrieffer (BCS) superfluid to the Bose-Einstein condensate (BEC) induced by the size quantization of the particle motion. Quantized perpendicular spectrum results in the formation of single-particle subbands (shells) so that the aggregate fermionic condensate becomes a coherent mixture of subband condensates. Each time when the lower edge of a subband crosses the chemical potential, the BCS-BEC crossover is approached in this subband, and the aggregate condensate contains both the BCS and BEC-like components.

preprint2012arXiv

Controllable spin singlet - spin triplet transition in three concentric quantum rings through magnetic field and confinement potential

We present a theoretical study of the spectrum of electrons confined in triple concentric rings. An unusual ordering and rich variety of angular momentum transitions are found that depend on the coupling between the rings and the confinement potential of the rings. Using the Configuration Interaction (CI) method, we calculated the two electron energy spectrum. Spin singlet to spin triplet transitions of the electron ground state are predicted and a fractional Aharonov-Bohm effect is found. We show that both the period and amplitude of the spin singlet - triplet energy gap depend strongly on the confinement potential and the external magnetic field. The spin singlet - triplet transition is found to depend on the spin Zeeman energy, especially for rings with weak confinement and in the presence of large magnetic field. The amplitude of the spin singlet - triplet energy gap depends on the Landé $g$-factor but the period of the transitions is independent of $g$.

preprint2012arXiv

Effect of grain boundary on the buckling of graphene nanoribbons

The buckling of graphene nano-ribbons containing a grain boundary is studied using atomistic simulations where free and supported boundary conditions are invoked. We found that when graphene contains a small angle grain boundary the buckling strains are larger when the ribbons with free (supported) boundary condition are subjected to compressive tension parallel (perpendicular) to the grain boundary. The shape of the deformations of the buckled graphene nanoribbons depends on the boundary conditions and the presence of the grain boundary and the direction of applied in-plane compressive tension. Large angle grain boundary results in smaller buckling strains as compared to perfect graphene or to a small angle grain boundary.

preprint2012arXiv

Electron-phonon bound state in graphene

The fine structure of the Dirac energy spectrum in graphene induced by electron-optical phonon coupling is investigated in the portion of the spectrum near the phonon emission threshold. The derived new dispersion equation in the immediate neighborhood below the phonon threshold corresponds to an electron-phonon bound state. We find that the singular vertex corrections beyond perturbation theory increase strongly the electron-phonon binding energy. The predicted enhancement of the effective electron-phonon coupling can be measured using angle-resolved spectroscopy.

preprint2012arXiv

Extended Ginzburg-Landau formalism: systematic expansion in small deviation from the critical temperature

Based on the Gor'kov formalism for a clean s-wave superconductor, we develop an extended version of the single-band Ginzburg-Landau (GL) theory by means of a systematic expansion in the deviation from the critical temperature T_c, i.e., tau=1-T/T_c. We calculate different contributions to the order parameter and the magnetic field: the leading contributions (~ tau^1/2 in the order parameter and ~ tau in the magnetic field) are controlled by the standard Ginzburg-Landau (GL) theory, while the next-to-leading terms (~ tau^3/2 in the gap and ~ tau^2 in the magnetic field) constitute the extended GL (EGL) approach. We derive the free-energy functional for the extended formalism and the corresponding expression for the current density. To illustrate the usefulness of our formalism, we calculate, in a semi-analytical form, the temperature-dependent correction to the GL parameter at which the surface energy becomes zero, and analytically, the temperature dependence of the thermodynamic critical field. We demonstrate that the EGL formalism is not just a mathematical extension to the theory - variations of both the gap and the thermodynamic critical field with temperature calculated within the EGL theory are found in very good agreement with the full BCS results down to low temperatures, which dramatically improves the applicability of the formalism compared to its standard predecessor.

preprint2012arXiv

Ginzburg-Landau theory of the zig-zag transition in quasi-one-dimensional classical Wigner crystals

We present a mean-field description of the zig-zag phase transition of a quasi-one-dimensional system of strongly interacting particles, with interaction potential $r^{-n}e^{-r/λ}$, that are confined by a power-law potential ($y^α$). The parameters of the resulting one-dimensional Ginzburg-Landau theory are determined analytically for different values of $α$ and $n$. Close to the transition point for the zig-zag phase transition, the scaling behavior of the order parameter is determined. For $α=2$ the zig-zag transition from a single to a double chain is of second order, while for $α>2$ the one chain configuration is always unstable and for $α<2$ the one chain ordered state becomes unstable at a certain critical density resulting in jumps of single particles out of the chain.

preprint2012arXiv

Induced polarization and electronic properties of carbon doped boron-nitride nanoribbons

The electronic properties of boron-nitride nanoribbons (BNNRs) doped with a line of carbon atoms are investigated by using density functional calculations. Three different configurations are possible: the carbon atoms may replace a line of boron or nitrogen atoms or a line of alternating B and N atoms which results in very different electronic properties. We found that: i) the NCB arrangement is strongly polarized with a large dipole moment having an unexpected direction, ii) the BCB and NCN arrangement are non-polar with zero dipole moment, iii) the doping by a carbon line reduces the band gap independent of the local arrangement of boron and nitrogen around the carbon line, iv) an electric field parallel to the carbon line polarizes the BN sheet and is found to be sensitive to the presence of carbon dopants, and v) the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital decreases linearly with increasing applied electric field directed parallel to the carbon line. We show that the polarization and energy gap of carbon doped BNNRs can be tuned by an electric field applied parallel along the carbon line.

preprint2012arXiv

Large Magnetoresistance Oscillations in Mesoscopic Superconductors Due to Current-Excited Moving Vortices

We show in the case of a superconducting Nb ladder that a mesoscopic superconductor typically exhibits magnetoresistance oscillations whose amplitude and temperature dependence are different from those stemming from the Little-Parks effect. We demonstrate that these large resistance oscillations (as well as the monotonic background on which they are superimposed) are due to {\it current-excited moving vortices}, where the applied current in competition with the oscillating Meissner currents imposes/removes the barriers for vortex motion in increasing magnetic field. Due to the ever present current in transport measurements, this effect should be considered in parallel with the Little-Parks effect in low-$T_c$ samples, as well as with recently proposed thermal activation of dissipative vortex-antivortex pairs in high-$T_c$ samples.

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

Piezoelectric surface acoustical phonon limited mobility of electrons in graphene on a GaAs substrate

We study the mobility of Dirac fermions in monolayer graphene on a GaAs substrate, restricted by the combined action of the extrinsic potential of piezoelectric surface acoustical phonons of GaAs (PA) and of the intrinsic deformation potential of acoustical eigen-phonons in graphene (DA). In the high temperature ($T$) regime the momentum relaxation rate exhibits the same linear dependence on $T$ but different dependences on the carrier density $n$, corresponding to the mobility $μ\propto 1/\sqrt{n}$ and $1/n$, respectively for the PA and DA scattering mechanisms. In the low $T$ Bloch-Grüneisen regime, the mobility shows the same square-root density dependence, $μ\propto \sqrt{n}$, but different temperature dependences, $μ\propto T^{-3}$ and $ T^{-4}$, respectively for PA and DA phonon scattering.

preprint2012arXiv

Plasmon and dielectric background inhomogeneity enhancement of Coulomb drag in graphene double-layer structures

The drag of massless fermions in graphene double-layer structures is investigated in a wide rage of temperatures and inter-layer separations. We show that the inhomogeneity of the dielectric background in such graphene structures for experimentally relevant parameters results in a significant enhancement of the drag resistivity. At intermediate temperatures the dynamical screening via plasmon-mediated drag enhances the drag resistivity and results in an upturn in its behavior at large inter-layer separations. In a range of inter-layer separations, corresponding to the strong-to-weak crossover coupling of graphene layers, we find that the drag resistivity decreases approximately quadratically with the inter-layer spacing. This dependence weakens with a decrease of the inter-layer spacing while for larger separations we recover the cubic (quartic) dependence at intermediate (low) temperatures.

preprint2012arXiv

Plasmonic excitations in Coulomb coupled N-layer graphene structures

We study Dirac plasmons and their damping in spatially separated $N$-layer graphene structures at finite doping and temperatures. The plasmon spectrum consists of one optical excitation with a square-root dispersion and $N-1$ acoustical excitations with linear dispersions, which are undamped at zero temperature within a triangular energy region outside the electron-hole continuum. For any finite number of graphene layers we have found that the energy and weight of the optical plasmon increase in the long wavelength limit, respectively, as square-root and linear functions of $N$. This is in agreement with recent experimental findings. With an increase of the number of multilayer acoustical plasmon modes, the energy and weight of the upper lying branches also exhibit an enhancement with $N$. This increase is strongest for the uppermost acoustical mode so that its energy can exceed at some value of momentum the plasmon energy in an individual graphene sheet. Meanwhile, the energy of the low lying acoustical branches decreases weakly with $N$ as compared with the single acoustical mode in double-layer graphene structures. Our numerical calculations provide a detailed understanding of the overall behavior of the wave vector dependence of the optical and acoustical multilayer plasmon modes and show how their dispersion and damping are modified as a function of temperature, interlayer spacing, and inlayer carrier density in (un)balanced graphene multilayer structures.

preprint2012arXiv

Predicted field-dependent increase of critical currents in asymmetric superconducting nanocircuits

The critical current of a thin superconducting strip of width $W$ much larger than the Ginzburg-Landau coherence length $ξ$ but much smaller than the Pearl length $Λ= 2 λ^2/d$ is maximized when the strip is straight with defect-free edges. When a perpendicular magnetic field is applied to a long straight strip, the critical current initially decreases linearly with $H$ but then decreases more slowly with $H$ when vortices or antivortices are forced into the strip. However, in a superconducting strip containing sharp 90-degree or 180-degree turns, the zero-field critical current at H=0 is reduced because vortices or antivortices are preferentially nucleated at the inner corners of the turns, where current crowding occurs. Using both analytic London-model calculations and time-dependent Ginzburg-Landau simulations, we predict that in such asymmetric strips the resulting critical current can be {\it increased} by applying a perpendicular magnetic field that induces a current-density contribution opposing the applied current density at the inner corners. This effect should apply to all turns that bend in the same direction.

preprint2012arXiv

Resonant harmonic generation and collective spin rotations in electrically driven quantum dots

Spin rotations induced by an AC electric field in a two-electron double quantum dot are studied by an exact numerical solution of the time dependent Schroedinger equation in the context of recent electric dipole spin resonance experiments based on the Pauli blockade. We demonstrate that the splitting of the main resonance line by the spin exchange coupling is accompanied by the appearance of fractional resonances and that both these effects are triggered by interdot tunnel coupling. We find that the AC driven system generates residual but distinct harmonics of the driving frequency which are amplified when tuned to the main transition frequency. The mechanism is universal for electron systems in electrically driven potentials and works also in the absence of electron-electron interaction or spin-orbit coupling.

preprint2012arXiv

Scattering of a Dirac electron on a mass barrier

The interaction of a wave packet (and in particular the wave front) with a mass barrier is investigated in one dimension. We discuss the main features of the wave packet that are inherent to two-dimensional wave packets, such as compression during reflection, penetration in the case when the energy is lower than the height of the barrier, waving tails, precursors, and the retardation of the reflected and penetrated wave packets. These features depend on the wave-packet envelope function which we demonstrate by considering the case of a rectangular wave packet with sharp front and trailing edges and a smooth Gaussian wave packet. The method of Fourier integral for obtaining the nonstationary solutions is used.

preprint2012arXiv

Spatially-dependent sensitivity of superconducting meanders as single-photon detectors

The photo-response of a thin current-carrying superconducting stripe with a 90-degree turn is studied within the time-dependent Ginzburg-Landau theory. We show that the photon acting near the inner corner (where the current density is maximal due to the current crowding [J. R. Clem and K. K. Berggren, Phys. Rev. B {\bf 84}, 174510 (2011)]) triggers the nucleation of superconducting vortices at currents much smaller than the expected critical one, but {\it does not} bring the system to a higher resistive state and thus remains undetected. The transition to the resistive state occurs only when the photon hits the stripe away from the corner due to there uniform current distribution across the sample, and dissipation is due to the nucleation of a kinematic vortex-antivortex pair near the photon incidence. We propose strategies to account for this problem in the measurements.

preprint2012arXiv

Spin Orbit Mediated Manipulation of Heavy Hole Spin Qubit in Gated Semiconductor Nanodevices

A novel spintronic nanodevice is proposed that is capable to manipulate the single heavy hole spin state in a coherent manner. It can act as a single quantum logic gate. The heavy hole spin transformations are realized by transporting the hole around closed loops defined by metal gates deposited on top of the nanodevice. The device exploits Dresselhaus spin orbit interaction which translates the spatial motion of the hole into a rotation of the spin. The proposed quantum gate operates on sub nanosecond time scales and requires only the application of a weak static voltage which allows for addressing heavy hole spin qubit individually. Our results are supported by quantum mechanical time dependent calculations within the four band Luttinger Kohn model.

preprint2012arXiv

Strain engineered graphene using a nanostructured substrate: I Deformations

Using atomistic simulations we investigate the morphological properties of graphene deposited on top of a nanostructured substrate. Sinusoidally corrugated surfaces, steps, elongated trenches, one dimensional and cubic barriers, spherical bubbles, Gaussian bump and Gaussian depression are considered as support structures for graphene. The graphene-substrate interaction is governed by van der Waals forces and the profile of the graphene layer is determined by minimizing the energy using molecular dynamics simulations. Based on the obtained optimum configurations, we found that: (i) for graphene placed over sinusoidally corrugated substrates with corrugation wave lengths longer than 2\,nm, the graphene sheet follows the substrate pattern while for supported graphene it is always suspended across the peaks of the substrate, (ii) the conformation of graphene to the substrate topography is enhanced when increasing the energy parameter in the van der Waals model, (iii) the adhesion of graphene into the trenches depends on the width of the trench and on graphene's orientation, i.e. in contrast to a small width (3 nm) nanoribbon with armchair edges, the one with zig-zag edges follows the substrate profile, (iv) atomic scale graphene follows a Gaussian bump substrate but not the substrate with a Gaussian depression, and (v) the adhesion energy due to van der Waals interaction varies in the range [0.1-0.4] J/m^2.

preprint2012arXiv

Strain engineered graphene using a nanostructured substrate: II Pseudo-magnetic fields

The strain induced pseudo-magnetic field in supported graphene deposited on top of a nanostructured substrate is investigated by using atomistic simulations. Step, elongated trench, one dimensional barrier, spherical bubbles, Gaussian bump and Gaussian depression are considered as support structures for graphene. From the obtained optimum configurations we found very strong induced pseudo-magnetic fields which can reach up to $\sim$ 1000\,T due to the strain-induced deformations in the supported graphene. Different magnetic confinements with controllable geometries are found by tuning the pattern of the substrate. The resulting induced magnetic fields for graphene on top of a step, barrier and trench are calculated. In contrast to the step and trench the middle part of graphene on top of a barrier has zero pseudo-magnetic field. This study provides a theoretical background for designing magnetic structures in graphene by nanostructuring substrates. We found that altering the radial symmetry of the deformation, changes the six-fold symmetry of the induced pseudo-magnetic field.

preprint2012arXiv

Structural transitions in vertically and horizontally coupled parabolic channels of Wigner crystals

Structural phase transitions in two vertically or horizontally coupled channels of strongly interacting particles are investigated. The particles are free to move in the $x$-direction but are confined by a parabolic potential in the $y$-direction. They interact with each other through a screened power-law potential ($r^{-n}e^{-r/λ}$). In vertically coupled systems the channels are stacked above each other in the direction perpendicular to the $(x,y)$-plane, while in horizontally coupled systems both channels are aligned in the confinement direction. Using Monte Carlo (MC) simulations we obtain the ground state configurations and the structural transitions as a function of the linear particle density and the separation between the channels. At zero temperature the vertically coupled system exhibits a rich phase diagram with continuous and discontinuous transitions. On the other hand the vertically coupled system exhibits only a very limited number of phase transitions due to its symmetry. Further we calculated the normal modes for the Wigner crystals in both cases. From MC simulations we found that in the case of vertically coupled systems the zigzag transition is only possible for low densities. A Ginzburg-Landau theory for the zigzag transition is presented, which predicts correctly the behavior of this transition from which we interpret the structural phase transition of the Wigner crystal through the reduction of the Brillouin zone.

preprint2012arXiv

Superconducting transition temperature of Pb nanofilms: Impact of the thickness-dependent oscillations of the phonon mediated electron-electron coupling

To date, several experimental groups reported measurements of the thickness dependence of T_c of atomically uniform single-crystalline Pb nanofilms. The reported amplitude of the T_c-oscillations varies significantly from one experiment to another. Here we propose that the reason for this unresolved issue is an interplay of the quantum-size variations in the single-electron density of states with thickness-dependent oscillations in the phonon mediated electron-electron coupling. Such oscillations in the coupling depend on the substrate material, the quality of the interface, the protection cover and other details of the fabrication process, changing from one experiment to another. This explains why the available data do not exhibit one-voice consistency about the amplitude of the T_c-oscillations. Our analyses are based on a numerical solution of the Bogoliubov-de Gennes equations for a superconducting slab.

preprint2012arXiv

Thermal rippling behavior of graphane

Thermal fluctuations of single layer hydrogenated graphene (graphane) are investigated using large scale atomistic simulations. By analyzing the mean square value of the height fluctuations $<h^2>$ and the height-height correlation function $H(q)$ for different system sizes and temperatures we show that hydrogenated graphene is an un-rippled system in contrast to graphene. The height fluctuations are bounded, which is confirmed by a $ H(q) $ tending to a constant in the long wavelength limit instead of showing the characteristic scaling law $ q^{4-η} (η\simeq 0.85)$ predicted by membrane theory. This unexpected behaviour persists up to temperatures of at least 900 K and is a consequence of the fact that in graphane the thermal energy can be accommodated by in-plane bending modes, i.e. modes involving C-C-C bond angles in the buckled carbon layer, instead of leading to significant out-of-plane fluctuations that occur in graphene.

preprint2012arXiv

Thermomechanical properties of graphene: valence force field model approach

Using the valence force field model of Perebeinos and Tersoff [Phys. Rev. B {\bf79}, 241409(R) (2009)], different energy modes of suspended graphene subjected to tensile or compressive strain are studied. By carrying out Monte Carlo simulations it is found that: i) only for small strains ($|\varepsilon| \lessapprox 0.02$) the total energy is symmetrical in the strain, while it behaves completely different beyond this threshold; ii) the important energy contributions in stretching experiments are stretching, angle bending, out-of-plane term and a term that provides repulsion against $π-π$ misalignment; iii) in compressing experiments the two latter terms increase rapidly and beyond the buckling transition stretching and bending energies are found to be constant; iv) from stretching-compressing simulations we calculated the Young modulus at room temperature 350$\pm3.15$\,N/m, which is in good agreement with experimental results (340$\pm50$\,N/m) and with ab-initio results [322-353]\,N/m; v) molar heat capacity is estimated to be 24.64\,J/mol$^{-1}$K$^{-1}$ which is comparable with the Dulong-Petit value, i.e. 24.94\,J/mol$^{-1}$K$^{-1}$ and is almost independent of the strain; vi) non-linear scaling properties are obtained from height-height correlations at finite temperature; vii) the used valence force field model results in a temperature independent bending modulus for graphene, and viii) the Gruneisen parameter is estimated to be 0.64.

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

Transition from single-file to two-dimensional diffusion of interacting particles in a quasi-one-dimensional channel

Diffusive properties of a monodisperse system of interacting particles confined to a \textit{quasi}-one-dimensional (Q1D) channel are studied using molecular dynamics (MD) simulations. We calculate numerically the mean-squared displacement (MSD) and investigate the influence of the width of the channel (or the strength of the confinement potential) on diffusion in finite-size channels of different shapes (i.e., straight and circular). The transition from single-file diffusion (SFD) to the two-dimensional diffusion regime is investigated. This transition (regarding the calculation of the scaling exponent ($α$) of the MSD $<Δx^{2}(t)>$ $\propto t^α$) as a function of the width of the channel, is shown to change depending on the channel's confinement profile. In particular the transition can be either smooth (i.e., for a parabolic confinement potential) or rather sharp/stepwise (i.e., for a hard-wall potential), as distinct from infinite channels where this transition is abrupt. This result can be explained by qualitatively different distributions of the particle density for the different confinement potentials.

preprint2012arXiv

Tunable optical Aharonov-Bohm effect in a semiconductor quantum ring

By applying an electric field perpendicular to a semiconductor quantum ring we show that it is possible to modify the single particle wave function between quantum dot (QD)-like to ring-like. The constraints on the geometrical parameters of the quantum ring to realize such a transition are derived. With such a perpendicular electric field we are able to tune the Aharanov-Bohm (AB) effect for both single particles and for excitons. The tunability is in both the strength of the AB-effect as well as in its periodicity. We also investigate the strain induce potential inside the self assembled quantum ring and the effect of the strain on the AB effect.

preprint2012arXiv

Two-band superconductors: Extended Ginzburg-Landau formalism by a systematic expansion in small deviation from the critical temperature

We derive the extended Ginzburg-Landau (GL) formalism for a clean s-wave two-band superconductor by employing a systematic expansion of the free-energy functional and the corresponding matrix gap equation in powers of the small deviation from the critical temperature tau = 1-T/T_c. The two lowest orders of this expansion produce the equation for T_c and the GL theory. It is shown that in agreement with previous studies, the two-band GL theory maps onto the single-band GL model and thus fails to describe the difference in the spatial profiles of the two band condensates. We prove that except for some very special cases, this difference appears already in the leading correction to the GL theory, which constitutes the extended GL formalism. We derive linear differential equations that determine the leading corrections to the band order parameters and magnetic field, discuss the validity of these equations, and consider examples of an important interplay between the band condensates. Finally, we present numerical results for the thermodynamic critical magnetic field and temperature-dependent band gaps (at zero field), which are in a very good agreement with those obtained from the full BCS approach in a wide temperature range. To this end, we emphasize the advantages of our extended GL theory in comparison with the often used two-component GL-like model based on an unreconstructed two-band generalization of the Gor'kov derivation.

preprint2012arXiv

Two-band superconductors: Hidden criticality deep in the superconducting state

We show that two-band superconductors harbor hidden criticality deep in the superconducting state, stemming from the critical temperature of the weaker band taken as an independent system. For sufficiently small interband coupling $γ$ the coherence length of the weaker band exhibits a remarkable deviation from the conventional monotonic increase with temperature, namely, a pronounced peak close to the hidden critical point. The magnitude of the peak scales proportionally to γ^(-μ), with the Landau critical exponent μ= 1/3, the same as found for the mean-field critical behavior with respect to the source field in ferromagnets and ferroelectrics. Here reported hidden criticality of multi-band superconductors can be experimentally observed by, e.g., imaging of the variations of the vortex core in a broader temperature range. Similar effects are expected for the superconducting multilayers.

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

Using self-driven microswimmers for particle separation

Microscopic self-propelled swimmers capable of autonomous navigation through complex environments provide appealing opportunities for localization, pick-up and delivery of micro-and nanoscopic objects. Inspired by motile cells and bacteria, man-made microswimmers have been fabricated, and their motion in patterned surroundings has been experimentally studied. We propose to use self-driven artificial microswimmers for separation of binary mixtures of colloids. We revealed different regimes of separation including one with a velocity inversion. Our finding could be of use for various biological and medical applications.

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

BCS-BEC crossover-like phenomena driven by quantum-size effects in quasi-one-dimensional fermionic condensates

Quantum confinement is known to influence fermionic condensates, resulting in quantum-size oscillations of superfluid/superconducting properties. Here we show that the impact of quantum-size effects is even more dramatic. Under realistic conditions, a significant phase-space reconfiguration induced by quantum-size effects opens a quasi-molecule channel in the fermionic pairing so that the condensed pairs exhibit features typical of a molecular state. As an illustration we consider a quasi-one-dimensional fermionic condensate, as realized, e.g., in cigar-shaped atomic Fermi gases or superconducting quantum wires. In this case the transverse quantization of the particle motion favors pairing through a coherent superposition of quantum channels that are formed due to the grouping of single-particle levels into a series of well distinguished subbands. Whenever the bottom of a subband approaches the Fermi level, the longitudinal spatial distribution of fermions in a condensed pair becomes strongly localized within the corresponding quantum channel. The fermionic pairs in this channel resemble molecules with bosonic character.

preprint2011arXiv

Bilayer graphene with single and multiple electrostatic barriers: band structure and transmission

We evaluate the electronic transmission and conductance in bilayer graphene through a finite number of potential barriers. Further, we evaluate the dispersion relation in a bilayer graphene superlattice with a periodic potential applied to both layers. As a model we use the tight-binding Hamiltonian in the continuum approxi- mation. For zero bias the dispersion relation shows a finite gap for carriers with zero momentum in the direction parallel to the barriers. This is in contrast to single-layer graphene where no such gap was found. A gap also appears for a finite bias. Numerical results for the energy spectrum, conductance, and the density of states are presented and contrasted with those pertaining to single-layer graphene.

preprint2011arXiv

Chiral states in bilayer graphene: magnetic field dependence and gap opening

At the interface of electrostatic potential kink profiles one dimensional chiral states are found in bilayer graphene (BLG). Such structures can be created by applying an asymmetric potential to the upper and the lower layer of BLG. We found that: i) due to the strong confinement by the single kink profile the uni-directional states are only weakly affected by a magnetic field, ii) increasing the smoothness of the kink potential results in additional bound states which are topologically different from those chiral states, and iii) in the presence of a kink-antikink potential the overlap between the oppositely moving chiral states results in the appearance of crossing and anti-crossing points in the energy spectrum. This leads to the opening of tunable minigaps in the spectrum of the uni-directional topological states.

preprint2011arXiv

Conditions for non-monotonic vortex interaction in two-band superconductors

We describe a semi-analytic approach to the two-band Ginzburg-Landau theory, which predicts the behavior of vortices in two-band superconductors. We show that the character of the short-range vortex-vortex interaction is determined by the sign of the normal domain - superconductor interface energy, in analogy with the conventional differentiation between type-I and type-II superconductors. However, we also show that the long-range interaction is determined by a modified Ginzburg-Landau parameter $κ^*$, different from the standard $κ$ of a bulk superconductor. This opens the possibility for non-monotonic vortex-vortex interaction, which is temperature-dependent, and can be further tuned by alterations of the material on the microscopic scale.

preprint2011arXiv

Different length-scales for order parameters in two-gap superconductors: the extended Ginzburg-Landau theory

Using the Ginzburg-Landau theory extended to the next-to-leading order we determine numerically the healing lengths of the two order parameters at the two-gap superconductor/normal metal interface. We demonstrate on several examples that those can be significantly different even in the strict domain of applicability of the Ginzburg-Landau theory. This justifies the use of this theory to describe relevant physics of two-gap superconductors, distinguishing them from their single-gap counterparts. The calculational degree of complexity increases only slightly with respect to the conventional Ginzburg-Landau expansion, thus the extended Ginzburg-Landau model remains numerically far less demanding compared to the full microscopic approaches.

preprint2011arXiv

Dirac and Klein-Gordon particles in one-dimensional periodic potentials

We evaluate the dispersion relation for massless fermions, described by the Dirac equation, and for zero-spin bosons, described by the Klein-Gordon equation, moving in two dimensions and in the presence of a one-dimensional periodic potential. For massless fermions the dispersion relation shows a zero gap for carriers with zero momentum in the direction parallel to the barriers in agreement with the well-known "Klein paradox". Numerical results for the energy spectrum and the density of states are presented. Those for fermions are appropriate to graphene in which carriers behave relativistically with the "light speed" replaced by the Fermi velocity. In addition, we evaluate the transmission through a finite number of barriers for fermions and zero-spin bosons and relate it with that through a superlattice.

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

Energy levels of triangular and hexagonal graphene quantum dots: a comparative study between the tight-binding and the Dirac approach

The Dirac equation is solved for triangular and hexagonal graphene quantum dots for different boundary conditions in the presence of a perpendicular magnetic field. We analyze the influence of the dot size and its geometry on their energy spectrum. A comparison between the results obtained for graphene dots with zigzag and armchair edges, as well as for infinite-mass boundary condition, is presented and our results show that the type of graphene dot edge and the choice of the appropriate boundary conditions have a very important influence on the energy spectrum. The single particle energy levels are calculated as function of an external perpendicular magnetic field which lifts degeneracies. Comparing the energy spectra obtained from the tight-binding approximation to those obtained from the continuum Dirac equation approach, we verify that the behavior of the energies as function of the dot size or the applied magnetic field are qualitatively similar, but in some cases quantitative differences can exist.

preprint2011arXiv

Enhancement of the retrapping current of superconducting microbridges of finite length

We theoretically find that the resistance of a superconducting microbridge/nanowire {\it decreases} while the retrapping current $I_r$ for the transition to the superconducting state {\it increases} when one suppresses the magnitude of the order parameter $|Δ|$ in the attached superconducting leads. This effect is a consequence of the increased energy interval for diffusion of the 'hot' nonequilibrium quasiparticles (induced by the oscillations of $|Δ|$ in the center of the microbridge) to the leads. The effect is absent in short microbridges (with length less than the coherence length) and it is relatively weak in long microbridges (with length larger than the inelastic relaxation length of the nonequilibrium distribution function). A nonmonotonous dependence of $I_r$ on the length of the microbridge is predicted. Our results are important for the explanation of the enhancement of the critical current and the appearance of negative magnetoresistance observed in many recent experiments on superconducting microbridges/nanowires.

preprint2011arXiv

Extended Ginzburg-Landau formalism for two-band superconductors

Recent observation of unusual vortex patterns in MgB_2 single crystals raised speculations about a new type of superconductivity (named "type-1.5"). However, the strict application of the standard two-band Ginzburg-Landau (GL) theory results in simply proportional order parameters of the two bands - and does not support the "type-1.5" behavior. Here we derive the extended GL formalism (accounting all terms of the next order over the small t= 1- T/T_c parameter) for a two-band clean s-wave superconductor and show that the two condensates generally have different spatial scales, with difference disappearing only in the limit T -> T_c. The extended version of the two-band GL formalism improves the validity of GL theory below T_c, and suggests revisiting of the earlier calculations based on the standard model.

preprint2011arXiv

Fano resonances and electron spin transport through a two-dimensional spin-orbit-coupled quantum ring

Electron transport through a spin-orbit-coupled quantum ring is investigated within linear response theory. We show that the finite width of the ring results in the appearance of Fano resonances in the conductance. This turns out to be a consequence of the spin-orbit interaction that leads to a breaking of the parity of the states localized in the ring. The resonances appear when the system is close to maxima of Aharonov-Casher conductance oscillations where spin transfer is heavily modified. When the spin-orbit coupling strength is detuned from the Aharonov-Casher maxima the resonances are broadened resulting in a dependence of the spin transport on the electron Fermi energy in contrast to predictions from one-dimensional models

preprint2011arXiv

Formation of multi-quanta vortices in superconductors: electronic, calorimetric and magnetic evidence

The ground state with vorticity larger than one in mesoscopic superconductors in applied magnetic field may manifest as a `giant'-vortex, where all vortices coalesce into a single singularity of the order parameter. Such a multi-quanta vortex may split into individual vortices (and vice versa) as a function of e.g. applied current, magnetic field or temperature. Here we show that such transitions can be identified by heat-capacity measurements, as the formation or splitting of a giant-vortex results in a clear jump in measured heat capacity vs. external drive. We attribute this phenomenon to an abrupt change in the density of states of the quasiparticle excitations in the vortex core(s), and further link it to a sharp change of the magnetic susceptibility at the transition - proving that formation of a giant-vortex can also be detected by conventional magnetometry.

preprint2011arXiv

Graphene nano ribbons subjected to axial stress

Atomistic simulations are used to study the bending of rectangular graphene nano ribbons subjected to axial stress both for free boundary and supported boundary conditions. The shape of the deformations of the buckled graphene nano ribbons, for small values of the stress, are sine waves where the number of nodal lines depend on the longitudinal size of the system and the applied boundary condition. The buckling strain for the supported boundary condition is found to be independent of the longitudinal size and estimated to be 0.86$%$. From a calculation of the free energy at finite temperature we find that the equilibrium projected two-dimensional area of the graphene nano ribbon is less than the area of a flat sheet. At the optimum length the boundary strain for the supported boundary condition is 0.48$%$.

preprint2011arXiv

Heating of quasiparticles driven by oscillations of the order parameter in short superconducting microbridges

We predict 'heating' of quasiparticles driven by order parameter oscillations in the resistive state of short superconducting microbridges. The finite relaxation time of the magnitude of the order parameter $|Δ|$ and the dependence of the spectral functions both on $|Δ|$ and the supervelocity $Q$ are the origin of this effect. Our result is opposite to those of Aslamazov and Larkin (Zh. Eks. Teor. Fiz. {\bf 70}, 1340 (1976)) and Schmid, Schon and Tinkham (Phys. Rev. B {\bf 21} 5076 (1980)) where 'cooling' of quasiparticles was found.

preprint2011arXiv

Kronig-Penney model on bilayer graphene: spectrum and transmission periodic in the strength of the barriers

We show that the transmission through single and double δ-function potential barriers of strength P in bilayer graphene is periodic in P with period π. For a certain range of P values we find states that are bound to the potential barrier and that run along the potential barrier. Similar periodic behaviour is found for the conductance. The spectrum of a periodic succession of δ-function barriers (Kronig-Penney model) in bilayer graphene is periodic in P with period 2π. For P smaller than a critical value, the spectrum exhibits two Dirac points while for P larger than this value an energy gap opens. These results are extended to the case of a superlattice of δ-function barriers with P alternating in sign between successive barriers; the corresponding spectrum is periodic in P with period π.

preprint2011arXiv

Landau levels in asymmetric graphene trilayers

The electronic spectrum of three coupled graphene layers (graphene trilayers) is investigated in the presence of an external magnetic field. We obtain analytical expressions for the Landau level spectrum for both the ABA and ABC - type of stacking, which exhibit very different dependence on the magnetic field. The effect of layer asymmetry and of external gate voltages can strongly influence the properties of the system.

preprint2011arXiv

Lattice thermal properties of Graphane: thermal contraction, roughness and heat capacity

Using atomistic simulations we determine the roughness and the thermal properties of a suspended graphane sheet. As compared to graphene we found that hydrogenated graphene has: 1) a larger thermal contraction, 2) the roughness exponent at room temperature is smaller, i.e. $\simeq$ 1.0 versus $\simeq$ 1.2 for graphene, 3) the wave lengths of the induced ripples in graphane cover a wide range corresponding to length scales in the range (30-125)\,Åat room temperature, and 4) the heat capacity of graphane is estimated to be 29.32$\pm$0.23\,J/molK which is 14.8% larger than the one for graphene, i.e. 24.98$\pm$0.14\,J/molK. Above 1500\,K we found that graphane buckles when its edges are supported in the $x-y$ plane.

preprint2011arXiv

Nanoindentation of a circular sheet of bilayer graphene

Nanoindentation of bilayer graphene is studied using molecular dynamics simulations. We compared our simulation results with those from elasticity theory as based on the nonlinear Föppl-Hencky equations with rigid boundary condition. The force deflection values of bilayer graphene are compered to those of monolayer graphene. Young's modulus of bilayer graphene is estimated to be 0.8 TPa which is close to the value for graphite. Moreover, an almost flat bilayer membrane at low temperature under central load has a 14$%$ smaller Young's modulus as compared to the one at room temperature.

preprint2011arXiv

Pattern formation and phase transitions in systems with non-monotonic interaction

We analyzed pattern formation and identified different phases in a system of particles interacting through a non-monotonic short-range repulsive (r<r_c) and long-range attractive (r>r_c) potential, using molecular-dynamics simulations. Depending on parameters, the interaction potential models the interparticle interaction in various physical systems ranging from atoms, molecules and colloids to neutron stars and vortices in low-kappa type-II superconductors and in recently discovered "type-1.5" superconductors. We constructed the phase diagram "critical radius r_c - density n" and proposed a new approach to characterization of the phases. Namely, we elaborated a set of quantitative criteria in order to identify the different phases using the Radial Distribution Function (RDF), the local density function, and the occupation factor.

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

Scattering of Dirac electrons by circular mass barriers: valley filter and resonant scattering

The scattering of two-dimensional (2D) massless Dirac electrons is investigated in the presence of a random array of circular mass barriers. The inverse momentum relaxation time and the Hall factor are calculated and used to obtain parallel and perpendicular resistivity components within linear transport theory. We found a non zero perpendicular resistivity component which has opposite sign for electrons in the different K and K' valleys. This property can be used for valley filter purposes. The total cross-section for scattering on penetrable barriers exhibit resonances due to the presence of quasi-bound states in the barriers that show up as sharp gaps in the cross-section while for Schrödinger electrons they appear as peaks.

preprint2011arXiv

Single-layer and bilayer graphene superlattices: collimation, additional Dirac points and Dirac lines

We review the energy spectrum and transport properties of several types of one- dimensional superlattices (SLs) on single-layer and bilayer graphene. In single-layer graphene, for certain SL parameters an electron beam incident on a SL is highly collimated. On the other hand there are extra Dirac points generated for other SL parameters. Using rectangular barriers allows us to find analytic expressions for the location of new Dirac points in the spectrum and for the renormalization of the electron velocities. The influence of these extra Dirac points on the conductivity is investigated. In the limit of δ-function barriers, the transmission T through, conductance G of a finite number of barriers as well as the energy spectra of SLs are periodic functions of the dimensionless strength P of the barriers, Pδ(x) ~ V (x). For a Kronig-Penney SL with alternating sign of the height of the barriers the Dirac point becomes a Dirac line for P = π/2 + nπ with n an integer. In bilayer graphene, with an appropriate bias applied to the barriers and wells, we show that several new types of SLs are produced and two of them are similar to type I and type II semiconductor SLs. Similar as in single-layer graphene extra "Dirac" points are found. Non-ballistic transport is also considered.

preprint2011arXiv

Strain-induced band gaps in bilayer graphene

We present a tight-binding investigation of strained bilayer graphene within linear elasticity theory, focusing on the different environments experienced by the A and B carbon atoms of the different sublattices. We find that the inequivalence of the A and B atoms is enhanced by the application of perpendicular strain $ε_{zz}$, which provides a physical mechanism for opening a band gap, most effectively obtained when pulling the two graphene layers apart. In addition, perpendicular strain introduces electron-hole asymmetry and can result in linear electronic dispersion near the K-point. When applying lateral strain to one layer and keeping the other layer fixed, we find the opening of an indirect band gap for small deformations. Our findings suggest experimental means for strain-engineered band gaps in bilayer graphene.

preprint2011arXiv

Superconducting nanofilms: molecule-like pairing induced by quantum confinement

Quantum confinement of the perpendicular motion of electrons in single-crystalline metallic superconducting nanofilms splits the conduction band into a series of single-electron subbands. A distinctive feature of such a nanoscale multi-band superconductor is that the energetic position of each subband can vary significantly with changing nanofilm thickness, substrate material, protection cover and other details of the fabrication process. It can occur that the bottom of one of the available subbands is situated in the vicinity of the Fermi level. We demonstrate that the character of the superconducting pairing in such a subband changes dramatically and exhibits a clear molecule-like trend, which is very similar to the well-known crossover from the Bardeen-Cooper-Schrieffer regime to Bose-Einstein condensation (BCS-BEC) observed in trapped ultracold fermions. For Pb nanofilms with thickness of 4 and 5 monolayers (ML) this will lead to a spectacular scenario: up to half of all the Cooper pairs nearly collapse, shrinking in the lateral size (parallel to the nanofilm) down to a few nanometers. As a result, the superconducting condensate will be a coherent mixture of almost molecule-like fermionic pairs with ordinary, extended Cooper pairs.

preprint2011arXiv

The electric field as a novel switch for uptake/release of hydrogen storage in nitrogen doped graphene

Nitrogen-doped graphene was recently synthesized and was reported to be a catalyst for hydrogen dissociative adsorption under a perpendicular applied electric field (F). In this work, the diffusion of H atoms on N-doped graphene, in the presence and absence of an applied perpendicular electric field, is studied using density functional theory. We demonstrate that the applied field can significantly facilitate the binding of hydrogen molecules on N-doped graphene through dissociative adsorption and diffusion on the surface. By removing the applied field the absorbed H atoms can be released efficiently. Our theoretical calculation indicates that N-doped graphene is a promising hydrogen storage material with reversible hydrogen adsorption/desorption where the applied electric field can act as a switch for the uptake/release processes.

preprint2011arXiv

Tuning of the spin-orbit interaction in a quantum dot by an in-plane magnetic field

Using an exact diagonalization approach we show that one- and two-electron InAs quantum dots exhibit avoided crossing in the energy spectra that are induced by the spin-orbit coupling in the presence of an in-plane external magnetic field. The width of the avoided crossings depends strongly on the orientation of the magnetic field which reveals the intrinsic anisotropy of the spin-orbit coupling interactions. We find that for specific orientations of the magnetic field avoided crossings vanish. Value of this orientation can be used to extract the ratio of the strength of Rashba and Dresselhaus interactions. The spin-orbit anisotropy effects for various geometries and orientations of the confinement potential are discussed. Our analysis explains the physics behind the recent measurements performed on a gated self-assembled quantum dot [S. Takahashi et al. Phys. Rev. Lett. 104, 246801 (2010)].

preprint2011arXiv

Valley-dependent Brewster angles and Goos-Hanchen effect in strained graphene

We demonstrate theoretically how local strains in graphene can be tailored to generate a valley polarized current. By suitable engineering of local strain profiles, we find that electrons in opposite valleys (K or K') show different Brewster-like angles and Goos-Hänchen shifts, exhibiting a close analogy with light propagating behavior. In a strain-induced waveguide, electrons in K and K' valleys have different group velocities, which can be used to construct a valley filter in graphene without the need for any external fields.

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

(Giant) Vortex - (anti) vortex interaction in bulk superconductors: The Ginzburg-Landau theory

The vortex-vortex interaction potential in bulk superconductors is calculated within the Ginzburg-Landau (GL) theory and is obtained from a numerical solution of a set of two coupled non-linear GL differential equations for the vector potential and the superconducting order parameter, where the merger of vortices into a giant vortex is allowed. Further, the interaction potentials between a vortex and a giant vortex and between a vortex and an antivortex are obtained for both type-I and type-II superconductors. Our numerical results agree asymptotically with the analytical expressions for large inter-vortex separations which are available in the literature. We propose new empirical expressions valid over the full interaction range, which are fitted to our numerical data for different values of the GL parameter.

preprint2010arXiv

Anticrossing-induced optical excitonic Aharonov-Bohm effect in strained type-I semiconductor nanorings

The exciton states in strained (In,Ga)As nanorings embedded in a GaAs matrix are computed. The strain distribution is extracted from the continuum mechanical model, and the exact diagonalization approach is employed to compute the exciton states. Weak oscillations of the ground exciton state energy with the magnetic field normal to the ring are an expression of the excitonic Aharonov-Bohm effect. Those oscillations arise from anticrossings between the ground and the second exciton state and can be enhanced by increasing the ring width. Simultaneously, the oscillator strength for exciton recombination exhibits oscillations, which are superposed on a linear increase with magnetic field. The obtained results are contrasted with previous theoretical results for 1D rings, and differences are explained to arise from different confinement potentials for the electron and the hole, and the large diamagnetic shift present in the analyzed type-I rings. Furthermore, our theory agrees qualitatively well with previous photoluminescence measurements on type-II InP/GaAs quantum dots.

preprint2010arXiv

Buckled circular monolayer graphene: a graphene nano-bowl

We investigate the stability of circular monolayer graphene subjected to a radial load using non-equilibrium molecular dynamics simulations. When monolayer graphene is radially stressed, after some small circular strain ($\sim 0.4%$) it buckles and bends into a new bowl like shape. Young's modulus is calculated from the linear relation between stress and strain before the buckling threshold, which is in agreement with experimental results. The prediction of elasticity theory for the buckling threshold of a radially stressed plate is presented and its results are compared to the one of our atomistic simulation. Jarzynski equality is used to estimate the difference between the free energy of the non-compressed states and the buckled states. From a calculation of the free energy we obtain the optimum radius for which the system feels the minimum boundary stress.

preprint2010arXiv

Carbon clusters: From ring structures to nanographene

The lowest energy configurations of Cn(n =< 55) clusters are obtained using the energy mini- mization technique with the conjugate gradient (CG) method where a modified Brenner potential is invoked to describe the carbon and hydrocarbon interaction. We found that the ground state configuration consists of a single ring for small number of C atoms and multi-ring structures are found with increasing n, which can be in planar, bowl-like or cap-like form. Contrary to previous predictions, the binding energy Eb does not show even-odd oscillations and only small jumps are found in the Eb(n) curve as a consequence of specific types of edges or equivalently the number of secondary atoms. We found that hydrogenation of the edge atoms may change the ground state configuration of the nanocluster. In both cases we determined the magic clusters. Special attention is paid to trigonal and hexagonal shaped carbon clusters and to clusters having a graphene-like configuration. Trigonal clusters are never the ground state, while hexagonal shaped clusters are only the ground state when they have zigzag edges.

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.

preprint2010arXiv

Condensate entanglement and multigap superconductivity in nanoscale superconductors

A Green functions approach is used to study superconductivity in nanofilms and nanowires. We show that the superconducting condensate results from the multimodal entanglement, or internal Josephson coupling, of the subcondensates associated with the manifold of Fermi surface subparts resulting from size-quantisation. This entanglement is of critical importance in these systems, since without it superconductivity would be extremely weak, if not completely negligible. Further, the multimodal character of the condensate generally results in multigap superconductivity, with great quantitative consequence for the value of the critical parameters. Our approach suggests that these are universal characteristics of confined superconductors.

preprint2010arXiv

D- shallow donor near a semiconductor-metal and a semiconductor-dielectric interface

The ground state energy and the extend of the wavefunction of a negatively charged donor (D-) located near a semiconductor-metal or a semiconductor-dielectric interface is obtained. We apply the effective mass approximation and use a variational two-electron wavefunction that takes into account the influence of all image charges that arise due to the presence of the interface, as well as the correlation between the two electrons bound to the donor. For a semiconductor-metal interface, the D- binding energy is enhanced for donor positions d>1.5a_B (a_B is the effective Bohr radius) due to the additional attraction of the electrons with their images. When the donor approaches the interface (i.e. d<1.5a_B) the D- binding energy drops and eventually it becomes unbound. For a semiconductor-dielectric (or a semiconductor-vacuum) interface the D- binding energy is reduced for any donor position as compared to the bulk case and the system becomes rapidly unbound when the donor approaches the interface.

preprint2010arXiv

Effect of correlated noise on quasi-1D diffusion

Single-file diffusion (SFD) of an infinite one-dimensional chain of interacting particles has a long-time mean-square displacement (MSD) ~t^1/2, independent of the type of inter-particle repulsive interaction. This behavior is also observed in finite-size chains, although only for certain intervals of time t depending on the chain length L, followed by the ~t for t->infinity, as we demonstrate for a closed circular chain of diffusing interacting particles. Here we show that spatial correlation of noise slows down SFD and can result, depending on the amount of correlated noise, in either subdiffusive behavior ~t^alpha, where 0<alpha<1/2, or even in a total suppression of diffusion (in the limit N-> infinity). Spatial correlation can explain the subdiffusive behavior in recent SFD experiments in circular channels.

preprint2010arXiv

Electronic and magnetic properties of superlattices of graphene/graphane nanoribbons with different edge hydrogenation

Zigzag graphene nanoribbons patterned on graphane are studied using spin-polarized ab initio calculations. We found that the electronic and magnetic properties of the graphene/graphane superlattice strongly depends on the degree of hydrogenation at the interfaces between the two materials. When both zigzag interfaces are fully hydrogenated, the superlattice behaves like a freestanding zigzag graphene nanoribbon, and the magnetic ground state is antiferromagnetic. When one of the interfaces is half hydrogenated, the magnetic ground state becomes ferromagnetic, and the system is very close to being a half metal with possible spintronics applications whereas the magnetic ground state of the superlattice with both interfaces half hydrogenated is again antiferromagnetic. In this last case, both edges of the graphane nanoribbon also contribute to the total magnetization of the system. All the spin-polarized ground states are semiconducting, independent of the degree of hydrogenation of the interfaces. The ab initio results are supplemented by a simple tight-binding analysis that captures the main qualitative features. Our ab initio results show that patterned hydrogenation of graphene is a promising way to obtain stable graphene nanoribbons with interesting technological applications.

preprint2010arXiv

Enhanced stability of hydrogen atoms at the graphene/graphane interface of nanoribbons

The thermal stability of graphene/graphane nanoribbons (GGNRs) is investigated using density functional theory. It is found that the energy barriers for the diffusion of hydrogen atoms on the zigzag and armchair interfaces of GGNRs are 2.86 and 3.17 eV, respectively, while the diffusion barrier of an isolated H atom on pristine graphene was only ~0.3 eV. These results unambiguously demonstrate that the thermal stability of GGNRs can be enhanced significantly by increasing the hydrogen diffusion barriers through graphene/graphane interface engineering. This may provide new insights for viable applications of GGNRs.

preprint2010arXiv

Enhancing the critical current in quasiperiodic pinning arrays below and above the matching magnetic flux

Quasiperiodic pinning arrays, as recently demonstrated theoretically and experimentally using a five-fold Penrose tiling, can lead to a significant enhancement of the critical current Ic as compared to "traditional" regular pinning arrays. However, while regular arrays showed only a sharp peak in Ic(Phi) at the matching flux Phi1 and quasiperiodic arrays provided a much broader maximum at Phi<Phi1, both types of pinning arrays turned out to be inefficient for fluxes larger than Phi1. We demonstrate theoretically and experimentally the enhancement of Ic(Phi) for Phi>Phi1 by using non-Penrose quasiperiodic pinning arrays. This result is based on a qualitatively different mechanism of flux pinning by quasiperiodic pinning arrays and could be potentially useful for applications in superconducting micro-electronic devices operating in a broad range of magnetic fields.

preprint2010arXiv

Extra Dirac points in the energy spectrum for superlattices on single-layer graphene

We investigate the emergence of extra Dirac points in the electronic structure of a periodically spaced barrier system, i.e., a superlattice, on single-layer graphene, using a Dirac-type Hamiltonian. Using square barriers allows us to find analytic expressions for the occurrence and location of these new Dirac points in k-space and for the renormalization of the electron velocity near them in the low-energy range. In the general case of unequal barrier and well widths the new Dirac points move away from the Fermi level and for given heights of the potential barriers there is a minimum and maximum barrier width outside of which the new Dirac points disappear. The effect of these extra Dirac points on the density of states and on the conductivity is investigated.

preprint2010arXiv

First-principles investigation of graphene fluoride and graphane

Different stoichiometric configurations of graphane and graphene fluoride are investigated within density functional theory. Their structural and electronic properties are compared, and we indicate the similarities and differences among the various configurations. Large differences between graphane and graphene fluoride are found that are caused by the presence of charges on the fluorine atoms. A new configuration that is more stable than the boat configuration is predicted for graphene fluoride. We also perform GW calculations for the electronic band gap of both graphene derivatives. These band gaps and also the calculated Young's moduli are at variance with available experimental data. This might indicate that the experimental samples contain a large number of defects or are only partially covered with H or F.

preprint2010arXiv

Geometry-induced localization of thermal fluctuations in ultrathin superconducting structures

Thermal fluctuations of the order parameter in an ultrathin triangular shaped superconducting structure are studied near $T_{c}$, in zero applied field. We find that the order parameter is prone to much larger fluctuations in the corners of the structure as compared to its interior. This geometry-induced localization of thermal fluctuations is attributed to the fact that condensate confinement in the corners is characterised by a lower effective dimensionality, which favors stronger fluctuations.

preprint2010arXiv

Guided nucleation of superconductivity on a graded magnetic substrate

We demonstrate the controlled spatial nucleation of superconductivity in a thin film deposited on periodic arrays of ferromagnetic dots with gradually increasing diameter. The perpendicular magnetization of the dots induces vortex-antivortex molecules in the sample, with the number of (anti)vortices increasing with magnet size. The resulting gradient of antivortex density between the dots predetermines local nucleation of superconductivity in the sample as a function of the applied external field and temperature. In addition, the compensation between the applied magnetic field and the antivortices results in an unprecedented enhancement of the critical temperature.

preprint2010arXiv

Metallic nanograins: spatially nonuniform pairing induced by quantum confinement

It is well-known that the formation of discrete electron levels strongly influences the pairing in metallic nanograins. Here we focus on another effect of quantum confinement in superconducting grains that was not studied previously, i.e., spatially nonuniform pairing. This effect is very significant when single-electron levels form bunches and/or a kind of shell structure: in highly symmetric grains the order parameter can exhibit variations with position by an order of magnitude. Nonuniform pairing is closely related to a quantum-confinement induced modification of the pairing-interaction matrix elements and size-dependent pinning of the chemical potential to groups of degenerate or nearly degenerate levels. For illustration we consider spherical metallic nanograins. We show that the relevant matrix elements are as a rule enhanced in the presence of quantum confinement, which favors spatial variations of the order parameter, compensating the corresponding energy cost. The size-dependent pinning of the chemical potential further increases the spatial variation of the pair condensate. The role of nonuniform pairing is smaller in less symmetric confining geometries and/or in the presence of disorder. However, it always remains of importance when the energy spacing between discrete electron levels $δ$ is approaching the scale of the bulk gap $Δ_B$, i.e., $δ> 0.1$-$0.2\,Δ_B$.

preprint2010arXiv

Proximity induced pseudogap in mesoscopic superconductor/normal-metal bilayers

Recent scanning tunneling microscopy measurements of the proximity effect in Au/La$_{2-x}$Sr$_{x}$CuO$_{4}$ and La$_{1.55}$Sr$_{0.45}$CuO$_{4}$/La$_{2-x}$Sr$_{x}$CuO$_{4}$ bilayers showed a proximity-induced pseudogap [Yuli et al., Phys. Rev. Lett. {\bf 103}, 197003 (2009)]. We describe the proximity effect in mesoscopic superconductor/normal-metal bilayers by using the Bogoliubov-de Gennes equations for a tight-binding Hamiltonian with competing antiferromagnetic and d-wave superconductivity orders . The temperature dependent local density of states is calculated as a function of the distance from the interface. Bound state due to both d-wave and spin density wave gaps are formed in the normal metal for energies less than the respective gaps. If there is a mismatch between the Fermi velocities in the two layers we observe that these states will shift in energy when spin density wave order is present, thus inducing a minigap at finite energy. We conclude that the STM measurement in the proximity structures is able to distinguish between the two scenarios proposed for the pseudogap (competing or precursor to superconductivity).

preprint2010arXiv

Simplified model for the energy levels of quantum rings in single layer and bilayer graphene

Within a minimal model, we present analytical expressions for the eigenstates and eigenvalues of carriers confined in quantum rings in monolayer and bilayer graphene. The calculations were performed in the context of the continuum model, by solving the Dirac equation for a zero width ring geometry, i.e. by freezing out the carrier radial motion. We include the effect of an external magnetic field and show the appearance of Aharonov-Bohm oscillations and of a non-zero gap in the spectrum. Our minimal model gives insight in the energy spectrum of graphene-based quantum rings and models different aspects of finite width rings.

preprint2010arXiv

Stacking Order dependent Electric Field tuning of the Band Gap in Graphene Multilayers

The effect of different stacking order of graphene multilayers on the electric field induced band gap is investigated. We considered a positively charged top and a negatively charged back gate in order to independently tune the band gap and the Fermi energy of three and four layer graphene systems. A tight-binding approach within a self-consistent Hartree approximation is used to calculate the induced charges on the different graphene layers. We found that the gap for trilayer graphene with the ABC stacking is much larger than the corresponding gap for the ABA trilayer. Also we predict that for four layers of graphene the energy gap strongly depends on the choice of stacking, and we found that the gap for the different types of stacking is much larger as compared to the case of Bernal stacking. Trigonal warping changes the size of the induced electronic gap by approximately 30% for intermediate and large values of the induced electron density.

preprint2010arXiv

Structural and dynamical properties of a quasi-one-dimensional classical binary system

The ground state configurations and the \lq{}\lq{}normal\rq{}\rq{} mode spectra of a $quasi$-one-dimensional (Q1D) binary system of charged particles interacting through a screened Coulomb potential are presented. The minimum energy configurations were obtained analytically and independently through molecular dynamic simulations. A rich variety of ordered structures were found as a function of the screening parameter, the particle density, and the ratio between the charges of the distinct types of particles. Continuous and discontinuous structural transitions, as well as an unexpected symmetry breaking in the charge distribution are observed when the density of the system is changed. For near equal charges we found a disordered phase where a mixing of the two types of particles occurs. The phonon dispersion curves were calculated within the harmonic approximation for the one- and two-chain structures.

preprint2010arXiv

The ground states of the two-component order parameter superconductor

We show that in presence of an applied external field the two-component order parameter superconductor falls in two categories of ground states, namely, in the traditional Abrikosov ground state or in a new ground state fitted to describe a superconducting layer with texture, that is, patched regions separated by a phase difference of $π$. The existence of these two kinds of ground states follows from the sole assumption that the total supercurrent is the sum of the two individual supercurrents and is independent of any consideration about the free energy expansion. Uniquely defined relations between the current density and the superfluid density hold for these two ground states, which also determine the magnetization in terms of average values of the order parameters. Because these ground state conditions are also Bogomolny equations we construct the free energy for the two-component superconductor which admits the Bogomolny solution at a special coupling value.

preprint2010arXiv

Topological confinement in graphene bilayer quantum rings

We demonstrate the existence of localized electron and hole states in a ring-shaped potential kink in biased bilayer graphene. Within the continuum description, we show that for sharp potential steps the Dirac equation describing carrier states close to the K (or K') point of the first Brillouin zone can be solved analytically for a circular kink/anti-kink dot. The solutions exhibit interfacial states which exhibit Aharonov-Bohm oscillations as functions of the height of the potential step and/or the radius of the ring.

preprint2010arXiv

Two-component mixture of charged particles confined in a channel: melting

The melting of a binary system of charged particles confined in a {\it quasi}-one-dimensional parabolic channel is studied through Monte Carlo simulations. At zero temperature the particles are ordered in parallel chains. The melting is anisotropic and different melting temperatures are obtained according to the spatial direction, and the different types of particles present in the system. Melting is very different for the single-, two- and four-chain configurations. A temperature induced structural phase transition is found between two different four chain ordered states which is absent in the mono-disperse system. In the mixed regime, where the two types of particles are only slightly different, melting is almost isotropic and a thermally induced homogeneous distribution of the distinct types of charges is observed.

preprint2010arXiv

Vibrational properties of graphene fluoride and graphane

The vibrational properties of graphene fluoride and graphane are studied using ab initio calculations. We find that both sp3 bonded derivatives of graphene have different phonon dispersion relations and phonon density of states as expected from the different masses associated with the attached atoms fluorine and hydrogen, respectively. These differences manifest themselves in the predicted temperature behavior of the constant-volume specific heat of both compounds.

preprint2010arXiv

Vortex manipulation in a superconducting matrix with view on applications

We show how a single flux quantum can be effectively manipulated in a superconducting film with a matrix of blind holes. Such a sample can serve as a basic memory element, where the position of the vortex in a [k x l] matrix of pinning sites defines the desired combination of n bits of information (2^n=k*l). Vortex placement is achieved by strategically applied current and the resulting position is read-out via generated voltage between metallic contacts on the sample. Such a device can also act as a controllable source of a nanoengineered local magnetic field for e.g. spintronics applications.

preprint2010arXiv

Vortex matter in mesoscopic two-gap superconducting disks: influence of Josephson and magnetic coupling

The effects of the coupling between two electronic condensates in two-gap mesoscopic superconductors are studied within the Ginzburg-Landau theory using a finite difference technique. In applied magnetic field, we derive the dependency of the size of the vortex on the sample size and the strength of the Josephson coupling. In addition, we elaborate on the dependence of the critical temperature and field on the parameters of the coupled condensates. We demonstrate further the existence and stability of fractional states, for which the two condensates comprise different vorticity. Moreover, we also found pronounced asymmetric fractional states and we show their experimentally observable magnetic response. Finally we introduce the magnetic coupling between condensates, and study in particular the case where one band is type II and the other is type I, i.e. the sample is effectively of I.x type. The calculated M(H) loops show a clear signature of the mixed type of superconductivity, which we find to be strongly affected by the ratio of the coherence lengths in the two condensates.

preprint2010arXiv

Vortex patterns in a superconducting-ferromagnetic rod

A superconducting rod with a magnetic moment on top develops vortices obtained here through 3D calculations of the Ginzburg-Landau theory. The inhomogeneity of the applied field brings new properties to the vortex patterns that vary according to the rod thickness. We find that for thin rods (disks) the vortex patterns are similar to those obtained in presence of a homogeneous magnetic field instead because they consist of giant vortex states. For thick rods novel patterns are obtained as vortices are curve lines in space that exit through the lateral surface.

preprint2010arXiv

Yukawa particles confined in a channel and subject to a periodic potential: ground state and normal modes

We consider a classical system of two-dimensional (2D) charged particles, which interact through a repulsive Yukawa potential $exp(-r/λ)/r$, confined in a parabolic channel which limits the motion of the particles in the $y$-direction. Along the $x$-direction, the particles are also subject to a periodic potential substrate. The ground state configurations and the normal mode spectra of the system are obtained as function of the periodicity and strength of the periodic potential ($V_0$), and density. An interesting set of tunable ground state configurations are found, with first and second order structural transitions between them. A magic configuration with particles aligned in each minimum of the periodic potential is obtained for V_0 larger than some critical value which has a power law dependence on the density. The phonon spectrum of different configurations were also calculated. A localization of the modes into a small frequency interval is observed for a sufficient strength of the periodic potential. A tunable band-gap is found as a function of $V_0$. This model system can be viewed as a generalization of the Frenkel and Kontorova model.

preprint2009arXiv

Confinement effects on intermediate state flux patterns in mesoscopic type-I superconductors

Intermediate state (IS) flux structures in mesoscopic type-I superconductors are investigated within the Ginzburg-Landau theory. In addition to well-established tubular and laminar structures, the strong confinement leads to the formation of (i) a phase of singly quantized vortices, which is typical for type-II superconductors and (ii) a ring of a normal domain at equilibrium. The stability region and the formation process of these IS flux structures are strongly influenced by the geometry of the sample.

preprint2009arXiv

Electrostatically confined Quantum Rings in bilayer Graphene

We propose a new system where electron and hole states are electrostatically confined into a quantum ring in bilayer graphene. These structures can be created by tuning the gap of the graphene bilayer using nanostructured gates or by position-dependent doping. The energy levels have a magnetic field ($B_{0}$) dependence that is strikingly distinct from that of usual semiconductor quantum rings. In particular, the eigenvalues are not invariant under a $B_0 \to -B_0$ transformation and, for a fixed total angular momentum index $m$, their field dependence is not parabolic, but displays two minima separated by a saddle point. The spectra also display several anti-crossings, which arise due to the overlap of gate-confined and magnetically-confined states.

preprint2009arXiv

Kink-antikink vortex transfer in periodic-plus-random pinning potential: Theoretical analysis and numerical experiments

The influence of random pinning on the vortex dynamics in a periodic square potential under an external drive is investigated. Using theoretical approach and numerical experiments, we found several dynamical phases of vortex motion that are different from the ones for a regular pinning potential. Vortex transfer is controlled by kinks and antikinks, which either preexist in the system or appear spontaneously in pairs and then propagate in groups. When kinks and antikinks collide, they annihilate.

preprint2009arXiv

Local current injection into mesoscopic superconductors for the manipulation of quantum states

We perform strategic current injection in a small mesoscopic superconductor and control the (non)equilibrium quantum states in an applied homogeneous magnetic field. In doing so, we realize a current-driven splitting of multi-quanta vortices, current-induced transitions between states with different angular momenta, and current-controlled switching between otherwise degenerate quantum states. These fundamental phenomena form the basis for discussed electronic and logic applications, and are confirmed in both theoretical simulations and multiple-small-tunnel-junction transport measurements.

preprint2009arXiv

Magnetic Kronig-Penney model for Dirac electrons in single-layer graphene

The properties of Dirac electrons in a magnetic superlattice (SL) on graphene consisting of very high and thin (delta-function) barriers are investigated. We obtain the energy spectrum analytically and study the transmission through a finite number of barriers. The results are contrasted with those for electrons described by the Schrodinger equation. In addition, a collimation of an incident beam of electrons is obtained along the direction perpendicular to that of the SL. We also highlight the analogy with optical media in which the refractive index varies in space.

preprint2009arXiv

Quantum rings with time dependent spin-orbit coupling: Rabi oscillations, spintronic Schrodinger-cat states, and conductance properties

The strength of the (Rashba-type) spin-orbit coupling in mesoscopic semiconductor rings can be tuned with external gate voltages. Here we consider the case of a periodically changing spin-orbit interaction strength as induced by sinusoidal voltages. In a closed one dimensional quantum ring with weak spin-orbit coupling, Rabi oscillations are shown to appear. We find that the time evolution of initially localized wave packets exhibits a series of collapse and revival phenomena. Partial revivals -- that are typical in nonlinear systems -- are shown to correspond to superpositions of states localized at different spatial positions along the ring. These "spintronic Schrodinger-cat sates" appear periodically, and similarly to their counterparts in other physical systems, they are found to be sensitive to environment induced disturbances. The time dependent spin transport problem, when leads are attached to the ring, is also solved. We show that the "sideband currents" induced by the oscillating spin-orbit interaction strength can become the dominant output channel, even in the presence of moderate thermal fluctuations and random scattering events.

preprint2009arXiv

Resonant Tunneling through S- and U-shaped Graphene Nanoribbons

We theoretically investigate resonant tunneling through S- and U-shaped nanostructured graphene nanoribbons. A rich structure of resonant tunneling peaks are found eminating from different quasi-bound states in the middle region. The tunneling current can be turned on and off by varying the Fermi energy. Tunability of resonant tunneling is realized by changing the width of the left and/or right leads and without the use of any external gates.

preprint2009arXiv

Second generation of vortex-antivortex states in mesoscopic superconductors: stabilization by artificial pinning

Antagonistic symmetries of superconducting polygons and their induced multi-vortex states in a homogeneous magnetic field may lead to appearance of antivortices in the vicinity of the superconducting/normal state boundary (where mesoscopic confinement is particularly strong). Resulting vortex-antivortex (V-Av) molecules match the sample symmetry, but are extremely sensitive to defects and fluctuations and remain undetected experimentally. Here we show that V-Av states can re-appear deep in the superconducting state due to an array of perforations in a polygonal setting, surrounding a central hole. Such states are no longer caused by the symmetry of the sample but rather by pinning itself, which prevents the vortex-antivortex annihilation. As a result, even micron-size, clearly spaced V-Av molecules can be stabilized in large mesoscopic samples.

preprint2009arXiv

Stability of spintronic devices based on quantum ring networks

Transport properties in mesoscopic networks are investigated, where the strength of the (Rashba-type) spin-orbit coupling is assumed to be tuned with external gate voltages. We analyze in detail to what extent the ideal behavior and functionality of some promising network-based devices are modified by random (spin-dependent) scattering events and by thermal fluctuations. It is found that although the functionality of these devices is obviously based on the quantum coherence of the transmitted electrons, there is a certain stability: moderate level of errors can be tolerated. For mesoscopic networks made of typical semiconductor materials, even cryogenic temperatures can smear out the desired transport properties. When the energy distribution of the input carriers is narrow enough, it turns out that the devices can operate close to their ideal limits even at relative high temperature. As an example, we present results for two different networks: one that realizes a Stern-Gerlach device and another that simulates a spin quantum walker. Finally we propose a simple network that can act as a narrow band energy filter even in the presence of random scatterers.

preprint2009arXiv

Superconducting nanowires: quantum confinement and spatially dependent Hartree-Fock potential

It is well-known that in bulk, the solution of the Bogoliubov-de Gennes equations is the same whether or not the Hartree-Fock term is included. Here the Hartree-Fock potential is position independent and, so, gives the same contribution to both the single-electron energies and the Fermi level (the chemical potential). Thus, the single-electron energies measured from the Fermi level (they control the solution) stay the same. It is not the case for nanostructured superconductors, where quantum confinement breaks the translational symmetry and results in a position dependent Hartree-Fock potential. In this case its contribution to the single-electron energies depends on the relevant quantum numbers. We numerically solved the Bogoliubov-de Gennes equations with the Hartree-Fock term for a clean superconducting nanocylinder and found a shift of the curve representing the thickness-dependent oscillations of the critical superconducting temperature to larger diameters.

preprint2009arXiv

Transverse instabilities of multiple vortex chains in superconductor-ferromagnet bilayers

Using scanning tunneling microscopy and Ginzburg-Landau simulations we explore vortex configurations in magnetically coupled NbSe$_2$-Permalloy superconductor-ferromagnet bilayer. The Permalloy film with stripe domain structure induces periodic local magnetic induction in the superconductor creating a series of pinning-antipinning channels for externally added magnetic flux quanta. Such laterally confined Abrikosov vortices form quasi-1D arrays (chains). The transitions between multichain states occur through propagation of kinks at the intermediate fields. At high fields we show that the system becomes non-linear due to a change in both the number of vortices and the confining potential. The longitudinal instabilities of the resulting vortex structures lead to vortices `levitating' in the anti-pinning channels.

preprint2007arXiv

Symmetric and asymmetric states in mesoscopic superconducting wire in voltage driven regime

The response of a mesoscopic homogeneous superconducting wire, connected with bulk normal metal reservoirs, is theoretically investigated as function of the applied voltage. The finite relaxation length of the nonequilibrium quasiparticle distribution function $\bar{L_E}$ is included where we assumed that our wire is in the dirty limit. We found that {\it both} symmetric and asymmetric states can exist which are characterized by a stationary symmetric and asymmetric distribution of the order parameter with respect to the center of the wire. Current voltage characteristics of the wire with length $L>\bar{L_E}$ being in the symmetric state show a pronounced S-behavior. The asymmetric state may exist only for voltages larger than some critical value and coexist with the symmetric state in a finite voltage interval. For wires with $L \sim \bar{L_E}$ the asymmetric state survives up to higher values of the voltage than the symmetric one and may exist both in the voltage and the current driven regimes. We propose an experiment to observe reversible switching between those stationary symmetric and asymmetric states.

preprint2006arXiv

Quantum dot size dependent influence of the substrate orientation on the electronic and optical properties of InAs/GaAs quantum dots

Using 3D k.p calculation including strain and piezoelectricity we predict variation of electronic and optical properties of InAs/GaAs quantum dots (QDs) with the substrate orientation. The QD transition energies are obtained for high index substrates [11k], where k = 1,2,3 and are compared with [001]. We find that the QD size in the growth direction determines the degree of influence of the substrate orientation: the flatter the dots, the larger the difference from the reference [001] case.

preprint2006arXiv

Superconducting films with antidot arrays - novel behavior of the critical current

Novel behavior of the critical current density $j_{c}$ of a regularly perforated superconducting film is found, as a function of applied magnetic field $H$. Previously pronounced peaks of $j_{c}$ at matching fields were always found to decrease with increasing $H$. Here we found a {\it reversal of this behavior} for particular geometrical parameters of the antidot lattice and/or temperature. This new phenomenon is due to a strong ``caging'' of interstitial vortices between the pinned ones. We show that this vortex-vortex interaction can be further tailored by an appropriate choice of the superconducting material, described by the Ginzburg-Landau parameter $κ$. In effective type-I samples we predict that the peaks in $j_{c}(H)$ at the matching fields are transformed into a {\it step-like behavior}.

preprint2006arXiv

Vortex configurations and critical parameters in superconducting thin films containing antidot arrays: Nonlinear Ginzburg-Landau theory

Using the non-linear Ginzburg-Landau (GL) theory, we obtain the possible vortex configurations in superconducting thin films containing a square lattice of antidots. The equilibrium structural phase diagram is constructed which gives the different ground-state vortex configurations as function of the size and periodicity of the antidots for a given effective GL parameter $κ^{*}$. Giant-vortex states, combination of giant- and multi-vortex states, as well as symmetry imposed vortex-antivortex states are found to be the ground state for particular geometrical parameters of the sample. The antidot occupation number $n_o$ is calculated as a function of related parameters and comparison with existing expressions for the saturation number $n_s$ and with experimental results is given. For a small radius of antidots a triangular vortex lattice is obtained, where some of the vortices are pinned by the antidots and some of them are located between them. Transition between the square pinned and triangular vortex lattices is given for different values of the applied field. The enhanced critical current at integer and rational matching fields is found, where the level of enhancement at given magnetic field directly depends on the vortex-occupation number of the antidots. For certain parameters of the antidot lattice and/or temperature the critical current is found to be larger for higher magnetic fields. Superconducting/normal $H-T$ phase boundary exhibits different regimes as antidots are made larger, and we transit from a plain superconducting film to a thin-wire superconducting network. Presented results are in good agreement with available experiments and suggest possible new experiments.

preprint2005arXiv

Vortex-Antivortex Lattices in Superconducting Films with Magnetic Pinning Arrays

Novel vortex structures are found when a thin superconducting film (SC) is covered with a lattice of out-of-plane magnetized magnetic dots (MDs). The stray magnetic field of the dots confines the vortices to the MD regions, surrounded by antivortices which "crystallize" into regular lattices. First and second order transitions are found as magnetic array is made sparser or MD-magnetization larger. For sparse MD-arrays fractional vortex-antivortex states are formed, where the crystal-symmetry is combined with a non-uniform "charge" distribution. We demonstrate that due to the (anti)vortices and the supercurrents induced by the MDs, the critical current of the sample actually increases if exposed to a homogeneous external magnetic field, contrary to conventional SC behavior.

preprint2003arXiv

Current distribution in a two dimensional electron gas exposed to a perpendicular nonhomogeneous magnetic field of a chess configuration

We have studied a finite two-dimensional electron system exposed to a normal nonhomogeneous magnetic field of a chess configuration. Using the conformal mapping method we obtain an exact analytical solution for the electric field distribution in terms of the Jacobi functions. The obtained formula is exploited to calculate the physical quantities of interest: the current density distribution, the linear density of charges accumulated along the magnetic interfaces, the magneto-resistance, and the Hall resistance.

preprint2002arXiv

Spin-current modulation and square-wave transmission through periodically stubbed electron waveguides

Ballistic spin transport through waveguides, with symmetric or asymmetric double stubs attached to them periodically, is studied systematically in the presence of a weak spin-orbit coupling that makes the electrons precess. By an appropriate choice of the waveguide length and of the stub parameters injected spin-polarized electrons can be blocked completely and the transmission shows a periodic and nearly square-type behavior, with values 1 and 0, with wide gaps when only one mode is allowed to propagate in the waveguide. A similar behavior is possible for a certain range of the stub parameters even when two-modes can propagate in the waveguide and the conductance is doubled. Such a structure is a good candidate for establishing a realistic spin transistor. A further modulation of the spin current can be achieved by inserting defects in a finite-number stub superlattice. Finite-temperature effects on the spin conductance are also considered.

preprint2001arXiv

Magneto-exciton in single and coupled type II quantum dots

We studied the exciton energy in a type II quantum disk as a function of the magnetic field, disk radius R and height d. We found angular momentum transitions for dots with d >> 2R. Application of an electric field perpendicular to the disk showed a non-linear Stark shift. In the case of three vertically coupled dots angular momentum transistions were found for small interdot distances which disappeared with increasing interdot distance.

preprint2001arXiv

Spin-engineered quantum dots

Spatially nonhomogeneously spin polarized nuclei are proposed as a new mechanism to monitor electron states in a nanostructure, or as a means to createn and, if necessary, reshape such nanostructures in the course of the experiment. We found that a polarization of nulear spins may lift the spin polarization of the electron states in a nanostructure and, if sufficiently strong, leads to a polarization of the electron spins. Polarized nuclear spins may form an energy landscape capable of binding electrons with energy up to several meV and the localization radius $ >$ 100Å.

preprint2000arXiv

Ground state of excitons and charged excitons in a quantum well

A variational calculation of the ground state of a neutral exciton and of positively and negatively charged excitons (trions) in single quantum well is presented. We study the dependance of the correlation energy and of the binding energy on the well width and on the hole mass. Our results are are compared with previous theoretical results and with avalaible experimental data.

preprint2000arXiv

Non-homogeneous magnetic field induced magnetic edge states and their transport in a quantum wire

The spectrum of magnetic edge states and their transport properties in the presence of a perpendicular non-homogeneous magnetic field in a quantum wire formed by a parabolic confining potential are obtained. Systems are studied where the magnetic field exhibits a discontinuous jump in the transverse direction and changes its sign, strength, and both sign and strength at the magnetic interface. The energy spectra and wave functions of these systems, the corresponding group velocities along the interface and the particle average positions normal to the interface are calculated. The resistance of the quantum wire in the presence of such a magnetic interface is obtained both in the ballistic and the diffusive regimes as a function of the Fermi energy and of the homogeneous background magnetic field. The results are compared with those for the case of a homogeneous field.

preprint2000arXiv

Non-Quantized Penetration of Magnetic Field in the Vortex State of Superconductors

As first pointed out by Bardeen and Ginzburg in the early sixties, the amount of magnetic flux carried by vortices depends on their distance to the sample edge and can be smaller than one flux quantum, f0 = h/2e. In bulk superconductors, this reduction of flux becomes negligible already at submicron distances from the edge but, in thin films, the effect may survive at much larger distances. In the absence of any experimental observation, such flux reduction is perceived to be an exotic or unimportant effect, and it is often assumed that magnetic field enters type-II superconductors in units of f0. Here we report the measurements of the amount of flux associated with the entrance of individual vortices in a superconducting film and show that the flux that they bring in, always differs substantially from f0. We have observed vortices that carry as little as 0.001 f0 as well as seemingly "negative vortices" whose penetration leads to the expulsion of magnetic field. We distinguish two phenomena responsible for non-quantized flux penetration: the finite-size effect and a non-linear screening of magnetic field in the presence of a surface barrier. The latter effect has previously not been considered but is likely to cause non-quantized penetration in many cases.

preprint1999arXiv

Enhanced stability of the square lattice of a classical bilayer Wigner crystal

The stability and melting transition of a single layer and a bilayer crystal consisting of charged particles interacting through a Coulomb or a screened Coulomb potential is studied using the Monte-Carlo technique. A new melting criterion is formulated which we show to be universal for bilayer as well as for single layer crystals in the case of (screened) Coulomb, Lennard--Jones and 1/r^{12} repulsive inter-particle interactions. The melting temperature for the five different lattice structures of the bilayer Wigner crystal is obtained, and a phase diagram is constructed as a function of the interlayer distance. We found the surprising result that the square lattice has a substantial larger melting temperature as compared to the other lattice structures. This is a consequence of the specific topology of the defects which are created with increasing temperature and which have a larger energy as compared to the defects in e.g. a hexagonal lattice.

preprint1999arXiv

Transitions between different superconducting states in mesoscopic disks

Using a linear analysis, we study the stability of giant-vortex states in very thin disks. The vortex expulsion and penetration fields are obtained for finite thickness disks from a numerical solution of the non-linear Ginzburg-Landau (GL) equations. Using an extension of the London approximation, in which the phase distribution of the order parameter is prescribed and the superconducting electron density is found numerically, we consider the free energy behavior for transitions between different superconducting states.

preprint1998arXiv

Binding energy of localized biexcitons in quantum wells

A variational calculation of the ground state energy of a biexciton in a GaAs/AlGaAs quantum well is presented. The well width fluctuations leading to trapping of the biexcitons are modeled by a parabolic potential. The results obtained for different well widths are compared with recent experimental data. Good agreement is obtained both for the biexciton binding energy and for the Haynes factor. We find that the structure of a biexciton is similar to the one of the H_2 molecule.

preprint1998arXiv

Vortex phase diagram for mesoscopic superconducting disks

Solving numerically the 3D non linear Ginzburg-Landau (GL) equations, we study equilibrium and nonequilibrium phase transitions between different superconducting states of mesoscopic disks which are thinner than the coherence length and the penetration depth. We have found a smooth transition from a multi-vortex superconducting state to a giant vortex state with increasing both the disk thickness and the magnetic field. A vortex phase diagram is obtained which shows, as function of the magnetic field, a re-entrant behavior between the multi-vortex and the giant vortex state.

F. M. Peeters

What is connected

Connect this record

See the researcher in context

Building this map preview

230 published item(s)

Alternating superconducting and charge density wave monolayers within bulk 6R-TaS2

Anisotropic and tunable optical conductivity of a two-dimensional semi-Dirac system in the presence of elliptically polarized radiation

Electronic Mach-Zehnder interference in a bipolar hybrid monolayer-bilayer graphene junction

Isolated and hybrid bilayer graphene rings

Strong anisotropic optical properties of 8-Pmmn borophene: a many-body perturbation study

Reversible ratchet effects in a narrow superconducting ring

Hematite at its thinnest limit

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

Three-dimensional electron-hole superfluidity in a superlattice close to room temperature

Tight-Binding Studio: A Technical Software Package to Find the Parameters of Tight-Binding Hamiltonian

All-strain based valley filter in graphene nanoribbons using snake states

Anomalous Dynamical Behavior of Freestanding Graphene Membranes

Bilayer SnS$_{2}$: Easy-tunable Stacking Sequence by Charging and Loading Pressure

Commensurability Effects in Viscosity of Nanoconfined Water

Comment on "Creating in-plane pseudomagnetic fields in excess of 1000 T by misoriented stacking in a graphene bilayer"

Computing Optical Properties of Ultra-thin Crystals

Electric field induced structural changes of water confined between two graphene layers

Large gap electron-hole superfluidity and shape resonances in coupled graphene nanoribbons

Mg(OH)2 -WS2 Heterobilayer: Electric Field Tunable Bandgap Crossover

Optical Properties of GaS-Ca(OH)$_2$ bilayer heterostructure

Piezoelectricity in asymmetrically strained bilayer graphene

Quantum Transport Characteristics of Lateral pn-Junction of Single Layer TiS3

Quantum transport in graphene Hall bars: Effects of vacancy disorder

Spin transitions driven by electric dipole spin resonance in fluorinated single- and bilayer-graphene quantum dots

Strain controlled valley filtering in multi-terminal graphene structures

Strain-induced topological phase transition in phosphorene and phosphorene nanoribbons

Strained graphene Hall bar

The split-operator technique for the study of spinorial wavepacket dynamics

Theoretical investigation of electron-hole complexes in anisotropic two-dimensional materials

Transport properties of bilayer graphene in a strong in-plane magnetic field

Valley filtering using electrostatic potentials in bilayer graphene

AA-stacked bilayer square ice between graphene layers?

Ag and Au Atoms Intercalated in Bilayer Heterostructures of Transition Metal Dichalcogenides and Graphene

Aharonov-Bohm oscillations in a mesoscopic ring with asymmetric arm-dependent injection

An efficient finite-difference scheme for computation of electron states in free-standing and core-shell quantum wires

Anisotropic exciton Stark shift in black phosphorus

Chain-like transitions in Wigner crystals: Sequential or non-sequential?

Distinct magnetic signatures of fractional vortex configurations in multiband superconductors

Electron polarization function and plasmons in metallic armchair graphene nanoribbons

Electronic and Magnetic Properties of 1T-TiSe2 Nanoribbons

Giant paramagnetic Meissner effect in multiband superconductors

GraphAne: From Synthesis to Applications

Graphene Ripples as a Realization of a Two-Dimensional Ising Model: A Scanning Tunneling Microscope Study

Hamiltonian of a many-electron system with single-electron and electron-pair states in a two-dimensional periodic potential

Hexagonal AlN: Dimensional-Crossover-Driven Bandgap Transition

High temperature conductance of a two-dimensional superlattice controlled by spin-orbit interaction

Magneto-optical transport properties of monolayer phosphorene

Oscillating spin-orbit interaction in two-dimensional superlattices: sharp transmission resonances and time-dependent spin polarized currents

Pentagonal Monolayer Crystals of Carbon, Boron Nitride, and Silver Azide

Spin- and valley-dependent transport through arrays of ferromagnetic silicene junctions

Spin-orbit interaction controlled properties of two-dimensional superlattices: Spintronic crystals

Stable Half-Metallic Monolayers of FeCl$_2$

Theory of anharmonic phonons in 2D crystals

Vacancy Formation and Oxidation Characteristics of Single Layer TiS3

Veselago lensing in graphene with a p-n junction: classical versus quantum effects

Vortical versus skyrmionic states in mesoscopic \emph{p}-wave superconductors

Disordered graphene Josephson junctions

Electron energy and temperature relaxation in graphene on a piezoelectric substrate

Electronic properties of graphene nano-flakes: Energy gap, permanent dipole, termination effect and Raman spectroscopy

Electronic structure of graphene hexagonal flake subjected to triaxial stress

Enhancement of electron-hole superfluidity in double few-layer graphene

Generic ordering of structural transitions in quasi-one-dimensional Wigner crystals

Graphene on boron-nitride: Moiré pattern in the van der Waals energy

Graphene on hexagonal lattice substrate: Stress and Pseudo-magnetic field

Interplay between snake and quantum edge states in a graphene Hall bar with a pn-junction

Magnetic particles confined in a modulated channel: structural transitions tunable by tilting a magnetic field

Membrane amplitude and triaxial stress in twisted bilayer graphene deciphered using first-principles directed elasticity theory and scanning tunneling microscopy

Multilayer graphene, Moiré patterns, grain boundaries and defects identified by scanning tunneling microscopy on the m-plane, non-polar surface of SiC

Non-linear response to electric field in extended Hubbard models

Oscillating spin-orbit interaction as a source of spin-polarized wave packets in two-terminal nanoscale devices

Plasmons and their interaction with electrons in trilayer graphene

Realization of free-standing silicene using bilayer graphene

Rippling, buckling and melting of single- and multi-layer MoS$_2$

Spin and valley polarization of plasmons in silicene due to external fields