Source author record

Gerardo G. Naumis

Gerardo G. Naumis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

Why the first magic-angle is different from others in twisted graphene bilayers: interlayer currents, kinetic and confinement energy and wavefunction localization

The chiral Hamiltonian for twisted graphene bilayers is analyzed in terms of its squared Hamiltonian which removes the particle-hole symmetry and thus one bipartite lattice, allowing to write the Hamiltonian in terms of a $2\times 2$ matrix. This brings to the front the three main physical actors of twisted systems: kinetic energy, confinement potential and an interlayer interaction operator which is divided in two parts: a non-Abelian interlayer operator and an operator which contains an interaction energy between layers. Here, each of these components is analyzed as a function of the angle of rotation, as well as in terms of the wave-function localization properties. In particular, it is proved that the non-Abelian operator represents interlayer currents between each layer triangular sublattices, i.e., a second-neighbor interlayer current between bipartite sublattices. A crossover is seen between such contributions and thus the first magic angle is different from other higher order magic angles. Such angles are determined by a balance between the negative energy contribution from interlayer currents and the positive contributions from the kinetic and confinement energies. A perturbative analysis performed around the first magic angle allows to explore analytically the details of such energy balance.

preprint2021arXiv

Electronic spectrum of Kekule patterned graphene considering second neighbor-interactions

The effects of second-neighbor interactions in Kekule patterned graphene electronic properties are studied starting from a tight-binding Hamiltonian. Thereafter, a low-energy effective Hamiltonian is obtained by projecting the high energy bands at the Gamma point into the subspace defined by the Kekule wave vector. The spectrum of the low energy Hamiltonian is in excellent agreement with the one obtained from a numerical diagonalization of the full tight-binding Hamiltonian. The main effect of the second-neighbour interaction is that a set of bands gains an effective mass and a shift in energy, thus lifting the degeneracy of the conduction bands at the Dirac point. This band structure is akin to a spin-one Dirac cone, a result expected for honeycomb lattices with a distinction between one third of the atoms in one sublattice. Finally, we present a study of Kekule patterned graphene nanoribbons. This shows that the previous effects are enhanced as the width decreases. Moreover, edge states become dispersive, as expected due to second neighbors interaction, but here the Kek-Y bond texture results in an hybridization of both edge states. The present study shows the importance of second neighbors in realistic models of Kekule patterned graphene, specially at surfaces.

preprint2020arXiv

Electron transitions for Dirac Hamiltonians with flat-bands under electromagnetic radiation and its application to the $α-\mathcal{T}_3$ graphene model

In a system with a Dirac-like linear dispersion there are always states that fulfill the resonance condition for electromagnetic radiation of arbitrary frequency $Ω$. When a flat band is present two kinds of resonant transitions are found. Considering the $α-\mathcal{T}_3$ graphene model as a minimal model with a flat band and Dirac cones, and describing the dynamics using the interaction picture, we study the band transitions induced by an external electromagnetic field. We found that transitions depend upon the relative angle between the electron momentum and the electromagnetic field wave vector. For parallel incidence, the transitions are found using Floquet theory while for other angles perturbation theory is used. In all cases, the transition probabilities and the frequencies are found. For some special values of the parameter $α$ or by charge doping, the system behaves as a three level or a two-level Rabi system. All these previous results were compared with numerical simulations. A good agreement was found between both. The obtained results are useful to provide a quantum control of the system.

preprint2020arXiv

Electronic and optical conductivity of kekulé-patterned graphene: Intravalley and intervalley transport

A Kubo formalism is used to calculate the electronic and optical conductivity of a graphene superlattice with Y-shaped kekulé bond texture, similar to that visualized in recent experiments. We show that new conduction channels between the valleys in graphene are opened by the kekulé distortion. This intervalley contribution to the electronic transport is not present in pristine graphene and here appears due to the folding of the Dirac cones $K$, $K'$ on top of each other. The contribution of intervalley transport to the conductivity of graphene, as well as the modification of the intravalley transport, are analyzed in detail for different frequency, temperature, chemical potential and scattering rate limits. We obtain analytical expressions for the conductivity that reproduce previous expressions used to fit experimental measurements in graphene and compare with direct numerical evaluations of the Kubo formula, finding great agreement. The optical absorption arising from intervalley transitions displays a maximum at a special frequency that can be tuned by doping. Our results show how the single parameter describing the valley coupling in this system could be obtained by measuring graphene's optical absorbance in the region where interband transitions are blocked by the impurities. Finally, we use Fermi's golden rule to independently verify some of the previous results.

