Abstract: The time evolution of a wave packet in strained graphene is studied within the tight-binding model and continuum model. The effect of an external magnetic field, as well as a strain-induced pseudomagnetic field, on the wave-packet trajectories and zitterbewegung are analyzed. Combining the effects of strain with those of an external magnetic field produces an effective magnetic field which is large in one of the Dirac cones, but can be practically zero in the other. We construct an efficient valley filter, where for a propagating incoming wave packet consisting of momenta around the K and K' Dirac points, the outgoing wave packet exhibits momenta in only one of these Dirac points while the components of the packet that belong to the other Dirac point are reflected due to the Lorentz force. We also found that the zitterbewegung is permanent in time in the presence of either external or strain-induced magnetic fields, but when both the external and strain-induced magnetic fields are present, the zitterbewegung is transient in one of the Dirac cones, whereas in the other cone the wave packet exhibits permanent spatial oscillations.

Abstract: We investigate the dynamics of a charged particle moving in a graphene layer and in a two-dimensional electron gas, where it obeys the Dirac and the Schrödinger equations, respectively. The charge carriers are described as Gaussian wavepackets. The dynamics of the wavepackets is studied numerically by solving both quantum-mechanical and relativistic equations of motion. The scattering of such wavepackets by step-like magnetic and potential barriers is analysed for different values of wavepacket energy and width. We find: (1) that the average position of the wavepacket does not coincide with the classical trajectory, and (2) that, for slanted incidence, the path of the centre of mass of the wavepacket does not have to penetrate the barrier during the scattering process. Trembling motion of the charged particle in graphene is observed in the absence of an external magnetic field and can be enhanced by a substrate-induced mass term.

Abstract: We consider a classical system of two-dimensional (2D) charged particles, interacting through a repulsive Yukawa potential exp(-r/λ)/r, and confined in a parabolic channel that limits the motion of the particles in the y direction. Along the x direction, the particles are subject to a periodic potential. The ground-state configurations and the normal-mode spectra of the system are obtained as a function of the periodicity and strength of the periodic potential (V0) and density. An interesting set of tunable ground-state configurations are found, with first- or second-order structural transitions between them. A configuration with particles aligned, perpendicular to the x direction, in each minimum of the periodic potential is obtained for V0 larger than some critical value that has a power-law dependence on the density. The phonon spectrum of different configurations was also calculated. A localization of the modes into a small frequency interval is observed for sufficiently large strength of the periodic potential, and a tunable gap in the phonon spectrum is found as a function of V0.

Abstract: The properties of a new class of phosphorene nanotubes (PNT) are investigated by performing first-principles calculations. We demonstrate that it is advantageous to use blue phosphorene in order to make small nanotubes and propose a way to create low-energy PNTs by the inclusion of defect lines. Five different types of defect lines are investigated and incorporated in various combinations. The resulting defect-induced faceted PNTs have negligible bending stresses which leads to a reduction in the formation energy with respect to round PNTs. Our armchair faceted PNTs have similar formation energies than the recently proposed multiphase faceted PNTs, but they have a larger variety of possible structures. Our zigzag faceted PNTs have lower formation energies than round tubes and multiphase faceted nanotubes. The electronic properties of the defect-induced faceted PNTs are determined by the defect lines which control the band gap and the shape of the electronic states at the band edges. These band gaps increase with the radius of the nanotubes and converge to those of isolated defect lines.

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: Water confined between two graphene layers with a separation of a few A forms a layered two-dimensional ice structure. Using large scale molecular dynamics simulations with the adoptable ReaxFF interatomic potential we found that flat monolayer ice with a rhombic-square structure nucleates between the graphene layers which is nonpolar and nonferroelectric. We provide different energetic considerations and H-bonding results that explain the interlayer and intralayer properties of two-dimensional ice. The controversial AA stacking found experimentally [Algara-Siller et al., Nature (London) 519, 443 (2015)] is consistent with our minimum-energy crystal structure of bilayer ice. Furthermore, we predict that an odd number of layers of ice has the same lattice structure as monolayer ice, while an even number of ice layers exhibits the square ice AA stacking of bilayer ice.

Abstract: The structural transitions of the ground state of a system of repulsively interacting particles confined in a quasi-one-dimensional channel, and the effect of the interparticle interaction as well as the functional form of the confinement potential on those transitions are investigated. Although the nonsequential ordering of transitions (non-SOT), i.e., the 1 – 2 – 4 – 3 – 4 – 5 – 6 – ... sequence of chain configurations with increasing density, is widely robust as predicted in a number of theoretical studies, the sequential ordering of transitions (SOT), i.e., the 1 – 2 – 3 – 4 – 5 – 6 – ... chain, is found as the ground state for long-ranged interparticle interaction and hard-wall-like confinement potentials. We found an energy barrier between every two different phases around its transition point, which plays an important role in the preference of the system to follow either a SOT or a non-SOT. However, that preferential transition requires also the stability of the phases during the transition. Additionally, we analyze the effect of a small structural disorder on the transition between the two phases around its transition point. Our results show that a small deformation of the triangular structure changes dramatically the picture of the transition between two phases, removing in a considerable region the non-SOT in the system. This feature could explain the fact that the non-SOT is, up to now, not observed in experimental systems, and suggests a more advanced experimental setup to detect the non-SOT.

