Researcher profile

Efthimios Kaxiras

Efthimios Kaxiras contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
23works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

23 published item(s)

preprint2022arXiv

Hydrogen Atoms on Zigzag Graphene Nanoribbons: Chemistry and Magnetism Meet at the Edge

Although the unconventional $π$-magnetism at the zigzag edges of graphene holds promise for a wide array of applications, whether and to what degree it plays a role in their chemistry remains poorly understood. Here, we investigate the addition of a hydrogen atom $-$ the simplest yet the most experimentally relevant adsorbate $-$ to zigzag graphene nanoribbons (ZGNRs). We show that the $π$-magnetism governs the chemistry of ZGNRs, giving rise to a site-dependent reactivity of the carbon atoms and driving the hydrogenation process to the nanoribbon edges. Conversely, the chemisorbed hydrogen atom governs the $π$-magnetism of ZGNRs, acting as a spin-$\frac{1}{2}$ paramagnetic center in the otherwise antiferromagnetic ground state and spin-polarizing the charge carriers at the band extrema. Our findings establish a comprehensive picture of the peculiar interplay between chemistry and magnetism that emerges at the zigzag edges of graphene.

preprint2022arXiv

Non-Abelian topological defects and strain mapping in 2D moiré materials

We present a general method to analyze the topological nature of the domain boundary connectivity that appeared in relaxed moiré superlattice patterns at the interface of 2-dimensional (2D) van der Waals (vdW) materials. At large enough moiré lengths, all moiré systems relax into commensurated 2D domains separated by networks of dislocation lines. The nodes of the 2D dislocation line network can be considered as vortex-like topological defects. We find that a simple analogy to common topological systems with an $S^1$ order parameter, such as a superconductor or planar ferromagnet, cannot correctly capture the topological nature of these defects. For example, in twisted bilayer graphene, the order parameter space for the relaxed moiré system is homotopy equivalent to a punctured torus. Here, the nodes of the 2D dislocation network can be characterized as elements of the fundamental group of the punctured torus, the free group on two generators, endowing these network nodes with non-Abelian properties. Extending this analysis to consider moiré patterns generated from any relative strain, we find that antivortices occur in the presence of anisotropic heterostrain, such as shear or anisotropic expansion, while arrays of vortices appear under twist or isotropic expansion between vdW materials. Experimentally, utilizing the dark field imaging capability of transmission electron microscopy (TEM), we demonstrate the existence of vortex and antivortex pair formation in a moiré system, caused by competition between different types of heterostrains in the vdW interfaces. We also present a methodology for mapping the underlying heterostrain of a moiré structure from experimental TEM data, which provides a quantitative relation between the various components of heterostrain and vortex-antivortex density in moiré systems.

preprint2022arXiv

Outstanding thermoelectric performance predicted for out-of-plane p-doped GeSe

The record-breaking thermoelectric performance of tin selenide (SnSe) has motivated the investigation of analogue compounds with the same structure. A promising candidate that emerged recently is germanium selenide (GeSe). Here, using extensive first-principles calculations of the hole-phonon and hole-impurity scattering, we investigate the thermoelectric transport properties of the orthorhombic phase of p-doped GeSe. We predict outstanding thermoelectric performance for GeSe over a broad range of temperatures due to its high Seebeck coefficients, extremely low Lorenz numbers, ultralow total thermal conductivity, and relatively large band gap. In particular, the out-of-plane direction in GeSe presents equivalent or even higher performance than SnSe for temperatures above 500 K. By extending the analysis to 900 K, we obtained an ultrahigh value for the thermoelectric figure of merit (zT = 3.2) at the optimal hole density of 4x10^19 cm^-3. Our work provides strong motivation for continued experimental work focusing on improving the GeSe doping efficiency in order to achieve this optimal hole density.

preprint2022arXiv

Twistronics of Janus transition metal dichalcogenide bilayers

