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: The strain-induced pseudomagnetic field in supported graphene deposited on top of a nanostructured substrate is investigated by using atomistic simulations. A step, an elongated trench, a one-dimensional barrier, a spherical bubble, a Gaussian bump, and a Gaussian depression are considered as support structures for graphene. From the obtained optimum configurations we found very strong induced pseudomagnetic fields which can reach up to similar to 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 pseudomagnetic 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 sixfold symmetry of the induced pseudomagnetic field.

Abstract: 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/lambda)). 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 horizontally 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.

Abstract: 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, protective 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 thicknesses of 4 and 5 monolayers (MLs) 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.

Abstract: We apply neutron diffraction, high-resolution synchrotron x-ray diffraction, magnetization measurements, electronic structure calculations, and quantum Monte-Carlo simulations to unravel the structure and magnetism of (CuCl)LaTa2O7. Despite the pseudo-tetragonal crystallographic unit cell, this compound features an orthorhombic superstructure, similar to the Nb-containing (CuX)LaNb2O7 with X = Cl and Br. The spin lattice entails dimers formed by the antiferromagnetic fourth-neighbor coupling J(4), as well as a large number of nonequivalent interdimer couplings quantified by an effective exchange parameter J(eff). In (CuCl)LaTa2O7, the interdimer couplings are sufficiently strong to induce the long-range magnetic order with the Neel temperature T-N similar or equal to 7 K and the ordered magnetic moment of 0.53 mu(B), as measured with neutron diffraction. This magnetic behavior can be accounted for by J(eff)/J(4) similar or equal to 1.6 and J(4) similar or equal to 16 K. We further propose a general magnetic phase diagram for the (CuCl)LaNb2O7-type compounds, and explain the transition from the gapped spin-singlet (dimer) ground state in (CuCl)LaNb2O7 to the long-range antiferromagnetic order in (CuCl)LaTa2O7 and (CuBr)LaNb2O7 by an increase in the magnitude of the interdimer couplings J(eff)/J(4), with the (CuCl)LaM2O7 (M = Nb, Ta) compounds lying on different sides of the quantum critical point that separates the singlet and long-range-ordered magnetic ground states.

Abstract: 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 unrippled 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-eta)(eta similar or equal to 0.85) predicted by membrane theory. This unexpected behavior 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.

Abstract: 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 regimes. Self-consistent results for the pair correlation and superconducting current resemble those reported previously for single-layer graphene except at the Dirac point, where remarkable differences in the proximity effect are found, as well as a suppression of the superconducting current in the 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 the junction length is larger than the Fermi length.

Abstract: 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 standard GL theory. It is shown that in agreement with previous studies, this 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 this difference appears already in the leading correction to the standard 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 for recent materials of interest, which are in 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.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: The electrical transport properties of a MoGe thin film with a honeycomb array of nanoscale holes are investigated. The critical current of the system shows nonmatching anomalies as a function of applied magnetic field, enabling us to distinguish between multiquanta vortices trapped in the holes and interstitial vortices located between the holes. The number of vortices trapped in each hole is found to be larger than the saturation number predicted for an isolated hole and shows a nonlinear field dependence, leading to the caging effect as predicted from the Ginzburg-Landau (GL) theory. Our experimental results are supplemented by numerical simulations based on the GL theory.

Abstract: The scattering of a Gaussian wave packet in armchair and zigzag graphene edges is theoretically investigated by numerically solving the time-dependent Schrodinger equation for the tight-binding model Hamiltonian. Our theory allows us 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-known skipping orbits are observed. However, our results demonstrate that in the case of a pseudomagnetic field, induced by nonuniform strain, the scattering by an armchair edge results in a nonpropagating edge state.

Abstract: 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 (~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. The 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.

Abstract: We describe a semianalytic 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 nonmonotonic vortex-vortex interaction, which is temperature dependent, and can be further tuned by alterations of the material on the microscopic scale.

Abstract: We report extensive measurements on a new compound (Yb0.24Sn0.76)Ru that crystallizes in the cubic CsCl structure. Valence-band photoemission (PES) and L3 x-ray absorption show no divalent component in the 4f configuration of Yb. Inelastic neutron scattering (INS) indicates that the eight-fold degenerate J-multiplet of Yb3+ is split by the crystalline electric field (CEF) into a Γ7-doublet ground state and a Γ8 quartet at an excitation energy 20 meV. The magnetic susceptibility can be fit very well by this CEF scheme under the assumption that a Γ6-excited state resides at 32 meV; however, the Γ8/Γ6 transition expected at 12 meV was not observed in the INS. The resistivity follows a Bloch-Grüneisen law shunted by a parallel resistor, as is typical of systems subject to phonon scattering with no apparent magnetic scattering. All of these properties can be understood as representing simple local moment behavior of the trivalent Yb ion. At 1 K there is a peak in specific heat that is too broad to represent a magnetic-phase transition, consistent with absence of magnetic reflections in neutron diffraction. On the other hand this peak also is too narrow to represent the Kondo effect in the Γ7-doublet ground state. On the basis of the field dependence of the specific heat, we argue that antiferromagnetic (AF) short-range order (SRO) (possibly coexisting with Kondo physics) occurs at low temperatures. The long-range magnetic order is suppressed because the Yb site occupancy is below the percolation threshold for this disordered compound.

Abstract: We study the electronic structure of diluted F atoms chemisorbed on graphene using density functional theory calculations. We show that the nature of the chemical bonding of a F atom adsorbed on top of a C atom in graphene strongly depends on carrier doping. In neutral samples the F impurities induce a sp(3)-like bonding of the C atom below, generating a local distortion of the hexagonal lattice. As the graphene is electron-doped, the C atom retracts back to the graphene plane and for high doping (10(14) cm(-2)) its electronic structure corresponds to a nearly pure sp(2) configuration. We interpret this sp(3)-sp(2) doping-induced crossover in terms of a simple tight-binding model and discuss the physical consequences of this change.