Vacancy clustering effect on the electronic and transport properties of bilayer graphene nanoribbons.

Experimental realizations of two-dimensional materials are hardly free of structural defects such as e.g. vacancies, which, in turn, modify drastically its pristine physical defect-free properties. In this work, we explore effects due to point defect clustering on the electronic and transport properties of bilayer graphene nanoribbons, for AA and AB stacking and zigzag and armchair boundaries, by means of the tight-binding approach and scattering matrix formalism. Evident vacancy concentration signatures exhibiting a maximum amplitude and an universality regardless of the system size, stacking and boundary types, in the density of states around the zero-energy level are observed. Our results are explained via the coalescence analysis of the strong sizeable vacancy clustering effect in the system and the breaking of the inversion symmetry at high vacancy densities, demonstrating a similar density of states for two equivalent degrees of concentration disorder, below and above the maximum value.

Zitterbewegung of Moiré Excitons in Twisted MoS_{2}/WSe_{2} Heterobilayers.

The moiré pattern observed in stacked noncommensurate crystal lattices, such as heterobilayers of transition metal dichalcogenides, produces a periodic modulation of their band gap. Excitons subjected to this potential landscape exhibit a band structure that gives rise to a quasiparticle dubbed the moiré exciton. In the case of MoS_{2}/WSe_{2} heterobilayers, the moiré trapping potential has honeycomb symmetry and, consequently, the moiré exciton band structure is the same as that of a Dirac-Weyl fermion, whose mass can be further tuned down to zero with a perpendicularly applied field. Here we show that, analogously to other Dirac-like particles, the moiré exciton exhibits a trembling motion, also known as Zitterbewegung, whose long timescales are compatible with current experimental techniques for exciton dynamics. This promotes the study of the dynamics of moiré excitons in van der Waals heterostructures as an advantageous solid-state platform to probe Zitterbewegung, broadly tunable by gating and interlayer twist angle.

Transition-metal adatoms on 2D-GaAs: a route to chiral magnetic 2D materials by design.

Using relativistic density-functional calculations, we examine the magneto-crystalline anisotropy and exchange properties of transition-metal atoms adsorbed on 2D-GaAs. We show that single Mn and Mo atom (Co and Os) strongly bind on 2D-GaAs, and induce local out-of-plane (in-plane) magnetic anisotropy. When a pair of TM atoms is adsorbed on 2D-GaAs in a close range from each other, magnetisation properties change (become tunable) with respect to concentrations and ordering of the adatoms. In all cases, we reveal presence of strong Dzyaloshinskii-Moriya interaction. These results indicate novel pathways towards two-dimensional chiral magnetic materials by design, tailored for desired applications in magneto-electronics.

Effect of zitterbewegung on the propagation of wave packets in ABC-stacked multilayer graphene: an analytical and computational approach.

The time evolution of a low-energy two-dimensional Gaussian wave packet in ABC-stacked n-layer graphene (ABC-NLG) is investigated. Expectation values of the position (x, y) of center-of-mass and the total probability densities of the wave packet are calculated analytically using the Green's function method. These results are confirmed using an alternative numerical method based on the split-operator technique within the Dirac approach for ABC-NLG, which additionally allows to include external fields and potentials. The main features of the zitterbewegung (trembling motion) of wave packets in graphene are demonstrated and are found to depend not only on the wave packet width and initial pseudospin polarization, but also on the number of layers. Moreover, the analytical and numerical methods proposed here allow to investigate wave packet dynamics in graphene systems with an arbitrary number of layers and arbitrary potential landscapes.

Magnetic field induced vortices in graphene quantum dots.

The energy spectrum and local current patterns in graphene quantum dots (QD) are investigated for different geometries in the presence of an external perpendicular magnetic field. Our results demonstrate that, for specific geometries and edge configurations, the QD exhibits vortex and anti-vortex patterns in the local current density, in close analogy to the vortex patterns observed in the probability density current of semiconductor QD, as well as in the order parameter of mesoscopic superconductors.

Structure and reentrant percolation in an inverse patchy colloidal system.

Two-dimensional systems of inverse patchy colloids modeled as disks with a central charge and having their surface decorated with oppositely pointlike charged patches are investigated using molecular dynamics simulations. The self-assembly of the patchy colloids leads to diverse ground state configurations ranging from crystalline arrangements of monomers to linear clusters, ramified linear clusters and to percolated configurations. Two structural phase diagrams are constructed: (1) as a function of the net charge and area fraction, and (2) as a function of the net charge and the range of the pair interaction potential. An interesting reentrant percolation transition is obtained as a function of the net charge of the colloids. We identify distinct mechanisms that lead to the percolation transition.

Giant paramagnetic Meissner effect in multiband superconductors.

Superconductors, ideally diamagnetic when in the Meissner state, can also exhibit paramagnetic behavior due to trapped magnetic flux. In the absence of pinning such paramagnetic response is weak, and ceases with increasing sample thickness. Here we show that in multiband superconductors paramagnetic response can be observed even in slab geometries, and can be far larger than any previous estimate - even multiply larger than the diamagnetic Meissner response for the same applied magnetic field. We link the appearance of this giant paramagnetic response to the broad crossover between conventional Type-I and Type-II superconductors, where Abrikosov vortices interact non-monotonically and multibody effects become important, causing unique flux configurations and their locking in the presence of surfaces.

