Abstract: Using the Ginzburg-Landau (GL) theory, we calculate the stability of sample symmetry-induced vortex-antivortex molecules in a mesoscopic superconducting bilayer exposed to a homogeneous magnetic field. We demonstrate the conditions under which the two condensates cooperatively broaden the field-temperature stability range of the composite (joint) vortex-antivortex state. In cases when such broadening is not achieved, a reentrance of the vortex-antivortex state is found at lower temperatures. In a large portion of the phase diagram noncomposite states are possible, in which the antivortex is present in only one of the layers. In this case, we demonstrate that the vortex-antivortex molecule in one of the layers can be pinned and enlarged by interaction with a vortex molecule in the other. Using analogies in the respective GL formalisms, we map our findings for the bilayer onto mesoscopic two-band superconductors.

Abstract: Recent advances in manufacturing arrays of artificial pinning sites, i.e., antidots, blind holes and magnetic dots, allowed an effective control of magnetic flux in superconductors. An array of magnetic bars deposited on top of a superconducting film was shown to display different pinning regimes depending on the direction of the in-plane magnetization of the bars. Changing the sign of their magnetization results in changes in the induced magnetic pinning potentials. By numerically solving the time-dependent Ginzburg-Landau equations in a superconducting film with periodic arrays of zigzag-arranged magnetic bars, we revealed various flux dynamics regimes. In particular, we demonstrate flux pinning and flux flow, depending on the direction of the magnetization of the magnetic bars. Remarkably, the revealed different flux-motion regimes are associated with different mechanisms of vortex-antivortex dynamics. For example, we found that for an 'antiparallel' configuration of magnetic bars this dynamics involves a repeating vortex-antivortex generation and annihilation. We show that the depinning transition and the onset of flux flow can be manipulated by the magnetization of the bars and the geometry of the array. This provides an effective control of the depinning critical current that can be useful for possible fluxonics applications.

Abstract: 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. DOI: 10.1103/PhysRevB.87.035131

Abstract: 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. DOI: 10.1103/PhysRevE.87.012307

Abstract: We propose a simple way of probing the number of modes contributing to the channeling in graphene waveguides which are formed by a gauge potential produced by mechanical strain. The energy mode structure for both homogeneous and non-homogeneous strain regimes is carefully studied using the continuum description of the Dirac equation. We found that high strain values privilege negative (instead of positive) group velocities throughout the guidance, sorting the types of modes flowing through it. We also show how the effect of a substrate-induced gap competes against the strain.

Abstract: The adsorption characteristics of alkali, alkaline-earth, and transition metal adatoms on silicene, a graphene-like monolayer structure of silicon are analyzed by means of first-principles calculations. In contrast to graphene, interaction between the metal atoms and the silicene surface is quite strong due to its highly reactive buckled hexagonal structure. In addition to structural properties, we also calculate the electronic band dispersion, net magnetic moment, charge transfer, work function, and dipole moment of the metal adsorbed silicene sheets. Alkali metals, Li, Na, and K, adsorb to hollow sites without any lattice distortion. As a consequence of the significant charge transfer from alkalis to silicene, metalization of silicene takes place. Trends directly related to atomic size, adsorption height, work function, and dipole moment of the silicene/alkali adatom system are also revealed. We found that the adsorption of alkaline-earth metals on silicene is entirely different from their adsorption on graphene. The adsorption of Be, Mg, and Ca turns silicene into a narrow gap semiconductor. Adsorption characteristics of eight transition metals Ti, V, Cr, Mn, Fe, Co, Mo, and W are also investigated. As a result of their partially occupied d orbital, transition metals show diverse structural, electronic, and magnetic properties. Upon the adsorption of transition metals, depending on the adatom type and atomic radius, the system can exhibit metal, half-metal, and semiconducting behavior. For all metal adsorbates, the direction of the charge transfer is from adsorbate to silicene, because of its high surface reactivity. Our results indicate that the reactive crystal structure of silicene provides a rich playground for functionalization at nanoscale. DOI: 10.1103/PhysRevB.87.085423

