Abstract: Using an effective one-dimensional cluster expansion in combination with first-principles electronic structure calculations we have studied the energetics and electronic properties of Sb-Te layered systems. For a Te concentration between 0 and 60 at. % an almost continuous series of metastable structures is obtained consisting of consecutive Sb bilayers next to consecutive Sb2Te3 units, with the general formula (Sb-2)(n)(Sb2Te3)(m) (n, m = 1,2, ... ). Between 60 and 100 at.% no stable structures are found. We account explicitly for the weak van derWaals bonding between Sb bilayers and Sb2Te3 units by using a recently developed functional, which strongly improves the interlayer bonding distances. At T = 0 K, no evidence is found for the existence of two separate single-phase regions delta and gamma and a two-phase region delta + gamma. Metastable compounds with a Te concentration between 0 and 40 at. % are semimetallic, whereas compounds with a Te concentration between 50 and 60 at. % are semiconducting. Compounds with an odd number of Sb layers are metallic and have a much higher formation energy than those with an even number of consecutive Sb layers, thereby favoring the formation of Sb bilayers.

Abstract: 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 bumps, and Gaussian depressions 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 wavelengths 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 the graphene's orientation, i. e., in contrast to a small-width (3 nm) nanoribbon with armchair edges, the one with zigzag 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).

Abstract: Using atomistic simulations, we determine the roughness and the thermal properties of a suspended graphane sheet. As compared to graphene, we found that (i) hydrogenated graphene has a larger thermal contraction, (ii) the roughness exponent at room temperature is smaller, i.e., ≃ 1.0 versus ≃ 1.2 for graphene, (iii) the wavelengths of the induced ripples in graphane cover a wide range corresponding to length scales in the range 30125 Å at room temperature, and (iv) the heat capacity of graphane is estimated to be 29.32±0.23 J/mol K, which is 14.8% larger than that for graphene, i.e., 24.98±0.14 J/mol K. Above 1500 K, we found that graphane buckles when its edges are supported in the x-y plane.

Abstract: Using first principles calculations we showed that the electronic and optical properties of single-layer black phosphorus (BP) depend strongly on the applied strain. Due to the strong anisotropic atomic structure of BP, its electronic conductivity and optical response are sensitive to the magnitude and the orientation of the applied strain. We found that the inclusion of many body effects is essential for the correct description of the electronic properties of monolayer BP; for example, while the electronic gap of strainless BP is found to be 0.90 eV by using semilocal functionals, it becomes 2.31 eV when many-body effects are taken into account within the G(0)W(0) scheme. Applied tensile strain was shown to significantly enhance electron transport along zigzag direction of BP. Furthermore, biaxial strain is able to tune the optical band gap of monolayer BP from 0.38 eV (at -8% strain) to 2.07 eV (at 5.5%). The exciton binding energy is also sensitive to the magnitude of the applied strain. It is found to be 0.40 eV for compressive biaxial strain of -8%, and it becomes 0.83 eV for tensile strain of 4%. Our calculations demonstrate that the optical response of BP can be significantly tuned using strain engineering which appears as a promising way to design novel photovoltaic devices that capture a broad range of solar spectrum.

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: Using multicomponent Ginzburg-Landau simulations, we show a plethora of vortex states possible in mesoscopic three-band superconductors. We find that mesoscopic confinement stabilizes chiral states, with nontrivial phase differences between the band condensates, as the ground state of the system. As a consequence, we report the broken-symmetry vortex states, the chiral states where vortex cores in different band condensates do not coincide (split-core vortices), as well as fractional-flux vortex states with broken time-reversal symmetry.

Abstract: 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 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.

Abstract: Using the tight-binding approach, we calculate the electronic structure of triangular and hexagonal MoS2 quantum dots. Due to the orbital asymmetry we show that it is possible to form quantum dots with the same shape but having different electronic properties. The electronic states of triangular and hexagonal quantum dots are explored, as well as the local and total density of states and the convergence towards the bulk spectrum with dot size is investigated. Our calculations show that: (1) edge states appear in the band gap, (2) that there are a larger number of electronic states in the conduction band as compared to the valence band, and (3) the relative number of edge states decreases with increasing dot size.

Abstract: Using the tight-binding formalism, we study spin and charge transport through a zigzag silicene ribbon subject to an external electric field E-z. The effect of an exchange field M-z is also taken into account and its consequences on the band structure as well as spin transport are evaluated. We show that the band structure lacks spin inversion symmetry in the presence of intrinsic spin-orbit interaction in combination of E-z and M-z fields. Our quantum transport calculations indicate that for certain energy ranges of the incoming electrons the silicene ribbon can act as a controllable high-efficiency spin polarizer. The polarization maxima occur simultaneously with the van Hove singularities of the local density of states. In this case, the combination of electric and exchange fields is the key to achieving nearly perfect spin polarization, which also leads to the appearance of additional narrow plateaus in the quantum conductance. Moreover, we demonstrate that the output current still remains completely spin-polarized for low-energy carriers even when a few edge vacancies are present.

