Refine
Year of publication
Document Type
- Preprint (37)
- Article (35)
- Contribution to a Periodical (1)
Has Fulltext
- yes (73)
Is part of the Bibliography
- no (73)
Keywords
- density functional theory (4)
- α-RuCl3 (3)
- electronic band structure (2)
- electronic transport (2)
- p-n junction (2)
- photoelectron spectroscopy (2)
- scanning tunneling microscopy (2)
- transition-metal oxides (2)
- (F)EBID (1)
- 2D vdW magnets (1)
Institute
The rapid spread of the Coronavirus (COVID-19) confronts policy makers with the problem of measuring the effectiveness of containment strategies, balancing public health considerations with the economic costs of social distancing measures. We introduce a modified epidemic model that we name the controlled-SIR model, in which the disease reproduction rate evolves dynamically in response to political and societal reactions. An analytic solution is presented. The model reproduces official COVID-19 cases counts of a large number of regions and countries that surpassed the first peak of the outbreak. A single unbiased feedback parameter is extracted from field data and used to formulate an index that measures the efficiency of containment strategies (the CEI index). CEI values for a range of countries are given. For two variants of the controlled-SIR model, detailed estimates of the total medical and socio-economic costs are evaluated over the entire course of the epidemic. Costs comprise medical care cost, the economic cost of social distancing, as well as the economic value of lives saved. Under plausible parameters, strict measures fare better than a hands-off policy. Strategies based on current case numbers lead to substantially higher total costs than strategies based on the overall history of the epidemic.
Predicting the cumulative medical load of COVID-19 outbreaks after the peak in daily fatalities
(2021)
The distinct ways the COVID-19 pandemic has been unfolding in different countries and regions suggest that local societal and governmental structures play an important role not only for the baseline infection rate, but also for short and long-term reactions to the outbreak. We propose to investigate the question of how societies as a whole, and governments in particular, modulate the dynamics of a novel epidemic using a generalization of the SIR model, the reactive SIR (short-term and long-term reaction) model. We posit that containment measures are equivalent to a feedback between the status of the outbreak and the reproduction factor. Short-term reaction to an outbreak corresponds in this framework to the reaction of governments and individuals to daily cases and fatalities. The reaction to the cumulative number of cases or deaths, and not to daily numbers, is captured in contrast by long-term reaction. We present the exact phase space solution of the controlled SIR model and use it to quantify containment policies for a large number of countries in terms of short and long-term control parameters. We find increased contributions of long-term control for countries and regions in which the outbreak was suppressed substantially together with a strong correlation between the strength of societal and governmental policies and the time needed to contain COVID-19 outbreaks. Furthermore, for numerous countries and regions we identified a predictive relation between the number of fatalities within a fixed period before and after the peak of daily fatality counts, which allows to gauge the cumulative medical load of COVID-19 outbreaks that should be expected after the peak. These results suggest that the proposed model is applicable not only for understanding the outbreak dynamics, but also for predicting future cases and fatalities once the effectiveness of outbreak suppression policies is established with sufficient certainty. Finally, we provide a web app (https://itp.uni-frankfurt.de/covid-19/) with tools for visualising the phase space representation of real-world COVID-19 data and for exporting the preprocessed data for further analysis.
Recent experimental findings have reported the presence of unconventional charge orders in the enlarged (2 × 2) unit-cell of kagome metals AV3Sb5 (A = K, Rb, Cs) and hinted towards specific topological signatures. Motivated by these discoveries, we investigate the types of topological phases that can be realized in such kagome superlattices. In this context, we employ a recently introduced statistical method capable of constructing topological models for any generic lattice. By analyzing large data sets generated from symmetry-guided distributions of randomized tight-binding parameters, and labeled with the corresponding topological index, we extract physically meaningful information. We illustrate the possible real-space manifestations of charge and bond modulations and associated flux patterns for different topological classes, and discuss their relation to present theoretical predictions and experimental signatures for the AV3Sb5 family. Simultaneously, we predict higher-order topological phases that may be realized by appropriately manipulating the currently known systems.
Electron beam-induced deposition with tungsten hexacarbonyl W(CO)6 as precursors leads to granular deposits with varying compositions of tungsten, carbon and oxygen. Depending on the deposition conditions, the deposits are insulating or metallic. We employ an evolutionary algorithm to predict the crystal structures starting from a series of chemical compositions that were determined experimentally. We show that this method leads to better structures than structural relaxation based on estimated initial structures. We approximate the expected amorphous structures by reasonably large unit cells that can accommodate local structural environments that resemble the true amorphous structure. Our predicted structures show an insulator-to-metal transition close to the experimental composition at which this transition is actually observed and they also allow comparison with experimental electron diffraction patterns.
The description of quantized collective excitations stands as a landmark in the quantum theory of condensed matter. A prominent example occurs in conventional magnets, which support bosonic magnons—quantized harmonic fluctuations of the ordered spins. In striking contrast is the recent discovery that strongly spin-orbital-coupled magnets, such as α-RuCl3, may display a broad excitation continuum inconsistent with conventional magnons. Due to incomplete knowledge of the underlying interactions unraveling the nature of this continuum remains challenging. The most discussed explanation refers to a coherent continuum of fractional excitations analogous to the celebrated Kitaev spin liquid. Here, we present a more general scenario. We propose that the observed continuum represents incoherent excitations originating from strong magnetic anharmonicity that naturally occurs in such materials. This scenario fully explains the observed inelastic magnetic response of α-RuCl3 and reveals the presence of nontrivial excitations in such materials extending well beyond the Kitaev state.
Motivated by recent experimental suggestions of charge-order-driven ferroelectricity in organic charge-transfer salts, such as κ-(BEDT-TTF)2Cu[N(CN)2]Cl, we investigate magnetic and charge-ordered phases that emerge in an extended two-orbital Hubbard model on the anisotropic triangular lattice at 3/4 filling. This model takes into account the presence of two organic BEDT-TTF molecules, which form a dimer on each site of the lattice, and includes short-range intramolecular and intermolecular interactions and hoppings. By using variational wave functions and quantum Monte Carlo techniques, we find two polar states with charge disproportionation inside the dimer, hinting to ferroelectricity. These charge-ordered insulating phases are stabilized in the strongly correlated limit and their actual charge pattern is determined by the relative strength of intradimer to interdimer couplings. Our results suggest that ferroelectricity is not driven by magnetism, since these polar phases can be stabilized also without antiferromagnetic order and provide a possible microscopic explanation of the experimental observations. In addition, a conventional dimer-Mott state (with uniform density and antiferromagnetic order) and a nonpolar charge-ordered state (with charge-rich and charge-poor dimers forming a checkerboard pattern) can be stabilized in the strong-coupling regime. Finally, when electron–electron interactions are weak, metallic states appear, with either uniform charge distribution or a peculiar 12-site periodicity that generates honeycomb-like charge order.
One of the most challenging problems in solid state systems is the microscopic analysis of electronic correlations. A paramount minimal model that encodes correlation effects is the Hubbard Hamiltonian, which—regardless of its simplicity—is exactly solvable only in a few limiting cases and approximate many-body methods are required for its solution. In this review, an overview on the non-perturbative two-particle self-consistent method (TPSC), which was originally introduced to describe the electronic properties of the single-band Hubbard model, is presented. A detailed derivation of the multi-orbital generalization of TPSC is introduced here and particular features of the method on exemplary interacting models in comparison to dynamical mean-field theory results are discussed.
Electronic states with non-trivial topology host a number of novel phenomena with potential for revolutionizing information technology. The quantum anomalous Hall effect provides spin-polarized dissipation-free transport of electrons, while the quantum spin Hall effect in combination with superconductivity has been proposed as the basis for realizing decoherence-free quantum computing. We introduce a new strategy for realizing these effects, namely by hole and electron doping kagome lattice Mott insulators through, for instance, chemical substitution. As an example, we apply this new approach to the natural mineral herbertsmithite. We prove the feasibility of the proposed modifications by performing ab-initio density functional theory calculations and demonstrate the occurrence of the predicted effects using realistic models. Our results herald a new family of quantum anomalous Hall and quantum spin Hall insulators at affordable energy/temperature scales based on kagome lattices of transition metal ions.
Type-II multiferroic materials, in which ferroelectric polarization is induced by inversion non-symmetric magnetic order, promise new and highly efficient multifunctional applications based on mutual control of magnetic and electric properties. However, to date this phenomenon is limited to low temperatures. Here we report giant pressure-dependence of the multiferroic critical temperature in CuBr2: at 4.5 GPa it is enhanced from 73.5 to 162 K, to our knowledge the highest TC ever reported for non-oxide type-II multiferroics. This growth shows no sign of saturating and the dielectric loss remains small under these high pressures. We establish the structure under pressure and demonstrate a 60\% increase in the two-magnon Raman energy scale up to 3.6 GPa. First-principles structural and magnetic energy calculations provide a quantitative explanation in terms of dramatically pressure-enhanced interactions between CuBr2 chains. These large, pressure-tuned magnetic interactions motivate structural control in cuprous halides as a route to applied high-temperature multiferroicity.
Formation of Hubbard-like bands as a fingerprint of strong electron-electron interactions in FeSe
(2017)
We use angle-resolved photo-emission spectroscopy (ARPES) to explore the electronic structure of single crystals of FeSe over a wide range of binding energies and study the effects of strong electron-electron correlations. We provide evidence for the existence of "Hubbard-like bands" at high binding energies consisting of incoherent many-body excitations originating from Fe 3d states in addition to the renormalized quasiparticle bands near the Fermi level. Many high energy features of the observed ARPES data can be accounted for when incorporating effects of strong local Coulomb interactions in calculations of the spectral function via dynamical mean-field theory, including the formation of a Hubbard-like band. This shows that over the energy scale of several eV, local correlations arising from the on-site Coulomb repulsion and Hund's coupling are essential for a proper understanding of the electronic structure of FeSe and other related iron based superconductors.
Controlling and understanding electron correlations in quantum matter is one of the most challenging tasks in materials engineering. In the past years a plethora of new puzzling correlated states have been found by carefully stacking and twisting two-dimensional van der Waals materials of different kind. Unique to these stacked structures is the emergence of correlated phases not foreseeable from the single layers alone. In Ta-dichalcogenide heterostructures made of a good metallic “1H”- and a Mott insulating “1T”-layer, recent reports have evidenced a cross-breed itinerant and localized nature of the electronic excitations, similar to what is typically found in heavy fermion systems. Here, we put forward a new interpretation based on first-principles calculations which indicates a sizeable charge transfer of electrons (0.4-0.6 e) from 1T to 1H layers at an elevated interlayer distance. We accurately quantify the strength of the interlayer hybridization which allows us to unambiguously determine that the system is much closer to a doped Mott insulator than to a heavy fermion scenario. Ta-based heterolayers provide therefore a new ground for quantum-materials engineering in the regime of heavily doped Mott insulators hybridized with metallic states at a van der Waals distance.
Critical spin liquid versus valence-bond glass in a triangular-lattice organic antiferromagnet
(2019)
In the quest for materials with unconventional quantum phases, the organic triangular-lattice antiferromagnet κ-(ET)2Cu2(CN)3 has been extensively discussed as a quantum spin liquid (QSL) candidate. The description of its low temperature properties has become, however, a particularly challenging task. Recently, an intriguing quantum critical behaviour was suggested from low-temperature magnetic torque experiments. Here we highlight significant deviations of the experimental observations from a quantum critical scenario by performing a microscopic analysis of all anisotropic contributions, including Dzyaloshinskii–Moriya and multi-spin scalar chiral interactions. Instead, we show that disorder-induced spin defects provide a comprehensive explanation of the low-temperature properties. These spins are attributed to valence bond defects that emerge spontaneously as the QSL enters a valence-bond glass phase at low temperature. This theoretical treatment is applicable to a general class of frustrated magnetic systems and has important implications for the interpretation of magnetic torque, nuclear magnetic resonance, thermal transport and thermodynamic experiments.
The phase diagram of the square lattice bilayer Hubbard model: a variational Monte Carlo study
(2014)
We investigate the phase diagram of the square lattice bilayer Hubbard model at half-filling with the variational Monte Carlo method for both the magnetic and the paramagnetic case as a function of the interlayer hopping and on-site Coulomb repulsion U. With this study we resolve some discrepancies in previous calculations based on the dynamical mean-field theory, and we are able to determine the nature of the phase transitions between metal, Mott insulator and band insulator. In the magnetic case we find only two phases: an antiferromagnetic Mott insulator at small for any value of U and a band insulator at large . At large U values we approach the Heisenberg limit. The paramagnetic phase diagram shows at small a metal to Mott insulator transition at moderate U values and a Mott to band insulator transition at larger U values. We also observe a re-entrant Mott insulator to metal transition and metal to band insulator transition for increasing in the range of . Finally, we discuss the phase diagrams obtained in relation to findings from previous studies based on different many-body approaches.
We use the topological heavy fermion (THF) model and its Kondo Lattice (KL) formulation to study the symmetric Kondo state in twisted bilayer graphene. Via a large-N approximation, we find a symmetric Kondo (SK) state in KL mode at fillings ν=0,±1,±2. In the SK state, all symmetries are preserved and the local moments are Kondo screened by the conduction electrons. At the mean-field level of the THF model at ν=0,±1,±2,±3, we also find a similar symmetric state. We study the stability of the symmetric state by comparing its energy with the ordered states and find the ordered states to have lower energy. However, moving away from integer fillings by doping holes to the light bands, we find the energy difference is reduced, which suggests the loss of ordering and a tendency towards Kondo screening. In order to include many-body effects beyond the mean-field approximation, we perform dynamical mean-field theory (DMFT) calculations on the THF model. We find the spin susceptibility follows a Curie behavior at ν=0,±1,±2 down to ∼2K where the onset of screening of the local moment becomes visible. This hints to very low Kondo temperatures at these fillings, in agreement with the outcome of our mean-field calculations. At non-integer filling ν=±0.5,±0.8,±1.2 DMFT shows deviations from a 1/T-susceptibility at much higher temperatures, suggesting a more effective screening of local moments with doping. Finally, we study the effect of a C3z-rotational-symmetry-breaking strain via mean-field approaches and find that a symmetric phase (that only breaks C3z symmetry) can be stabilized at sufficiently large strain at ν=0,±1,±2. Our results suggest that a symmetric Kondo phase is strongly suppressed at integer fillings, but could be stabilized either at non-integer fillings or by applying strain.
Under temperature or pressure tuning, tetragonal EuPd2Si2 is known to undergo a valence transition from nearly divalent to nearly trivalent Eu accompanied by a volume reduction. Albeit intensive work, its microscopic origin is still being discussed. Here, we investigate the mechanism of the valence transition under volume compression by ab initio density functional theory (DFT) calculations. Our analysis of the electronic and magnetic properties of EuPd2Si2 when approaching the valence transition shows an enhanced c-f hybridization between localized Eu 4f states and itinerant conduction states (Eu 5d, Pd 4d, and Si 3p) where an electronic charge redistribution takes place. We observe that the change in the electronic structure is intimately related to the volume reduction where Eu-Pd(Si) bond lengths shorten and, for the transition to happen, we trace the delicate balance between electronic bandwidth, crystal field splitting, Coulomb repulsion, Hund's coupling and spin-orbit coupling. In a next step we compare and benchmark our DFT results to surface-sensitive photoemission data in which the mixed-valent properties of EuPd2Si2 are reflected in a simultaneous observation of divalent and trivalent signals from the Eu 4f shell. The study serves as well to explore the limits of density functional theory and the choice of exchange correlation functionals to describe such a phenomenon as a valence transition.
The maximum recoverable strain of most crystalline solids is less than 1% because plastic deformation or fracture usually occurs at a small strain. In this work, we show that a SrNi2P2 micropillar exhibits pseudoelasticity with a large maximum recoverable strain of ~14% under uniaxial compression via unique reversible structural transformation, double lattice collapse-expansion that is repeatable under cyclic loading. Its high yield strength (~3.8±0.5 GPa) and large maximum recoverable strain bring out the ultrahigh modulus of resilience (~146±19MJ/m3) a few orders of magnitude higher than that of most engineering materials. The double lattice collapse-expansion mechanism shows stress-strain behaviors similar with that of conventional shape memory alloys, such as hysteresis and thermo-mechanical actuation, even though the structural changes involved are completely different. Our work suggests that the discovery of a new class of high performance ThCr2Si2-structured materials will open new research opportunities in the field of pseudoelasticity
We demonstrate ultra-sharp (≲10 nm) lateral p-n junctions in graphene using electronic transport, scanning tunneling microscopy, and first principles calculations. The p-n junction lies at the boundary between differentially-doped regions of a graphene sheet, where one side is intrinsic and the other is charge-doped by proximity to a flake of α-RuCl3 across a thin insulating barrier. We extract the p-n junction contribution to the device resistance to place bounds on the junction width. We achieve an ultra-sharp junction when the boundary between the intrinsic and doped regions is defined by a cleaved crystalline edge of α-RuCl3 located 2 nm from the graphene. Scanning tunneling spectroscopy in heterostructures of graphene, hexagonal boron nitride, and α-RuCl3 shows potential variations on a sub-10 nm length scale. First principles calculations reveal the charge-doping of graphene decays sharply over just nanometers from the edge of the α-RuCl3 flake.
In the novel stoichiometric iron-based material RbEuFe4As4 superconductivity coexists with a peculiar long-range magnetic order of Eu 4f states; their coexistance is puzzling and represents a challenge for both experiment and theory. Using angle-resolved photoemission spectroscopy, resonant photoemission spectroscopy, Andreev reflection spectroscopy and scanning tunneling spectroscopy we have addressed this puzzle and unambigously shown that Fe- and Eu-derived states are largely decoupled and that superconducting and a long range magnetic orders exist almost independently from each other.
We present results of hard X-ray angle-resolved photoemission spectroscopy and photoemission diffraction measurements performed on high-quality single crystals of the valence transition compound EuPd2Si2 for temperatures 25~K ≤ T ≤ 300~K. At low temperatures we observe a Eu 4f valence v=2.5, % occupation number n=6.5, which decreases to v=2.1 for temperatures above the valence transition around TV≈160~K. The experimental valence numbers resulting from an evaluation of the Eu(III)/Eu(II) 3d core levels, are used for calculating band structures using density functional theory. The valence transition significantly changes the band structure as determined by angle-resolved photoemission spectroscopy. In particular, the Eu 5d valence bands are shifted to lower binding energies with increasing Eu 4f occupancy. To a lesser extent, bands derived from the Si 3p and Pd 4d orbitals are also affected. This observation suggests a partial charge transfer between Eu and Pd/Si sites. Comparison with {\it ab-initio} theory shows a good agreement with experiment, in particular concerning the unequal band shift with increasing Eu 4f occupancy.
The discovery of the 1144-phase, e.g. CaKFe4As4, creates opportunities to build novel intermetallics with alternative stacking of two parent compounds. Here we formalize the idea by defining a class of bulk crystalline solids with A-B stacking (including 1144-phases and beyond), which is a generalization of hetero-structures from few-layer or thin-film semi-conductors to bulk intermetallics. Theoretically, four families of phosphides \textit{AB}(TM)4P4 (TM=Fe, Ru, Co, Ni) are investigated by first-principles calculations, wherein configurational, vibrational and electronic degrees of freedom are considered. It predicts a variety of stable 1144-phases (especially Ru- and Fe-phosphides). Stability rules are found and structural/electronic properties are discussed. Experimentally, we synthesize high-purity CaKRu4P4 as a proof of principle example. The synthetic method is simple and easily applied. Moreover, it alludes to a strategy to explore complex multi-component compounds, facilitated by a phase diagram coordinated by collective descriptors.
Rich functionalities of transition-metal oxides and their interfaces bear an enormous technological potential. Its realization in practical devices requires, however, a significant improvement of yet relatively low electron mobility in oxide materials. Recently, a mobility boost of about two orders of magnitude has been demonstrated at the spinel/perovskite {\gamma}-Al2O3/SrTiO3 interface compared to the paradigm perovskite/perovskite LaAlO3/SrTiO3. We explore the fundamental physics behind this phenomenon from direct measurements of the momentum-resolved electronic structure of this interface using resonant soft-X-ray angle-resolved photoemission. We find an anomaly in orbital ordering of the mobile electrons in {\gamma}-Al2O3/SrTiO3 which depopulates electron states in the top STO layer. This rearrangement of the mobile electron system pushes the electron density away from the interface that reduces its overlap with the interfacial defects and weakens the electron-phonon interaction, both effects contributing to the mobility boost. A crystal-field analysis shows that the band order alters owing to the symmetry breaking between the spinel {\gamma}-Al2O3 and perovskite SrTiO3. The band-order engineering exploiting the fundamental symmetry properties emerges as another route to boost the performance of oxide devices.
In magic angle twisted bilayer graphene, transport, thermodynamic and spectroscopic experiments pinpoint at a competition between distinct low-energy states with and without electronic order, as well as a competition between localized and delocalized charge carriers. In this study, we utilize Dynamical Mean Field Theory (DMFT) on the topological heavy Fermion (THF) model of twisted bilayer graphene to investigate the emergence of electronic correlations and long-range order in the absence of strain. We explain the nature of emergent insulating and correlated metallic states, as well as transitions between them driven by three central phenomena: (i) the formation of local spin and valley isospin moments around 100K, (ii) the ordering of the local isospin moments around 10K, and (iii) a cascadic redistribution of charge between localized and delocalized electronic states upon doping. At integer fillings, we find that low energy spectral weight is depleted in the symmetric phase, while we find insulating states with gaps enhanced by exchange coupling in the zero-strain ordered phases. Doping away from integer filling results in distinct metallic states: a "bad metal" above the ordering temperature, where coherence of the low-energy electronic excitations is suppressed by scattering off the disordered local moments, and a "good metal" in the ordered states with coherence of quasiparticles facilitated by isospin order. Upon doping, there is charge transfer between the localized and delocalized orbitals of the THF model such that they get periodically filled and emptied in between integer fillings. This charge reshuffling manifests itself in cascades of doping-induced Lifshitz transitions, local spectral weight redistributions and periodic variations of the electronic compressibility ranging from nearly incompressible to negative.
Strontium ruthenate Sr2RuO4 is an unconventional superconductor whose pairing symmetry has not been fully clarified, despite more than two decades of intensive research. Recent NMR Knight shift experiments have rekindled the Sr2RuO4 pairing debate by giving strong evidence against all odd-parity pairing states, including chiral p-wave pairing that was for a long time the leading pairing candidate. Here, we exclude additional pairing states by analyzing recent elastocaloric measurements [YS. Li et al., Nature 607, 276--280 (2022)]. To be able to explain the elastocaloric experiment, we find that unconventional even-parity pairings must include either large dx2−y2-wave or large {dxz∣dyz}-wave admixtures, where the latter possibility arises because of the body-centered point group symmetry. These {dxz∣dyz}-wave admixtures take the form of distinctively body-centered-periodic harmonics that have horizontal line nodes. Hence gxy(x2−y2)-wave and dxy-wave pairings are excluded as possible dominant even pairing states.
The existence of bound states induced by local impurities coupled to an insulating host depends decisively on the global topological properties of the host's electronic structure. In this context, we consider magnetic impurities modelled as classical unit-length spins that are exchange-coupled to the spinful Haldane model on the honeycomb lattice. We investigate the spectral flow of bound states with the coupling strength J in both the topologically trivial and Chern-insulating phases. In addition to conventional k-space topology, an additional, spatially local topological feature is available, based on the space of impurity-spin configurations forming, in case of R impurities, an R-fold direct product of two-dimensional spheres. Global k-space and local S-space topology are represented by different topological invariants, the first (k-space) Chern number and the R-th (S-space) spin-Chern number. We demonstrate that there is a local S-space topological transition as a function of J associated with a change in the spin Chern number and work out the implications of this for the J-dependent local electronic structure close to the impurities and, in particular, for in-gap bound states. The critical exchange couplings' dependence on the parameters of the Haldane model, and thus on the k-space topological state, is obtained numerically to construct local topological phase diagrams for systems with R=1 and R=2 impurity spins.
Although iron-based catalysts are regarded as a promising alternative to precious metal catalysts, their precise electronic structures during catalysis still pose challenges for computational descriptions. A particularly urgent question is the influence of the environment on the electronic structure, and how to describe this properly with computational methods. Here, we study an iron porphyrin chloride complex adsorbed on a graphene sheet using density functional theory calculations to detail how much the electronic structure is influenced by the presence of a graphene layer. Our results indicate that weak interactions due to van der Waals forces dominate between the porphyrin complex and graphene, and only a small amount of charge is transferred between the two entities. Furthermore, the interplay of the ligand field environment, strong p − d hybridization, and correlation effects within the complex are strongly involved in determining the spin state of the iron ion. By bridging molecular chemistry and solid state physics, this study provides first steps towards a joint analysis of the properties of iron-based catalysts from first principles.
Although iron-based catalysts are regarded as a promising alternative to precious metal catalysts, their precise electronic structures during catalysis still pose challenges for computational descriptions. A particularly urgent question is the influence of the environment on the electronic structure, and how to describe this properly with computational methods. Here, we study an iron porphyrin chloride complex adsorbed on a graphene sheet using density functional theory calculations to detail how much the electronic structure is influenced by the presence of a graphene layer. Our results indicate that weak interactions due to van der Waals forces dominate between the porphyrin complex and graphene, and only a small amount of charge is transferred between the two entities. Furthermore, the interplay of the ligand field environment, strong p − d hybridization, and correlation effects within the complex are strongly involved in determining the spin state of the iron ion. By bridging molecular chemistry and solid state physics, this study provides first steps towards a joint analysis of the properties of iron-based catalysts from first principles.
The discovery of the 1144-phase, e.g. CaKFe4As4, creates opportunities to build novel intermetallics with alternative stacking of two parent compounds. Here we formalize the idea by defining a class of bulk crystalline solids with A-B stacking (including 1144-phases and beyond), which is a generalization of hetero-structures from few-layer or thin-film semi-conductors to bulk intermetallics. Theoretically, four families of phosphides \textit{AB}(TM)4P4 (TM=Fe, Ru, Co, Ni) are investigated by first-principles calculations, wherein configurational, vibrational and electronic degrees of freedom are considered. It predicts a variety of stable 1144-phases (especially Ru- and Fe-phosphides). Stability rules are found and structural/electronic properties are discussed. Experimentally, we synthesize high-purity CaKRu4P4 as a proof of principle example. The synthetic method is simple and easily applied. Moreover, it alludes to a strategy to explore complex multi-component compounds, facilitated by a phase diagram coordinated by collective descriptors.
The discovery of the 1144-phase, e.g. CaKFe4As4, creates opportunities to build novel intermetallics with alternative stacking of two parent compounds. Here we formalize the idea and introduce a concept, namely hetero-crystals (HC), to describe a unique class of bulk crystalline solids with such A-B stacking (including 1144-phases and beyond). HC generalizes hetero-structures from few-layer or thin-film semi-conductors to bulk intermetallics. We illustrate the HC concept with the example of 1144-phases. Theoretically, four families of phosphides \textit{AB}(TM)4P4 (TM=Fe, Ru, Co, Ni) are investigated by first-principles calculations, wherein configurational, vibrational and electronic degrees of freedom are considered. It predicts a variety of stable 1144-phases (especially Ru- and Fe-phosphides). Stability rules are found and structural/electronic properties are discussed. Experimentally, we synthesize high-purity CaKRu4P4 as a proof of principle example of such a HC. The synthetic method is simple and can be applied to other HC. Moreover, HC alludes to a new strategy to explore complex multi-component compounds, facilitated by a new phase diagram coordinated by collective descriptors.
Twisted heterostructures of van der Waals materials have received much attention for their many remarkable properties. Here, we present a comprehensive theory of the long-range ordered magnetic phases of twisted bilayer α-RuCl3 via a combination of first-principles calculations and atomistic simulations. While a monolayer exhibits zigzag antiferromagnetic order with three possible ordering wave vectors, a rich phase diagram is obtained for moiré superlattices as a function of interlayer exchange and twist angle. For large twist angles, each layer spontaneously picks a single zigzag ordering wave vector, whereas, for small twist angles, the ground state involves a combination of all three wave vectors in a complex hexagonal domain structure. This multi-domain order minimizes the interlayer energy while enduring the energy cost due to the domain wall formation. Our results indicate that magnetic frustration due to stacking-dependent interlayer exchange in moiré superlattices can be used to tune the magnetic ground state and enhance quantum fluctuations in α-RuCl3.
Evolution of nematic fluctuations in CaK(Fe1−xNix)4As4 with spin-vortex crystal magnetic order
(2020)
The CaK(Fe1−xNix)4As4 superconductors resemble the archetypal 122-type iron-based materials but have a crystal structure with distinctly lower symmetry. This family hosts one of the few examples of the so-called spin-vortex crystal magnetic order, a non-collinear magnetic configuration that preserves tetragonal symmetry, in contrast to the orthorhombic collinear stripe-type magnetic configuration common to the 122-type systems. Thus, nematic order is completely absent from its phase diagram. To investigate the evolution of nematic fluctuations in such a case, we present elastoresistance and elastic modulus measurements in CaK(Fe1−xNix)4As4 (x=0−0.05) combined with phenomenological modeling and density functional theory. We find clear experimental signatures of considerable nematic fluctuations, including softening of the Young's modulus Y[110] and a Curie-Weiss type divergence of the B2g elastoresistance coefficient in CaK(Fe0.951Ni0.049)4As4. Overall, nematic fluctuations within this series bear strong similarities to the hole-doped Ba1−xKxFe2As2 series, including a substitution-induced sign change. Our theoretical analysis addresses the effect of the specific crystal symmetry of the 1144-type structure in determining its magnetic ground state and on the nematic fluctuations.
We investigate the magnetism of a previously unexplored distorted spin-1/2 kagome model consisting of three symmetry-inequivalent nearest-neighbor antiferromagnetic Heisenberg couplings and uncover a rich ground state phase diagram even at the classical level. Using analytical arguments and numerical techniques we identify a collinear Q⃗ =0 magnetic phase, two unusual non-collinear coplanar Q⃗ =(1/3,1/3) phases and a classical spin liquid phase with a degenerate manifold of non-coplanar ground states, resembling the jammed spin liquid phase found in the context of a bond-disordered kagome antiferromagnet. We further show with density functional theory calculations that the recently synthesized Y-kapellasite Y3Cu9(OH)19Cl8 is a realization of this model and predict its ground state to lie in the region of Q⃗ =(1/3,1/3) order, which remains stable even after inclusion of quantum fluctuation effects within variational Monte Carlo and pseudofermion functional renormalization group. Interestingly, the excitation spectrum of Y-kapellasite lies between that of an underlying triangular lattice of hexagons and a kagome lattice of trimers. The presented model opens a new direction in the study of kagome antiferromagnets.
We investigate the magnetism of a previously unexplored distorted spin-1/2 kagome model consisting of three symmetry-inequivalent nearest-neighbor antiferromagnetic Heisenberg couplings Jhexagon, J and J', and uncover a rich ground state phase diagram even at the classical level. Using analytical arguments and numerical techniques we identify a collinear Q = 0 magnetic phase, two unusual non-collinear coplanar Q = (1/3,1/3) phases and a classical spin liquid phase with a degenerate manifold of non-coplanar ground states, resembling the jammed spin liquid phase found in the context of a bond-disordered kagome antiferromagnet. We further show with density functional theory calculations that the recently synthesized Y-kapellasite Y3Cu9(OH)19Cl8 is a realization of this model and predict its ground state to lie in the region of Q = (1/3,1/3) order, which remains stable even after inclusion of quantum fluctuation effects within variational Monte Carlo and pseudofermion functional renormalization group. The presented model opens a new direction in the study of kagome antiferromagnets.
The discovery of superconductivity in layered vanadium-based kagome metals AV3Sb5 (A: K, Rb, Cs) has added a new family of materials to the growing class of possible unconventional superconductors. However, the nature of the superconducting pairing in these materials remains elusive. We present a microscopic theoretical study of the leading superconducting instabilities on the kagome lattice based on spin- and charge-fluctuation mediated Cooper pairing. The applied methodology includes effects of both on-site and nearest-neighbor repulsive Coulomb interactions. Near the upper van Hove filling -- relevant for the AV3Sb5 materials -- we find a rich phase diagram with several pairing symmetries being nearly degenerate. In particular, while a substantial fraction of the phase diagram is occupied by a spin-singlet order parameter transforming as a two-dimensional irreducible representation of the point group, several nodal spin-triplet pairing states remain competitive. We compute the band and interaction parameter-dependence of the hierarchy of the leading superconducting instabilities, and determine the detailed momentum dependence of the resulting preferred gap structures. Crucially, for moderate values of the interaction parameters, the individual pairing states depend strongly on momentum and exhibit multiple nodes on the Fermi surface. We discuss the properties of these superconducting gap structures in light of recent experimental developments of the AV3Sb5 materials.
The Kitaev material α-RuCl3 is among the most prominent candidates to host a quantum spin-liquid state endowed with fractionalized excitations. Recent experimental and theoretical investigations have separately revealed the importance of both the magnetoelastic coupling and the magnetic anisotropy, in dependence of the applied magnetic field direction. In this combined theoretical and experimental research, we investigate the anisotropic magnetic and magnetoelastic properties for magnetic fields applied along the main crystallographic axes as well as for fields canted out of the honeycomb plane. We found that the magnetostriction anisotropy is unusually large compared to the anisotropy of the magnetization, which is related to the strong magnetoelastic Γ′˜-type coupling in our \textit{ab-initio} derived model. We observed large, non-symmetric magnetic anisotropy for magnetic fields canted out of the honeycomb ab-plane in opposite directions, namely towards the +c∗ or −c∗ axes, respectively. The observed directional anisotropy is explained by considering the relative orientation of the magnetic field with respect to the co-aligned RuCl6 octahedra. Magnetostriction measurements in canted fields support this non-symmetric magnetic anisotropy, however these experiments are affected by magnetic torque effects. Comparison of theoretical predictions with experimental findings allow us to recognize the significant contribution of torque effects in experimental setups where α-RuCl3 is placed in canted magnetic fields.
Layered {\alpha}-RuCl3 is a promising material to potentially realize the long-sought Kitaev quantum spin liquid with fractionalized excitations. While evidence of this exotic state has been reported under a modest in-plane magnetic field, such behavior is largely inconsistent with theoretical expectations of Kitaev phases emerging only in out-of-plane fields. These predicted field-induced states have been mostly out of reach due to the strong easy-plane anisotropy of bulk crystals, however. We use a combination of tunneling spectroscopy, magnetotransport, electron diffraction, and ab initio calculations to study the layer-dependent magnons, anisotropy, structure, and exchange coupling in atomically thin samples. Due to structural distortions, the sign of the average off-diagonal exchange changes in monolayer {\alpha}-RuCl3, leading to a reversal of magnetic anisotropy to easy-axis. Our work provides a new avenue to tune the magnetic interactions in {\alpha}-RuCl3 and allows theoretically predicted quantum spin liquid phases for out-of-plane fields to be more experimentally accessible.
Under temperature or pressure tuning, tetragonal EuPd2Si2 is known to undergo a valence transition from nearly divalent to nearly trivalent Eu accompanied by a volume reduction. Albeit intensive work, its origin is not yet completely understood. Here, we investigate the mechanism of the valence transition under volume compression by density functional theory calculations (DFT). Our analysis suggests that the transition is a consequence of an enhanced c−f hybridization between localized Eu 4f states and itinerant conduction states (Eu 5d, Pd 4d, and Si 3p) where the interplay of the electronic bandwidth, crystal field environment, Coulomb repulsion, Hund's coupling and spin-orbit coupling plays a crucial role for the transition to happen. The change in the electronic structure is intimately related to the volume reduction where Eu-Pd(Si) bond lengths shorten. In a next step we compare our DFT results to surface-sensitive photoemission data in which the mixed-valent properties of EuPd2Si2 are reflected in a simultaneous observation of divalent and trivalent signals from the Eu 4f shell.
We investigate the magnetism of a previously unexplored distorted spin-1/2 kagome model consisting of three symmetry-inequivalent nearest-neighbor antiferromagnetic Heisenberg couplings Jhexagon, J and J', and uncover a rich ground state phase diagram even at the classical level. Using analytical arguments and numerical techniques we identify a collinear Q = 0 magnetic phase, two unusual non-collinear coplanar Q = (1/3,1/3) phases and a classical spin liquid phase with a degenerate manifold of non-coplanar ground states, resembling the jammed spin liquid phase found in the context of a bond-disordered kagome antiferromagnet. We further show with density functional theory calculations that the recently synthesized Y-kapellasite Y3Cu9(OH)19Cl8 is a realization of this model and predict its ground state to lie in the region of Q = (1/3,1/3) order, which remains stable even after inclusion of quantum fluctuation effects within variational Monte Carlo and pseudofermion functional renormalization group. The presented model opens a new direction in the study of kagome antiferromagnets.
Due to the small photon momentum, optical spectroscopy commonly probes magnetic excitations only at the center of the Brillouin zone; however, there are ways to override this restriction. In case of the distorted kagome quantum magnet Y-kapellasite, Y3Cu9(OH)19Cl8, under scrutiny here, the spin (magnon) density of states (SDOS) can be accessed over the entire Brillouin zone through three-center magnon excitations. This mechanism is aided by the three different magnetic sublattices and strong short-range correlations in the distorted kagome lattice. The results of THz time-domain experiments agree remarkably well with linear spin-wave theory (LSWT). Relaxing the conventional zone-center constraint of photons gives a new aspect to probe magnetism in matter.
The rich functionalities of transition-metal oxides and their interfaces bear an enormous technological potential. Its realization in practical devices requires, however, a significant improvement of yet relatively low electron mobility in oxide materials. Recently, a mobility boost of about 2 orders of magnitude has been demonstrated at the spinel–perovskite γ-Al2O3/SrTiO3 interface compared to the paradigm perovskite–perovskite LaAlO3/SrTiO3 interface. We explore the fundamental physics behind this phenomenon from direct measurements of the momentum-resolved electronic structure of this interface using resonant soft-X-ray angle-resolved photoemission. We find an anomaly in orbital ordering of the mobile electrons in γ-Al2O3/SrTiO3 which depopulates electron states in the top SrTiO3 layer. This rearrangement of the mobile electron system pushes the electron density away from the interface, which reduces its overlap with the interfacial defects and weakens the electron–phonon interaction, both effects contributing to the mobility boost. A crystal-field analysis shows that the band order alters owing to the symmetry breaking between the spinel γ-Al2O3 and perovskite SrTiO3. Band-order engineering, exploiting the fundamental symmetry properties, emerges as another route to boost the performance of oxide devices.
The families of organic charge-transfer salts 𝜅-(BEDT-TTF)2𝑋 and 𝜅-(BETS)2𝑋, where BEDT-TTF and BETS stand for the organic donor molecules C10H8S8 and C10H8S4Se4, respectively, and X for an inorganic electron acceptor, have been proven to serve as a powerful playground for the investigation of the physics of frustrated Mott insulators. These materials have been ascribed a model character, since the dimerization of the organic molecules allows to map these materials onto a single band Hubbard model, in which the dimers reside on an anisotropic triangular lattice. By changing the inorganic unit X or applying physical pressure, the correlation strength and anisotropy of the triangular lattice can be varied. This has led to the discovery of a variety of exotic phenomena, including quantum-spin liquid states, a plethora of long-range magnetic orders in proximity to a Mott metal-insulator transition, and unconventional superconductivity. While many of these phenomena can be described within this effective one-band Hubbard model on a triangular lattice, it has become evident in recent years that this simplified description is insufficient to capture all observed magnetic and electronic properties. The ingredients for generalized models that are relevant include, but are not limited to, spin-orbit coupling, intra-dimer charge and spin degrees of freedom, electron-lattice coupling, as well as disorder effects. Here, we review selected theoretical and experimental discoveries that clearly demonstrate the relevance thereof. At the same time, we outline that these aspects are not only relevant to this class of organic charge-transfer salts, but are also receiving increasing attention in other classes of inorganic strongly correlated electron systems. This reinforces the model character that the 𝜅-phase organic charge-transfer salts have for understanding and discovering novel phenomena in strongly correlated electron systems from a theoretical and experimental point of view.
The search for materials with topological properties is an ongoing effort. In this article we propose a systematic statistical method, supported by machine learning techniques, that is capable of constructing topological models for a generic lattice without prior knowledge of the phase diagram. By sampling tight-binding parameter vectors from a random distribution, we obtain data sets that we label with the corresponding topological index. This labeled data is then analyzed to extract those parameters most relevant for the topological classification and to find their most likely values. We find that the marginal distributions of the parameters already define a topological model. Additional information is hidden in correlations between parameters. Here we present as a proof of concept the prediction of the Haldane model as the prototypical topological insulator for the honeycomb lattice in Altland-Zirnbauer (AZ) class A. The algorithm is straightforwardly applicable to any other AZ class or lattice, and could be generalized to interacting systems.
Recent density functional theory (DFT) calculations for KFe2As2 have been shown to be insufficient to satisfactorily describe angle-resolved photoemission (ARPES) measurements as well as observed de Haas–van Alphen (dHvA) frequencies. In the present work, we extend DFT calculations based on the full-potential linear augmented plane-wave method by dynamical mean field theory (DFT+DMFT) to include correlation effects beyond the local density approximation. We present results for two sets of reported crystal structures. Our calculations indicate that KFe2As2 is a moderately correlated metal with a mass renormalization factor of the Fe $3{\rm d}$ orbitals between 1.6 and 2.7. Furthermore, the obtained shape and size of the Fermi surface are in good agreement with ARPES measurements and we observe some topological changes with respect to DFT calculations such as the opening of an inner hole cylinder at the Z point. As a result, our calculated dHvA frequencies differ greatly from existing DFT results and qualitatively agree with experimental data. On this basis, we argue that correlation effects are important to understand the -presently under debate- nature of the superconducting state in KFe2As2.
We demonstrate ultra-sharp (≲10 nm) lateral p-n junctions in graphene using electronic transport, scanning tunneling microscopy, and first principles calculations. The p-n junction lies at the boundary between differentially-doped regions of a graphene sheet, where one side is intrinsic and the other is charge-doped by proximity to a flake of α-RuCl3 across a thin insulating barrier. We extract the p-n junction contribution to the device resistance to place bounds on the junction width. We achieve an ultra-sharp junction when the boundary between the intrinsic and doped regions is defined by a cleaved crystalline edge of α-RuCl3 located 2 nm from the graphene. Scanning tunneling spectroscopy in heterostructures of graphene, hexagonal boron nitride, and α-RuCl3 shows potential variations on a sub-10 nm length scale. First principles calculations reveal the charge-doping of graphene decays sharply over just nanometers from the edge of the α-RuCl3 flake.
Strontium ruthenate Sr2RuO4 is an unconventional superconductor whose pairing symmetry has not been fully clarified, despite more than two decades of intensive research. Recent NMR Knight shift experiments have rekindled the Sr2RuO4 pairing debate by giving strong evidence against all odd-parity pairing states, including chiral p-wave pairing that was for a long time the leading pairing candidate. Here, we exclude additional pairing states by analyzing recent elastocaloric measurements [YS. Li et al., Nature 607, 276--280 (2022)]. To be able to explain the elastocaloric experiment, we find that unconventional even-parity pairings must include either large dx2−y2-wave or large {dxz∣dyz}-wave admixtures, where the latter possibility arises because of the body-centered point group symmetry. These {dxz∣dyz}-wave admixtures take the form of distinctively body-centered-periodic harmonics that have horizontal line nodes. Hence gxy(x2−y2)-wave and dxy-wave pairings are excluded as possible dominant even pairing states.
Electronic systems living on Archimedean lattices such as kagome and square–octagon networks are presently being intensively discussed for the possible realization of topological insulating phases. Coining the most interesting electronic topological states in an unbiased way is however not straightforward due to the large parameter space of possible Hamiltonians. A possible approach to tackle this problem is provided by a recently developed statistical learning method (Mertz and Valentí in Phys Rev Res 3:013132, 2021. https://doi.org/10.1103/PhysRevResearch.3.013132), based on the analysis of a large data sets of randomized tight-binding Hamiltonians labeled with a topological index. In this work, we complement this technique by introducing a feature engineering approach which helps identifying polynomial combinations of Hamiltonian parameters that are associated with non-trivial topological states. As a showcase, we employ this method to investigate the possible topological phases that can manifest on the square–octagon lattice, focusing on the case in which the Fermi level of the system lies at a high-order van Hove singularity, in analogy to recent studies of topological phases on the kagome lattice at the van Hove filling.
Continued advances in quantum technologies rely on producing nanometer-scale wires. Although several state-of-the-art nanolithographic technologies and bottom-up synthesis processes have been used to engineer these wires, critical challenges remain in growing uniform atomic-scale crystalline wires and constructing their network structures. Here, we discover a simple method to fabricate atomic-scale wires with various arrangements, including stripes, X-junctions, Y-junctions, and nanorings. Single-crystalline atomic-scale wires of a Mott insulator, whose bandgap is comparable to those of wide-gap semiconductors, are spontaneously grown on graphite substrates by pulsed-laser deposition. These wires are one unit cell thick and have an exact width of two and four unit cells (1.4 and 2.8 nm) and lengths up to a few micrometers. We show that the nonequilibrium reaction-diffusion processes may play an essential role in atomic pattern formation. Our findings offer a previously unknown perspective on the nonequilibrium self-organization phenomena on an atomic scale, paving a unique way for the quantum architecture of nano-network.
In the search for novel organic charge transfer salts with variable degrees of charge transfer we have studied the effects of two modifications of the recently synthesized donor–acceptor system [tetramethoxypyrene (TMP)]–[tetracyanoquinodimethane (TCNQ)]. One is of chemical nature by substituting the acceptor TCNQ molecules by F4TCNQ molecules. The second consists in simulating the application of uniaxial pressure along the stacking axis of the system. In order to test the chemical substitution, we have grown single crystals of the TMP–F4TCNQ complex and analyzed its electronic structure via electronic transport measurements, ab initio density functional theory (DFT) calculations and UV/VIS/IR absorption spectroscopy. This system shows an almost ideal geometrical overlap of nearly planar molecules stacked alternately (mixed stack) and this arrangement is echoed by a semiconductor-like transport behavior with an increased conductivity along the stacking direction. This is in contrast to TMP–TCNQ which shows a less pronounced anisotropy and a smaller conductivity response. Our band structure calculations confirm the one-dimensional behavior of TMP–F4TCNQ with pronounced dispersion only along the stacking axis. Infrared measurements illustrating the C[triple bond, length as m-dash]N vibration frequency shift in F4TCNQ suggest however no improvement in the degree of charge transfer in TMP–F4TCNQ with respect to TMP–TCNQ. In both complexes about 0.1e is transferred from TMP to the acceptor. Concerning the pressure effect, our DFT calculations on the designed TMP–TCNQ and TMP–F4TCNQ structures under different pressure conditions show that application of uniaxial pressure along the stacking axis of TMP–TCNQ may be the route to follow in order to obtain a much more pronounced charge transfer.
RuO₂: a puzzle to be solved
(2023)
Altermagnetism is a topic that has lately been gaining attention and the RuO2 compound is among one of the most studied altermagnetic candidates. However, the survey of available literature on RuO2 properties suggests that there is no consensus about the magnetism of this material. By performing density functional theory calculations, we show that the electronic properties of stoichiometric RuO2 are described in terms of a smaller Hubbard U within DFT+U than the value required to have magnetism. We further argue that Ru vacancies can actually aid the formation of a magnetic state in RuO2. This in turn suggests that a characterization of the amount of Ru vacancies in experimental samples might help the resolution of the controversy between the different experimental results.
Topological semimetal antiferromagnets provide a rich source of exotic topological states which can be controlled by manipulating the orientation of the Néel vector, or by modulating the lattice parameters through strain. We investigate via ab initio density functional theory calculations, the effects of shear strain on the bulk and surface states n two antiferromagnetic EuCd2As2 phases with out-of-plane and in-plane spin configurations. When magnetic moments are along the c-axis, a 3% longitudinal or diagonal shear strain can tune the Dirac semimetal phase to an axion insulator phase, characterized by the parity-based invariant η4I=2. For an in-plane magnetic order, the axion insulator phase remains robust under all shear strains. We further find that for both magnetic orders, the bulk gap increases and a surface gap opens on the (001) surface up to 16 meV. Because of a nonzero η4I index and gapped states on the (001) surface, hinge modes are expected to happen on the side surface states between those gapped surface states. This result can provide a valuable insight in the realization of the long-sought axion states.