Abstract: We present a theoretical study on plasmon-polariton modes in graphene placed inside an optical cavity and driven by a source-to-drain electric field. The electron velocity and electron temperature are determined by solving self-consistently the momentum-and energy-balance equations in which electron interactions with impurities, acoustic-, and optic-phonons are included. Based on many-body self-consistent field theory, we develop a tractable approach to study plasmon-polariton in an electron gas system. We find that when graphene is placed inside a Fabry-Perot cavity, two branches of the plasmon-polariton modes can be observed and these modes are very much optic-or plasmon-like. The frequencies of these modes depend markedly on driving electric field especially at higher resonant frequency regime. Moreover, the plasmon-polariton frequency in graphene is in terahertz (THz) bandwidth and can be tuned by changing the cavity length, gate voltage, and driving electric field. This work is pertinent to the application of graphene-based structures as tunable THz plasmonic devices. (C) 2015 AIP Publishing LLC.

Abstract: We study the influence of confinement on the atomic collapse due to a Coulomb impurity placed at the center of a graphene quantum dot of radius R. We apply the zigzag or infinite-mass boundary condition and consider both a point-size and a finite-size impurity. As a function of the impurity strength Za, the energy spectra are discrete. In the case of the zigzag boundary condition, the degenerate (with respect to the angular momentum m) zero-energy levels are pulled down in energy as Z alpha increases, and they remain below epsilon = – Z alpha. Our results show that the energy levels exhibit a 1/R dependence in the subcritical regime [Z alpha < |km + 1/2|, k = 1 (-1) for the K (K') valley]. In the supercritical regime (Z alpha > |km + 1/2|) we find a qualitatively very different behavior where the levels decrease as a function of R in a nonmonotonic manner. While the valley symmetry is preserved in the presence of the impurity, we find that the impurity breaks electron-hole symmetry. We further study the energy spectrum of zigzag quantum dots in gapped graphene. Our results show that as the gap increases, the lowest electron states are pushed into the gap by the impurity.

Abstract: Optical and electronic properties of black phosphorus strongly depend on the number of layers and type of stacking. Using first-principles calculations within the framework of density functional theory, we investigate the electronic properties of bilayer black phosphorus with an interlayer twist angle of 90 degrees. These calculations are complemented with a simple (k) over right arrow . (p) over right arrow model which is able to capture most of the low energy features and is valid for arbitrary twist angles. The electronic spectrum of 90 degrees twisted bilayer black phosphorus is found to be x-y isotropic in contrast to the monolayer. However x-y anisotropy, and a partial return to monolayer-like behavior, particularly in the valence band, can be induced by an external out-of-plane electric field. Moreover, the preferred hole effective mass can be rotated by 90 degrees simply by changing the direction of the applied electric field. In particular, a +0.4 (-0.4) V angstrom(1) out-of-plane electric field results in a similar to 60% increase in the hole effective mass along the y (x) axis and enhances the m(y)*/m(x)* (m(x)*/m(y)*) ratio as much as by a factor of 40. Our DFT and (k) over right arrow . (p) over right arrow simulations clearly indicate that the twist angle in combination with an appropriate gate voltage is a novel way to tune the electronic and optical properties of bilayer phosphorus and it gives us a new degree of freedom to engineer the properties of black phosphorus based devices.

Abstract: <script type='text/javascript'>document.write(unpmarked('Quantum electrodynamics predicts that heavy atoms (Z \u003E Z(c) approximate to 170) will undergo the process of atomic collapse where electrons sink into the positron continuum and a new family of so-called collapsing states emerges. The relativistic electrons in graphene exhibit the same physics but at a much lower critical charge (Z(c) approximate to 1) which has made it possible to confirm this phenomenon experimentally. However, there exist conflicting predictions on the effect of a magnetic field on atomic collapse. These theoretical predictions are based on the continuum Dirac-Weyl equation, which does not have an exact analytical solution for the interplay of a supercritical Coulomb potential and the magnetic field. Approximative solutions have been proposed, but because the two effects compete on similar energy scales, the theoretical treatment varies depending on the regime which is being considered. These limitations are overcome here by starting from a tight-binding approach and computing exact numerical results. By avoiding special limit cases, we found a smooth evolution between the different regimes. We predict that the atomic collapse effect persists even after the magnetic field is activated and that the critical charge remains unchanged. We show that the atomic collapse regime is characterized: (1) by a series of Landau level anticrossings and (2) by the absence of root B scaling of the Landau levels with regard to magnetic field strength.'));

Abstract: Dielectric properties of nano-confined water are important in several areas of science, i.e., it is relevant in the dielectric double layer that exists in practically all heterogeneous fluid-based systems. Molecular dynamics simulations are used to predict the in-plane dielectric properties of confined water in planar channels of width ranging from sub-nanometer to bulk. Because of suppressed rotational degrees of freedom near the confining walls, the dipole of the water molecules tends to be aligned parallel to the walls, which results in a strongly enhanced in-plane dielectric constant (epsilon (parallel to)) reaching values of about 120 for channels with height 8 angstrom < h < 10 angstrom. With the increase in the width of the channel, we predict that epsilon (parallel to) decreases nonlinearly and reaches the bulk value for h > 70 angstrom. A stratified continuum model is proposed that reproduces the h > 10 angstrom dependence of epsilon (parallel to). For sub-nanometer height channels, abnormal behavior of epsilon (parallel to) is found with two orders of magnitude reduction of epsilon (parallel to) around h similar to 7.5 angstrom, which is attributed to the formation of a particular ice phase that exhibits long-time (similar to mu s) stable ferroelectricity. This is of particular importance for the understanding of the influence of confined water on the functioning of biological systems.

Abstract: Few-electron systems confined in two-dimensional parabolic quantum dots at high magnetic fields are studied by the Hartree-Fock (HF) and exact diagoiialization methods. A generalized multicenter Gaussian basis is proposed in the HF method. A comparison of the HF and exact, results allows as to discuss the relevance of the symmetry of the charge density distribution for the accuracy of the HF method. It is shown that the energy estimates obtained with the broken-symmetry HF wave functions become exact in the infinite magnetic-field limit. In this limit the charge density of the broken-symmetry solution call be identified with the classical charge distribution.

Abstract: We investigate the ground and excited states of a bipolar artificial molecule composed of two vertically coupled quantum dots containing different type of carriers-electrons and holes-in equilibrium. The approach based on exact diagonalization is used and reveals an intricate pattern of ground-state angular momentum switching and a rearrangement of approximate single-particle levels as a function of the interdot coupling strength.

Abstract: Using Brownian dynamics simulations, we investigate the dynamics of colloids confined in two-dimensional narrow channels driven by a nonuniform force Fdr(y). We considered linear-gradient, parabolic, and deltalike driving-force profiles. This driving force induces melting of the colloidal solid (i.e., shear-induced melting), and the colloidal motion experiences a transition from elastic to plastic regime with increasing Fdr. For intermediate Fdr (i.e., in the transition region) the response of the system, i.e., the distribution of the velocities of the colloidal chains õi(y), in general does not coincide with the profile of the driving force Fdr(y), and depends on the magnitude of Fdr, the width of the channel, and the density of colloids. For example, we show that the onset of plasticity is first observed near the boundaries while the motion in the central region is elastic. This is explained by: (i) (in)commensurability between the chains due to the larger density of colloids near the boundaries, and (ii) the gradient in Fdr. Our study provides a deeper understanding of the dynamics of colloids in channels and could be accessed in experiments on colloids (or in dusty plasma) with, e.g., asymmetric channels or in the presence of a gradient potential field.

Abstract: We study nonlinear flux dynamics in a current-carrying type-I superconductor. The stray magnetic field of the current induces the intermediate state, where nucleation of flux domains is discretized to a single fluxoid at a time, while their final shape (tubular or laminar), size, and nucleation rate depend on applied current and edge conditions. The current induces opposite flux domains on opposite sides of the sample, and subsequently drives them to annihilation-which is also discretized, as a sequence of vortex-antivortex pairs. The discretization of both nucleation and annihilation leaves measurable traces in the voltage across the sample and in locally probed magnetization. The reported dynamic phenomena thus provide an unambiguous proof of a flux quantum being the smallest building block of the intermediate state in type-I superconductors.

Abstract: Moire patterns in the pseudo-magnetic field and in the strain profile of graphene (GE) when put on top of a hexagonal lattice substrate are predicted from elasticity theory. The van der Waals interaction between GE and the substrate induces out-of-plane deformations in graphene which results in a strain field, and consequently in a pseudo-magnetic field. When the misorientation angle is about 0.5 degrees, a three-fold symmetric strain field is realized that results in a pseudo-magnetic field very similar to the one proposed by F. Guinea, M. I. Katsnelson, and A. K. Geim [Nature Phys. 6, 30 (2010)]. Our results show that the periodicity and length of the pseudo-magnetic field can be tuned in GE by changing the misorientation angle and substrate adhesion parameters and a considerable energy gap (23 meV) can be obtained due to out-of-plane deformation of graphene which is in the range of recent experimental measurements (20-30 meV). (C) 2014 AIP Publishing LLC.

Abstract: A combination of infrared spectroscopy and magnetotransport is used to investigate the impurity band and the magnetic-field-induced metal-insulator transition in n-type GaAs/AlxGa1-xAs superlattices. The dropping of the Fermi level from the conduction band into the impurity band upon increasing magnetic field is observed in a sample doped to n=4n(c), where n(c) is the critical density according to the Mott criterion. The metal-insulator transition takes place while the Fermi level is in the impurity band, with no qualitative change from the metallic to the insulating side. Due to the anisotropy of the superlattice band structure, the metal-insulator transition is shifted to higher magnetic field, when the magnetic field is tilted away from the growth axis towards the layer planes.

Abstract: Epitaxial graphene is grown on a non-polar n(+) 6H-SiC m-plane substrate and studied using atomic scale scanning tunneling microscopy. Multilayer graphene is found throughout the surface and exhibits rotational disorder. Moire patterns of different spatial periodicities are found, and we found that as the wavelength increases, so does the amplitude of the modulations. This relationship reveals information about the interplay between the energy required to bend graphene and the interaction energy, i.e. van der Waals energy, with the graphene layer below. Our experiments are supported by theoretical calculations which predict that the membrane topographical amplitude scales with the Moire pattern wavelength, L as L-1 + alpha L-2. (C) 2014 Elsevier Ltd. All rights reserved.

Abstract: The instability of a free-standing one-sided hydrogenated/fluorinated graphene nanoribbon, i.e., graphone/fluorographene, is studied using ab initio, semiempirical, and large-scale molecular dynamics simulations. Free-standing semi-infinite armchairlike hydrogenated/fluorinated graphene (AC-GH/AC-GF) and boatlike hydrogenated/fluorinated graphene (B-GH/B-GF) (nanoribbons which are periodic along the zigzag direction) are unstable and spontaneously transform into spiral structures. We find that rolled, spiral B-GH and B-GF are energetically more favorable than spiral AC-GH and AC-GF which is opposite to the double-sided flat hydrogenated/fluorinated graphene, i.e., graphane/fluorographene. We found that the packed, spiral structures exhibit an unexpected localized highest occupied molecular orbital and lowest occupied molecular orbital at the edges with increasing energy gap during rolling. These rolled hydrocarbon structures are stable beyond room temperature up to at least T = 1000 K within our simulation time of 1 ns. DOI: 10.1103/PhysRevB.87.075448