Researcher profile

David Vanderbilt

David Vanderbilt contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

19 published item(s)

preprint2022arXiv

Direct visualization of surface spin-flip transition in MnBi4Te7

We report direct visualization of spin-flip transition of the surface layer in antiferromagnet MnBi4Te7, a natural superlattice of alternating MnBi2Te4 and Bi2Te3 layers, using cryogenic magnetic force microscopy (MFM). The observation of magnetic contrast across domain walls and step edges confirms that the antiferromagnetic order persists to the surface layers. The magnetic field dependence of the MFM images reveals that the surface magnetic layer undergoes a first-order spin-flip transition at a magnetic field that is lower than the bulk transition, in excellent agreement with a revised Mills' model. Our analysis indicates an enhancement of the order parameter in the surface magnetic layer, implying robust ferromagnetism in the single-layer limit. The direct visualization of surface spin-flip transition not only opens up exploration of surface metamagnetic transitions in layered antiferromagnets, but also provides experimental support for realizing quantized transport in ultra-thin films of MnBi4Te7 and other natural superlattice topological magnets.

preprint2022arXiv

Electronic correlation in nearly free electron metals with beyond-DFT methods

For more than three decades, nearly free electron elemental metals have been a topic of debate because the computed bandwidths are significantly wider in the local density approximation to density-functional theory (DFT) than indicated by angle-resolved photoemission (ARPES) experiments. Here, we systematically investigate this using first-principles calculations for alkali and alkaline-earth metals using DFT and various beyond-DFT methods such as meta-GGA, G$_0$W$_0$, hybrid functionals (YS-PBE0, B3LYP), and LDA+eDMFT. We find that the static non-local exchange, as partly included in the hybrid functionals, significantly increase the bandwidths even compared to LDA, while the G$_0$W$_0$ bands are only slightly narrower than in LDA. The agreement with the ARPES is best when the local approximation to the self-energy is used in the LDA+eDMFT method. We infer that even moderately correlated systems with partially occupied s-orbitals, which were assumed to approximate the uniform electron gas, are very well described in terms of short-range dynamical correlations that are only local to an atom.

preprint2022arXiv

Engineering magnetic topological insulators in Eu$_5M_2X_6$ Zintls

Magnetic topological insulator provide a prominent material platform for quantum anomalous Hall physics and axion electrodynamics. However, the lack of material realizations with cleanly gapped surfaces hinders technological utilization of these exotic quantum phenomena. Here, using the Zintl concept and the properties of non-symmorphic space groups, we computationally engineer magnetic topological insulators. Specifically, we explore Eu$_5M_2X_6$ ($M$=metal, $X$=pnictide) Zintl compounds and find that Eu$_5$Ga$_2$Sb$_6$, Eu$_5$Tl$_2$Sb$_6$ and Eu$_5$In$_2$Bi$_6$ form stable structures with non-trivial $\mathbb{Z}_2$ indices. We also show that epitaxial and uniaxial strain can be used to control the $\mathbb{Z}_2$ index and the bulk energy gap. Finally, we discuss experimental progress towards the synthesis of the proposed candidates and provide insights that can be used in the search for robust magnetic topological insulators in Zintl compounds.

preprint2022arXiv

Vibrational fingerprints of ferroelectric hafnia

Hafnia (HfO2) is a promising material for emerging chip applications due to its high-k dielectric behaviour, suitability for negative capacitance heterostructures, scalable ferroelectricity, and silicon compatibility. The lattice dynamics along with phononic properties such as thermal conductivity, contraction, and heat capacity are under-explored, primarily due to the absence of high quality single crystals. Herein, we report the vibrational properties of a series of HfO2 crystals stabilized with yttrium (chemical formula HfO2:xY, where x = 20, 12, 11, 8, and 0%) and compare our findings with a symmetry analysis and lattice dynamics calculations. We untangle the effects of Y by testing our calculations against the measured Raman and infrared spectra of the cubic, antipolar orthorhombic, and monoclinic phases and then proceed to reveal the signature modes of polar orthorhombic hafnia. This work provides a spectroscopic fingerprint for several different phases of HfO2 and paves the way for an analysis of mode contributions to high-k dielectric and ferroelectric properties for chip technologies.