Twisted multilayers of two-dimensional (2D) materials are an increasingly important platform for investigating quantum phases of matter, and in particular, strongly correlated electrons. The moiré pattern introduced by the relative twist between layers creates effective potentials of long-wavelength, leading to electron localization. However, in contrast to the abundance of 2D materials, few twisted heterostructures have been studied until now. Here we develop a first-principle continuum theory to study the electronic bands introduced by moire patterns of twisted Janus transition metal dichalcogenides (TMD) homo- and hetero-bilayers. The model includes lattice relaxation, stacking-dependent effective mass, and Rashba spin-orbit coupling. We then perform a high-throughput generation and characterization of DFT-extracted continuum models for more than a hundred possible combinations of materials and stackings. Our model predicts that the moiré physics and emergent symmetries depend on chemical composition, vertical layer orientation, and twist angle, so that the minibands wavefunctions can form triangular, honeycomb, and Kagome networks. Rashba spin-orbit effects, peculiar of these systems, can dominate the moiré bandwidth at small angles. Our work enables the detailed investigation of Janus twisted heterostructures, allowing the discovery and control of novel electronic phenomena.

preprint2021arXiv

Ferromagnetic helical nodal line and Kane-Mele spin-orbit coupling in kagome metal Fe3Sn2

The two-dimensional kagome lattice hosts Dirac fermions at its Brillouin zone corners K and K', analogous to the honeycomb lattice. In the density functional theory electronic structure of ferromagnetic kagome metal Fe$_3$Sn$_2$, without spin-orbit coupling we identify two energetically split helical nodal lines winding along $z$ in the vicinity of K and K' resulting from the trigonal stacking of the kagome layers. We find that hopping across A-A stacking introduces a layer splitting in energy while that across A-B stacking controls the momentum space amplitude of the helical nodal lines. The effect of spin-orbit coupling is found to resemble that of a Kane-Mele term, where the nodal lines can either be fully gapped to quasi-two-dimensional massive Dirac fermions, or remain gapless at discrete Weyl points depending on the ferromagnetic moment orientation. Aside from numerically establishing Fe$_3$Sn$_2$ as a model Dirac kagome metal, our results provide insights into materials design of topological phases from the lattice point of view, where paradigmatic low dimensional lattice models often find realizations in crystalline materials with three-dimensional stacking.

preprint2021arXiv

Imprinting Tunable $π$-Magnetism in Graphene Nanoribbons via Edge Extensions

Magnetic carbon nanostructures are currently under scrutiny for a wide spectrum of applications. Here, we theoretically investigate armchair graphene nanoribbons patterned with asymmetric edge extensions consisting of laterally fused naphtho groups, as recently fabricated via on-surface synthesis. We show that an individual edge extension acts as a spin-$\frac{1}{2}$ center and develops a sizable spin-polarization of the conductance around the band edges. The Heisenberg exchange coupling between a pair of edge extensions is dictated by the position of the second naphtho group in the carbon backbone, thus enabling ferromagnetic, antiferromagnetic, or non-magnetic states. The periodic arrangement of edge extensions yields full spin-polarization at the band extrema, and the accompanying ferromagnetic ground state can be driven into non-magnetic or antiferromagnetic phases through external stimuli. Overall, our work reveals precise tunability of the $π$-magnetism in graphene nanoribbons induced by naphtho groups, thereby establishing these one-dimensional architectures as suitable platforms for logic spintronics.

preprint2021arXiv

Modeling the effect of the vaccination campaign on the Covid-19 pandemic

Population-wide vaccination is critical for containing the SARS-CoV-2 (Covid-19) pandemic when combined with restrictive and prevention measures. In this study, we introduce SAIVR, a mathematical model able to forecast the Covid-19 epidemic evolution during the vaccination campaign. SAIVR extends the widely used Susceptible-Infectious-Removed (SIR) model by considering the Asymptomatic (A) and Vaccinated (V) compartments. The model contains several parameters and initial conditions that are estimated by employing a semi-supervised machine learning procedure. After training an unsupervised neural network to solve the SAIVR differential equations, a supervised framework then estimates the optimal conditions and parameters that best fit recent infectious curves of 27 countries. Instructed by these results, we performed an extensive study on the temporal evolution of the pandemic under varying values of roll-out daily rates, vaccine efficacy, and a broad range of societal vaccine hesitancy/denial levels. The concept of herd immunity is questioned by studying future scenarios which involve different vaccination efforts and more infectious Covid-19 variants.

preprint2021arXiv

Twistons in a Sea of Magic

