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 single crystal of Cu74.08Al23.13Be2.79 undergoes a martensitic phase transition at 246 and 232 K under heating and cooling, respectively. The phase fronts between the austenite and martensite regions of the sample are weakly mobile with a power-law resonance under external stress fields. Surprisingly, the martensite phase is elastically much harder than the austenite phase showing that interfaces between various crystallographic variants are strongly pinned and cannot be moved by external stress while the phase boundary between the austenite and martensite regions in the sample remains mobile. This unusual behavior was studied by dynamical mechanical analysis (DMA) and resonant ultrasound spectroscopy. The remnant strain, storage modulus, and internal friction were recorded simultaneously for different applied forces in DMA. With increasing forces, the remnant strain increases monotonously while the internal friction peak height shows a minimum at 300 mN. Transmission electron microscopy shows that the pinning is generated by dislocations which are inherited from the austenite phase.

Abstract: Using scanning tunneling microscopy and Ginzburg-Landau simulations, we explore vortex configurations in magnetically coupled NbSe2/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-one-dimensional 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 nonlinear 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 antipinning channels.

Abstract: Anharmonic effects in an atomic monolayer thin crystal with honeycomb lattice structure are investigated by analytical and numerical lattice dynamical methods. Starting from a semiempirical model for anharmonic couplings of third and fourth orders, we study the in-plane and out-of-plane (flexural) mode components of the generalized wave vector dependent Gruneisen parameters, the thermal tension and the thermal expansion coefficients as a 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 a 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 T-alpha from negative to positive thermal expansion is lowered with smaller system size. Renormalization of the flexural mode has the opposite effect and leads to values of T-alpha approximate to 300 K for systems of macroscopic size. 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.

Abstract: Magnesium hydroxide [Mg(OH)(2)] has a layered brucitelike structure in its bulk form and was recently isolated as a new member of two-dimensional 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 WS2 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 separates 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 G(0)W(0) calculations, we show that the low-energy spectrum of the heterobilayer is dominated by the intralyer excitons of the WS2 monolayer. 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.

Abstract: Transition metal dichalcogenides (TMDs) are considered as promising candidates for next generation of electronic and optoelectronic devices. To make use of these materials, for instance in field effect transistor applications, it is mandatory to know the detailed properties of contacts of such TMDs with metal electrodes. Here, we investigate the role of the contact structure on the electronic and transport properties of metal-MoSe2 interfaces. Two different contact types, namely face and edge contacts, are studied. We consider both low (Sc) and high (Au) work function metals in order to thoroughly elucidate the role of the metal work function and the type of metal. First principles plane wave calculations and transport calculations based on nonequilibrium Green's function formalism reveal that the contact type has a large impact on the electronic and transport properties of metal-MoSe2 interfaces. For the Sc electrode, the Schottky barrier heights are around 0.25 eV for face contact and bigger than 0.6 eV for edge contact. For the Au case, we calculate very similar barrier heights for both contact types with an average value of 0.5 eV. Furthermore, while the face contact is found to be highly advantageous as compared to the edge contact for the Sc electrode, the latter contact becomes a better choice for the Au electrode. Our findings provide guidelines for the fabrication of TMD-based devices.

Abstract: The dissociation of 60degrees and screw dislocations in diamond is modeled in an approach combining isotropic elasticity theory with ab initio-based tight-binding total-energy calculations. Both dislocations are found to dissociate with a substantial lowering of their line energies. For the 60degrees dislocation, however, an energy barrier to dissociation is found. We investigate the core structure of a screw dislocation distinguishing “shuffle,” “mixed,” and “glide” cores. The latter is found to be the most stable undissociated screw dislocation. Further, the glide motion of 90degrees and 30degrees partials is discussed in terms of a process involving the thermal formation and subsequent migration of kinks along the dislocation line. The calculated activation barriers to dislocation motion show that the 30degrees partial is less mobile than the 90degrees partial. Finally, high-resolution electron microscopy is performed on high-temperature, high-pressure annealed natural brown diamond, allowing the core regions of 60degrees dislocations to be imaged. The majority of dislocations are found to be dissociated. However, in some cases, undissociated 60degrees dislocations were also observed.

Abstract: By using density functional theory (DFT) calculations, we predict a puckered, dynamically stable Janus single-layer black arsenic-phosphorus (b-AsP), which is composed of two different atomic sublayers, arsenic and phosphorus atoms. The calculated phonon spectrum reveals that Janus single-layer b-AsP is dynamically stable with either pure or coupled optical phonon branches arising from As and P atoms. The calculated Raman spectrum indicates that due to the relatively strong P-P bonds, As atoms have no contribution to the highfrequency optical vibrations. In addition, the orientation-dependent isovolume heat capacity reveals anisotropic contributions of LA and TA phonon branches to the low-temperature thermal properties. Unlike pristine single layers of b-As and b-P, Janus single-layer b-AsP exhibits additional out-of-plane asymmetry which leads to important consequences for its electronic, optical, and elastic properties. In contrast to single-layer b-As, Janus single-layer b-AsP is found to possess a direct band gap dominated by the P atoms. Moreover, real and imaginary parts of the dynamical dielectric function, including excitonic effects, reveal the highly anisotropic optical feature of the Janus single-layer. A tight-binding (TB) model is also presented for Janus single-layer b-AsP, and it is shown that, with up to seven nearest hoppings, the TB model reproduces well the DFT band structure in the low-energy region around the band gap. This TB model can be used in combination with the Green's function approach to study, e.g., quantum transport in finite systems based on Janus single-layer b-AsP. Furthermore, the linear-elastic properties of Janus single-layer b-AsP are investigated, and the orientation-dependent in-plane stiffness and Poisson ratio are calculated. It is found that the Janus single layer exhibits strong in-plane anisotropy in its Poisson ratio much larger than that of single-layer b-P. This Janus single layer is relevant for promising applications in optical dichroism and anisotropic nanoelasticity.

Abstract: The critical current of a thin superconducting strip of width W much larger than the Ginzburg-Landau coherence length xi but much smaller than the Pearl length Lambda = 2 lambda(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 degrees or 180 degrees 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 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.

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.