Abstract: The effect of the number of stacking layers and the type of stacking on the electronic and optical properties of bilayer and trilayer black phosphorus are investigated by using first-principles calculations within the framework of density functional theory. We find that inclusion of many-body effects (i.e., electron-electron and electron-hole interactions) modifies strongly both the electronic and optical properties of black phosphorus. While trilayer black phosphorus with a particular stacking type is found to be a metal by using semilocal functionals, it is predicted to have an electronic band gap of 0.82 eV when many-body effects are taken into account within the G(0)W(0) scheme. Though different stacking types result in similar energetics, the size of the band gap and the optical response of bilayer and trilayer phosphorene are very sensitive to the number of layers and the stacking type. Regardless of the number of layers and the type of stacking, bilayer and trilayer black phosphorus are direct band gap semiconductors whose band gaps vary within a range of 0.3 eV. Stacking arrangements that are different from the ground state structure in both bilayer and trilayer black phosphorus exhibit significant modified valence bands along the zigzag direction and result in larger hole effective masses. The optical gap of bilayer (trilayer) black phosphorus varies by 0.4 (0.6) eV when changing the stacking type. The calculated binding energy of the bound exciton hardly changes with the type of stacking and is found to be 0.44 (0.30) eV for bilayer (trilayer) phosphorous.

Abstract: Nitrogenated holey graphene (NHG) is a recently synthesized two-dimensional material. In this paper the structural and electronic properties of heterostructures of graphene and NHG are investigated using first-principles and tight-binding calculations. Due to the lattice mismatch between NHG and graphene, the formation of a moire pattern is preferred in the graphene/NHG heterostructure, instead of a lattice-coherent structure. In moire-patterned graphene/NHG, the band gap opening at the K point is negligible, and the linear band dispersion of graphene survives. Applying an electric field modifies the coupling strength between the two atomic layers. The Fermi velocity upsilon(F) is reduced as compared to the one of pristine graphene, and its magnitude depends on the twist angle theta between graphene and NHG: For theta = 0 degrees, upsilon(F) is 30% of that of graphene, and it increases rapidly to a value of 80% with increasing theta. The heterostructure exhibits electron-hole asymmetry in upsilon(F), which is large for small theta. In NHG encapsulated between two graphene layers, a “Dirac ring” appears around the K point. Its presence is robust with respect to the relative stacking of the two graphene layers. These findings can be useful for future applications of graphene/NHG heterostructures.

Abstract: Experiments have confirmed that double monolayer graphene does not generate finite-temperature electron-hole superfluidity, because of very strong screening of the pairing attraction. The linear dispersing energy bands in monolayer graphene block any attempt to reduce the strength of the screening. We propose a hybrid device with two sheets of monolayer graphene in a modulated periodic perpendicular magnetic field. The field preserves the isotropic Dirac cones of the original monolayers but reduces the slope of the cones, making the monolayer Fermi velocity v(F) smaller. We demonstrate that with current experimental techniques, the reduction in vF can weaken the screening sufficiently to allow electron-hole superfluidity at measurable temperatures.

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: 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 electron distribution function (EDF) in an electron cyclotron resonance (ECR) discharge is far from Maxwellian. The self-consistent simulation of ECR discharges requires a calculation of the EDF on every magnetic line for various ion density profiles. The straightforward self-consistent simulation of ECR discharges using the Monte Carlo technique for the EDF calculation is very computer time expensive, since the electron and ion time scales are very different. An electron Boltzmann kinetic equation averaged over the fast electron bouncing and pitch-angle scattering was derived in order to develop an effective and operative tool for the fast modeling (FM) of low-pressure ECR discharges. An analytical solution for the EDF in a loss cone was derived. To check the validity of the FM, one-dimensional (in coordinate) and two-dimensional (in velocity) Monte Carlo simulation codes were developed. The validity of the fast modeling method is proved by comparison with the Monte Carlo simulations. The complete system of equations for FM is presented and ready for use in a comprehensive study of ECR discharges. The variations of plasma density and of wall and sheath potentials are analyzed by solving a self-consistent set of equations for the EDF.