Magic angle twisted trilayer graphene (TTG) has recently emerged as a new platform to engineer strongly correlated flat bands. Here, we reveal the structural and electronic properties of TTG using low temperature scanning tunneling microscopy at twist angles for which superconductivity has been observed. Real trilayer samples deviate from their idealized structure due to a strong reconstruction of the moiré lattice, which locks layers into near-magic angle, mirror symmetric domains comparable in size to the superconducting coherence length. The price for this magic relaxation is the introduction of an array of localized twist angle faults, termed twistons. These novel, gate-tunable moiré defects offer a natural explanation for the superconducting dome observed in transport and provide an avenue to probe superconducting pairing mechanisms through disorder tuning.

preprint2020arXiv

Boosting the efficiency of ab initio electron-phonon coupling calculations through dual interpolation

The coupling between electrons and phonons in solids plays a central role in describing many phenomena, including superconductivity and thermoelecric transport. Calculations of this coupling are exceedingly demanding as they necessitate integrations over both the electron and phonon momenta, both of which span the Brillouin zone of the crystal, independently. We present here an ab initio method for efficiently calculating electron-phonon mediated transport properties by dramatically accelerating the computation of the double integrals with a dual interpolation technique that combines maximally localized Wannier functions with symmetry-adapted plane waves. The performance gain in relation to the current state-of-the-art Wannier-Fourier interpolation is approximately 2n_s \times M, where n_s is the number of crystal symmetry operations and M, a number in the range 5 - 60, governs the expansion in star functions. We demonstrate with several examples how our method performs some ab initio calculations involving electron-phonon interactions.

preprint2020arXiv

Duality between atomic configurations and Bloch states in twistronic materials

The relative orientation (twist) of successive layers of stacked two-dimensional (2D) materials creates variations in the interlayer atomic registry. The variations often form a super lattice, called a moiré pattern, which can alter electronic properties. In this work we introduce a classification of the single-particle electronic structures that can occur in twisted stacks of 2D layers by characterizing them as "moiré molecules" or "moiré crystals". The molecules generate localized electronic states and moiré flat bands, while the crystals are sometimes unconventional and produce electronic banding in the configuration basis. The underpinning of this classification is the duality between interlayer configuration and monolayer Bloch momentum in moiré Hamiltonians. We apply this understanding to diagrams of local electron density in untwisted geometries to produce intuitive and quantitative predictions of twistronic properties. We provide a conceptual introduction to this framework through a one-dimensional model, and then apply it to 2D twisted bilayers of the semi-metal graphene, and of MoS$_2$, a representative material of the transition metal dichalcogenide (TMDC) family of semiconductors. This level of thorough understanding of twistronic phenomena is vital in the search for new material platforms for localized moiré electrons.

preprint2020arXiv

Effects of Structural Distortions on the Electronic Structure of T-type Transition Metal Dichalcogenides

Single-layer transition metal dichalcogenides (TMDCs) can adopt two distinct structures corresponding to different coordination of the metal atoms. TMDCs adopting the T-type structure exhibit a rich and diverse set of phenomena, including charge density waves (CDW) in a $\sqrt{13}\times\sqrt{13}$ supercell pattern in TaS$_2$ and TaSe$_2$, and a possible excitonic insulating phase in TiSe$_2$. These properties make the T-TMDCs desirable components of layered heterostructure devices. In order to predict the emergent properties of combinations of different layered materials, one needs simple and accurate models for the constituent layers which can take into account potential effects of lattice mismatch, relaxation, strain, and structural distortion. Previous studies have developed ab initio tight-binding Hamiltonians for H-type TMDCs [arXiv:1709.07510]. Here we extend this work to include T-type TMDCs. We demonstrate the capabilities of our model using three example systems: a 1-dimensional sinusoidal ripple, the 2$\times$2 CDW in TiSe$_2$, and the $\sqrt{13}\times\sqrt{13}$ CDW in TaS$_2$. Using the technique of band unfolding we compare the electronic structure of the distorted crystals to the pristine band structure and find excellent agreement with direct DFT calculations, provided the magnitude of the distortions remains in the linear regime.

preprint2020arXiv

Electronic structure calculations of twisted multi-layer graphene superlattices