preprint2022arXiv

Vibrational properties of CuInP2S6 across the ferroelectric transition

In order to explore the properties of a two-sublattice ferroelectric, we measured the infrared and Raman scattering response of CuInP2S6 across the ferroelectric and glassy transitions and compared our findings to a symmetry analysis, calculations of phase stability, and lattice dynamics. In addition to uncovering displacive character and a large hysteresis region surrounding the ferroelectric transition temperature T_C, we identify the vibrational modes that stabilize the polar phase and confirm the presence of two ferroelectric variants with opposite polarizations. Below TC, a poorly understood relaxational or glassy transition at Tg is characterized by local structure changes in the form of subtle peak shifting and activation of low frequency out-of-plane Cu- and In-containing modes. The latter are due to changes in the Cu/In coordination environments and associated order-disorder processes. Moreover, Tg takes place in two steps with another large hysteresis region and significant underlying scattering. Combined with imaging of the room temperature phase separation, this effort lays the groundwork for studying CuInP2S6 under external stimuli and in the ultra-thin limit.

preprint2021arXiv

Humble planar defects in SiGe nanopillars

We report a new \{001\} planar defect found in SiGe nanopillars. The defect structure, determined by atomic resolution electron microscopy, matches the Humble defect model proposed for diamond. We also investigated several possible variants of the Humble structure using first principles calculations and found that the one lowest in energy was also in best agreement with the STEM images. The pillar composition has been analyzed with electron energy loss spectroscopy, which hints at how the defect is formed. Our results show that the structure and formation process of defects in nanostructured group IV semiconductors can be different from their bulk counterparts.

preprint2021arXiv

Lattice dynamics and structural transition of the hyperhoneycomb iridate $β$-Li$_2$IrO$_3$ investigated by high-pressure Raman scattering

We report a polarized Raman scattering study of the lattice dynamics of $β$-Li$_2$IrO$_3$ under hydrostatic pressures up to 7.62 GPa. At ambient pressure, $β$-Li$_2$IrO$_3$ exhibits the hyperhoneycomb crystal structure and a magnetically ordered state of spin-orbit entangled Jeff = 1/2 moments that is strongly influenced by bond-directional (Kitaev) exchange interactions. At a critical pressure of ~ 4.1 GPa, the phonon spectrum changes abruptly consistent with the reported structural transition into a monoclinic, dimerized phase. A comparison to the phonon spectra obtained from density functional calculations shows reasonable overall agreement. The calculations also indicate that the high-pressure phase is a nonmagnetic insulator driven by the formation of Ir-Ir dimer bonds. Our results thus indicate a strong sensitivity of the electronic properties of $β$-Li$_2$IrO$_3$ to the pressure-induced structural transition.

preprint2021arXiv

Orbital-selective Mott phase and non-Fermi liquid in FePS$_3$

The layered metal phosphorous trisulfide FePS$_3$ is reported to be a Mott insulator at ambient conditions and to undergo structural and insulator-metal phase transitions under pressure. However, the character of the resulting metallic states has not been understood clearly so far. Here, we theoretically study the phase transitions of FePS$_3$ using first-principles methods based on density functional theory and embedded dynamical mean field theory. We find that the Mott transition in FePS$_3$ can be orbital-selective, with $t_{2g}$ states undergoing a correlation-induced insulator-to-metal transition while $e_g$ states remain gapped. We show that this orbital-selective Mott phase, which occurs only when non-hydrostatic pressure is used, is a bad metal (or non-Fermi liquid) with large fluctuating moments due to Hund's coupling. Further application of pressure increases the crystal-field splitting and converts the system to a conventional Fermi liquid with low-spin configurations dominant. Our results show that FePS$_3$ is a novel example of a system that realizes an orbital-selective Mott phase, allowing tuning between correlated and uncorrelated metallic properties in an accessible pressure range ($\leq$ 18 GPa).