preprint2020arXiv

Floquet spectrum and electronic transitions of tilted anisotropic Dirac materials under electromagnetic radiation: monodromy matrix approach

We analyze the quasienergy-spectrum and the valence to conduction-band transition probabilities of a tilted anisotropic Dirac material subject to linearly and circularly polarized electromagnetic fields. The quasienergy-spectrum is numerically calculated from the monodromy matrix of the Schrödinger equation via the Floquet theorem for arbitrarily intense electromagnetic fields. To asses the valence to conduction-band transition times we deduced a Rabi-like formula in the rotating wave approximation. In the strong-field regime the spectrum as a function of the momentum components divides into two very distinctive regions. In the first, located around the Dirac point, the quasi-spectrum is significantly distorted by the field as the electronic parameters are renormalized by electronic-dressing. In the second, all the characteristics of the free carrier spectrum are retained. Linearly polarized light anisotropically deforms the spectrum according to the field polarization direction. Dirac-like points form around the original Dirac point. The quasi spectrum of circularly polarized light, instead, exhibits a gap formation in the Dirac point and has elliptical symmetry. We show that, in contrast to the single-photon resonant transitions that characterize the weak-field regime, the strong-field regime is dominated by multiphoton resonances.

preprint2020arXiv

Floquet spectrum for anisotropic and tilted Dirac materials under linearly polarized light at all field intensities

The Floquet spectrum in an anisotropic tilted Dirac semimetal modulated by linearly polarized light is addressed through the solution of the time-dependent Schrödinger equation for the two-dimensional Dirac Hamiltonian via the Floquet theorem. The time-dependent wave functions and the quasienergy spectrum of the two-dimensional Dirac Hamiltonian under the normal incidence of linearly polarized waves are obtained for an arbitrarily intense electromagnetic radiation. We applied a set of unitary transformations to reduce the Schrödinger equation to an ordinary second-order differential Hill equation with complex coefficients. Through the stability analysis of this differential equation, the weak and strong field regimes are clearly distinguished in the quasi-spectrum. In the weak electric field regime, above a certain threshold given by the field parameters, the spectrum mostly resembles that of free electrons in graphene. Below this threshold, in the strong electric field regime, the spectrum abruptly becomes highly anisotropic and a gap opens up. As an example, we apply the results to the particular case of borophene.

preprint2019arXiv

How the overlap of excluded volume determines the configurational energy landscape and "thermodynamics" in the "one to five hard disks in a box" system

In this work, the effects of excluded volume are studied in the "one to five hard disks in a box" system. For one and two disks in different types of cages, the attractive and repulsive forces are calculated analytically. Attractive forces are due to excluded volume overlap, and thus allows to understand that in hard-core systems, second neighbors have an effective interaction, a fact that is usually neglected when considering only collisions with first-neighbors. The same effects are observed for numerical computations with five-disks in a box, as the distributions of distances between disks and disks to walls suggest. The results indicate that first and second-neighbor excluded volume interactions are essential to determine the energy landscape topology. This work supports previous observations that suggests that second-neighbor excluded volume interactions in hard-core systems are related with phase transitions.

preprint2015arXiv

Analytical calculation of electron's group velocity surfaces in uniform strained graphene

Electron group velocity for graphene under uniform strain is obtained analitically by using the Tight-Binding approx- imation. Such closed analytical expressions are useful in order to calculate electronic, thermal and optical properties of strained graphene. These results allow to understand the behavior of electrons when graphene is subjected to strong strain and nonlinear corrections, for which the usual Dirac approach is not longer valid. Some particular cases of uni- axial and shear strain were analized. The evolution of the electron group velocity indicates a break up of the trigonal warping symmetry, which is replaced by a warping consistent with the symmetry of the strained reciprocal lattice. The Fermi velocity becomes strongly anisotropic, i.e, for a strong pure shear-strain (20% of the lattice parameter), the two inequivalent Dirac cones merge and the Fermi velocity is zero in one of the principal axis of deformation. We found that non-linear terms are essential to describe the effects of deformation for electrons near or at the Fermi energy.

preprint2015arXiv

Comment on "Penrose Tilings as Jammed Solids"