Quantum confinement endows two-dimensional (2D) layered materials with exceptional physics and novel properties compared to their bulk counterparts. Although certain two- and few-layer configurations of graphene have been realized and studied, a systematic investigation of the properties of arbitrarily layered graphene assemblies is still lacking. We introduce theoretical concepts and methods for the processing of materials information, and as a case study, apply them to investigate the electronic structure of multi-layer graphene-based assemblies in a high-throughput fashion. We provide a critical discussion of patterns and trends in tight binding band structures and we identify specific layered assemblies using low-dispersion electronic bands as indicators of potentially interesting physics like strongly correlated behavior. A combination of data-driven models for visualization and prediction is used to intelligently explore the materials space. This work more generally aims to increase confidence in the combined use of physics-based and data-driven modeling for the systematic refinement of knowledge about 2D layered materials, with implications for the development of novel quantum devices.

preprint2020arXiv

How carbon vacancies can affect the properties of group IV color centers in diamond: A study of thermodynamics and kinetics

Recently there has been much interest in using Group IV elements from the Periodic Table to fabricate and study X$V$ color centers in diamond where X = Si, Ge, Sn, or Pb and $V$ is a carbon vacancy. These Group IV color centers have a number of interesting spin and optical properties which could potentially make them better candidates than N$V^-$ centers for important applications in quantum computing and quantum information processing. Unfortunately, the very same ion implantation process that is required to create these X$V$ color centers in diamond necessarily also produces many carbon vacancies ($V_{\rm C}$) which can form complexes with these color centers ($V_{\rm C}-$X$V$) that can dramatically affect the properties of the isolated X$V$ color centers. The main focus of this work is to use density-functional theory (DFT) to study the thermodynamics and kinetics of the formation of these $V_{\rm C}-$X$V$ complexes and to suggest experimental ways to impede this process such as varying the Fermi level of the host diamond material through chemical doping or applying an external electrical bias. We also include a discussion of how the simple presence of many $V_{\rm C}$ can negatively impact the spin coherence times ($T_2$) of Group IV color centers through the presence of acoustic phonons.

preprint2020arXiv

LAN -- A materials notation for 2D layered assemblies

Two-dimensional (2D) layered materials offer intriguing possibilities for novel physics and applications. Before any attempt at exploring the materials space in a systematic fashion, or combining insights from theory, computation and experiment, a formal description of information about an assembly of arbitrary composition is required. Here, we introduce a domain-generic notation that is used to describe the space of 2D layered materials from monolayers to twisted assemblies of arbitrary composition, existent or not-yet-fabricated. The notation corresponds to a theoretical materials concept of stepwise assembly of layered structures using a sequence of rotation, vertical stacking, and other operations on individual 2D layers. Its scope is demonstrated with a number of example structures using common single-layer materials as building blocks. This work overall aims to contribute to the systematic codification, capture and transfer of materials knowledge in the area of 2D layered materials.

preprint2020arXiv

Moiré metrology of energy landscapes in van der Waals heterostructures

The emerging field of twistronics, which harnesses the twist angle between two-dimensional materials, represents a promising route for the design of quantum materials, as the twist-angle-induced superlattices offer means to control topology and strong correlations. At the small twist limit, and particularly under strain, as atomic relaxation prevails, the emergent moiré superlattice encodes elusive insights into the local interlayer interaction. Here we introduce moiré metrology as a combined experiment-theory framework to probe the stacking energy landscape of bilayer structures at the 0.1 meV/atom scale, outperforming the gold-standard of quantum chemistry. Through studying the shapes of moiré domains with numerous nano-imaging techniques, and correlating with multi-scale modelling, we assess and refine first-principle models for the interlayer interaction. We document the prowess of moiré metrology for three representative twisted systems: bilayer graphene, double bilayer graphene and H-stacked $MoSe_2/WSe_2$. Moiré metrology establishes sought after experimental benchmarks for interlayer interaction, thus enabling accurate modelling of twisted multilayers.

preprint2020arXiv

Moiré of Moiré: Modeling Mechanical Relaxation in Incommensurate Trilayer van der Waals Heterostructures

