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 effect of an applied electric field and the effect of charging are investigated on themagnetic anisotropy (MA) of various stable two-dimensional (2D) crystals such as graphene, FeCl2, graphone, fluorographene, and MoTe2 using first-principles calculations. We found that themagnetocrystalline anisotropy energy of Co-on-graphene and Os-doped-MoTe2 systems change linearly with electric field, opening the possibility of electric field tuningMAof these compounds. In addition, charging can rotate the easy-axis direction ofCo-on-graphene andOs-doped-MoTe2 systems from the out-of-plane (in-plane) to in-plane (out-of-plane) direction. The tunable MA of the studied materials is crucial for nanoscale electronic technologies such as data storage and spintronics devices. Our results show that controlling the MA of the mentioned 2D crystal structures can be realized in various ways, and this can lead to the emergence of a wide range of potential applications where the tuning and switching of magnetic functionalities are important.

Abstract: Employing density functional theory-based methods, we investigate monolayer and bilayer structures of hexagonal SnS2, which is a recently synthesized monolayer metal dichalcogenide. Comparison of the 1H and 1T phases of monolayer SnS2 confirms the ground state to be the 1T phase. In its bilayer structure we examine different stacking configurations of the two layers. It is found that the interlayer coupling in bilayer SnS2 is weaker than that of typical transition-metal dichalcogenides so that alternative stacking orders have similar structural parameters and they are separated with low energy barriers. A possible signature of the stacking order in the SnS2 bilayer has been sought in the calculated absorbance and reflectivity spectra. We also study the effects of the external electric field, charging, and loading pressure on the characteristic properties of bilayer SnS2. It is found that (i) the electric field increases the coupling between the layers at its preferred stacking order, so the barrier height increases, (ii) the bang gap value can be tuned by the external E field and under sufficient E field, the bilayer SnS2 can become a semimetal, (iii) the most favorable stacking order can be switched by charging, and (iv) a loading pressure exceeding 3 GPa changes the stacking order. The E-field tunable band gap and easily tunable stacking sequence of SnS2 layers make this 2D crystal structure a good candidate for field effect transistor and nanoscale lubricant applications.

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: The structural and mechanical properties of N-doped graphene (NG) are investigated using reactive force field (ReaxFF) potentials in large-scale molecular dynamics simulations. We found that ripples, which are induced by the dopants, change the roughness of NG, which depends on the number of dopants and their local arrangement. For any doping ratio N/C, the NG becomes ferroelectric with a net dipole moment. The formation energy increases nonlinearly with N/C ratio, while the Young's modulus, tensile strength, and intrinsic strain decrease with the number of dopants. Our results for the structural deformation and the thermoelectricity of the NG sheet are in good agreement with recent experiments and ab initio calculations.

Abstract: Finding novel atomically thin heterostructures and understanding their characteristic properties are critical for developing better nanoscale optoelectronic devices. In this study, we investigate the electronic and optical properties of a GaS-Ca(OH)(2) heterostructure using first-principle calculations. The band gap of the GaS-Ca(OH)(2) heterostructure is significantly reduced when compared to those of the isolated constituent layers. Our calculations showthat the GaS-Ca(OH)(2) heterostructure is a type-II heterojunction which can be used to separate photoinduced charge carriers where electrons are localized in GaS and holes in the Ca(OH)(2) layer. This leads to spatially indirect excitons which are important for solar energy and optoelectronic applications due to their long lifetime. By solving the Bethe-Salpeter equation on top of a single shot GW calculation (G(0)W(0)), the dielectric function and optical oscillator strength of the constituent monolayers and the heterostructure are obtained. The oscillator strength of the optical transition for the GaS monolayer is an order of magnitude larger than the Ca(OH)(2) monolayer. We also found that the calculated optical spectra of different stacking types of the heterostructure show dissimilarities, although their electronic structures are rather similar. This prediction can be used to determine the stacking type of ultrathin heterostructures.

Abstract: We investigate silicene superlattices in the presence of a tunable barrier potential U, an exchange field M, and a perpendicular electric field E-z. The resulting miniband structure depends on the spin and valley indices and on the fields M and E-z. These fields determine the minigaps and also affect the additional Dirac points brought about by the periodic potential U. In addition, we consider diffusive transport and assess its dependence on the spin and valley indices as well as on temperature. The corresponding spin and valley polarizations strongly depend on the potential U and can be made almost 100% at very low temperatures at particular values of the Fermi energy.