In a recent letter, Stenull and Lubensky claim that periodic approximants of Penrose tilings, which are generically isostatic, have a nonzero bulk modulus B when disordered, and, therefore, Penrose tilings are good models of jammed packings. The claim of a nonzero B, which is made on the basis of a normal mode analysis of periodic Penrose approximants for a single value of the disorder epsilon, is the central point of their letter: other properties of Penrose tilings, such as the vanishing of the shear modulus, and a flat density of vibrational states, are already shared by most geometrically disordered isostatic networks studied so far. In this comment, Conjugate Gradient is used to solve the elastic equations on approximants with up to 8x10^4 sites for several values of epsilon, to show beyond reasonable doubt that Stenull and Lubensky's claim is incorrect. The bulk modulus of generic Penrose tilings is zero asymptotically. According to our results, B grows as (epsilon^2 L^3) when (epsilon^2 L^3) << 10^2, then saturates, and finally decays as (epsilon^2 L^3)^{-2/3} ~ 1/L^2 for epsilon^2 L^3 >> 10^2. Stenull and Lubensky seem to have only analyzed one value of epsilon for which saturation is reached at the largest size studied. This led them to a wrong conclusion. We support our results by also considering generic Penrose approximants with fixed boundaries, whose bulk modulus constitutes a strict upper bound for that of periodic systems, finding that these have a vanishing B as well for large L. We conclude that the main point in Stenull and Lubensky letter is unjustified. Penrose tilings are no better models of jammed packings than any of the previously studied isostatic networks with geometric disorder.

preprint2015arXiv

Generalizing the Fermi velocity of strained graphene from uniform to nonuniform strain

The relevance of the strain-induced Dirac point shift to obtain the appropriate anisotropic Fermi velocity of strained graphene is demonstrated. Then a critical revision of the available effective Dirac Hamiltonians is made by studying in detail the limiting case of a uniform strain. An effective Dirac Hamiltonian for nonuniform strain is thus reported, which takes into account all strain-induced effects: changes in the nearest-neighbor hopping parameters, the reciprocal lattice deformation and the true shift of the Dirac point. Pseudomagnetic fields are thus explained by means of position-dependent Dirac cones, whereas complex gauge fields appear as a consequence of a position-dependent Fermi velocity. Also, position-dependent Fermi velocity effects on the spinor wavefunction are considered for interesting cases of deformations such as flexural modes.

preprint2015arXiv

Heider balance, asymmetric ties, and gender segregation

To remove a cognitive dissonance in interpersonal relations, people tend to divide our acquaintances into friendly and hostile parts, both groups internally friendly and mutually hostile. This process is modeled as an evolution towards the Heider balance. A set of differential equations have been proposed and validated (Kulakowski {\it et al}, IJMPC 16 (2005) 707) to model the Heider dynamics of this social and psychological process. Here we generalize the model by including the initial asymmetry of the interprersonal relations and the direct reciprocity effect which removes this asymmetry. Our model is applied to the data on enmity and friendship in 37 school classes and 4 groups of teachers in México. For each class, a stable balanced partition is obtained into two groups. The gender structure of the groups reveals stronger gender segregation in younger classes, i.e. of age below 12 years, a fact consistent with other general empirical results.

preprint2015arXiv

Minimal cooling speed for glass transition in a simple solvable energy landscape model

The minimal cooling speed required to form a glass is obtained for a simple solvable energy landscape model. The model, made from a two-level system modified to include the topology of the energy landscape, is able to capture either a glass transition or a crystallization depending on cooling rate. In this setup, the minimal cooling speed to achieve glass formation is then found to be related with the relaxation time and with the thermal history. In particular, we obtain that the thermal history encodes small fluctuations around the equilibrium population which are exponentially amplified near the glass transition, which mathematically corresponds to the boundary layer of the master equation. Finally, to verify our analytical results, a kinetic Monte-Carlo simulation was implemented.

preprint2014arXiv

Anisotropic AC conductivity of strained graphene

The density of states and the AC conductivity of graphene under uniform strain are calculated using a new Dirac Hamiltonian that takes into account the main three ingredients that change the electronic properties of strained graphene: the real displacement of the Fermi energy, the reciprocal lattice strain and the changes in the overlap of atomic orbitals. Our simple analytical expressions of the density of states and the AC conductivity generalizes previous expressions only available for uniaxial strain. The results suggest a way to measure the Gruneisen parameter that appears in any calculation of strained graphene, as well as the emergence of a sort of Hall effect due to shear strain.

preprint2014arXiv

Mapping of strained graphene into one-dimensional Hamiltonians: quasicrystals and modulated crystals