The incommensurate stacking of multi-layered two-dimensional materials is a challenging problem from a theoretical perspective and an intriguing avenue for manipulating their physical properties. Here we present a multi-scale model to obtain the mechanical relaxation pattern of twisted trilayer van der Waals (vdW) heterostructures with two independent twist angles, a generally incommensurate system without a supercell description. We adopt the configuration space as a natural description of such incommensurate layered materials, based on the local environment of atomic positions, bypassing the need for commensurate approximations. To obtain the relaxation pattern, we perform energy minimization with respect to the relaxation displacement vectors. We use a continuum model in combination with the Generalized Stacking Fault energy to describe the interlayer coupling, obtained from first-principles calculations based on Density Functional Theory. We show that the relaxation patterns of twisted trilayer graphene and $\mathrm{WSe_2}$ are "moiré of moiré", as a result of the incommensurate coupling two bilayer moiré patterns. We also show that, in contrast to the symmetry-preserving in-plane relaxation in twisted bilayers, trilayer relaxation can break the two-fold rotational symmetry about the xy-plane when the two twist angles are equal.

preprint2020arXiv

The first 100 days: modeling the evolution of the COVID-19 pandemic

A simple analytical model for modeling the evolution of the 2020 COVID-19 pandemic is presented. The model is based on the numerical solution of the widely used Susceptible-Infectious-Removed (SIR) populations model for describing epidemics. We consider an expanded version of the original Kermack-McKendrick model, which includes a decaying value of the parameter $β$ (the effective contact rate) due to externally imposed conditions, to which we refer as the forced-SIR (FSIR) model. We introduce an approximate analytical solution to the differential equations that represent the FSIR model which gives very reasonable fits to real data for a number of countries over a period of 100 days (from the first onset of exponential increase, in China). The proposed model contains 3 adjustable parameters which are obtained by fitting actual data (up to April 28, 2020). We analyze these results to infer the physical meaning of the parameters involved. We use the model to make predictions about the total expected number of infections in each country as well as the date when the number of infections will have reached 99% of this total. We also compare key findings of the model with recently reported results on the high contagiousness and rapid spread of the disease.

preprint2020arXiv

Topological flat bands in frustrated kagome lattice CoSn

Electronic flat bands in momentum space, arising from strong localization of electrons in real space, are an ideal stage to realize strong correlation phenomena. In certain lattices with built-in geometrical frustration, electronic confinement and flat bands can naturally arise from the destructive interference of electronic hopping pathways. Such lattice-borne flat bands are often endowed with nontrivial topology if combined with spin-orbit coupling, while their experimental realization in condensed matter system has been elusive so far. Here, we report the direct observation of topological flat bands in the vicinity of the Fermi level in frustrated kagome system CoSn, using angle-resolved photoemission spectroscopy and band structure calculations. The flat band manifests itself as a dispersionless electronic excitation along the G-M high symmetry direction, with an order of magnitude lower bandwidth (below 150 meV) compared to the Dirac bands originating from the same orbitals. The frustration-driven nature of the flat band is directly confirmed by the real-space chiral d-orbital texture of the corresponding effective Wannier wave functions. Spin-orbit coupling opens a large gap of 80 meV at the quadratic band touching point between the Dirac and flat bands, endowing a nonzero Z2 topological invariant to the flat band in the two-dimensional Brillouin zone. Our observation of lattice-driven topological flat band opens a promising route to engineer novel emergent phases of matter at the crossroad between strong correlation physics and electronic topology.

preprint2020arXiv

Unifying framework for strong and fragile liquids via machine learning: a study of liquid silica

The fragility of a glassforming liquid characterizes how rapidly its relaxation dynamics slow down with cooling. The viscosity of strong liquids follows an Arrhenius law with a temperature-independent barrier height to rearrangements responsible for relaxation, whereas fragile liquids experience a much faster increase in their dynamics, suggesting a barrier height that increases with decreasing temperature. Strong glassformers are typically network glasses, while fragile glassformers are typically molecular or hard-sphere-like. As a result of these differences at the microscopic level, strong and fragile glassformers are usually treated separately from a theoretical point of view. Silica is the archetypal strong glassformer at low temperatures, but also exhibits a mysterious strong-to-fragile crossover at higher temperatures. Here we show that softness, a structure-based machine learned parameter that has previously been applied to fragile glassformers provides a useful description of model liquid silica in the strong and fragile regimes, and through the strong-to-fragile crossover. Just as for fragile glassformers, the relationship between softness and dynamics is invariant and Arrhenius in all regimes, but the average softness changes with temperature. The strong-to-fragile crossover in silica is not due to a sudden, qualitative change in structure, but can be explained by a simple Arrhenius form with a continuously and linearly changing local structure. Our results unify the study of liquid silica under a single simple conceptual picture.