Tuning the polarized quantum phonon transmission in graphene nanoribbons.

We propose systems that allow a tuning of the phonon transmission function T(ω) in graphene nanoribbons by using C(13) isotope barriers, antidot structures, and distinct boundary conditions. Phonon modes are obtained by an interatomic fifth-nearest neighbor force-constant model (5NNFCM) and T(ω) is calculated using the non-equilibrium Green's function formalism. We show that by imposing partial fixed boundary conditions it is possible to restrict contributions of the in-plane phonon modes to T(ω) at low energy. On the contrary, the transmission functions of out-of-plane phonon modes can be diminished by proper antidot or isotope arrangements. In particular, we show that a periodic array of them leads to sharp dips in the transmission function at certain frequencies [Formula: see text] which can be pre-defined as desired by controlling their relative distance and size. With this, we demonstrated that by adequate engineering it is possible to govern the magnitude of the ballistic transmission functions [Formula: see text] in graphene nanoribbons. We discuss the implications of these results in the design of controlled thermal transport at the nanoscale as well as in the enhancement of thermo-electric features of graphene-based materials.

Dynamical properties and melting of binary two-dimensional colloidal alloys.

A two-dimensional (2D) binary colloidal system consisting of interacting dipoles is investigated using an analytical approach. Within the harmonic approximation we obtain the phonon spectrum of the system as a function of the composition, dipole-moment ratio, and mass ratio between the small and big particles. Through a systematic analysis of the phonon spectra we are able to determine the stability region of the different lattice structures of the colloidal alloys. The gaps in the phonon frequency spectrum, the optical frequencies in the long-wavelength limit, and the sound velocity are discussed as well. Using the modified Lindemann criterion and within the harmonic approximation we estimate the melting temperature of the sublattice generated by the big particles.

Angular melting scenarios in binary dusty-plasma Coulomb balls: magic versus normal clusters.

Molecular-dynamic simulations were performed in order to investigate the melting processes of isotropically confined binary systems. We considered two species of particles, which differ by their amount of electric charge. A Lindemann type of criterion was used to determine the angular melting temperature. We demonstrate that the magic-to-normal cluster transition can evolve in two distinct ways, that is, through a structural phase transition of the first order or via a smooth transition where an increase of the shells' width leads to a continuous decreasing mechanical stability of the system. Moreover, for large systems, we demonstrate that the internal cluster exerts a minor effect on the mechanical stability of the external shell. Furthermore, we show that highly symmetric configurations, such as those found for multiple ring structures, have large mechanical stability, i.e., high angular melting temperature.

Tunable diffusion of magnetic particles in a quasi-one-dimensional channel.

The diffusion of a system of ferromagnetic dipoles confined in a quasi-one-dimensional parabolic trap is studied using Brownian dynamics simulations. We show that the dynamics of the system is tunable by an in-plane external homogeneous magnetic field. For a strong applied magnetic field, we find that the mobility of the system, the exponent of diffusion, and the crossover time among different diffusion regimes can be tuned by the orientation of the magnetic field. For weak magnetic fields, the exponent of diffusion in the subdiffusive regime is independent of the orientation of the external field.

Single-file diffusion in periodic energy landscapes: the role of hydrodynamic interactions.

We report on the dynamical properties of interacting colloids confined to one dimension and subjected to external periodic energy landscapes. We particularly focus on the influence of hydrodynamic interactions on the mean-square displacement. Using Brownian dynamics simulations, we study colloidal systems with two types of repulsive interparticle interactions, namely, Yukawa and superparamagnetic potentials. We find that in the homogeneous case, hydrodynamic interactions lead to an enhancement of the particle mobility and the mean-square displacement at long times scales as t(α), with α=1/2+ε and ε being a small correction. This correction, however, becomes much more important in the presence of an external field, which breaks the homogeneity of the particle distribution along the line and, therefore, promotes a richer dynamical scenario due to the hydrodynamical coupling among particles. We provide here the complete dynamical scenario in terms of the external potential parameters: amplitude and commensurability.

Bilayer crystals of charged magnetic dipoles: structure and phonon spectrum.

We study the structure and phonon spectrum of a two-dimensional bilayer system of classical charged dipoles oriented perpendicular to the plane of the layers for equal density in each layer. This system can be tuned through six different crystalline phases by changing the interlayer separation or the charge and/or dipole moment of the particle. The presence of the charge on the dipole particles is responsible for the nucleation of five staggered phases and a disordered phase which are not found in the magnetic dipole bilayer system. These extra phases are a consequence of the competition between the repulsive Coulomb and the attractive dipole interlayer interaction. We present the phase diagram and determine the order of the phase transitions. The phonon spectrum of the system was calculated within the harmonic approximation, and a nonmonotonic behavior of the phonon spectrum is found as a function of the effective strength of the interparticle interaction. The stability of the different phases is determined.

