Abstract: Similar to the spontaneous formation of colonies of bacteria, flocks of birds, or schools of fish, “living crystals” can be formed by artificial self-propelled particles such as Janus colloids. Unlike usual solids, these “crystals” are far from thermodynamic equilibrium. They fluctuate in time forming a crystalline structure, breaking apart and re-forming again. We propose a method to stabilize living crystals by applying a weak confinement potential that does not suppress the ability of the particles to perform self-propelled motion, but it stabilizes the structure and shape of the dynamical clusters. This gives rise to such configurations of living crystals as “living shells” formed by Janus colloids. Moreover, the shape of the stable living clusters can be controlled by tuning the potential strength. Our proposal can be verified experimentally with either artificial microswimmers such as Janus colloids, or with living active matter.

Abstract: The autoregressive neural networks are emerging as a powerful computational tool to solve relevant problems in classical and quantum mechanics. One of their appealing functionalities is that, after they have learned a probability distribution from a dataset, they allow exact and efficient sampling of typical system configurations. Here we employ a neural autoregressive distribution estimator (NADE) to boost Markov chain Monte Carlo (MCMC) simulations of a paradigmatic classical model of spin-glass theory, namely, the two-dimensional Edwards-Anderson Hamiltonian. We show that a NADE can be trained to accurately mimic the Boltzmann distribution using unsupervised learning from system configurations generated using standard MCMC algorithms. The trained NADE is then employed as smart proposal distribution for the Metropolis-Hastings algorithm. This allows us to perform efficient MCMC simulations, which provide unbiased results even if the expectation value corresponding to the probability distribution learned by the NADE is not exact. Notably, we implement a sequential tempering procedure, whereby a NADE trained at a higher temperature is iteratively employed as proposal distribution in a MCMC simulation run at a slightly lower temperature. This allows one to efficiently simulate the spin-glass model even in the low-temperature regime, avoiding the divergent correlation times that plague MCMC simulations driven by local-update algorithms. Furthermore, we show that the NADE-driven simulations quickly sample ground-state configurations, paving the way to their future utilization to tackle binary optimization problems.

Abstract: Due to the potential application of DNA for biophysics and optoelectronics, the electronic energy states and transitions of this genetic material have attracted a great deal of attention recently. However, the fluorescence and corresponding physical process of DNA under optical excitation with photon energies below ultraviolet are still not fully clear. In this work, we experimentally investigate the photoluminescence (PL) properties of single-stranded DNA (ssDNA) samples under near-ultraviolet (NUV) and visible excitations (270 similar to 440 nm). Based on the dependence of the PL peak wavelength (lem) upon the excitation wavelength (lex), the PL behaviors of ssDNA can be approximately classified into two categories. In the relatively short excitation wavelength regime, lem is nearly constant due to exciton-like transitions associated with delocalized excitonic states and excimer states. In the relatively long excitation wavelength range, a linear relation of lem = Alex + B with A 0 or A < 0 can be observed, which comes from electronic transitions related to coupled vibrational-electronic levels. Moreover, the transition channels in different excitation wavelength regimes and the effects of strand length and base type can be analyzed on the basis of these results. These important findings not only can give a general description of the electronic energy states and transitional behaviors of ssDNA samples under NUV and visible excitations, but also can be the basis for the application of DNA in nanoelectronics and optoelectronics.

Abstract: 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.

Abstract: We present a theoretical study of the electronic structure of a heavily Si delta-doped layer in a GaAs/AlxGa1-xAs/GaAs quantum barrier. In this class of structures the effect of DX centers on the electronic properties can be tuned by changing the AlxGa1-xAs barrier width and/or the Al concentration, which leads to a lowering of the DX level with respect to the Fermi energy without disturbing the wave functions much. A self-consistent approach is developed in which the effective confinement potential and the Fermi energy of the system, the energies, the wave functions, and the electron densities of the discrete subbands have been obtained as a function of both the material parameters of the samples and the experimental conditions. The effect of DX centers on such structures at nonzero temperature and under an external pressure is investigated for three different models: (1) the DX(nc)(0) model with no correlation effects, (2) the d(+)/DX(0) model, and (3) the d(+)/DX(-) model with inclusion of correlation effects. In the actual calculation, influences of the background accepters, the discontinuity of the effective mass of the electrons at the interfaces of the different materials, band nonparabolicity, and the exchange-correlation energy of the electrons have been taken into account. We have found that (1) introducing a quantum barrier into delta-doped GaAs makes it possible to control the energy gaps between different electronic; subbands; (2) the electron wave functions are mon spread out when the repellent effect of the barriers is increased as compared to those in delta-doped GaAs; (3) increasing the quantum-barrier height and/or the application of hydrostatic pressure are helpful to experimentally observe the effect of the DX centers through a decrease of the total free-electron density; and (4) the correlation effects of the charged impurities are important for the systems under study.