Abstract: We use the anisotropic time-dependent Ginzburg-Landau theory to investigate the effect of a square array of out-of-plane magnetic dots on the dynamics of Josephson vortices (fluxons) in artificial stacks of superconducting-normal-superconducting (SNS) Josephson junctions in the presence of external DC and AC currents. Periodic pinning due to the magnetic dots distorts the triangular lattice of fluxons and results in the appearance of commensurability features in the current-voltage characteristics of the system. For the larger values of the magnetization, additional peaks appear in the voltage-time characteristics of the system due to the creation and annihilation of vortex-antivortex pairs. Peculiar changes in the response of the system to the applied current is found resulting in a “superradiant” vortex-flow state at large current values, where a rectangular lattice of moving vortices is formed. Synchronizing the motion of fluxons by adding a small ac component to the biasing dc current is realized. However, we found that synchronization becomes difficult for large magnetization of the dots due to the formation of vortex-antivortex pairs.

Abstract: Within a tight binding approach we investigate the energy levels of hexagonal and triangular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We study AA- and AB-(Bernal) stacked BLG QDs and obtain the energy levels in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). Our results show that the size dependence of the energy levels is different from that of monolayer graphene QDs. The energy spectrum of AB-stacked BLG QDs with zigzag edges exhibits edge states which spread out into the opened energy gap in the presence of a perpendicular electric field. We found that the behavior of these edges states is different for the hexagonal and triangular geometries. In the case of AA-stacked BLG QDs, the electron and hole energy levels cross each other in both cases of armchair and zigzag edges as the dot size or the applied bias increases.

Abstract: Using density functional theory, we study the electronic properties of a graphene-silicene bilayer (GSB). A single layer of silicene binds to the graphene layer with adhesion energy of about 25 meV/atom. This adhesion energy between the two layers follows accurately the well-known -1/z(2) dispersion energy as found between two infinite parallel plates. In small flakes of GSB with hydrogenated edges, negative charge is transferred from the graphene layer to the silicene layer, producing a permanent and a switchable polar bilayer, while in an infinite GSB, the negative charge is transferred from the silicene layer to the graphene layer. The graphene-silicene bilayer is a good candidate for a nanocapacitor with piezoelectric capabilities. We found that the permanent dipole of the bilayer can be tuned by an external perpendicular electric field.

Abstract: The Majorana zero modes and persistent spin current in mesoscopic d-wave-superconducting loops with spin-orbit (SO) interaction are investigated by numerically solving the spin-generalized Bogoliubov-de Gennes equations self-consistently. For some appropriate strength of the SO coupling, Majorana zero-energy states and sharp jumps of the spin-polarized currents can be observed when the highest energy levels cross the Fermi energy in the spectrum, leading to spin currents with opposite chirality flowing near the inner and outer edges of the sample. When the threaded magnetic flux turns on, four flux-dependent patterns of the persistent spin current with step-like features show up, accompanied by Majorana edge modes at flux values where the energy gap closes. Moreover, the Majorana zero mode is highly influenced by the direction of the Zeeman field. A finite in-plane field can lead to the gap opening since the inversion symmetry is broken. Remarkably, multiple Majorana zero-energy states occur in the presence of an out-of-plane field h(z), and the number of steps in the spin current evolution can be effectively tuned by the field strength due to the shift of Majorana zero modes. Finally, when the loop sample contains surface indentation defects, zero-energy modes can always show up in the presence of an appropriate h(z). Interestingly, multiple Majorana states may be present in the system with a corner defect even if h(z) = 0.

Abstract: In order to study quantum effects in a two-dimensional crystal lattice of a free-standing monolayer graphene, we have performed both path-integral Monte Carlo (PIMC) and classical Monte Carlo (MC) simulations for temperatures up to 2000 K. The REBO potential is used for the interatomic interaction. The total energy, interatomic distance, root-mean-square displacement of the atom vibrations, and the free energy of the graphene layer are calculated. The obtained lattice vibrational energy per atom from the classical MC simulation is very close to the energy of a three-dimensional harmonic oscillator 3k(B)T. The PIMC simulation shows that quantum effects due to zero-point vibrations are significant for temperatures T < 1000 K. The quantum contribution to the lattice vibrational energy becomes larger than that of the classical lattice for T < 400 K. The lattice expansion due to the zero-point motion causes an increase of 0.53% in the lattice parameter. A minimum in the lattice parameter appears at T similar or equal to 500 K. Quantum effects on the atomic vibration amplitude of the graphene lattice and its free energy are investigated.