Abstract: Neutral excitons in strained axially symmetric In(Ga)As/GaAs quantum dots with a ringlike shape are investigated. Similar to experimental self-assembled quantum rings, the analyzed quantum dots have volcano-like shapes. The continuum mechanical model is employed to determine the strain distribution, and the single-band envelope function approach is adopted to compute the electron states. The hole states are determined by the axially symmetric multiband Luttinger-Kohn Hamiltonian, and the exciton states are obtained from an exact diagonalization. We found that the presence of the inner layer covering the ring opening enhances the excitonic Aharonov-Bohm (AB) oscillations. The reason is that the hole becomes mainly localized in the inner part of the quantum dot due to strain, whereas the electron resides mainly inside the ring-shaped rim. Interestingly, larger AB oscillations are found in the analyzed quantum dot than in a fully opened quantum ring of the same width. Comparison with the unstrained ringlike quantum dot shows that the amplitude of the excitonic Aharonov-Bohm oscillations are almost doubled in the presence of strain. The computed oscillations of the exciton energy levels are comparable in magnitude to the oscillations measured in recent experiments. DOI: 10.1103/PhysRevB.87.085314

Abstract: We study the mobility of Dirac fermions in monolayer graphene on a GaAs substrate, limited by the combined action of the extrinsic potential of piezoelectric surface acoustical phonons of GaAs (PA) and of the intrinsic deformation potential of acoustical phonons in graphene (DA). In the high-temperature (T) regime, the momentum relaxation rate exhibits the same linear dependence on T but different dependencies on the carrier density n, corresponding to the mobility mu proportional to 1 root n and 1/n, respectively for the PA and DA scattering mechanisms. In the low-T Bloch-Gruneisen regime, the mobility shows the same square-root density dependence mu proportional to root n, but different temperature dependencies mu proportional to T-3 and T-4, respectively for PA and DA phonon scattering. DOI: 10.1103/PhysRevB.87.075443

Abstract: Ab initio calculations within the density-functional theory formalism are performed to investigate the chemical functionalization of a graphene-like monolayer of siliconsilicenewith B, N, Al, or P atoms. The structural, electronic, magnetic, and vibrational properties are reported. The most preferable adsorption sites are found to be valley, bridge, valley and hill sites for B, N, Al, and P adatoms, respectively. All the relaxed systems with adsorbed/substituted atoms exhibit metallic behavior with strongly bonded B, N, Al, and P atoms accompanied by an appreciable electron transfer from silicene to the B, N, and P adatom/substituent. The Al atoms exhibit opposite charge transfer, with n-type doping of silicene and weaker bonding. The adatoms/substituents induce characteristic branches in the phonon spectrum of silicene, which can be probed by Raman measurements. Using molecular dynamics, we found that the systems under study are stable up to at least T=500 K. Our results demonstrate that silicene has a very reactive and functionalizable surface.

Abstract: Large-scale atomistic simulations using the reactive force field approach are implemented to investigate the thermomechanical properties of fluorinated graphene (FG). A set of parameters for the reactive force field potential optimized to reproduce key quantum mechanical properties of relevant carbon-fluorine cluster systems are presented. Molecular dynamics simulations are used to investigate the thermal rippling behavior of FG and its mechanical properties and compare them with graphene, graphane and a sheet of boron nitride. 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 unrippled system in contrast to the thermal rippling behavior of graphene. 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 armchair and zigzag directions, respectively. DOI: 10.1103/PhysRevB.87.104114

Abstract: 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 nonhomogeneous 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. DOI: 10.1103/PhysRevB.87.134509

Abstract: 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. DOI: 10.1103/PhysRevB.87.134510

Abstract: 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 Gruneisen parameter are estimated to be 25.2 J mol(-1) K-1, 7.2 x 10(-6) K-1, and 0.89, respectively.

Abstract: The electronic tunneling properties of the two stable forms of trilayer graphene (TLG), rhombohedral ABC and Bernal ABA, are examined for p-n and p-n-p 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 p-n-p 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.

Abstract: The tunneling current has been computed for a cylindrical nanowire tunneling field-effect transistor (TFET) with an all-round gate that covers the source region. Being the underlying mechanism, band-to-band tunneling, mediated by electron-phonon interaction, is pronouncedly affected by carrier confinement in the radial direction and, therefore, involves the self-consistent solution of the Schrodinger and Poisson equations. The latter has been accomplished by exploiting a non-linear variational principle within the framework of the modified local density approximation taking into account the nonparabolicity of both the valence band and conduction band in relatively thick wires. Moreover, while the effective-mass approximation might still provide a reasonable description of the conduction band in relatively thick wires, we have found that the nonparabolicity of the valence band needs to be included. As a major conclusion, it is observed that confinement effects in nanowire tunneling field-effect transistors have a stronger impact on the onset voltage of the tunneling current in comparison with planar TFETs. On the other hand, the value of the onset voltage is found to be overestimated when the valence band nonparabolicity is ignored. (C) 2013 AIP Publishing LLC.