Abstract: We have studied the temperature and field dependencies of the critical current I(C) in the Nb-Fe(0.1)Si(0.9)-Nb Josephson junction with a tunneling barrier formed by a paramagnetic insulator. We demonstrate that in these junctions coexistence of both the 0 and the pi states within one tunnel junction occurs, and leads to the appearance of a sharp cusp in the temperature dependence I(C)(T), similar to the I(C)(T) cusp found for the 0-pi transition in metallic pi junctions. This cusp is not related to the 0-pi temperature-induced transition itself, but is caused by the different temperature dependencies of the opposing 0 and pi supercurrents through the barrier.

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 address the nature of the disordered state that results from the adsorption of adatoms in graphene. For adatoms that sit at the center of the honeycomb plaquette, as in the case of most transition metals, we show that the ones that form a zero-energy resonant state lead to Anderson localization in the vicinity of the Dirac point. Among those, we show that there is a symmetry class of adatoms where Anderson localization is suppressed, leading to an exotic metallic state with large and rare charge droplets, that localizes only at the Dirac point. We identify the experimental conditions for the observation of the Anderson transition for adatoms in graphene.

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

Abstract: Typical Raman spectra of transition-metal dichalcogenides (TMDs) display two prominent peaks, E-2g and A(1g), that are well separated from each other. We find that these modes are degenerate in bulk WSe2 yielding one single Raman peak in contrast to other TMDs. As the dimensionality is lowered, the observed peak splits in two. In contrast, our ab initio calculations predict that the degeneracy is retained even for WSe2 monolayers. Interestingly, for minuscule biaxial strain, the degeneracy is preserved, but once the crystal symmetry is broken by a small uniaxial strain, the degeneracy is lifted. Our calculated phonon dispersion for uniaxially strained WSe2 shows a good 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 the interaction between WSe2 and the substrate. Furthermore, we find that WSe2 undergoes an indirect-to-direct band-gap transition from bulk to monolayers, which is ubiquitous for semiconducting TMDs. These results not only allow us to understand the vibrational and electronic properties of WSe2, but also point to effects of the interaction between the monolayer TMDs and the substrate on the vibrational and electronic properties. DOI: 10.1103/PhysRevB.87.165409

Abstract: Themean-field Hubbard model is used to investigate the formation of the antiferromagnetic phase in hexagonal graphene rings with inner zigzag edges. The outer edge of the ring was taken to be either zigzag or armchair, and we found that both types of structures can have a larger antiferromagnetic interaction as compared with hexagonal dots. This difference could be partially ascribed to the larger number of zigzag edges per unit area in rings than in dots. Furthermore, edge states localized on the inner ring edge are found to hybridize differently than the edge states of dots, which results in important differences in the magnetism of graphene rings and dots. The largest staggered magnetization is found when the outer edge has a zigzag shape. However, narrow rings with armchair outer edge are found to have larger staggered magnetization than zigzag hexagons. The edge defects are shown to have the least effect on magnetization when the outer ring edge is armchair shaped. DOI: 10.1103/PhysRevB.87.085434

Abstract: The technique of beam optics is applied to the description of the wave function of Dirac electrons. This approach is illustrated by considering electron transmission through simple nonhomogeneous structures, such as flat and bent p-n junctions and superlattices. We found that a convex p-n junction compresses the beam waist, while a concave interface widens it without loosing its focusing properties. At a flat p-n junction the waist of the transmitted Gaussian beam can be narrowed or widened, depending on the angle of incidence. A general condition is derived for the occurrence of beam collimation in a superlattice which is less stringent than previous discussed.

Abstract: The hydrogen interstitial and the substitutional AlZn, GaZn, and InZn are all shallow donors in ZnO and lead to n-type conductivity. Although shallow donors are expected to repel each other, we show by first-principles calculations that in ZnO these shallow donor impurities attract and form a complex, leading to a donor level deep in the band gap. This puts a limit on the n-type conductivity of (Al,Ga,In)-doped ZnO in the presence of hydrogen.

Abstract: Biexcitons in fractional dimensional spaces are studied using variational quantum Monte Carlo. We investigate the biexciton binding energy as a function of the electron-hole mass fraction sigma as well as study the dimensional dependence of biexcitons for sigma = 0 and sigma = 1. As our first application of this model we treat the H(2) molecule in two and three dimensions. Next we investigate biexcitons in carbon nanotubes within the fractional dimensional model. To this end we find a relation between the nanotube radius and the effective dimension. The results of both applications are compared with results obtained using different models and we find a reasonable agreement. Within the fractional dimensional model we find that the biexciton binding energy in carbon nanotubes accurately scales as E(B)(r,epsilon) = 1280 meV angstrom/(r epsilon), as a function of radius r and the dielectric screening epsilon.