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: Plasmon excitations in metallic armchair graphene nanoribbons are investigated using the random phase approximation. An exact analytical expression for the polarization function of Dirac fermions is obtained, valid for arbitrary temperature and doping. We find that at finite temperatures, due to the phase space redistribution among inter-band and intra-band electronic transitions in the conduction and valence bands, the full polarization function becomes independent of temperature and position of the chemical potential. It is shown that for a given width of nanoribbon there exists a single plasmon mode whose energy dispersion is determined by the graphene's fine structure constant. In the case of two Coulomb-coupled nanoribbons, this plasmon splits into in-phase and out-of-phase plasmon modes with splitting energy determined by the inter-ribbon spacing.

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 propose that low-strain Si/O superlattices can be constructed by connecting reconstructed Si{001} surfaces by Si-O-Si bridges. Ab initio calculations show that our models are energetically more favorable than all the models proposed so far. The part of our Si/O superlattice model is experimentally accessible just by oxidizing a Si( 001) substrate. To complete our Si/O superlattice model, we propose a three-step method. We also explore the potential of our Si/O superlattice models for new materials used in future Si electronics. We find that the location of the channel where the carriers travel can be controlled between the interfaces and the Si layers by the insertion of O atoms into the Si-Si dimers. By revealing the origins of the interface electron and hole states, we find that similar interface states should be easily achieved for Si slabs and Si substrates. Interestingly, the interface electrons and holes have small effective masses in the direction parallel to the channel and large effective masses in the direction normal to the channel, which makes the Si/O superlattices attractive to be used for channel materials. We also find that the valley splitting of Si is enhanced by the formation of the Si/O/Si interfaces, which is ideal for developing Si-based qubits. Our findings open new perspectives to design and control the electronic properties of Si.

Abstract: A tight-binding approach based on the Chebyshev-Bogoliubov-de Gennes method is used to describe disordered single-layer graphene Josephson junctions. Scattering by vacancies, ripples, or charged impurities is included. We compute the Josephson current and investigate the nature of multiple Andreev reflections, which induce bound states appearing as peaks in the density of states for energies below the superconducting gap. In the presence of single-atom vacancies, we observe a strong suppression of the supercurrent, which is a consequence of strong intervalley scattering. Although lattice deformations should not induce intervalley scattering, we find that the supercurrent is still suppressed, which is due to the presence of pseudomagnetic barriers. For charged impurities, we consider two cases depending on whether the average doping is zero, i.e., existence of electron-hole puddles, or finite. In both cases, short-range impurities strongly affect the supercurrent, similar to the vacancies scenario.

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.

Abstract: Large-scale atomistic simulations using the reactive empirical bond order force field approach is implemented to investigate thermal and mechanical properties of single-layer (SL) and multilayer (ML) molybdenum disulfide (MoS2). The amplitude of the intrinsic ripples of SL MoS2 are found to be smaller than those exhibited by graphene (GE). Furthermore, because of the van der Waals interaction between layers, the out-of-plane thermal fluctuations of ML MoS2 decreases rapidly with increasing number of layers. This trend is confirmed by the buckling transition due to uniaxial stress which occurs for a significantly larger applied tension as compared to graphene. For SL MoS2, the melting temperature is estimated to be 3700 K which occurs through dimerization followed by the formation of small molecules consisting of two to five atoms. When different types of vacancies are inserted in the SL MoS2 it results in a decrease of both the melting temperature as well as the stiffness.

Abstract: The Tomasch effect (TE) is due to quasiparticle interference (QPI) as induced by a nonuniform superconducting order parameter, which results in oscillations in the density of states (DOS) at energies above the superconducting gap. Quantum confinement in nanoscale superconductors leads to an inhomogenerous distribution of the Cooperpair condensate, which, as we found, triggers the manifestation of a new TE. We investigate the electronic structure of nanoscale superconductors by solving the Bogoliubov-de Gennes (BdG) equations self-consistently and describe the TE determined by two types of processes, involving two-or three-subband QPIs. Both types of QPIs result in additional BCS-like Bogoliubov-quasiparticles and BCS-like energy gaps leading to oscillations in the DOS and modulated wave patterns in the local density of states. These effects are strongly related to the symmetries of the system. A reduced 4 x 4 inter-subband BdG Hamiltonian is established in order to describe analytically the TE of two-subband QPIs. Our study is relevant to nanoscale superconductors, either nanowires or thin films, Bose-Einsten condensates, and confined systems such as two-dimensional electron gas interface superconductivity.

Abstract: We investigate mesoscopic Josephson-junction arrays created by patterning superconducting disks on monolayer graphene, concentrating on the high-T/T-c regime of these devices and the phenomena which contribute to the superconducting glass state in diffusive arrays. We observe features in the magnetoconductance at rational fractions of flux quanta per array unit cell, which we attribute to the formation of flux-quantized vortices. The applied fields at which the features occur are well described by Ginzburg-Landau simulations that take into account the number of unit cells in the array. We find that the mean conductance and universal conductance fluctuations are both enhanced below the critical temperature and field of the superconductor, with greater enhancement away from the graphene Dirac point.

Abstract: Using first-principles calculations and reactive force field molecular dynamics simulations, we study the structural properties and dynamics of a fluorine (F) atom, either adsorbed on the surface of single layer graphene (F/GE) or between the layers of AB stacked bilayer graphene (F@ bilayer graphene). It is found that the diffusion of the F atom is very different in those cases, and that the mobility of the F atom increases by about an order of magnitude when inserted between two graphene layers. The obtained diffusion constant for F/GE is twice larger than that experimentally found for gold adatom and theoretically found for C-60 molecule on graphene. Our study provides important physical insights into the dynamics of fluorine atoms between and on graphene layers and explains the mechanism behind the separation of graphite layers due to intercalation of F atoms.