Researcher profile

Arash A. Mostofi

Arash A. Mostofi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
5topics
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

11 published item(s)

preprint2022arXiv

Atomistic hartree theory and crystal field of twisted double bilayer graphene near the magic angle

Twisted double bilayer graphene (tDBLG) is a moiré material that has recently generated significant interest because of the observation of correlated phases near the magic angle. We carry out atomistic Hartree theory calculations to study the role of electron-electron interactions in the normal state. In contrast to twisted bilayer graphene (tBLG), we find that such interactions do not result in significant doping-dependent deformations of the electronic band structure. However, interactions play an important role for the electronic structure in the presence of a perpendicular electric field as they screen the external field. Finally, we analyze the contribution of the Hartree potential to the crystal field, i.e. the on-site energy difference between the inner and outer layers. We find that the on-site energy obtained from Hartree theory has the same sign, but a smaller magnitude compared to previous studies in which the on-site energy was determined by fitting tight-binding results to ab initio density-functional theory (DFT) band structures. To understand this quantitative difference, we analyze the ab initio Kohn-Sham potential obtained from DFT and find that a subtle interplay of electron-electron and electron-ion interactions determines the magnitude of the on-site potential.

preprint2022arXiv

Chiral valley phonons and flat phonon bands in moiré materials

We investigate the chirality of phonon modes in twisted bilayer WSe2 and demonstrate distinct chiral behavior of the $K/K^\prime$ valley phonons for twist angles close to $0^{\circ}$ and close to $60^{\circ}$. In particular, multiple chiral non-degenerate $K/K^\prime$ valley phonons are found for twist angles near $60^{\circ}$ whereas no non-degenerate chiral modes are found for twist angles close to $0^\circ$. Moreover, we discover two sets of emergent chiral valley modes that originate from an inversion symmetry breaking at the moiré scale and find similar modes in moiré patterns of strain-engineered bilayers WSe2 and MoSe2/WSe2 heterostructures. At the energy gap between acoustic and optical modes, the formation of flat phonon bands for a broad range of twist angles is observed in the twisted bilayer WSe2. Our findings are relevant for understanding electron-phonon and exciton-phonon scattering in moiré materials and also for the design of phononic analogues of flat band electrons.

preprint2022arXiv

Unconventional Superconductivity in Magic-Angle Twisted Trilayer Graphene

Magic-angle twisted trilayer graphene (MATTG) recently emerged as a highly tunable platform for studying correlated phases of matter, such as correlated insulators and superconductivity. Superconductivity occurs in a range of doping levels that is bounded by van Hove singularities which stimulates the debate of the origin and nature of superconductivity in this material. In this work, we discuss the role of spin-fluctuations arising from atomic-scale correlations in MATTG for the superconducting state. We show that in a phase diagram as function of doping ($ν$) and temperature, nematic superconducting regions are surrounded by ferromagnetic states and that a superconducting dome with $T_c \approx 2\,\mathrm{K}$ appears between the integer fillings $ν=-2$ and $ν= -3$. Applying a perpendicular electric field enhances superconductivity on the electron-doped side which we relate to changes in the spin-fluctuation spectrum. We show that the nematic unconventional superconductivity leads to pronounced signatures in the local density of states detectable by scanning tunneling spectroscopy measurements.

preprint2021arXiv

Structural origins of the infamous "Low Temperature Orthorhombic" to "Low Temperature Tetragonal" phase transition in high-Tc cuprates