Abstract: The purpose of this Brief Report is to point out the mistake in a formula for anisotropic electron scattering, previously published in Phys. Rev. A 41, 1112 (1990), which is widely used in Monte Carlo models of gas discharges. Anisotropic electron scattering is investigated based on the screened Coulomb potential between electrons and neutral atoms. The approach is also applied for electron scattering by nonpolar neutral molecules. Differential cross sections for electron scattering by Ar, N2, and CH4 are constructed on the basis of momentum and integrated cross sections. The formula derived in this paper is useful for Monte Carlo simulations of gas discharges.

Abstract: The ordering of N equally charged particles (-e) moving in two dimensions and confined by a Coulomb potential, resulting from a displaced positive charge Ze is discussed. This is a classical model system for atoms. We obtain the configurations of charged particles which, depending on the value of N and Z, may result in ring structures, hexagonal-type configurations, and for N/Z approximate to 1 in an inner structure of particles which is separated by an outer ring of particles. For N/Z << 1, the Hamiltonian of the parabolic confinement case is recovered. For N/Z approximate to 1, the configurations are very different from those found in the case of a parabolic confinement potential. A hydrodynamic analysis is presented in order to highlight the correlations effects.

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

Abstract: The diffusion of charged particles interacting through a repulsive Yukawa potential, exp(-r/lambda)/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.

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

Abstract: 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. DOI: 10.1103/PhysRevE.87.012307

Abstract: Janus particles present an important class of building blocks for directional assembly. These are compartmentalized colloids with two different hemispheres. Here, we consider a two-dimensional model of Janus disks consisting of a hydrophobic semicircle and an electro-negatively charged one. Placed in a solution, the hydrophobic sides will attract each other while the charged sides will give rise to a repulsive force. Using molecular dynamics simulations, we study the morphology of these particles when confined in a channel-like environment using a one dimensional harmonic confinement potential. The interest to this system is first of all due to the fact that it could serve as a simple model for membrane formation. Indeed, the recently synthesized new class of artificial amphiphiles, known as Janus dendrimers, were shown to self-assemble in bilayer structures mimicking biological membranes. In turn, Janus particles that combine the amphiphilicity and colloidal rigidity serve as a good model for Janus dendrimers. A variety of ordered membrane-like morphologies are found consisting of single and multiple chain configurations with different orientations of the particles with respect to each other that we summarize in a phase diagram.

Abstract: The ground state of colloidal magnetic particles in a modulated channel are investigated as a function of the tilt angle of an applied magnetic field. The particles are confined by a parabolic potential in the transversal direction while in the axial direction a periodic substrate potential is present. By using Monte Carlo simulations, we construct a phase diagram for the different crystal structures as a function of the magnetic field orientation, strength of the modulated potential, and the commensurability factor of the system. Interestingly, we found first-and second-order phase transitions between different crystal structures, which can be manipulated by the orientation of the external magnetic field. A reentrant behavior is found between two-and four-chain configurations, with continuous second-order transitions. Novel configurations are found consisting of frozen solitons of defects. By changing the orientation and/or strength of the magnetic field and/or the strength and periodicity of the substrate potential, the system transits through different phases.

Abstract: We expand on a recent study of a lattice model of interacting particles [Phys. Rev. Lett. 111, 110601 (2013)]. The adsorption isotherm and equilibrium fluctuations in particle number are discussed as a function of the interaction. Their behavior is similar to that of interacting particles in porous materials. Different expressions for the particle jump rates are derived from transition-state theory. Which expression should be used depends on the strength of the interparticle interactions. Analytical expressions for the self-and transport diffusion are derived when correlations, caused by memory effects in the environment, are neglected. The diffusive behavior is studied numerically with kinetic Monte Carlo (kMC) simulations, which reproduces the diffusion including correlations. The effect of correlations is studied by comparing the analytical expressions with the kMC simulations. It is found that the Maxwell-Stefan diffusion can exceed the self-diffusion. To our knowledge, this is the first time this is observed. The diffusive behavior in one-dimensional and higher-dimensional systems is qualitatively the same, with the effect of correlations decreasing for increasing dimension. The length dependence of both the self-and transport diffusion is studied for one-dimensional systems. For long lengths the self-diffusion shows a 1/L dependence. Finally, we discuss when agreement with experiments and simulations can be expected. The assumption that particles in different cavities do not interact is expected to hold quantitatively at low and medium particle concentrations if the particles are not strongly interacting.