The electronic properties of graphene under any arbitrary uniaxial strain field are obtained by an exact mapping of the corresponding tight-binding Hamiltonian into an effective one-dimensional modulated chain. For a periodic modulation, the system displays a rich behavior, including quasicrystals and modulated crystals with a complex spectrum, gaps at the Fermi energy and interesting localization properties. These features are explained by the incommensurate or commensurate nature of the potential, which leads to a dense filling of the reciprocal space in the former case. Thus, the usual perturbation theory approach breaks down in some cases, as is proved by analyzing a special momentum that uncouples the model into dimers.

preprint2014arXiv

Spectral butterfly, mixed Dirac-Schrödinger fermion behavior and topological states in armchair uniaxial strained graphene

An exact mapping of the tight-binding Hamiltonian for a graphene's nanoribbon under any armchair uniaxial strain into an effective one-dimensional system is presented. As an application, for a periodic modulation we have found a gap opening at the Fermi level and a complex fractal spectrum, akin to the Hofstadter butterfly resulting from the Harper model. The latter can be explained by the commensurability or incommensurability nature of the resulting effective potential. When compared with the zig-zag uniaxial periodic strain, the spectrum shows much bigger gaps, although in general the states have a more extended nature. For a special critical value of the strain amplitude and wavelength, a gap is open. At this critical point, the electrons behave as relativistic Dirac femions in one direction, while in the other, a non-relativistic Schrödinger behavior is observed. Also, some topological states were observed which have the particularity of not being completly edge states since they present some amplitude in the bulk. However, these are edge states of the effective system due to a reduced dimensionality through decoupling. These states also present the fractal Chern beating observed recently in quasiperiodic systems.

preprint2013arXiv

Chern and Majorana Modes of Quasiperiodic Systems

New types of self-similar states are found in quasiperiodic systems characterized by topological invariants-- the Chern numbers. We show that the topology introduces a competing length in the self-similar band edge states transforming peaks into doublets of size equal to the Chern number. This length intertwines with the quasiperiodicity and introduces an intrinsic scale, producing Chern-beats and nested regions where the fractal structure becomes smooth. Cherns also influence the zero-energy mode, that for quasiperiodic systems which exhibit exponential localization, is related to the ghost of the Majorana; the delocalized state at the onset to topological transition. The Chern and the Majorana, two distinct types of topological edge modes, exist in quasiperiodic superconducting wires.

preprint2013arXiv

Electron localization in disordered graphene: multifractal properties of the wavefunctions

An analysis of the electron localization properties in doped graphene is performed by doing a numerical multifractal analysis. By obtaining the singularity spectrum of a tight-binding model, it is found that the electron wave functions present a multifractal behavior. Such multifractality is preserved even for second neighbor interaction, which needs to be taken into account if a comparison is desired with experimental results. States close to the Dirac point have a wider multifractal character than those far from this point as the impurity concentration is increased. The analysis of the results allows to conclude that in the split-band limit, where impurities act as vacancies, the system can be well described by a chiral orthogonal symmetry class, with a singularity spectrum transition approaching freezing as disorder increases. This also suggests that in doped graphene, localization is in contrast with the conventional picture of Anderson localization in two dimensions, showing also that the common belief of the absence of quantum percolation in two dimensional systems needs to be revised.

preprint2013arXiv

Uniform pseudomagnetic fields in graphene: corrections due to lattice deformations

Very recently, the standard description of electrons in strained graphene has been completed by the explicit inclusion of the lattice deformation. Here, the effect of these lattice corrections is taken into account to find the mechanical deformations that generate a uniform pseudomagnetic field inside strained graphene. The main difference when compared with previous proposals in a rectangular geometry, is the need to stretch both zigzag and armchair edges.

preprint2012arXiv

A simple solvable energy landscape model that shows a thermodynamic phase transition and a glass transition

When a liquid melt is cooled, a glass or phase transition can be obtained depending on the cooling rate. Yet, this behavior has not been clearly captured in energy landscape models. Here a model is provided in which two key ingredients are considered based in the landscape, metastable states and their multiplicity. Metastable states are considered as in two level system models. However, their multiplicity and topology allows a phase transition in the thermodynamic limit, while a transition to the glass is obtained for fast cooling. By solving the corresponding master equation, the minimal speed of cooling required to produce the glass is obtained as a function of the distribution of metastable and stable states. This allows to understand cooling trends due to rigidity considerations in chalcogenide glasses.

preprint2012arXiv

Critical wave functions in disordered graphene

In order to elucidate the presence of non-localized states in doped graphene, an scaling analysis of the wave function moments known as inverse participation ratios is performed. The model used is a tight- binding hamiltonian considering nearest and next-nearest neighbors with random substitutional impurities. Our findings indicate the presence of non-normalizable wave functions that follow a critical (power-law) decay, which are between a metallic and insulating behavior. The power-law exponent distribution is robust against the inclusion of next-nearest neighbors and on growing the system size.