We undertake a detailed high-resolution diffraction study of a novel plain band insulator, La$_2$MgO$_4$, which may be viewed as a structural surrogate system of the undoped end-member of the high-T$_c$ superconductors, La$_{2-x-y}$A$^{2+}_x$RE$^{3+}_y$CuO$_{4}$ (A = Ba, Sr, RE= Rare Earth). We find that La$_2$MgO$_4$ exhibits the infamous low-temperature orthorhombic (LTO) to low-temperature tetragonal (LTT) phase transition that has been linked to the suppression of superconductivity in a variety of underdoped cuprates, including the well known La$_{2-x}$Ba$_{x}$CuO$_4$ ($x=0.125$). Furthermore, we find that the LTO-to-LTT phase transition in La$_2$MgO$_4$ occurs for an octahedral tilt angle in the 4 $^{\circ}$ to 5 $^{\circ}$ range, similar to that which has previously been identified as a critical tipping point for superconductivity in these systems. We show that this phase transition, occurring in a system lacking spin correlations and competing electronic states such as charge-density waves and superconductivity, can be understood by simply navigating the density-functional theory ground-state energy landscape as a function of the order parameter amplitude. This result calls for a careful re-investigation of the origins of the phase transitions in high-T$_c$ superconductors based on the hole-doped, $n = 1$ Ruddelsden-Popper lanthanum cuprates.

preprint2020arXiv

Automated high-throughput Wannierisation

Maximally-localised Wannier functions (MLWFs) are routinely used to compute from first-principles advanced materials properties that require very dense Brillouin zone integration and to build accurate tight-binding models for scale-bridging simulations. At the same time, high-throughput (HT) computational materials design is an emergent field that promises to accelerate the reliable and cost-effective design and optimisation of new materials with target properties. The use of MLWFs in HT workflows has been hampered by the fact that generating MLWFs automatically and robustly without any user intervention and for arbitrary materials is, in general, very challenging. We address this problem directly by proposing a procedure for automatically generating MLWFs for HT frameworks. Our approach is based on the selected columns of the density matrix method (SCDM) and we present the details of its implementation in an AiiDA workflow. We apply our approach to a dataset of 200 bulk crystalline materials that span a wide structural and chemical space. We assess the quality of our MLWFs in terms of the accuracy of the band-structure interpolation that they provide as compared to the band-structure obtained via full first-principles calculations. Finally, we provide a downloadable virtual machine that can be used to reproduce the results of this paper, including all first-principles and atomistic simulations as well as the computational workflows.

preprint2020arXiv

From first- to second-order phase transitions in hybrid improper ferroelectrics through entropy stabilisation

Hybrid improper ferroelectrics (HIFs) have been intensely studied over the last few years to gain understanding of their temperature induced phase transitions and ferroelectric switching pathways. Here we report a switching from first- to second-order phase transition pathway for topical HIFs $Ca_{3-x}Sr_xTi_2O_7$, which is driven by the differing entropies of the phases that we identify as being associated with the dynamic motion of octahedral tilts and rotations. A greater understanding of the transition pathways in this class of layered perovskites, which host many physical properties that are coupled to specific symmetries and octahedral rotation and tilt distortions -- such as superconductivity, negative thermal expansion, fast ion conductivity, ferroelectricity, among others -- is a crucial step in creating novel functional materials by design.

preprint2020arXiv

Hartree theory calculations of quasiparticle properties in twisted bilayer graphene

A detailed understanding of interacting electrons in twisted bilayer graphene (tBLG) near the magic angle is required to gain insights into the physical origin of the observed broken symmetry phases. Here, we present extensive atomistic Hartree theory calculations of the electronic properties of tBLG in the (semi-)metallic phase as function of doping and twist angle. Specifically, we calculate quasiparticle properties, such as the band structure, density of states (DOS) and local density of states (LDOS), which are directly accessible in photoemission and tunnelling spectroscopy experiments. We find that quasiparticle properties change significantly upon doping - an effect which is not captured by tight-binding theory. In particular, we observe that the partially occupied bands flatten significantly which enhances the density of states at the Fermi level. We predict a clear signature of this band flattening in the LDOS in the AB/BA regions of tBLG which can be tested in scanning tunneling experiments. We also study the dependence of quasiparticle properties on the dielectric environment of tBLG and discover that these properties are surprisingly robust as a consequence of the strong internal screening. Finally, we present a simple analytical expression for the Hartree potential which enables the determination of quasiparticle properties without the need for self-consistent calculations.

preprint2020arXiv

