Abstract: Large-scale atomistic simulations using the reactive force field approach are implemented to investigate the thermomechanical properties of fluorinated graphene (FG). A set of parameters for the reactive force field potential optimized to reproduce key quantum mechanical properties of relevant carbon-fluorine cluster systems are presented. Molecular dynamics simulations are used to investigate the thermal rippling behavior of FG and its mechanical properties and compare them with graphene, graphane and a sheet of boron nitride. The mean square value of the height fluctuations < h(2)> and the height-height correlation function H(q) for different system sizes and temperatures show that FG is an unrippled system in contrast to the thermal rippling behavior of graphene. The effective Young's modulus of a flake of fluorinated graphene is obtained to be 273 N/m and 250 N/m for a flake of FG under uniaxial strain along armchair and zigzag directions, respectively. DOI: 10.1103/PhysRevB.87.104114

Abstract: Layered structures of molybdenum disulfide (MoS2) belong to a new class of two-dimensional (2D) semiconductor materials in which monolayers exhibit a direct band gap in their electronic spectrum. This band gap has recently been shown to vanish due to the presence of metallic edge modes when MoS2 monolayers are terminated by zigzag edges on both sides. Here, we demonstrate that a zigzag nanoribbon of MoS2, when exposed to an external exchange field in combination with a transverse electric field, has the potential to exhibit a peculiar half-metallic nature and thereby allows electrons of only one spin direction to move. The peculiarity of such spin-selective conductors originates from a spin switch near the gap-closing region, so the allowed spin orientation can be controlled by means of an external gate voltage. It is shown that the induced half-metallic phase is resistant to random fluctuations of the exchange field as well as the presence of edge vacancies.

Abstract: Magnesium hydroxide [Mg(OH)(2)] has a layered brucitelike structure in its bulk form and was recently isolated as a new member of two-dimensional monolayer materials. We investigated the electronic and optical properties of monolayer crystals of Mg(OH)(2) and WS2 and their possible heterobilayer structure by means of first-principles calculations. It was found that both monolayers of Mg(OH)(2) and WS2 are direct-gap semiconductors and these two monolayers form a typical van der Waals heterostructure with a weak interlayer interaction and a type-II band alignment with a staggered gap that spatially separates electrons and holes. We also showed that an out-of-plane electric field induces a transition from a staggered to a straddling-type heterojunction. Moreover, by solving the Bethe-Salpeter equation on top of single-shot G(0)W(0) calculations, we show that the low-energy spectrum of the heterobilayer is dominated by the intralyer excitons of the WS2 monolayer. Because of the staggered interfacial gap and the field-tunable energy-band structure, the Mg(OH)(2)-WS2 heterobilayer can become an important candidate for various optoelectronic device applications in nanoscale.

Abstract: Magnetic nanoparticles are critical building blocks for future technologies ranging from nanomedicine to spintronics. Many related applications require nanoparticles with tailored magnetic properties. However, despite significant efforts undertaken towards this goal, a broad and poorly understood dispersion of magnetic properties is reported, even within monodisperse samples of the canonical ferromagnetic 3d transition metals. We address this issue by investigating the magnetism of a large number of size- and shape-selected, individual nanoparticles of Fe, Co, and Ni using a unique set of complementary characterization techniques. At room temperature, only superparamagnetic behavior is observed in our experiments for all Ni nanoparticles within the investigated sizes, which range from 8 to 20 nm. However, Fe and Co nanoparticles can exist in two distinct magnetic states at any size in this range: (i) a superparamagnetic state, as expected from the bulk and surface anisotropies known for the respective materials and as observed for Ni, and (ii) a state with unexpected stable magnetization at room temperature. This striking state is assigned to significant modifications of the magnetic properties arising from metastable lattice defects in the core of the nanoparticles, as concluded by calculations and atomic structural characterization. Also related with the structural defects, we find that the magnetic state of Fe and Co nanoparticles can be tuned by thermal treatment enabling one to tailor their magnetic properties for applications. This paper demonstrates the importance of complementary single particle investigations for a better understanding of nanoparticle magnetism and for full exploration of their potential for applications.