preprint2007arXiv

Tail universalities in rank distributions as an algebraic problem: the beta-like function

Although power laws of the Zipf type have been used by many workers to fit rank distributions in different fields like in economy, geophysics, genetics, soft-matter, networks etc., these fits usually fail at the tails. Some distributions have been proposed to solve the problem, but unfortunately they do not fit at the same time both ending tails. We show that many different data in rank laws, like in granular materials, codons, author impact in scientific journal, etc. are very well fitted by a beta-like function. Then we propose that such universality is due to the fact that a system made from many subsystems or choices, imply stretched exponential frequency-rank functions which qualitatively and quantitatively can be fitted with the proposed beta-like function distribution in the limit of many random variables. We prove this by transforming the problem into an algebraic one: finding the rank of successive products of a given set of numbers.

preprint2006arXiv

From the Hofstadter to the Fibonacci butterfly

We show that the electronic spectrum of a tight-binding Hamiltonian defined in a quasiperiodic chain with an on-site potential given by a Fibonacci sequence, can be obtained as a superposition of Harper potentials. The electronic spectrum of the Harper equation is a fractal set, known as Hofstadter butterfly. Here we show that is possible to construct a similar butterfly for the Fibonacci potential just by adding harmonics to the Harper potential. As a result, the equations in reciprocal space for the Fibonacci case have the form of a chain with a long range interaction between Fourier components. Then we explore the transformation between both spectra, and specifically the origin of energy gaps due to the analytical calculation of the components in reciprocal space of the potentials. We also calculate some localization properties by finding the correlator of each potential.

preprint2006arXiv

Minimal Multifractality in the Spectrum of a Quasiperiodic Hamiltonian

In many systems, the electronic energy spectrum is a continuous or singular continuous multifractal set with a distribution of scaling exponents. Here, we show that for a quasiperiodic potential, the multifractal energy spectrum can have a minimal dispersion of the scaling exponents. This is made by tuning the ratio between self-energies in a tigth-binding Hamiltonian defined in a quasiperiodic Fibonacci chain. The tuning is calculated numerically from a trace map, and coincides with the place where the scaling exponents of the map, obtained from cycles with period six and two, are equal. The diffusion of electronic wave packets also reflects the minimal multifractal fractal nature of the spectrum, by reducing intermittency. The present result can help to simplify the task of studying several problems in quasiperiodic systems, since the effects of multiscaling are better isolated.

preprint2005arXiv

Monte-Carlo rejection as a tool for measuring the energy landscape scaling of simple fluids

A simple modification of the Monte-Carlo algortihm is proposed to explore the topography and the scaling of the energy landscape. We apply this idea to a simple hard-core fluid. The results for different packing fractions show a power law scaling of the landscape boundary, with a characteristic scale that separates the values of the scaling exponents. Finally, it is shown how the topology determines the freezing point of the system due to the increasing importance and complexity of the boundary.

preprint2005arXiv

Variation of the glass transition temperature with rigidity and chemical composition

The effects of flexibility and chemical composition in the variation of the glass transition temperature are obtained by using the Lindemann criteria, that relates melting temperature with atomic vibrations. Using this criteria and that floppy modes at low frequencies enhance in a considerable way the average cuadratic displacement, we show that the consequence is a modified glass transition temperature. This approach allows to obtain in a simple way the empirically modified Gibbs-DiMarzio law, which has been widely used in chalcogenide glasses to fit the changes in the glass transition temperature with the chemical composition . The method predicts that the constant that appears in the law depends upon the ratio of two characteristic frequencies (or temperatures). Then, the constant for the Se-Ge-As glass is estimated by using the experimental density of vibrational states, and the result shows a very good agreement with the experimental fit from glass transition temperature variation.

preprint2004arXiv

Energy landscape and rigidity

The effects of floppy modes in the thermodynamical properties of a system are studied. From thermodynamical arguments, we deduce that floppy modes are not at zero frequency and thus a modified Debye model is used to take into account this effect. The model predicts a deviation from the Debye law at low temperatures. Then, the connection between the topography of the energy landscape, the topology of the phase space and the rigidity of a glass is explored. As a result, we relate the number of constraints and floppy modes with the statistics of the landscape. We apply these ideas to a simple model for which we provide an approximate expression for the number of energy basins as a function of the rigidity. This allows to understand certains features of the glass transition, like the jump in the specific heat or the reversible window observed in chalcogenide glasses.