The CECAM Electronic Structure Library and the modular software development paradigm

First-principles electronic structure calculations are very widely used thanks to the many successful software packages available. Their traditional coding paradigm is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, from the compiler up, with the exception of linear-algebra and message-passing libraries. This model has been quite successful for decades. The rapid progress in methodology, however, has resulted in an ever increasing complexity of those programs, which implies a growing amount of replication in coding and in the recurrent re-engineering needed to adapt to evolving hardware architecture. The Electronic Structure Library (\esl) was initiated by CECAM (European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure programs and redesign them as free, open-source libraries. They include "heavy-duty" ones with a high degree of parallelisation, and potential for adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by scientists when implementing new ideas. It is a community effort, undertaken by developers of various successful codes, now facing the challenges arising in the new model. This modular paradigm will improve overall coding efficiency and enable specialists (computer scientists or computational scientists) to use their skills more effectively. It will lead to a more sustainable and dynamic evolution of software as well as lower barriers to entry for new developers.

preprint2019arXiv

Critical role of device geometry for the phase diagram of twisted bilayer graphene

The effective interaction between electrons in two-dimensional materials can be modified by their environment, enabling control of electronic correlations and phases. Here, we study the dependence of electronic correlations in twisted bilayer graphene (tBLG) on the separation to the metallic gate(s) in two device configurations. Using an atomistic tight-binding model, we determine the Hubbard parameters of the flat bands as a function of gate separation, taking into account the screening from the metallic gate(s), the dielectric spacer layers and the tBLG itself. We determine the critical gate separation at which the Hubbard parameters become smaller than the critical value required for a transition from a correlated insulator state to a (semi-)metallic phase. We show how this critical gate separation depends on twist angle, doping and the device configuration. These calculations may help rationalise the reported differences between recent measurements of tBLG's phase diagram and suggests that correlated insulator states can be screened out in devices with thin dielectric layers.

preprint2019arXiv

Resonant and Bound States of Charged Defects in Two-Dimensional Semiconductors

A detailed understanding of charged defects in two-dimensional semiconductors is needed for the development of ultrathin electronic devices. Here, we study negatively charged acceptor impurities in monolayer WS$_2$ using a combination of scanning tunnelling spectroscopy and large-scale atomistic electronic structure calculations. We observe several localized defect states of hydrogenic wave function character in the vicinity of the valence band edge. Some of these defect states are bound, while others are resonant. The resonant states result from the multi-valley valence band structure of WS$_2$, whereby localized states originating from the secondary valence band maximum at $Γ$ hybridize with continuum states from the primary valence band maximum at K/K$^{\prime}$. Resonant states have important consequences for electron transport as they can trap mobile carriers for several tens of picoseconds.

preprint2019arXiv

Wannier90 as a community code: new features and applications

Wannier90 is an open-source computer program for calculating maximally-localised Wannier functions (MLWFs) from a set of Bloch states. It is interfaced to many widely used electronic-structure codes thanks to its independence from the basis sets representing these Bloch states. In the past few years the development of Wannier90 has transitioned to a community-driven model; this has resulted in a number of new developments that have been recently released in Wannier90 v3.0. In this article we describe these new functionalities, that include the implementation of new features for wannierisation and disentanglement (symmetry-adapted Wannier functions, selectively-localised Wannier functions, selected columns of the density matrix) and the ability to calculate new properties (shift currents and Berry-curvature dipole, and a new interface to many-body perturbation theory); performance improvements, including parallelisation of the core code; enhancements in functionality (support for spinor-valued Wannier functions, more accurate methods to interpolate quantities in the Brillouin zone); improved usability (improved plotting routines, integration with high-throughput automation frameworks), as well as the implementation of modern software engineering practices (unit testing, continuous integration, and automatic source-code documentation). These new features, capabilities, and code development model aim to further sustain and expand the community uptake and range of applicability, that nowadays spans complex and accurate dielectric, electronic, magnetic, optical, topological and transport properties of materials.