Abstract: Magnetization reversal is a well-studied problem with obvious applicability in computer hard drives. One can accomplish a magnetization reversal in at least one of two ways: application of a magnetic field or through a spin current. The latter is more amenable to a fully quantum-mechanical analysis. We formulate and solve the problem whereby a spin current interacts with a ferromagnetic Heisenberg spin chain, to eventually reverse the magnetization of the chain. Spin flips are accomplished through both elastic and inelastic scattering. A consequence of the inelastic-scattering channel, when it is no longer energetically possible, is the occurrence of a nonequilibrium bound state, which is an emergent property of the coupled local plus itinerant spin system. For certain definite parameter values the itinerant spin lingers near the local spins for some time, before eventually leaking out as an outwardly diffusing state. This phenomenon results in spin-flip dynamics and filtering properties for this type of system.

Abstract: Magnetotransport in bilayer graphene in the presence of a weak and periodic potential is investigated in the presence of a perpendicular magnetic field B. The modulation broadens the Landau levels into bands and for weak magnetic fields leads to the well-known Weiss oscillations in their bandwidth and their transport coefficients at very low B and to the Shubnikov-de Haas oscillations at larger B. The amplitude of the Weiss oscillations is severely reduced if the periodic potentials applied to the two layers oscillate out of phase. We also contrast some results with those corresponding to single-layer graphene. Relative to them the flat-band condition and the oscillation amplitude differ substantially, due to the interlayer coupling, and agree only when this coupling is extremely weak. We further show that the Hall conductivity exhibits the well-known steps at half-integer and integer multiples of 4e(2)/h in single-layer and bilayer graphene, respectively, even for very weak magnetic fields. The results are pertinent to weak and periodic corrugations when the potential modulation dominates the strain-induced magnetic modulation.

Abstract: Massless Dirac fermions (MDFs) emerge as quasiparticles in various novel materials such as graphene and topological insulators, and they exhibit several intriguing properties, of which Veselago focusing is an outstanding example with a lot of possible applications. However, up to now Veselago focusing merely occurred in p-n junction devices based on the isotropic MDF, which lacks the tunability needed for realistic applications. Here, motivated by the emergence of novel Dirac materials, we investigate the propagation behaviors of anisotropic MDFs in such a p-n junction structure. By projecting the Hamiltonian of the anisotropic MDF to that of the isotropic MDF and deriving an exact analytical expression for the propagator, precise Veselago focusing is demonstrated without the need for mirror symmetry of the electron source and its focusing image. We show a tunable focusing position that can be used in a device to probe masked atom-scale defects. This study provides an innovative concept to realize Veselago focusing relevant for potential applications, and it paves the way for the design of novel electron optics devices by exploiting the anisotropic MDF.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: Motivated by recent experimental observations of Tongay et al. [Nano Lett. 12, 5576 (2012)] we show how the electronic properties and Raman characteristics of single layer MoSe2 are affected by elastic biaxial strain. We found that with increasing strain: (1) the E' and E '' Raman peaks (E-2g and E-1g in bulk) exhibit significant redshifts (up to similar to 30 cm(-1)), (2) the position of the A'(1) peak remains at similar to 180 cm(-1) (A(1g) in bulk) and does not change considerably with further strain, (3) the dispersion of low energy flexural phonons crosses over from quadratic to linear, and (4) the electronic band structure undergoes a direct to indirect band gap crossover under similar to 3% biaxial tensile strain. Thus the application of strain appears to be a promising approach for a rapid and reversible tuning of the electronic, vibrational, and optical properties of single layer MoSe2 and similar MX2 dichalcogenides. DOI:10.1103/PhysRevB.87.125415

Abstract: Motivated by two-band superconductivity in, e.g., borides and pnictides, starting from the two-band Ginzburg-Landau energy functional, we discuss how the presence of two order parameters and the coupling between them influence a superconducting ring in the fluctuative regime. Our method is an extension of the von OppenRiedel formalism for rings; it is exact, but requires numerical implementation. We also study approximations for which analytic expressions can be obtained, and check their ranges of validity. We provide estimates for the temperature ranges where fluctuations are important, calculate the persistent current in MgB2 rings as a function of temperature and enclosed flux, and point out its additional dependence on the cross-section area of the wire from which the ring is made. We find temperature regions in which fluctuations enhance the persistent currents and regions where they inhibit the persistent current. The presence of two order parameters that can fluctuate independently always leads to larger averages of the order parameters at Tc, but yields larger persistent current only for appropriate parameters. In cases of very different material parameters for the two coupled condensates, the persistent current is inhibited.

