Abstract: Within the Ginzburg-Landau theory, we study the vortex structures in three-dimensional anisotropic mesoscopic superconductors in the presence of a uniform magnetic field. Anisotropy is included through varied Tc in different layers of the sample and leads to distinct differences in the vortex states and their free energy. Several unconventional states are found, some comprising vortex clusters or exhibiting asymmetry. In a tilted magnetic field, we found second-order transitions between different vortex states, although vortex entry is generally a first-order transition in mesoscopic samples. In multilayered samples the kinked vortex strings are formed owing to the competing interactions of vortices with Meissner currents and the weak-link boundaries. The length and deformation of vortex fragments are determined solely by the inclination and strength of applied magnetic field, and this lock-in does not depend on the degree of anisotropy between the superconducting layers.

Abstract: The superconducting state. of a thin superconducting disk with a hole is studied within the, nonlinear Ginzburg-Landau theory in which the demagnetization effect is accurately taken into account. We find that the flux through the hole is not quantized, the superconducting state is stabilized with increasing size of the hole for fixed radius of the disk, and a transition to a multivortex state is found if the disk is sufficiently large. Breaking the circular symmetry through a non-central-location of the hole in the disk favors the multivortex state.

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 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: Using the tight-binding approach, we calculate the electronic structure of triangular and hexagonal MoS2 quantum dots. Due to the orbital asymmetry we show that it is possible to form quantum dots with the same shape but having different electronic properties. The electronic states of triangular and hexagonal quantum dots are explored, as well as the local and total density of states and the convergence towards the bulk spectrum with dot size is investigated. Our calculations show that: (1) edge states appear in the band gap, (2) that there are a larger number of electronic states in the conduction band as compared to the valence band, and (3) the relative number of edge states decreases with increasing dot size.

Abstract: Stable monolayer materials based on existing, well known and stable two-dimensional crystal fluorographene and molybdenum disulfide are predicted to exhibit a huge magnetocrystalline anisotropy when functionalized with adsorbed transition metal atoms at vacant sides. Ab initio calculations within the density-functional theory formalism were performed to investigate the adsorption of the transitional metals in a single S (or F) vacancy of monolayer molybdenum disulfide (or fluorographene). We found strong bonding of the transitional metal atoms to the vacant sites with binding energies ranging from 2.5 to 5.2 eV. Our calculations revealed that these systems with adsorbed metal atoms exhibit a magnetic anisotropy, specifically the structures including Os and Ir show a giant magnetocrystalline anisotropy energy of 31-101 meV. Our results demonstrate the possibility of obtaining stable monolayer materials with huge magnetocrystalline anisotropy based on preexisting, well known and stable two-dimensional crystals: fluorographene and molybdenum disulfide. We believe that the results obtained here are useful not only for deeper understanding of the origin of magnetocrystalline anisotropy but also for the design of monolayer optoelectronic devices with novel functionalities.

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