Abstract: In order to study quantum effects in a two-dimensional crystal lattice of a free-standing monolayer graphene, we have performed both path-integral Monte Carlo (PIMC) and classical Monte Carlo (MC) simulations for temperatures up to 2000 K. The REBO potential is used for the interatomic interaction. The total energy, interatomic distance, root-mean-square displacement of the atom vibrations, and the free energy of the graphene layer are calculated. The obtained lattice vibrational energy per atom from the classical MC simulation is very close to the energy of a three-dimensional harmonic oscillator 3k(B)T. The PIMC simulation shows that quantum effects due to zero-point vibrations are significant for temperatures T < 1000 K. The quantum contribution to the lattice vibrational energy becomes larger than that of the classical lattice for T < 400 K. The lattice expansion due to the zero-point motion causes an increase of 0.53% in the lattice parameter. A minimum in the lattice parameter appears at T similar or equal to 500 K. Quantum effects on the atomic vibration amplitude of the graphene lattice and its free energy are investigated.

Abstract: The effect of the number of stacking layers and the type of stacking on the electronic and optical properties of bilayer and trilayer black phosphorus are investigated by using first-principles calculations within the framework of density functional theory. We find that inclusion of many-body effects (i.e., electron-electron and electron-hole interactions) modifies strongly both the electronic and optical properties of black phosphorus. While trilayer black phosphorus with a particular stacking type is found to be a metal by using semilocal functionals, it is predicted to have an electronic band gap of 0.82 eV when many-body effects are taken into account within the G(0)W(0) scheme. Though different stacking types result in similar energetics, the size of the band gap and the optical response of bilayer and trilayer phosphorene are very sensitive to the number of layers and the stacking type. Regardless of the number of layers and the type of stacking, bilayer and trilayer black phosphorus are direct band gap semiconductors whose band gaps vary within a range of 0.3 eV. Stacking arrangements that are different from the ground state structure in both bilayer and trilayer black phosphorus exhibit significant modified valence bands along the zigzag direction and result in larger hole effective masses. The optical gap of bilayer (trilayer) black phosphorus varies by 0.4 (0.6) eV when changing the stacking type. The calculated binding energy of the bound exciton hardly changes with the type of stacking and is found to be 0.44 (0.30) eV for bilayer (trilayer) phosphorous.

Abstract: Recent experimental studies showed that formation of a two-dimensional crystal structure of nitrogenated holey graphene (NHG) is possible. Similar to graphene, NHGs have an atomically thin and strong crystal structure. Using first-principles calculations, we investigate the structural, phononic, and thermal properties of monolayer NHG crystal. Our charge analysis reveals that the charged holey sites of NHG provide a reactive ground for further functionalization by adatoms or molecules. We also found that similar to graphene, the NHG structure has quite high-frequency phonon modes and the presence of nitrogen atoms leads to the emergence of additional vibrational modes. Our phonon analysis reveals the presence of three characteristic Raman-active modes of NHG. Furthermore, the analysis of constant-volume heat capacity showed that the NHG structure has a linear temperature dependence in the low-temperature region. The strong lattice structure and unique thermal properties of the NHG crystal structure are desirable in nanoscale device applications.

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: We examine theoretically the dynamics of two vortex shells in pinning-free superconducting thin disks in the Corbino geometry. In the first considered case, the inner shell is composed of vortices and the outer one of antivortices, corresponding to a state induced by the stray field of an off-plane magnetic dipole placed on top of the superconductor. In the second considered case, both shells comprise vortices induced by a homogeneous external field. We derive the equation of motion for each shell within the Bardeen-Stephen model and study the dynamics analytically by assuming both shells are rigid and commensurate. In both cases, two distinct regimes for vortex shell motion are identified: For low applied currents the entire configuration rotates rigidly, while above a threshold current the shells decouple from each other and rotate at different angular velocities. Analytical expressions for the decoupling current, the recombination time in the decoupled phases, as well as the voltage-current characteristics are presented. Our analytical results are in excellent agreement with numerical molecular dynamics simulations of the full many-vortex problem.

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.

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