preprint2019arXiv

Coexistence of ultraheavy and ultrarelativistic Dirac quasiparticles in sandwiched trilayer graphene

Electrons in quantum materials exhibiting coexistence of dispersionless (flat) bands piercing dispersive (steep) bands can give rise to strongly correlated phenomena, and are associated with unconventional superconductivity. It is known that in twisted trilayer graphene steep Dirac cones can coexist with band flattening, but the phenomenon is not stable under layer misalignments. Here we show that such a twisted sandwiched graphene (TSWG) -- a three-layer van der Waals heterostructure with a twisted middle layer -- can have very stable flat bands coexisting with Dirac cones near the Fermi energy when twisted to 1.5 degrees. These flat bands require a specific high-symmetry stacking order, and our atomistic calculations predict that TSWG always relaxes to it. Additionally, with external fields, we can control the relative energy offset between the Dirac cone vertex and the flat bands. Our work establishes twisted sandwiched graphene as a new platform for research into strongly interacting phases, and topological transport beyond Dirac and Weyl semimetals.

preprint2019arXiv

Creating Weyl nodes and controlling their energy by magnetization rotation

As they do not rely on the presence of any crystal symmetry, Weyl nodes are robust topological features of an electronic structure that can occur at any momentum and energy. Acting as sinks and sources of Berry curvature, Weyl nodes have been predicted to strongly affect the transverse electronic response, like in the anomalous Hall or Nernst effects. However, to observe large anomalous effects the Weyl nodes need to be close to or at the Fermi-level, which implies the band structure must be tuned by an external parameter, e.g. chemical doping or pressure. Here we show that in a ferromagnetic metal tuning of the Weyl node energy and momentum can be achieved by rotation of the magnetization. Taking Co$_3$Sn$_2$S$_2$ as an example, we use electronic structure calculations based on density-functional theory to show that not only new Weyl fermions can be created by canting the magnetization away from the easy axis, but also that the Weyl nodes can be driven exactly to the Fermi surface. We also show that the dynamics in energy and momentum of the Weyl nodes strongly affect the calculated anomalous Hall and Nernst conductivities.

preprint2019arXiv

Effects of Lithium Intercalation in Twisted Bilayer Graphene

We investigate the effects of lithium intercalation in twisted bilayers of graphene, using first-principles electronic structure calculations. To model this system we employ commensurate supercells that correspond to twist angles of 7.34$^\circ$ and 2.45$^\circ$. From the energetics of lithium absorption we demonstrate that for low Li concentration the intercalants cluster in the AA regions with double the density of a uniform distribution. The charge donated by the Li atoms to the graphene layers results in modifications to the band structure that can be qualitatively captured using a continuum model with modified interlayer couplings in a region of parameter space that has yet to be explored either experimentally or theoretically. Thus, the combination of intercalation and twisted layers simultaneously provides the means for spatial control over material properties and an additional knob with which to tune moiré physics in twisted bilayers of graphene, with potential applications ranging from energy storage and conversion to quantum information.

preprint2019arXiv

Imaginary-time time-dependent density functional theory and its application for robust convergence of electronic states

Reliable and robust convergence to the electronic ground state within density functional theory (DFT) Kohn-Sham (KS) calculations remains a thorny issue in many systems of interest. In such cases, charge sloshing can delay or completely hinder the convergence. Here, we use an approach based on transforming the time-dependent DFT equations to imaginary time, followed by imaginary-time evolution, as a reliable alternative to the self-consistent field (SCF) procedure for determining the KS ground state. We discuss the theoretical and technical aspects of this approach and show that the KS ground state should be expected to be the long-imaginary-time output of the evolution, independent of the exchange-correlation functional or the level of theory used to simulate the system. By maintaining self-consistency between the single-particle wavefunctions and the electronic density throughout the determination of the stationary state, our method avoids the typical difficulties encountered in SCF. To demonstrate dependability of our approach, we apply it to selected systems which struggle to converge with SCF schemes. In addition, through the van Leeuwen theorem, we affirm the physical meaningfulness of imaginary time TDDFT, justifying its use in certain topics of statistical mechanics such as in computing imaginary time path integrals.