preprint2021arXiv

Phonon spectrum of Pr$_2$Zr$_2$O$_7$ and Pr$_2$Ir$_2$O$_7$ as an evidence of coupling of the lattice with electronic and magnetic degrees of freedom

Magnetic materials with pyrochlore crystal structure form exotic magnetic states due to the high lattice frustration. In this work we follow the effects of coupling of the lattice and electronic and magnetic degrees of freedom in two Praseodymium-based pyrochlores Pr$_2$Zr$_2$O$_7$ and Pr$_2$Ir$_2$O$_7$. In both materials the presence of magnetic interactions does not lead to magnetically ordered low temperature states, however their electronic properties are different. A comparison of Raman phonon spectra of Pr$_2$Zr$_2$O$_7$ and Pr$_2$Ir$_2$O$_7$ allows us to identify magneto-elastic coupling in Pr$_2$Zr$_2$O$_7$ that elucidates its magnetic properties at intermediate temperatures, and allows us to characterize phonon-electron coupling in the semimetallic Pr$_2$Ir$_2$O$_7$. We also show that the effects of random disorder on the Raman phonon spectra is small.

preprint2021arXiv

Quadrupole moments, edge polarizations, and corner charges in the Wannier representation

The modern theory of polarization allows for the determination of the macroscopic end charge of a truncated one-dimensional insulator, modulo the charge quantum $e$, from a knowledge of bulk properties alone. A more subtle problem is the determination of the corner charge of a two-dimensional insulator, modulo $e$, from a knowledge of bulk and edge properties alone. While previous works have tended to focus on the quantization of corner charge in the presence of symmetries, here we focus on the case that the only bulk symmetry is inversion, so that the corner charge can take arbitrary values. We develop a Wannier-based formalism that allows the corner charge to be predicted, modulo $e$, only from calculations on ribbon geometries of two different orientations. We elucidate the dependence of the interior quadrupole and edge dipole contributions upon the gauge used to construct the Wannier functions, finding that while these are individually gauge-dependent, their sum is gauge-independent. From this we conclude that the edge polarization is not by itself a physical observable, and that any Wannier-based method for computing the corner charge requires the use of a common gauge throughout the calculation. We satisfy this constraint using two Wannier construction procedures, one based on projection and another based on a gauge-consistent nested Wannier construction. We validate our theory by demonstrating the correct prediction of corner charge for several tight-binding models. We comment on the relations between our approach and previous ones that have appeared in the literature.

preprint2020arXiv

Axion coupling in the hybrid Wannier representation

Many magnetic point-group symmetries induce a topological classification on crystalline insulators, dividing them into those that have a nonzero quantized Chern-Simons magnetoelectric coupling ("axion-odd" or "topological"), and those that do not ("axion-even" or "trivial"). For time-reversal or inversion symmetries, the resulting topological state is usually denoted as a "strong topological insulator" or an "axion insulator" respectively, but many other symmetries can also protect this "axion $Z_2$" index. Topological states are often insightfully characterized by choosing one crystallographic direction of interest, and inspecting the hybrid Wannier (or equivalently, the non-Abelian Wilson-loop) band structure, considered as a function of the two-dimensional Brillouin zone in the orthogonal directions. Here, we systematically classify the axion-quantizing symmetries, and explore the implications of such symmetries on the Wannier band structure. Conversely, we clarify the conditions under which the axion$Z_2$ index can be deduced from the Wannier band structure. In particular, we identify cases in which a counting of Dirac touchings between Wannier bands, or a calculation of the Chern number of certain Wannier bands, or both, allows for a direct determination of the axion $Z_2$ index. We also discuss when such symmetries impose a "flow" on the Wannier bands, such that they are always glued to higher and lower bands by degeneracies somewhere in the projected Brillouin zone, and the related question of when the corresponding surfaces can remain gapped, thus exhibiting a half-quantized surface anomalous Hall conductivity. Our formal arguments are confirmed and illustrated in the context of tight-binding models for several paradigmatic axion-odd symmetries including time reversal, inversion, simple mirror, and glide mirror symmetries.

