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: 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: Janus particles present an important class of building blocks for directional assembly. These are compartmentalized colloids with two different hemispheres. In this work we consider a three-dimensional model of Janus spheres that contain one hydrophobic and one charged hemisphere. Using molecular dynamics simulations, we study the morphology of these particles when confined in a channel-like environment. The interplay between the attractive and repulsive forces on each particle gives rise to a rich phase space where the relative orientation of each particle plays a dominant role in the formation of large-scale clusters. The interest in this system is primarily due to the fact that it could give a better understanding of the mechanisms of the formation of polar membranes. A variety of ordered membranelike morphologies is found consisting of single and multiple connected chain configurations. The helicity of these chains can be chosen by simply changing the salt concentration of the solution. Special attention is given to the formation of Bernal spirals. These helices are composed of regular tetrahedra and are known to exhibit nontrivial translational and rotational symmetry.

Abstract: We study ballistic transport of Dirac fermions in silicene through arrays of barriers, of width d, in the presence of an exchange field M and a tunable potential of height U or depth-U. The spin-and valley-resolved conductances as functions of U or M, exhibit resonances away from the Dirac point (DP) and close to it a pronounced dip that becomes a gap when a critical electric field E-z is applied. This gap widens by increasing the number of barriers and can be used to realize electric field-controlled switching of the current. The spin p(s) and valley p(v) polarizations of the current near the DP increase with Ez or M and can reach 100% for certain of their values. These field ranges widen significantly by increasing the number of barriers. Also, ps and pv oscillate nearly periodically with the separation between barriers or wells and can be inverted by reversing M.

Abstract: The electronic properties, carrier mobility, and strain response of TiS3 nanoribbons (TiS3 NRs) are investigated by first-principles calculations. We found that the electronic properties of TiS3 NRs strongly depend on the edge type (a or b). All a-TiS3 NRs are metallic with a magnetic ground state, while b-TiS3 NRs are direct band gap semiconductors. Interestingly, the size of the band gap and the band edge position are almost independent of the ribbon width. This feature promises a constant band gap in a b-TiS3 NR with rough edges, where the ribbon width differs in different regions. The maximum carrier mobility of b-TiS3 NRs is calculated by using the deformation potential theory combined with the effective mass approximation and is found to be of the order 10(3) cm(2) V-1 s(-1). The hole mobility of the b-TiS3 NRs is one order of magnitude lower, but it is enhanced compared to the monolayer case due to the reduction in hole effective mass. The band gap and the band edge position of b-TiS3 NRs are quite sensitive to applied strain. In addition we investigate the termination of ribbon edges by hydrogen atoms. Upon edge passivation, the metallic and magnetic features of a-TiS3 NRs remain unchanged, while the band gap of b-TiS3 NRs is increased significantly. The robust metallic and ferromagnetic nature of a-TiS3 NRs is an essential feature for spintronic device applications. The direct, width-independent, and strain-tunable band gap, as well as the high carrier mobility, of b-TiS3 NRs is of potential importance in many fields of nanoelectronics, such as field-effect devices, optoelectronic applications, and strain sensors.

Abstract: We numerically study the vortex-vortex interaction in multicomponent homogeneous Bose-Einstein condensates within the realm of the Gross-Pitaevskii theory. We provide strong evidence that pairwise vortex interaction captures the underlying mechanisms which determine the geometric configuration of the vortices, such as different lattices in many-vortex states, as well as the bound vortex states with two (dimer) or three (trimer) vortices. Specifically, we discuss and apply our theoretical approach to investigate intra- and intercomponent vortex-vortex interactions in two- and three-component Bose-Einstein condensates, thereby shedding light on the formation of the exotic vortex configurations. These results correlate with current experimental efforts in multicomponent Bose-Einstein condensates and the understanding of the role of vortex interactions in multiband superconductors.

Abstract: The control of side gates on the quantum electron transport in narrow graphene ribbons of different widths and edge types (armchair and zigzag) is investigated. The conductance exhibits Fano resonances with varying side gate potential. Resonant and antiresonant peaks in the conductance can be associated with the eigenstates of a closed system, and these peaks can be accurately fitted with a Fano line shape. The local density of states (LDOS) and the electron current show a specific behavior at these resonances, which depends on the ribbon edge type. In zigzag ribbons, transport is dominated by intervalley scattering, which is reflected in the transmission functions of individual modes. The side gates induce p-n interfaces near the edges at which the LDOS exhibits peaks. Near the resonance points, the electron current flows uniformly through the constriction, while near the antiresonances it creates vortices. In the armchair ribbons the LDOS spreads in areas of high potential, with current flowing near the edges.