Abstract: A strong in-plane magnetic field drastically alters the low-energy spectrum of bilayer graphene by separating the parabolic energy dispersion into two linear Dirac cones. The effect of this dramatic change on the transport properties strongly depends on the orientation of the in-plane magnetic field with respect to the propagation direction of the charge carriers and the angle at which they impinge on the electrostatic potentials. For magnetic fields oriented parallel to the potential boundaries an additional propagating mode that results from the splitting into Dirac cones enhances the transmission probability for charge carriers tunneling through the potentials and increases the corresponding conductance. Our results show that the chiral suppression of transmission at normal incidence, reminiscent of bilayer graphene's 2 pi Berry phase, is turned into a chiral enhancement when the magnetic field increases, thus indicating a transition from a bilayer to a monolayer-like system at normal incidence. Further, we find that the typical transmission resonances stemming from confinement in a potential barrier are shifted to higher energy and are eventually transformed into antiresonances with increasing magnetic field. For magnetic fields oriented perpendicular to the potential boundaries we find a very pronounced transition from a bilayer system to two separated monolayer-like systems with Klein tunneling emerging at certain incident angles symmetric around 0, which also leaves a signature in the conductance. For both orientations of the magnetic field, the transmission probability is still correctly described by pseudospin conservation. Finally, to motivate the large in-plane magnetic field, we show that its energy spectrum can be mimicked by specific lattice deformations such as a relative shift of one of the layers. With this equivalence we introduce the notion of an in-plane pseudomagnetic field.

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: The dynamics of driven pancake vortices in layered superconductors is studied using molecular-dynamics simulations. We found that, with increasing driving force, for strong interlayer coupling, the preexisted vortex lines either directly depin or first transform to two-dimensional (2D) pinned states before they are depinned, depending on the pinning strength. In a narrow region of pinning strengths, we found an interesting repinning process, which results in a negative differential resistance. For weak interlayer coupling, individually pinned pancake vortices first form disordered 2D flow and then transform to ordered three-dimensional (3D) flow with increasing driving force. However, for extremely strong pinning, the random pinning-induced thermal-like Langevin forces melt 3D vortex lines, which results in a persistent 2D flow in the fast-sliding regime. In the intermediate regime, the peak effect is found: With increasing driving force, the moving pancake vortices first crystallize to moving 3D vortex lines, and then these 3D vortex lines are melted, leading to the appearance of a reentrant 2D flow state. Our results are summarized in a dynamical phase diagram.

Abstract: Using the continuum model, we investigate the confined states and the corresponding wave functions of ABC-stacked trilayer graphene (TLG) quantum dots (QDs). First, a general infinite-mass boundary condition is derived and applied to calculate the electron and hole energy levels of a circular QD in both the absence and presence of a perpendicular magnetic field. Our analytical results for the energy spectra agree with those obtained by using the tight-binding model, where a TLG QD is surrounded by a staggered potential. Our findings show that (i) the energy spectrum exhibits intervalley symmetry E-K(e)(m) = -E-K'(h)(m) for the electron (e) and hole (h) states, where m is the angular momentum quantum number, (ii) the zero-energy Landau level (LL) is formed by the magnetic states with m <= 0 for both Dirac valleys, that is different from monolayer and bilayer graphene QD with infinite-mass potential in which only one of the cones contributes, and (iii) groups of three quantum Hall edge states in the tight-binding magnetic spectrum approach the zero LL, which results from the layer symmetry in TLG QDs.

Abstract: The mechanical properties of monolayer GaS and GaSe crystals are investigated in terms of their elastic constants: in-plane stiffness (C), Poisson ratio (nu), and ultimate strength (sigma(U)) by means of first-principles calculations. The calculated elastic constants are compared with those of graphene and monolayer MoS2. Our results indicate that monolayer GaS is a stiffer material than monolayer GaSe crystals due to the more ionic character of the Ga-S bonds than the Ga-Se bonds. Although their Poisson ratio values are very close to each other, 0.26 and 0.25 for GaS and GaSe, respectively, monolayer GaS is a stronger material than monolayer GaSe due to its slightly higher sU value. However, GaS and GaSe crystals are found to be more ductile and flexible materials than graphene and MoS2. We have also analyzed the band-gap response of GaS and GaSe monolayers to biaxial tensile strain and predicted a semiconductor-metal crossover after 17% and 14% applied strain, respectively, for monolayer GaS and GaSe. In addition, we investigated how the mechanical properties are affected by charging. We found that the flexibility of single layer GaS and GaSe displays a sharp increase under 0.1e/cell charging due to the repulsive interactions between extra charges located on chalcogen atoms. These charging-controllable mechanical properties of single layers of GaS and GaSe can be of potential use for electromechanical applications.