Abstract: Nanoindentation of bilayer graphene is studied using molecular-dynamics simulations. We compared our simulation results with those from elasticity theory as based on the nonlinear Föppl-Hencky equations with rigid boundary condition. The force-deflection values of bilayer graphene are compared to those of monolayer graphene. Youngs modulus of bilayer graphene is estimated to be 0.8 TPa which is close to the value for graphite. Moreover, an almost flat bilayer membrane at low temperature under central load has a 14% smaller Youngs modulus as compared to the one at room temperature.

Abstract: Neutral excitons in strained axially symmetric In(Ga)As/GaAs quantum dots with a ringlike shape are investigated. Similar to experimental self-assembled quantum rings, the analyzed quantum dots have volcano-like shapes. The continuum mechanical model is employed to determine the strain distribution, and the single-band envelope function approach is adopted to compute the electron states. The hole states are determined by the axially symmetric multiband Luttinger-Kohn Hamiltonian, and the exciton states are obtained from an exact diagonalization. We found that the presence of the inner layer covering the ring opening enhances the excitonic Aharonov-Bohm (AB) oscillations. The reason is that the hole becomes mainly localized in the inner part of the quantum dot due to strain, whereas the electron resides mainly inside the ring-shaped rim. Interestingly, larger AB oscillations are found in the analyzed quantum dot than in a fully opened quantum ring of the same width. Comparison with the unstrained ringlike quantum dot shows that the amplitude of the excitonic Aharonov-Bohm oscillations are almost doubled in the presence of strain. The computed oscillations of the exciton energy levels are comparable in magnitude to the oscillations measured in recent experiments. DOI: 10.1103/PhysRevB.87.085314

Abstract: Nonlinear dynamics of Josephson vortices (fluxons) in artificial stacks of superconducting-normal-superconducting Josephson junctions under simultaneously applied time-periodic ac and constant biasing dc currents is studied using the time dependent Ginzburg-Landau formalism with a Lawrence-Doniach extension. At zero external magnetic field and dc biasing current the resistive state of the system is characterized by periodic nucleation and annihilation of fluxon-antifluxon pairs, relative positions of which are determined by the state of neighboring junctions. Due to the mutual repulsive interaction, fluxons in different junctions move out of phase. Their collective motion can be synchronized by adding a small ac component to the biasing dc current. Coherent motion of fluxons is observed for a broad frequency range of the applied drive. In the coherent state the maximal output voltage, which is proportional to the number of junctions in the stack, is observed near the characteristic frequency of the system determined by the crossing of the fluxons across the sample. However, in this frequency range the dynamically synchronized state has an alternative-a less ordered state with smaller amplitude of the output voltage. Collective behavior of the junctions is strongly affected by the sloped sidewalls of the stack. Synchronization is observed only for weakly trapezoidal cross sections, whereas irregular motion of fluxons is observed for larger slopes of the sample edge.

Abstract: On the basis of the momentum-balance equation derived from the Boltzmann equation in which electron interactions with impurities and acoustic and optic phonons are included, we examine the dependence of the resistivity in graphene on temperature and electron density. Simple analytical expressions for the different contributions to the resistivity are obtained. Our results reproduce recent experimental findings and we are able to understand the different temperature dependence of the resistivity for low and high density samples.

Abstract: Phase-slip centers/lines and hot spots are the main mechanisms for dissipation in current-carrying superconducting thin films. The pulsed-current method has recently been shown to be an effective tool in studying the dynamics of phase-slip centers and their evolution to hot spots. We use the time-dependent Ginzburg-Landau theory in the study of the dynamics of the superconducting condensate in superconducting strips under external current and zero external magnetic field. We show that both the flux-flow state (i.e., slow-moving vortices) and the phase-slip line state (i.e., fast-moving vortices) are dynamically stable dissipative units with temperature smaller than the critical one, whereas hot spots, which are localized normal regions where the local temperature exceeds the critical value, expand in time, resulting ultimately in a complete destruction of the condensate. The response time of the system to abrupt switching on of the overcritical current decreases with increasing both the value of the current (at all temperatures) and temperature (for a given value of the applied current). Our results are in good qualitative agreement with experiments we have conducted on Nb thin strips.