Abstract: Ripples in pristine freestanding graphene naturally orient themselves in an array that is alternately curved-up and curved-down; maintaining an average height of zero. Using scanning tunneling microscopy (STM) to apply a local force, the graphene sheet will reversibly rise and fall in height until the height reaches 60%-70% of its maximum at which point a sudden, permanent jump occurs. We successfully model the ripples as a spin-half Ising magnetic system, where the height of the graphene plays the role of the spin. The permanent jump in height, controlled by the tunneling current, is found to be equivalent to an antiferromagnetic-to-ferromagnetic phase transition. The thermal load underneath the STM tip alters the local tension and is identified as the responsible mechanism for the phase transition. Four universal critical exponents are measured from our STM data, and the model provides insight into the statistical role of graphene's unusual negative thermal expansion coefficient.

Abstract: Motivated by a recent experiment that reported the successful synthesis of hexagonal (h) AlN [Tsipas et al., Appl. Phys. Lett. 103, 251605 (2013)], we investigate structural, electronic, and vibrational properties of bulk, bilayer, and monolayer structures of h-AlN by using first-principles calculations. We show that the hexagonal phase of the bulk h-AlN is a stable direct-band-gap semiconductor. The calculated phonon spectrum displays a rigid-layer shear mode at 274 cm(-1) and an E-g mode at 703 cm(-1), which are observable by Raman measurements. In addition, single-layer h-AlN is an indirect-band-gap semiconductor with a nonmagnetic ground state. For the bilayer structure, AA'-type stacking is found to be the most favorable one, and interlayer interaction is strong. While N-layered h-AlN is an indirect-band-gap semiconductor for N = 1 – 9, we predict that thicker structures (N >= 10) have a direct band gap at the Gamma point. The number-of-layer-dependent band-gap transitions in h-AlN is interesting in that it is significantly different from the indirect-to-direct crossover obtained in the transition-metal dichalcogenides.

Abstract: We numerically investigate the mixed pairing symmetries inmesoscopic superconducting loops in the presence of spin correlations by solving the Bogoliubov-de Gennes equations self-consistently. The spatial variations of the superconducting order parameters and the spontaneous magnetization are determined by the band structure. When the threaded magnetic flux turns on, the charge and spin currents both emerge and depict periodic evolution. In the case of a mesoscopic loop with dominant triplet p(x) +/- ip(y)-wave symmetry, a slight change of the chemical potential may lead to novel flux-dependent evolution patterns of the ground-state energy and the magnetization. The spin-polarized currents show pronounced quantum oscillations with fractional periods due to the appearance of energy jumps in flux, accompanied with a steplike feature of the enhanced spin current. Particularly, at some appropriate flux, the peaks of the zero-energy local density of states clearly indicate the occurrence of the odd-frequency pairing. In the case of a superconducting loop with dominant singlet d(x2-y2)-wave symmetry, the spatial profiles of the zero-energy local density of states and the magnetization show spin-dependent features on different sample diagonals. Moreover, the evolution of the flux-induced spin current always exhibits an hc/e periodicity.

Abstract: Ca(OH)(2) crystals, well known as portlandite, are grown in layered form, and we found that they can be exfoliated on different substrates. We performed first principles calculations to investigate the structural, electronic, vibrational, and mechanical properties of bulk, bilayer, and monolayer structures of this material. Different from other lamellar structures such as graphite and transition-metal dichalcogenides, intralayer bonding in Ca(OH)(2) is mainly ionic, while the interlayer interaction remains a weak dispersion-type force. Unlike well-known transition-metal dichalcogenides that exhibit an indirect-to-direct band gap crossover when going from bulk to a single layer, Ca(OH)(2) is a direct band gap semiconductor independent of the number layers. The in-plane Young's modulus and the in-plane shear modulus of monolayer Ca(OH)(2) are predicted to be quite low while the in-plane Poisson ratio is larger in comparison to those in the monolayer of ionic crystal BN. We measured the Raman spectrum of bulk Ca(OH)(2) and identified the high-frequency OH stretching mode A(1g) at 3620 cm(-1). In this study, bilayer and monolayer portlandite [Ca(OH)(2)] are predicted to be stable and their characteristics are analyzed in detail. Our results can guide further research on ultrathin hydroxites.