The electronic properties of graphene and graphene ribbons under simple shear strain.

A tight-binding model is used to study the energy band of graphene and graphene ribbon under simple shear strain. The ribbon consists of lines of carbon atoms in an armchair or zigzag orientation where a simple shear strain is applied in the x-direction keeping the atomic distances in the y-direction unchanged. Such modification in the lattice gives an energy band that differs in several aspects from the one without any shear and with pure shear. The changes in the spectrum depend on the line displacement of the ribbon, and also on the modified hopping parameter. It is also shown that this simple shear strain tunes the electronic properties of both graphene and graphene ribbon, opening and closing energy gaps for different displacements of the system. The modified density of states is also shown.

Transition from single-file to two-dimensional diffusion of interacting particles in a quasi-one-dimensional channel.

Diffusive properties of a monodisperse system of interacting particles confined to a quasi-one-dimensional channel are studied using molecular dynamics simulations. We calculate numerically the mean-squared displacement (MSD) and investigate the influence of the width of the channel (or the strength of the confinement potential) on diffusion in finite-size channels of different shapes (i.e., straight and circular). The transition from single-file diffusion to the two-dimensional diffusion regime is investigated. This transition [regarding the calculation of the scaling exponent (α) of the MSD (Δx(2)(t) ∝ t(α)] as a function of the width of the channel is shown to change depending on the channel's confinement profile. In particular, the transition can be either smooth (i.e., for a parabolic confinement potential) or rather sharp (i.e., for a hard-wall potential), as distinct from infinite channels where this transition is abrupt. This result can be explained by qualitatively different distributions of the particle density for the different confinement potentials.

Diffusion in a quasi-one-dimensional system on a periodic substrate.

The diffusion of charged particles interacting through a repulsive Yukawa potential, exp(-r/λ)/r, confined by a parabolic potential in the y direction and subjected to a periodic substrate potential in the x direction is investigated. Langevin dynamic simulations are used to investigate the effect of the particle density, the amplitude of the periodic substrate, and the range of the interparticle interaction potential on the diffusive behavior of the particles. We found that in general the diffusion is suppressed with increasing the amplitude of the periodic potential, but for specific values of the strength of the substrate potential a remarkable increase of the diffusion is found with increasing the periodic potential amplitude. In addition, we found a strong dependence of the diffusion on the specific arrangement of the particles, e.g., single-chain versus multichain configuration. For certain particle configurations, a reentrant behavior of the diffusion is found as a function of the substrate strength due to structural transitions in the ordering of the particles.

Binary dusty plasma Coulomb balls.

We investigated the mixing and segregation of a system consisting of two different species of particles, having different charges, interacting through a pure Coulomb potential, and confined in a three-dimensional parabolic trap. The structure of the cluster and its normal mode spectrum are analyzed as a function of the relative charge and the relative number of different types of particles. We found that (a) the system can be in a mixed or segregated state depending on the relative charge ratio parameter and (b) the segregation process is mediated by a first or second order structural phase transition which strongly influences the magic cluster properties of the system.

A two-component mixture of charged particles confined in a channel: melting.

The melting of a binary system of charged particles confined in a quasi-one-dimensional parabolic channel is studied through Monte Carlo simulations. At zero temperature the particles are ordered in parallel chains. The melting is anisotropic and different melting temperatures are obtained according to the spatial direction, and the different kinds of particles present in the system. Melting is very different for the single-, two- and four-chain configurations. A temperature induced structural phase transition is found between two different four-chain ordered states which is absent in the mono-disperse system. In the mixed regime, where the two kinds of particles are only slightly different, melting is almost isotropic and a thermally induced homogeneous distribution of the distinct kinds of charges is observed.

Melting and evaporation in classical two-dimensional clusters confined by a Coulomb potential.

The thermal properties of a two-dimensional classical cluster of negatively charged particles bound by a punctual positive charge are presented. The melting phenomenon is analyzed and the features which characterize such a solid-liquid transition are highlighted. We found that the presence of metastable states strongly modifies the melting scenario, and that the melting temperature of the system is determined by the height of the saddle point energy separating the ground state and the metastable state. Due to the particular type of confinement potential considered in this paper, we also found that, at sufficiently large temperature, the cluster can become thermally ionized.

Structure, normal mode spectra, and mixing of a binary system of charged particles confined in a parabolic trap.

We study the mixing of two different kinds of particles, having different charge and/or mass, interacting through a pure Coulomb potential, and confined in a parabolic trap. The structure of the cluster and its normal mode spectrum are analyzed as a function of the ratio of the charges (mass ratio) of the two types of particles. We show that particles are not always arranged in a shell structure. Mixing of the particles goes hand in hand with a large number of metastable states. The normal modes of the system are obtained, and we find that some of the special modes can be tuned by varying the ratio between the charges (masses) of the two species. The degree of mixing of the two type of particles is summarized in a phase diagram, and an order parameter that describes quantitatively the mixing between particles is defined.