Abstract: Magnetic skyrmions are topological spin configurations in materials with chiral Dzyaloshinskii-Moriya interaction (DMI), that are potentially useful for storing or processing information. To date, DMI has been found in few bulk materials, but can also be induced in atomically thin magnetic films in contact with surfaces with large spin-orbit interactions. Recent experiments have reported that isolated magnetic skyrmions can be stabilized even near room temperature in few-atom-thick magnetic layers sandwiched between materials that provide asymmetric spin-orbit coupling. Here we present the minimum-energy path analysis of three distinct mechanisms for the skyrmion collapse, based on ab initio input and the performed atomic-spin simulations. We focus on the stability of a skyrmion in three atomic layers of Co, either epitaxial on the Pt(111) surface or within a hybrid multilayer where DMI nontrivially varies per monolayer due to competition between different symmetry breaking from two sides of the Co film. In laterally finite systems, their constrained geometry causes poor thermal stability of the skyrmion toward collapse at the boundary, which we show to be resolved by designing the high-DMI structure within an extended film with lower or no DMI.

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 recently discovered structural transition in polymerized KC60 at about 50 K results in a doubling of the unit cell volume and accompanies the metal-insulator transition. Here we show that the ((a) over right arrow+(c) over right arrow,(b) over right arrow,(a) over right arrow-(c) over right arrow) superstructure results from small orientational charge density waves along the polymer chains and concomitant displacements of the surrounding K+ ions. The effect is specific for the space group Pmnn of KC60 and is absent in RbC60 and CsC60 (space group I2/m). The mechanism is relevant for the metal-insulator transition.

Abstract: Using the tight-binding approach, we study the influence of a perpendicular magnetic field on the energy levels of hexagonal, triangular, and circular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We obtain the energy levels for AB (Bernal)-stacked BLG QDs in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). We find different regions in the spectrum of biased QDs with respect to the crossing point between the lowest-electron and -hole Landau levels of a biased BLG sheet. Those different regions correspond to electron states that are localized at the center, edge, or corner of the BLG QD. Quantum Hall corner states are found to be absent in circular BLG QDs. The spatial symmetry of the carrier density distribution is related to the symmetry of the confinement potential, the position of zigzag edges, and the presence or absence of interlayer inversion symmetry.

Abstract: In a recent paper [Phys. Rev. B 89, 125418 (2014)], the authors argue that it is possible to map the electronic properties of twisted bilayer graphene to those of bilayer graphene in an in-plane magnetic field. However, their description of the low-energy dynamics of twisted bilayer graphene is restricted to the extended zone scheme and therefore neglects the effects of the superperiodic structure. If the energy spectrum is studied in the supercell Brillouin zone, we find that the comparison with an in-plane magnetic field fails because (i) the energy spectra of the two situations exhibit different symmetries and (ii) the low-energy spectra are very different.

Abstract: We reveal the emergence of metallicKondo clouds in an atomistic carbon quantum dot, realized as a single-atom junction in a suitably patterned graphene nanoflake. Using density functional dynamical mean-field theory (DFDMFT) we show how correlation effects lead to striking features in the electronic structure of our device, and how those are enhanced by the electron-electron interactions when graphene is patterned at the atomistic scale. Our setup provides a well-controlled environment to understand the principles behind the orbital-selective Kondo physics and the interplay between orbital and spin degrees of freedom in carbon-based nanomaterials, which indicate new pathways for spintronics in atomically patterned graphene.

Abstract: Using the tight-binding approach, we investigate confined states in two different hybrid monolayer-bilayer systems: (i) a hexagonal monolayer area surrounded by bilayer graphene in the presence of a perpendicularly applied electric field and (ii) a hexagonal bilayer graphene dot surrounded by monolayer graphene. The dependence of the energy levels on dot size and external magnetic field is calculated. We find that the energy spectrum for quantum dots with zigzag edges consists of states inside the gap which range from dot-localized states, edge states, to mixed states coexisting together, whereas for dots with armchair edges, only dot-localized states are observed.

Abstract: The application of graphene to electronic and optoelectronic devices is limited by the absence of reliable semiconducting variants of this material. A promising candidate in this respect is graphene oxide, with a band gap on the order of similar to 5 eV, however, this has a finite density of states at the Fermi level. Here, we examine the electronic structure of three variants of half-fluorinated carbon on Sic(0001), i.e., the (6 root 3 x 6 root 3) R30 degrees C/SiC “buffer layer,” graphene on this (6 root 3 x 6 root 3) R30 degrees C/SiC buffer layer, and graphene decoupled from the SiC substrate by hydrogen intercalation. Using angle-resolved photoemission, core level photoemission, and x-ray absorption, we show that the electronic, chemical, and physical structure of all three variants is remarkably similar, exhibiting a large band gap and a vanishing density of states at the Fermi level. These results are explained in terms of first-principles calculations. This material thus appears very suitable for applications, even more so since it is prepared on a processing-friendly substrate. We also investigate two separate UV photon-induced modifications of the electronic structure that transform the insulating samples (6.2-eV band gap) into semiconducting (similar to 2.5-eV band gap) and metallic regions, respectively.

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.