preprint2020arXiv

Berry Flux Diagonalization: Application to Electric Polarization

The switching polarization of a ferroelectric is a characterization of the current that flows due to changes in polarization when the system is switched between two states. Computation of this change in polarization in crystal systems has been enabled by the modern theory of polarization, where it is expressed in terms of a change in Berry phase as the material switches. It is straightforward to compute this change of phase, but only modulo $2π$, requiring a branch choice from among a lattice of values separated by $2π$. The measured switching polarization depends on the actual path along which the material switches, which in general involves nucleation and growth of domains and is therefore quite complex. In this work, we present a physically motivated approach for predicting the experimentally measured switching polarization that involves separating the change in phase between two states into as many gauge-invariant smaller phase changes as possible. As long as the magnitudes of these smaller phase changes remain smaller than $π$, their sum forms a phase change which corresponds to the change one would find along any path involving minimal evolution of the atomic and electronic structure. We show that for typical ferroelectrics, including those that would have otherwise required a densely sampled path, this technique allows the switching polarization to be computed without any need for intermediate sampling between oppositely polarized states.

preprint2020arXiv

Engineering Weyl phases and nonlinear Hall effects in T$_d$-MoTe$_2$

MoTe$_2$ has recently attracted much attention due to the observation of pressure-induced superconductivity, exotic topological phase transitions, and nonlinear quantum effects. However, there has been debate on the intriguing structural phase transitions among various observed phases of MoTe$_2$, and their connection to the underlying topological electronic properties. In this work, by means of density-functional theory (DFT+U) calculations, we investigate the structural phase transition between the polar T$_d$ and nonpolar 1T$'$ phases of MoTe$_2$ in reference to a hypothetical high-symmetry T$_0$ phase that exhibits higher-order topological features. In the T$_d$ phase we obtain a total of 12 Weyl points, which can be created/annihilated, dynamically manipulated, and switched by tuning a polar phonon mode. We also report the existence of a tunable nonlinear Hall effect in T$_d$-MoTe$_2$, and propose the use of this effect as a probe for the detection of polarity orientation in polar (semi)metals. By studying the role of dimensionality, we identify a configuration in which a nonlinear surface response current emerges. The potential technological applications of the tunable Weyl phase and the nonlinear Hall effect are discussed.

preprint2020arXiv

Gapless hinge states from adiabatic pumping of axion coupling

We demonstrate that chiral hinge modes naturally emerge in insulating crystals undergoing a slow cyclic evolution that changes the Chern-Simons axion angle $θ$ by $2π$. This happens when the surface (not just the bulk) returns to its initial state at the end of the cycle, in which case it must pass through a metallic state to dispose of the excess quantum of surface anomalous Hall conductivity pumped from the bulk. If two adjacent surfaces become metallic at different points along the cycle, there is an interval in which they are in topologically distinct insulating states, with chiral modes propagating along the connecting hinge. We illustrate these ideas for a tight-binding model consisting of coupled layers of the Haldane model with alternating parameters. The surface topology is determined in a slab geometry using two different markers, surface anomalous Hall conductivity and surface-localized charge pumping (flow of surface-localized Wannier bands), and we find that both correctly predict the appearance of gapless hinge modes in a rod geometry. When viewing the axion pump as a four-dimensional (4D) crystal with one synthetic dimension, the hinge modes trace Fermi arcs in the Brillouin zone of the 2D hinge connecting a pair of 3D surfaces of the 4D crystal.

preprint2020arXiv

Magnetic phase transitions and spin density distribution in the molecular multiferroic GaV$_4$S$_8$ system

