Abstract: Nitrogenated holey graphene (NHG) is a recently synthesized two-dimensional material. In this paper the structural and electronic properties of heterostructures of graphene and NHG are investigated using first-principles and tight-binding calculations. Due to the lattice mismatch between NHG and graphene, the formation of a moire pattern is preferred in the graphene/NHG heterostructure, instead of a lattice-coherent structure. In moire-patterned graphene/NHG, the band gap opening at the K point is negligible, and the linear band dispersion of graphene survives. Applying an electric field modifies the coupling strength between the two atomic layers. The Fermi velocity upsilon(F) is reduced as compared to the one of pristine graphene, and its magnitude depends on the twist angle theta between graphene and NHG: For theta = 0 degrees, upsilon(F) is 30% of that of graphene, and it increases rapidly to a value of 80% with increasing theta. The heterostructure exhibits electron-hole asymmetry in upsilon(F), which is large for small theta. In NHG encapsulated between two graphene layers, a “Dirac ring” appears around the K point. Its presence is robust with respect to the relative stacking of the two graphene layers. These findings can be useful for future applications of graphene/NHG heterostructures.

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: Recently, core-shell Ag-Co, Ag-Cu, and “onionlike” Cu-Co equilibrium configurations were predicted in the case of isolated face centered cubic (fcc) bimetallic clusters, and three shell onionlike configurations were predicted in the case of ternary metallic clusters with spherical and truncated octahedral morphologies. In the present paper, immiscible binary CuCo and ternary AgCuCo clusters with icosahedral structures are studied as functions of their size and composition. Clusters studied are formed by 13, 55, 147, 309, and 561 atoms corresponding to the five smallest possible closed shell icosahedral structures. An embedded atom model potential is used to describe their cohesion. Equilibrium configurations are investigated by means of Metropolis Monte Carlo free energy minimization in the (NPT) canonical ensemble. Most simulations are achieved at 10 and 300 K. The effect of temperature on segregation ordering is systematically investigated. Selected cases are used to identify the effect of size and composition on melting. In contrast with fcc clusters, homogeneous onionlike configurations of binary clusters are not predicted. When it is allowed by the composition, a complete outer shell is formed by Cu in binary Cu-Co clusters and by Ag in ternary Ag-Cu-Co clusters. Depending on temperature, Co may precipitate into decahedral groups under the Cu vertices of the icosahedra in binary clusters, while the Co-Cu configuration in ternary clusters drastically depends on the Ag coating. Despite the multicomponent character of the clusters and the immiscibility of the species forming them, for most compositions and sizes, equilibrium structures remain close to perfectly icosahedral at 10 K as well as at 300 K.

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.