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: The structural transitions of the ground state of a system of repulsively interacting particles confined in a quasi-one-dimensional channel, and the effect of the interparticle interaction as well as the functional form of the confinement potential on those transitions are investigated. Although the nonsequential ordering of transitions (non-SOT), i.e., the 1 – 2 – 4 – 3 – 4 – 5 – 6 – ... sequence of chain configurations with increasing density, is widely robust as predicted in a number of theoretical studies, the sequential ordering of transitions (SOT), i.e., the 1 – 2 – 3 – 4 – 5 – 6 – ... chain, is found as the ground state for long-ranged interparticle interaction and hard-wall-like confinement potentials. We found an energy barrier between every two different phases around its transition point, which plays an important role in the preference of the system to follow either a SOT or a non-SOT. However, that preferential transition requires also the stability of the phases during the transition. Additionally, we analyze the effect of a small structural disorder on the transition between the two phases around its transition point. Our results show that a small deformation of the triangular structure changes dramatically the picture of the transition between two phases, removing in a considerable region the non-SOT in the system. This feature could explain the fact that the non-SOT is, up to now, not observed in experimental systems, and suggests a more advanced experimental setup to detect the non-SOT.

Abstract: Simulations are carried out based on the dynamical mean-field theory (DMFT) in order to investigate the properties of correlated thin films for various values of the chemical potential, temperature, interaction strength, and applied transverse electric field. Application of a sufficiently strong field to a thin film at half filling leads to the appearance of conducting regions near the surfaces of the film, whereas in doped slabs the application of a field leads to a conductivity enhancement on one side of the film and a gradual transition to the insulating state on the opposite side. In addition to the inhomogeneous DMFT, a local density approximation (LDA) is considered in which the particle density n, quasiparticle residue Z, and spectral weight at the Fermi level A(ω=0) of each layer are approximated by a homogeneous bulk environment. A systematic comparison between the two approaches reveals that the less expensive LDA results are in good agreement with the DMFT approach, except close to the metal-to-insulator transition points and in the layers immediately at the film surfaces. LDA values for n are overall more reliable than those for Z and A(ω=0). The hysteretic behavior (memory effect) characteristic of the bulk doping driven Mott transition persists in the slab.

Abstract: Within a tight binding approach we investigate the energy levels of hexagonal and triangular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We study AA- and AB-(Bernal) stacked BLG QDs and obtain the energy levels in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). Our results show that the size dependence of the energy levels is different from that of monolayer graphene QDs. The energy spectrum of AB-stacked BLG QDs with zigzag edges exhibits edge states which spread out into the opened energy gap in the presence of a perpendicular electric field. We found that the behavior of these edges states is different for the hexagonal and triangular geometries. In the case of AA-stacked BLG QDs, the electron and hole energy levels cross each other in both cases of armchair and zigzag edges as the dot size or the applied bias increases.

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: The structural properties of partially hydrogenated and fluorinated graphene with different percentages of H/F atoms are investigated using molecular dynamics simulations based on reactive force field (ReaxFF) potentials. We found that the roughness of graphene varies with the percentage (p) of H or F and in both cases is maximal around p = 50%. Similar results were obtained for partially oxidized graphene. The two-dimensional area size of partially fluorinated and hydrogenated graphene exhibits a local minimum around p = 35% coverage. The lattice thermal contraction in partially functionalized graphene is found to be one order of magnitude larger than that of fully covered graphene. We also show that the armchair structure for graphene oxide (similar to the structure of fully hydrogenated and fluorinated graphene) is unstable. Our results show that the structure of partially functionalized graphene changes nontrivially with the C : H and C : F ratio as well as with temperature.

Abstract: We study transport properties of a MoS2 monolayer in the presence of a perpendicular magnetic field B. We derive and discuss its band structure and take into account spin and valley Zeeman effects. Compared to a conventional two-dimensional electron gas, these effects lead to new quantum Hall plateaus and new peaks in the longitudinal resistivity as functions of the magnetic field. The field B leads to a significant enhancement of the spin splitting in the conduction band, to a beating of the Shubnikov-de Haas (SdH) oscillations in the low-field regime, and to their splitting in the high-field regime. The Zeeman fields suppress significantly the beating of the SdH oscillations in the low-field regime and strongly enhance their splitting at high fields. The spin and valley polarizations show a similar beating pattern at low fields and are clearly separated at high fields in which they attain a value higher than 90%.

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: Recently discovered series of high Tc superconductors, characterized by the existence of two types of cations within the same layer, are presented. The first family concerns the mercury based cuprates, Hg(1-x)M(x)A(2)Ca(m-1)Cu(m)O(2m+2+delta), with A = Ba and/or Sr, which exhibit structures closely related to that of the thallium cuprates TlBa2Cam-1CumO2m+3. They differ from the thallium cuprates by a high oxygen deficiency at the level of the mercury layer. It is shown that cations such as M = Cu, Pb, Tl, Bi, Ce, Pr, Cr, V, Mo, W, Ti, Sr, Ca,... can partially substitute for mercury ions, stabilizing the structures. The cationic composition of the layer depends indeed on the nature of the M cation but also on that of the alkaline earth A. For given A and M cations, the a: value remains unchanged even when the number of copper layers varies. M and Hg cations are either statistically distributed over the same site or ordered. Different types of ordering have been detected. Another way of generating mixed layers is to shear periodically the structure, leading to the formation of the so called ''collapsed phase''. In the collapsed bismuth cuprates, bismuth and copper segments, a few octahedra long, alternate in strongly waving layers. In the collapsed oxycarbonates, carbonate groups and M cations are ordered within the intermediate layer so that they can be simply described from a partial and ordered substitution of carbon for Hg,TI, Bi and other M cations building the intermediate layer. The oxycarbonitrates (Y1-xCax)(n)Ba2nCu3n-1(C,N)O3O7n-3 can also be described as an ordered substitution of carbon for copper in the 123 matrix. The different families of superconducting materials which are generated by such mechanisms are described as well as the way the different species are distributed within the mixed layers. Their influence on the physical properties are discussed.