We have carried out neutron diffraction and small angle neutron scattering measurements on a high quality single crystal of the cubic lacunar spinel multiferroic, GaV$_4$S$_8$, as a function of magnetic field and temperature to determine the magnetic properties for the single electron that is located on the tetrahedrally coordinated V$_4$ molecular unit. Our results are in good agreement with the structural transition at 44 K from cubic to rhombohedral symmetry where the system becomes a robust ferroelectric, while long range magnetic order develops below 13 K in the form of an incommensurate cycloidal magnetic structure, which can transform into a Néel-type skyrmion phase in a modest applied magnetic field. Below 5.9(3) K, the crystal enters a ferromagnetic phase, and we find the magnetic order parameter indicates a long range ordered ground state with an ordered moment of 0.23(1) $μ_\mathrm{B}$ per V ion. Both polarized and unpolarized neutron data in the ferroelectric-paramagnetic phase have been measured to determine the magnetic form factor. The data are consistent with a model of the single spin being uniformly distributed across the V$_4$ molecular unit, rather than residing on the single apical V ion, in substantial agreement with the results of first-principles theory. In the magnetically ordered state, polarized neutron measurements are important since both the cycloidal and ferromagnetic order parameters are clearly coupled to the ferroelectricity, causing the structural peaks to be temperature and field dependent. For the ferromagnetic ground state, the spins are locked along the $[1,1,1]$ direction by a surprisingly large anisotropy.

preprint2020arXiv

Molecular Mott state in the deficient spinel GaV$_4$S$_8$

In this study, we investigated theoretically the Mott-insulating phase of a deficient spinel chalcogenide GaV$_4$S$_8$, which is known to form a tetrahedral V$_4$S$_4$ cluster unit that results in molecular orbitals (MOs) with a narrow bandwidth in the noninteracting limit. We used a cluster extension of charge self-consistent embedded dynamical mean-field theory to study the impact of strong intra-cluster correlations on the spectral properties as well as the structural degrees of freedom of the system. We found that the strong tetrahedral clustering renders the atomic Mott picture ineffective, and that the resulting MO picture is essential to describe the Mott phase. It was also found that, while the spectral properties can be qualitatively described by the truncation of the Hilbert space down to the lowest-energy MO, a proper description of the structural degrees of freedom requires the inclusion of multi-MO correlations that span a larger energy window. Specifically, we found that the lowest-energy MO description overemphasizes the clustering tendency, while the inclusion of the Hund's coupling between the lower- and higher-energy MOs corrects this tendency, bringing the theoretically predicted crystal structure into good agreement with the experiment.

preprint2020arXiv

Robust A-type order and spin-flop transition on the surface of the antiferromagnetic topological insulator MnBi$_2$Te$_4$

Here we present microscopic evidence of the persistence of uniaxial A-type antiferromagnetic order to the surface layers of MnBi$_2$Te$_4$ single crystals using magnetic force microscopy. Our results reveal termination-dependent magnetic contrast across both surface step edges and domain walls, which can be screened by thin layers of soft magnetism. The robust surface A-type order is further corroborated by the observation of termination-dependent surface spin-flop transitions, which have been theoretically proposed decades ago. Our results not only provide key ingredients for understanding the electronic properties of the antiferromagnetic topological insulator MnBi$_2$Te$_4$, but also open a new paradigm for exploring intrinsic surface metamagnetic transitions in natural antiferromagnets.

preprint2019arXiv

Systematic beyond-DFT study of binary transition metal oxides

Various methods going beyond density-functional theory (DFT), such as DFT+U, hybrid functionals, meta-GGAs, GW, and DFT-embedded dynamical mean field theory (eDMFT), have been developed to describe the electronic structure of correlated materials, but it is unclear how accurate these methods can be expected to be when applied to a given strongly correlated solid. It is thus of pressing interest to compare their accuracy as they apply to different categories of materials. Here, we introduce a novel paradigm in which a chosen set of beyond-DFT methods is systematically and uniformly tested on a chosen class of materials. For a first application, we choose the target materials to be the binary transition-metal oxides FeO, CoO, MnO, and NiO in their antiferromagnetic phase and present a head-to-head comparison of spectral properties as computed using the various methods. We also compare with available experimental angle-resolved photoemission spectroscopy (ARPES), inverse-photoemission spectroscopy, and with optical absorption. We find both B3LYP and eDMFT can reproduce the experiments quite well, with eDMFT doing best in particular when comparing with the ARPES data.

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.