Abstract: Two-dimensional materials composed of transition metal carbides and nitrides (MXenes) are poised to revolutionize energy conversion and storage. In this work, we used density functional theory (DFT) to investigate the adsorption of Mg and Na adatoms on five M2CS2 monolayers (where M = Mo, Nb, Ti, V, and Zr) for battery applications. We assessed the stability of the adatom (i.e. Na and Mg)-monolayer systems by calculating adsorption and formation energies, as well as voltages as a function of surface coverage. For instance, we found that Mo2CS2 cannot support a full layer of Na nor even a single Mg atom. Na and Mg exhibit the strongest binding on Zr2CS2, followed by Ti2CS2, Nb2CS2 and V2CS2. Using the nudged elastic band method (NEB), we computed promising diffusion barriers for both dilute and nearly full ion surface coverage cases. In the dilute ion adsorption case, a single Mg and Na atom on Ti2CS2 experience similar to 0.47 eV and similar to 0.10 eV diffusion barriers between the lowest energy sites, respectively. For a nearly full surface coverage, a Na ion moving on Ti2CS2 experiences a similar to 0.33 eV energy barrier, implying a concentration-dependent diffusion barrier. Our molecular dynamics results indicate that the three (one) layers (layer) of the Mg (Na) ion on both surfaces of Ti2CS2 remain stable at T = 300 K. While, according to voltage calculations, Zr2CS2 can store Na up to three atomic layers, our MD simulations predict that the outermost layers detach from the Zr2CS2 monolayer due to the weak interaction between Na ions and the monolayer. This suggests that MD simulations are essential to confirm the stability of an ion-electrode system – an insight that is mostly absent in previous studies.

Abstract: We report the formation of nanobubbles on graphene with a radius of the order of 1 nm, using ultralow energy implantation of noble gas ions (He, Ne, Ar) into graphene grown on a Pt(111) surface. We show that the universal scaling of the aspect ratio, which has previously been established for larger bubbles, breaks down when the bubble radius approaches 1 nm, resulting in much larger aspect ratios. Moreover, we observe that the bubble stability and aspect ratio depend on the substrate onto which the graphene is grown (bubbles are stable for Pt but not for Cu) and trapped element. We interpret these dependencies in terms of the atomic compressibility of the noble gas as well as of the adhesion energies between graphene, the substrate, and trapped atoms.

Abstract: The precise control of skyrmionics and magnonics in magnetic materials is key to the development of novel spin-based technology and information transport applications. Essentially, the inherent stability of magnetic skyrmions (provided by their topological charge) together with their extremely small size (down to a few nanometers) and the ultralow threshold current necessary to move them in nanostructures are the main advantages of skyrmionics. Not least, magnonics offers lower power consumption compared to electronics and the excitation of high frequency (sub-100~nm wavelength) magnons makes it possible for the creation of nanometric devices for ultrafast information transport. Even though extensive research has been carried out in recent years, the precise manipulation of skyrmions and spin waves (magnons) in nanostructures is not fully mastered and needs to be addressed before making functional skyrmionic and magnonic devices. In this thesis, we reveal multiple alternatives for the manipulation of skyrmions and spin-waves in different materials, such as bulk chiral magnets, heterochiral structures, magnet-supperconductor hybrids and two-dimentional magnetic materials. We make use of a multiscale model to numerically simulate the magnetic states at each considered material, from micromagnetic to atomistic control. We first explore the different nucleation mechanisms, activation energy, and the time evolution of the skyrmion formation in chiral magnetic films, crucial for the realization of skyrmion-based devices. We show that the skyrmion lattice is formed from the conical phase progressively, most probably by the formation of chiral bobbres, followed by the cylindrical growth of individual skyrmions from the film surface. That reflects a rod-like (one-dimensional) nucleation of the skyrmion phase, with an activation barrier of several electronvolts per skyrmion for the case of MnSi (Manganese monosilicide). In addition, we reveal the interesting blinking (creation-annihilation) behavior of skyrmions close to the phase boundary between the conical and skyrmion phases, where we recall that such switching between topologically distinct states has been proposed as a bit operation for information storage. Next, we discuss the motion of ferromagnetic and antiferromagnetic skyrmions in heterochiral magnets. We report the characteristic deflection of ferromagnetic skyrmions when moving across a heterochiral interface, where the extent of such deflection is tuned by the applied spin-polarized current and the magnitude of Dzyaloshinskii-Moriya interaction. Following, we show that the antiferromagnetic skyrmion achieves much higher velocities than its ferromagnetic counterpart, yet experiences far stronger confinement in nanoengineered heterochiral tracks, which reinforces antiferromagnetic skyrmions as a favorable choice for skyrmion-based devices. After that, we study the interesting coupling of magnetic skyrmions and superconducting vortices in magnet-superconductor heterostructures. We perform numerical simulations, based on experimental observations, to demonstrate that the stray field of magnetic skyrmions can nucleate antivortices in an adjacent superconducting film, giving rise to a hybrid topological object, the skyrmion-vortex pair, which harbor promising features for skyrmionics and quantum computing applications. We then explore the manipulation of a single skyrmion-vortex pair when currents are applied into both superconducting and magnetic parts of the heterostructure, which is of importance for the facilitated skyrmion guidance in racetrack applications. Afterwards, we make use of the high tunability of magnetic parameters in two-dimensional magnetic materials to reveal the rich phase diagram of exotic magnetic configurations in magnetic monolayers with suppressed nearest-neighbour exchange, where we show that several unique cycloidal, checkerboard, row-wise and spin-ice states are stabilized by the competition between the second-nearest-neighbor exchange, Dzyaloshinskii-Moriya, and dipolar interactions. Additionally, we show the coexistence of ferromagnetic and antiferromagnetic spin-cycloids, as well as novel types of skyrmions and chiral domain walls. Finally, in the last part of the thesis, we present the spin wave properties in the two-dimensional magnetic materials CrBr$3$ and CrI$3$. Using spin-dynamics simulations parametrized from first principles, we reveal that the spin wave dispersion in such materials can be tuned in a broad range of frequencies by strain-engineering, and that a designed pattern of strain, as well as structural defects (halide vacancies) can be turned useful in the design of spin-wave guides. Lastly, we discuss the realization of magnonic crystals by moiré-periodic modulation of magnetic parameters in van der Waals heterostructures, where we show that the several nanometer small periodicities in such samples are ideal for the interference of terahertz spin waves. Recalling the wide range of possibilities for manipulating spin waves in such two-dimensional materials, we therefore suggest these systems as a front-runner for prospective terahertz magnonic applications.