Abstract: Using time-differential perturbed angular correlation spectroscopy we have measured the electric field gradient (EFG) at 111Cd probe nuclei in solid Ce in a pressure range up to 8 GPa. Covering various allotropic phases of Ce, we find that the value of the EFG in the cubic α phase is almost four times larger than in the cubic γ phase and close to values in the noncubic phases α′ and α″. These results together with the differences in time modulation of the spectra are interpreted as evidence for quadrupolar electronic charge-density ordering and symmetry lowering at the γ→α transition while the lattice remains face-centered cubic

Abstract: Water confined between two graphene layers with a separation of a few A forms a 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 the graphene layers which is nonpolar and nonferroelectric. We provide different energetic considerations and H-bonding results that explain the interlayer and intralayer properties of two-dimensional ice. The controversial AA stacking found experimentally [Algara-Siller et al., Nature (London) 519, 443 (2015)] is consistent with our minimum-energy crystal structure of bilayer ice. Furthermore, we predict that an odd number of layers of ice has the same lattice structure as monolayer ice, while an even number of ice layers exhibits the square ice AA stacking of bilayer ice.

Abstract: We analyze theoretically and experimentally vortex configurations in mesoscopic superconducting squares. Our theoretical approach is based on the analytical solution of the London equation using Green's-function method. The potential-energy landscape found for each vortex configuration is then used in Langevin-type molecular-dynamics simulations to obtain stable vortex configurations. Metastable states and transitions between them and the ground state are analyzed. We present our results of the first direct visualization of vortex patterns in micrometer-sized Nb squares, using the Bitter decoration technique. We show that the filling rules for vortices in squares with increasing applied magnetic field can be formulated, although in a different manner than in disks, in terms of formation of vortex “shells”.

Abstract: We calculate the equilibrium properties and the dynamic response of two vertically coupled circular quantum dots populated by particles of different electrical charge sign, i.e., electrons and holes. The equilibrium density profiles are obtained and used to compute the frequencies and oscillator strengths of magnetoplasma excitations. We find a strong coupling between the modes derived from the center-of-mass modes of the individual dots which leads to an anticrossing with a pronounced oscillator strength transfer from the “acoustic” to the “optical” branch. Also, due to the breaking of the generalized Kohn theorem a number of other than center-of-mass modes are excited whose oscillator strengths, however, are rather weak.

Abstract: We calculate the fluorescence of electron spins confined to a plane and driven into resonance by a magnetic field gradient and a constant magnetic field applied at right angles to each other. We solve the equation of motion of two-dimensional electrons in the magnetic field gradient to derive the dispersion curve of spin oscillators, the amplitude of electron oscillations, the effective magnetic field sensed by the electron spin, and the rate at which electrons are injected from an electrode into spin oscillators. We then switch on the interaction between the spin magnetic dipole and the electromagnetic field to find the fluorescence power radiated by the individual spin oscillators. The rate of radiative decay is first derived, followed by the probability of sequential photon emission whereby a series of spontaneous decays occurs at random times separated by intervals during which the spin performs Rabi oscillations. The quantum correlations between random radiative decays manifest as bursts of emission at regular intervals along the wire. We integrate all multiphoton processes to obtain an exact analytical expression for the radiated electromagnetic power. The present theory obtains all parameters of the problem including magnetodipole coupling, the particle dwell time in the magnetic field gradient, and the spin polarization of the incoming current. The output power contains a fine structure arising from the anharmonicity of electron oscillations and from nonlinear optical effects which both give satellite emission peaks at odd multiples of the fundamental frequency.

Abstract: We calculate the spin-polarized electronic transport through a molecular bilayer spin valve from first principles, and establish the link between the magnetoresistance and the spin-dependent interactions at the metal-molecule interfaces. The magnetoresistance of a Fe vertical bar bilayer-C-70 vertical bar Fe spin valve attains a high value of 70% in the linearresponse regime, but it drops sharply as a function of the applied bias. The current polarization has a value of 80% in linear response and also decreases as a function of bias. Both these trends can be modeled in terms of prominent spin-dependent Fe vertical bar C-70 interface states close to the Fermi level, unfolding the potential of spinterface science to control and optimize spin currents.