Abstract: The interface-dependent electronic, vibrational, piezoelectric, and optical properties of van der Waals heterobilayers, formed by buckled GeO (b-GeO) and Janus MoSO structures, are investigated by means of first-principles calculations. The electronic band dispersions show that O/Ge and S/O interface formations result in a type-II band alignment with direct and indirect band gaps, respectively. In contrast, O/O and S/Ge interfaces give rise to the formation of a type-I band alignment with an indirect band gap. By considering the Bethe-Salpeter equation (BSE) on top of G(0)W(0) approximation, it is shown that different interfaces can be distinguished from each other by means of the optical absorption spectra as a consequence of the band alignments. Additionally, the low- and high-frequency regimes of the Raman spectra are also different for each interface type. The alignment of the individual dipoles, which is interface-dependent, either weakens or strengthens the net dipole of the heterobilayers and results in tunable piezoelectric coefficients. The results indicate that the possible heterobilayers of b-GeO/MoSO asymmetric structures possess various electronic, optical, and piezoelectric properties arising from the different interface formations and can be distinguished by means of various spectroscopic techniques.

Abstract: By means of density functional theory-based first-principle calculations, the structural, vibrational and electronic properties of single-layer Ge3N4 are investigated. Structural optimizations and phonon band dispersions reveal that single-layer ultrathin form of Ge3N4 possesses a dynamically stable buckled structure with large hexagonal holes. Predicted Raman spectrum of single-layer Ge3N4 indicates that the buckled holey structure of the material exhibits distinctive vibrational features. Electronic band dispersion calculations indicate the indirect band gap semiconducting nature of single-layer Ge3N4. It is also proposed that single-layer Ge3N4 forms type-II vertical heterostructures with various planar and puckered 2D materials except for single-layer GeSe which gives rise to a type-I band alignment. Moreover, the electronic properties of single-layer Ge3N4 are investigated under applied external in-plane strain. It is shown that while the indirect gap behavior of Ge3N4 is unchanged by the applied strain, the energy band gap increases (decreases) with tensile (compressive) strain.

Abstract: Two-dimensional crystals with angstrom-scale pores are widely considered as candidates for a next generation of molecular separation technologies aiming to provide extreme, exponentially large selectivity combined with high flow rates. No such pores have been demonstrated experimentally. Here we study gas transport through individual graphene pores created by low intensity exposure to low kV electrons. Helium and hydrogen permeate easily through these pores whereas larger species such as xenon and methane are practically blocked. Permeating gases experience activation barriers that increase quadratically with molecules' kinetic diameter, and the effective diameter of the created pores is estimated as similar to 2 angstroms, about one missing carbon ring. Our work reveals stringent conditions for achieving the long sought-after exponential selectivity using porous two-dimensional membranes and suggests limits on their possible performance. Two-dimensional membranes with angstrom-sized pores are predicted to combine high permeability with exceptional selectivity, but experimental demonstration has been challenging. Here the authors realize angstrom-sized pores in monolayer graphene and demonstrate gas transport with activation barriers increasing quadratically with the molecular kinetic diameter.

Abstract: The modification of the quantum states in a Dice lattice due to a Coulomb impurity are investigated. The energy-band structure of a pristine Dice lattice consists of a Dirac cone and a flat band at the Dirac point. We use the tight-binding formalism and find that the flat band states transform into a set of discrete bound states whose electron density is localized on a ring around the impurity mainly on two of the three sublattices. Its energy is proportional to the strength of the Coulomb impurity. Beyond a critical strength of the Coulomb potential atomic collapse states appear that have some similarity with those found in graphene with the difference that the flat band states contribute with an additional ringlike electron density that is spatially decoupled from the atomic collapse part. At large value of the strength of the Coulomb impurity the flat band bound states anticross with the atomic collapse states.

Abstract: We investigate the torque field and skyrmion motion at an interface between a ferromagnet hosting a skyrmion and a material with a strong spin-orbit interaction. We analyze both semiconductor materials and topological insulators using a Hamiltonian model that includes a linear term. The spin torque-inducing current is considered to flow in the single band limit; therefore, a quantum model of current is used. Skyrmion motion due to spin transfer torque proves to be more difficult in the presence of a spin-orbit interaction in the case where only interface in-plane currents are present. However, edge effects in narrow nanowires can be used to drive the skyrmion motion and to exert a limited control on its motion direction. We also show the differences and similarities between torque fields due to electric current in the many and single band limits. Published under an exclusive license by AIP Publishing.

Abstract: A twin boundary in one of the layers of a twisted van der Waals heterostructure separates regions with near opposite inter-layer twist angles. In a MoS₂/WSe₂bilayer, the regions with<inline-formula><tex-math><?CDATA $Rh^h$?></tex-math><math overflow=“scroll”><msubsup><mi>R</mi><mi>h</mi><mi>h</mi></msubsup></math><inline-graphic href=“tdmac529dieqn1.gif” type=“simple” /></inline-formula>and<inline-formula><tex-math><?CDATA $Rh^X$?></tex-math><math overflow=“scroll”><msubsup><mi>R</mi><mi>h</mi><mi>X</mi></msubsup></math><inline-graphic href=“tdmac529dieqn2.gif” type=“simple” /></inline-formula>stacking registry that defined the sub-lattices of the moiré honeycomb pattern would be mirror-reflected across such a twist boundary. In that case, we demonstrate that topologically protected chiral moiré exciton states are confined at the twist boundary. These are one-dimensional and uni-directional excitons with opposite velocities for excitons composed by electronic states with opposite valley/spin character, enabling intrinsic, guided, and far reaching valley-polarized exciton currents.

Abstract: Bare chemical vapor deposition (CVD) grown graphene (GRP) was anisotropically etched with various etching parameters. The morphological and structural characterizations were carried out by optical microscopy and the vibrational properties substrates were obtained by Raman spectroscopy. The ammonia adsorption and desorption behavior of graphene-based sensors were recorded via quartz crystal microbalance (QCM) measurements at room temperature. The etched samples for ambient NH3 exhibited nearly 35% improvement and showed high resistance to humidity molecules when compared to bare graphene. Besides exhibiting promising sensitivity to NH3 molecules, the etched graphene-based sensors were less affected by humidity. The experimental results were collaborated by Density Functional Theory (DFT) calculations and it was shown that while water molecules fragmented into H and O, NH3 interacts weakly with EGPR2 sample which reveals the enhanced sensing ability of EGPR2. Apparently, it would be more suitable to use EGRP2 in sensing applications due to its sensitivity to NH3 molecules, its stability, and its resistance to H2O molecules in humid ambient.

Abstract: Molecular dynamics simulations are used to investigate the effect of an AFM tip when indenting graphene nano bubbles filled by a noble gas (i.e. He, Ne and Ar) up to the breaking point. The failure points resemble those of viral shells as described by the Foppl-von Karman (FvK) dimensionless number defined in the context of elasticity theory of thin shells. At room temperature, He gas inside the bubbles is found to be in the liquid state while Ne and Ar atoms are in the solid state although the pressure inside the nano bubble is below the melting pressure of the bulk. The trapped gases are under higher hydrostatic pressure at low temperatures than at room temperature.

Abstract: The alpha-T3 lattice model interpolates a honeycomb (graphene-like) lattice and a T3 (also known as dice) lattice via the parameter alpha. These lattices are made up of three atoms per unit cell. This gives rise to an additional dispersionless flat band touching the conduction and valence bands. Electrons in this model are analogous to Dirac fermions with an enlarged pseudospin, which provides unusual tunneling features like omnidirectional Klein tunneling, also called super-Klein tunneling (SKT). However, it is unknown how small deviations in the equivalence between the atomic sites, i.e., variations in the alpha parameter, and the number of tunnel barriers changes the transmission properties. Moreover, it is interesting to learn how tunneling occurs through regions where the energy spectrum changes from linear with a middle flat band to a hyperbolic dispersion. In this paper we investigate these properties, its dependence on the number of square barriers and the alpha parameter for either gapped and gapless cases. Furthermore, we compare these results to the case where electrons tunnel from a region with linear dispersion to a region with a bandgap. In the latter case, contrary to tunneling through a potential barrier, the SKT is no longer observed. Finally, we find specific cases where transmission is allowed due to a symmetry breaking of sublattice equivalence.

Abstract: Using the continuum model, we investigate the electronic properties of two types of bilayer graphene (BLG) quantum ring (QR) geometries: (i) An isolated BLG QR and (ii) a monolayer graphene (MLG) with a QR put on top of an infinite graphene sheet (hybrid BLG QR). Solving the Dirac-Weyl equation in the presence of a perpendicular magnetic field and applying the infinite mass boundary condition at the ring boundaries, we obtain analytical results for the energy levels and corresponding wave spinors for both structures. In the case of isolated BLG QR, we observe a sizable and magnetically tunable band gap which agrees with the tight-binding transport simulations. Our analytical results also show the intervalley symmetry EeK (m) = ???EK??? h (m) between the electron (e) and the hole (h) states (m is the angular momentum quantum number) for the energy spectrum of the isolated BLG QR. The presence of interface boundary in a hybrid BLG QR modifies drastically the energy levels as compared with that of an isolated BLG QR. Its energy levels are tunable from MLG dot to isolated BLG QR and to MLG Landau energy levels as the magnetic field is varied. Our predictions can be verified experimentally using different techniques such as by magnetotransport measurements.

Abstract: This piece of work is twofold. First, the time evolution of wave-packets in 2D systems is analyzed by the Split-Operator technique in three different scenarios: in multilayer phosphorene, the transient oscillations in the time-dependent average of position and momentum were observed due to the zitterbewegung effect, and the wave packet propagates non-uniformly along the space deforming itself into an elliptical shape. These results were corroborated by the Green’s function formalism except for large values of the wave-vector and long times; in 2D semiconductor quantum wires (QWs) with anisotropic effective masses and different angle orientations with respect to the anisotropic axis. We have shown that the greater this angle, the smaller is the energy levels spacing implying in an increase of the accessible electronic states. Additionally, for non-null magnetic field, the quantum Hall edge states are significantly affected by the edge orientation. In the anisotropic case damped oscillations in the average values of velocity in both x and y directions where obtained. Theses oscillations are originated by the QW geometry but also from subwavepackets with different momentum orientations, whereas for isotropic QWs the wavepacket disperses without splitting; in the third scenario the split-operator technique was used to study the Landau levels, the wave packet trajectories and velocities of electrons in graphene at low-energy regime described by a modified Dirac equation where the momentum-operator is written in a generalized form as result of applying the position-dependent translation operator formalism (PDTO). In the second part of this thesis, the electronic and tunneling properties of α − T3 lattices were studied. Electrons in these lattices behave analogous to integer-spin Dirac Fermions. The presence of a third atomic site in the unit cell leads to a flat band in the energy spectrum, providing unique electronic and tunneling properties. The presence of a super-periodic potential and the inclusion of symmetry-breaking terms results in deviations of the atomic equivalence between the atomic sites affecting the Dirac points and the band-gap. Small deviations in the equivalence between the atomic sites and the number of barriers change the transmission properties in these lattices. Additionally, new tunneling regions are possible by adjusting the symmetry between the atomic sites and affect the omnidirectional total transmission called super-Klein tunneling observed in these lattices. We compare those results to the tunneling probabilities through regions where the energy spectrum changes from linear with a middle flat band to a hyperbolic dispersion.