Search

The Raman active G mode in graphene exhibits strong coupling to electrons, yet the comprehensive treatment of this interaction in the calculation of its temperature-dependent Raman spectrum remains incomplete. In this study, we calculate the temperature dependence of the G mode frequency and linewidth, and successfully explain the experimental trend, by accounting for the contributions arising from the first-order electron-phonon coupling, electron-mediated phonon-phonon coupling, and standard lattice anharmonicity. The generality of our approach enables its broad applicability to study phonon dynamics in materials where both electron-phonon coupling and anharmonicity are important., Comment: 14 pages, 7 figures, SI included in "anc" folder

Heat-to-charge conversion efficiency of thermoelectric materials is closely linked to the entropy per charge carrier. Thus, magnetic materials are promising building blocks for highly efficient energy harvesters, as their carrier entropy is boosted by a spin degree of freedom. In this work, we investigate how this spin entropy impacts heat-to-charge conversion in A-type antiferromagnet CrSBr. We perform simultaneous measurements of electrical conductance and thermocurrent while changing magnetic order using temperature and magnetic field as tuning parameters. We find a strong enhancement of the thermoelectric power factor around the N\'eel temperature. We further reveal that the power factor at low temperature can be increased by up to 600% upon applying a magnetic field. Our results demonstrate that the thermoelectric properties of 2D magnets can be optimized by exploiting the sizeable impact of spin entropy and confirm thermoelectric measurements as a sensitive tool to investigate subtle magnetic phase transitions in low-dimensional magnets., Comment: 22 pages, 4 figures

The van der Waals (vdW) type-II multiferroic NiI$_2$ has emerged as a candidate for exploring non-collinear magnetism and magnetoelectric effects in the 2D limit. Frustrated intralayer exchange interactions on a triangular lattice result in a helimagnetic ground state, with spin-induced improper ferroelectricity stabilized by the interlayer interactions. Here we investigate the magnetic and structural phase transitions in bulk NiI$_2$, using high-pressure Raman spectroscopy, optical linear dichroism, and x-ray diffraction. We obtain evidence for a significant pressure enhancement of the antiferromagnetic and helimagnetic transition temperatures, at rates of $\sim15.3/14.4$ K/GPa, respectively. These enhancements are attributed to a cooperative effect of pressure-enhanced interlayer and third-nearest-neighbor intralayer exchange. These results reveal a general path for obtaining high-temperature type-II multiferroicity via high pressures in vdW materials.

Transition metal dihalides have recently garnered interest in the context of two-dimensional van der Waals magnets as their underlying geometrically frustrated triangular lattice leads to interesting competing exchange interactions. In particular, NiI$_{2}$ is a magnetic semiconductor that has been long known for its exotic helimagnetism in the bulk. Recent experiments have shown that the helimagnetic state survives down to the monolayer limit with a layer-dependent magnetic transition temperature that suggests a relevant role of the interlayer coupling. Here, we explore the effects of hydrostatic pressure as a means to enhance this interlayer exchange and ultimately tune the electronic and magnetic response of NiI$_{2}$. We study first the evolution of the structural parameters as a function of external pressure using first-principles calculations combined with x-ray diffraction measurements. We then examine the evolution of the electronic structure and magnetic exchange interactions via first-principles calculations and Monte Carlo simulations. We find that the leading interlayer coupling is an antiferromagnetic second-nearest neighbor interaction that increases monotonically with pressure. The ratio between isotropic third- and first-nearest neighbor intralayer exchanges, which controls the magnetic frustration and determines the magnetic propagation vector $\mathbf{q}$ of the helimagnetic ground state, is also enhanced by pressure. As a consequence, our Monte Carlo simulations show a monotonic increase in the magnetic transition temperature, indicating that pressure is an effective means to tune the magnetic response of NiI$_{2}$.

The dynamical properties of nuclei, carried by the concept of phonon quasiparticles (QP), are central to the field of condensed matter. While the harmonic approximation can reproduce a number of properties observed in real crystals, the inclusion of anharmonicity in lattice dynamics is essential to accurately predict properties such as heat transport or thermal expansion. For highly anharmonic systems, non perturbative approaches are needed, which result in renormalized theories of lattice dynamics. In this article, we apply the Mori-Zwanzig projector formalism to derive an exact generalized Langevin equation describing the quantum dynamics of nuclei in a crystal. By projecting this equation on quasiparticles in reciprocal space, and with results from linear response theory, we obtain a formulation of vibrational spectra that fully accounts for the anharmonicity. Using a mode-coupling approach, we construct a systematic perturbative expansion in which each new order is built to minimize the following ones. With a truncation to the lowest order, we show how to obtain a set of self-consistent equations that can describe the lineshapes of quasiparticles. The only inputs needed for the resulting set of equations are the static Kubo correlation functions, which can be computed using (fully quantum) path-integral molecular dynamics or approximated with (classical or ab initio) molecular dynamics. We illustrate the theory with an application on fcc 4He, an archetypal quantum crystal with very strong anharmonicity.

Magnetic rare-earth orthoferrites $R$FeO$_3$ host a variety of functional properties from multiferroicity and strong magnetostriction, to spin-reorientation transitions and ultrafast light-driven manipulation of magnetism, which can be exploited in spintronics and next-generation devices. Among these systems, SmFeO$_3$ is attracting a particular interest for its rich phase diagram and the high temperature Fe-spin magnetic transitions, which combines with a very low temperature and as yet unclear Sm-spin ordering. Various experiments suggest that the interaction between the Sm and Fe magnetic moments (further supported by the magnetic anisotropy), is at the origin of the complex cascade of transitions, but a conclusive and clear picture has not yet been reached. In this work, by means of comprehensive first-principles calculations, we unravel the role of the magnetic Sm ions in the Fe-spin reorientation transition and in the detected anomalies in the lattice vibrational spectrum, which are a signature of a relevant spin-phonon coupling. By including both Sm-$f$ electrons and non-collinear magnetism, we find frustrated and anisotropic Sm interactions, and a large magnetocrystalline anisotropy mediated by the SOC of the Sm-$4f$ electrons, which drive the complex magnetic properties and phase diagram of SmFeO$_3$.

Diffusion is one of the most ubiquitous transport phenomena in nature. Experimentally, it can be tracked by following point spreading in space and time. Here, we introduce a spatiotemporal pump-probe microscopy technique that exploits the residual spatial temperature profile obtained through the transient reflectivity when probe pulses arrive before pump pulses. This corresponds to an effective pump-probe time delay of 13 ns, determined by the repetition rate of our laser system (76 MHz). This pre-time-zero technique enables probing the diffusion of long-lived excitations created by previous pump pulses with nanometer accuracy, and is particularly powerful for following in-plane heat diffusion in thin films. In contrast to existing techniques for quantifying thermal transport it does not require any material input parameters or strong heating. We demonstrate the direct determination of the thermal diffusivities of the layered materials MoSe$_2$ (0.18 cm$^2$/s), WSe$_2$ (0.20 cm$^2$/s), MoS$_2$ (0.35 cm$^2$/s), and WS$_2$ (0.59 cm$^2$/s). This technique paves the way for observing novel nanoscale thermal transport phenomena and tracking diffusion of a broad range of species.

We evidence the possibility for coherent electrical manipulation of the spin orientation of topologically protected edge states in a low-symmetry quantum spin Hall insulator. By using a combination of ab-initio simulations, symmetry-based modeling, and large-scale calculations of the spin Hall conductivity, it is shown that small electric fields can efficiently vary the spin textures of edge currents in monolayer 1T'-WTe2 by up to a 90-degree spin rotation, without jeopardizing their topological character. These findings suggest a new kind of gate-controllable spin-based device, topologically protected against disorder and of relevance for the development of topological spintronics.

The electronic structure of condensed matter can be significantly affected by the electron-phonon interaction, leading to important phenomena such as electrical resistance, superconductivity or the formation of polarons. This interaction is often neglected in band structure calculations but can have a strong impact on band gaps or optical spectra. Commonly used frameworks for electron-phonon energy corrections are the Allen-Heine-Cardona theory and the Fr\"ohlich model. While the latter shows qualitative agreement with experiment for many polar materials, its simplicity should bring hard limits to its applicability in real materials. Improvements can be made by introducing a generalized version of the model, which considers anisotropic and degenerate electronic bands, and multiple phonon branches. In this work, we search for trends and outliers on over a thousand materials in existing databases of phonon and electron band structures. We use our results to identify the limits of applicability of the standard Fr\"olich model by comparing to the generalized version, and by testing its basic hypothesis of a large radius for the polaronic wavefunction and the corresponding atomic displacement cloud. Among our extended set of materials, most exhibit large polaron behavior as well as validity of the perturbative treatment. For the valence band, there is also a significant fraction of the materials for which the perturbative treatment cannot be applied and/or for which the size of the self-trapping region is close to the atomic repetition distance. We find a large variety of behaviors, and employ much more accurate, fully ab initio Allen-Heine-Cardona calculations to understand extreme cases, where the Fr\"ohlich model should fail and unusually large zero-point renormalization energies occur.

Polarons are quasi-particles made from electrons interacting with vibrations in crystal lattices. They derive their name from the strong electron-vibration polar interaction in ionic systems, that induces associated spectroscopic and optical signatures of such quasi-particles in these materials. In this paper, we focus on diamond, a non-polar crystal with inversion symmetry which nevertheless shows characteristic signatures of polarons, better denoted "nonpolarons" in this case. The polaronic effects are produced by short-range crystal fields with only a small influence of long-range quadrupoles. The many-body spectral function has a characteristic energy dependence, showing a plateau structure that is similar to but distinct from the satellites observed in the polar Fr\"{o}hlich case. The temperature-dependent spectral function of diamond is determined by two methods: the standard Dyson-Migdal approach, which calculates electron-phonon interactions within the lowest-order expansion of the self-energy, and the cumulant expansion, which includes higher orders of electron-phonon interactions. The latter corrects the nonpolaron energies and broadening, providing a more realistic spectral function, which we examine in detail for both conduction and valence band edges., Comment: 10 pages 16 Figures + Subfigures 2 Tables 82 References Supplementary Information (SI) at the end: SI 10 pages SI 14 Figures + Subfigures SI 1 Table

Elastic knock-on is the main electron irradiation damage mechanism in metals including graphene. Atomic vibrations influence its cross-section, but only the out-of-plane direction has been considered so far in the literature. Here, we present a full three-dimensional theory of knock-on damage including the effect of temperature and vibrations to describe ejection into arbitrary directions. We thus establish a general quantitative description of electron irradiation effects through elastic scattering. Applying our methodology to in-plane jumps of pyridinic nitrogen atoms, we show their observed rates imply much stronger inelastic effects than in pristine graphene., Comment: 5 pages, 2 figures, 1 table (+ supplement)

Polarons, that is, charge carriers correlated with lattice deformations, are ubiquitous quasiparticles in semiconductors, and play an important role in electrical conductivity. To date most theoretical studies of so-called large polarons, in which the lattice can be considered as a continuum, have focused on the original Fr\"ohlich model: a simple (non-degenerate) parabolic isotropic electronic band coupled to one dispersionless longitudinal optical phonon branch. The Fr\"ohlich model allows one to understand characteristics such as polaron formation energy, radius, effective mass and mobility. Real cubic materials, instead, have electronic band extrema that are often degenerate or anisotropic and present several phonon modes. In the present work, we address such issues. We keep the continuum hypothesis inherent to the large polaron Fr\"ohlich model, but waive the isotropic and non-degeneracy hypotheses, and also include multiple phonon branches. For polaron effective masses, working at the lowest order of perturbation theory, we provide analytical results for the case of anisotropic electronic energy dispersion, with two distinct effective masses (uniaxial) and numerical simulations for the degenerate 3-band case, typical of III-V and II-VI semiconductor valence bands. We also deal with the strong-coupling limit, using a variational treatment: we propose trial wavefunctions for the above-mentioned cases, providing polaron radii and energies. Then, we evaluate the polaron formation energies, effective masses and localisation lengths using parameters representative of a dozen II-VI, III-V and oxide semiconductors, for both electron and hole polarons...In the non-degenerate case, we compare the perturbative approach with the Feynman path integral approach in characterisizing polarons in the weak coupling limit..., Comment: 18 pages, 10 figures in main text; 11 pages, 1 figure in supplementary material

Interference effects in finite sections of one-dimensional moir\'e crystals are investigated using a Landauer-B\"uttiker formalism within the tight-binding approximation. We explain interlayer transport in double-wall carbon nanotubes and design a predictive model. Wave function interference is visible at the mesoscale: in the strong coupling regime, as a periodic modulation of quantum conductance and emergent localized states; in the localized-insulating regime, as a suppression of interlayer transport, and oscillations of the density of states. These results could be exploited to design quantum electronic devices.

Many promising optoelectronic devices, such as broadband photodetectors, nonlinear frequency converters, and building blocks for data communication systems, exploit photoexcited charge carriers in graphene. For these systems, it is essential to understand, and eventually control, the cooling dynamics of the photoinduced hot-carrier distribution. There is, however, still an active debate on the different mechanisms that contribute to hot-carrier cooling. In particular, the intrinsic cooling mechanism that ultimately limits the cooling dynamics remains an open question. Here, we address this question by studying two technologically relevant systems, consisting of high-quality graphene with a mobility >10,000 cm$^2$V$^{-1}$s$^{-1}$ and environments that do not efficiently take up electronic heat from graphene: WSe$_2$-encapsulated graphene and suspended graphene. We study the cooling dynamics of these two high-quality graphene systems using ultrafast pump-probe spectroscopy at room temperature. Cooling via disorder-assisted acoustic phonon scattering and out-of-plane heat transfer to the environment is relatively inefficient in these systems, predicting a cooling time of tens of picoseconds. However, we observe much faster cooling, on a timescale of a few picoseconds. We attribute this to an intrinsic cooling mechanism, where carriers in the hot-carrier distribution with enough kinetic energy emit optical phonons. During phonon emission, the electronic system continuously re-thermalizes, re-creating carriers with enough energy to emit optical phonons. We develop an analytical model that explains the observed dynamics, where cooling is eventually limited by optical-to-acoustic phonon coupling. These fundamental insights into the intrinsic cooling mechanism of hot carriers in graphene will play a key role in guiding the development of graphene-based optoelectronic devices.

Interest in layered two dimensional materials, particularly stacked heterostructures of transition metal dichalcogenides, has led to the need for a better understanding of the structural and electronic changes induced by stacking. Here, we investigate the effects of idealized heterostructuring, with periodic commensurate stacking, on the structural, electronic, and vibrational properties, when compared to the counterpart bulk transition metal dichalcogenide. We find that in heterostructures with dissimilar chalcogen species there is a strong compression of the inter-layer spacing, compared to the bulk compounds. This compression of the heterostructure is caused by an increase in the strength of the induced polarization interaction between the layers, but not a full charge transfer. We argue that this effect is real, not due to the imposed commensurability, and should be observable in heterostructures combining different chalcogens. Interestingly, we find that incommensurate stacking of Ti-based dichalcogenides leads to the stabilization of the charge density wave phonon mode, which is unstable in the 1T phase at low temperature. Mixed Ti- and Zr-heterostructures are still unstable, but with a charge density wave or a ferroelectric instability., Comment: 10 pages, 3 tables, 4 figures

We present TB2J, a Python package for the automatic computation of magnetic interactions, including exchange and Dzyaloshinskii-Moriya interactions, between atoms of magnetic crystals from the results of density functional calculations. The program is based on the Green's function method with the local rigid spin rotation treated as a perturbation. As input,the package uses the output of either Wannier90, which is interfaced with many density functional theory packages,or of codes based on localized orbitals. A minimal user input is needed, which allows for easy integration into high-throughput workflows. The package is open source under BSD 2-Clause license, available at https://github.com/mailhexu/TB2J.

We include the treatment of quadrupolar fields beyond the Fr\"ohlich interaction in the first-principles electron-phonon vertex in semiconductors. Such quadrupolar fields induce long-range interactions that have to be taken into account for accurate physical results. We apply our formalism to Si (nonpolar), GaAs, and GaP (polar) and demonstrate that electron mobilities show large errors if dynamical quadrupoles are not properly treated.

We describe a new approach to compute the electron-phonon self-energy and carrier mobilities in semiconductors. Our implementation does not require a localized basis set to interpolate the electron-phonon matrix elements, with the advantage that computations can be easily automated. Scattering potentials are interpolated on dense $\mathbf{q}$ meshes using Fourier transforms and ab initio models to describe the long-range potentials generated by dipoles and quadrupoles. To reduce significantly the computational cost, we take advantage of crystal symmetries and employ the linear tetrahedron method and double-grid integration schemes, in conjunction with filtering techniques in the Brillouin zone. We report results for the electron mobility in Si, GaAs, and GaP obtained with this new methodology.

Thermoelectricity is a promising avenue for harvesting energy but large-scale applications are still hampered by the lack of highly-efficient low-cost materials. Recently, Fe$_2YZ$ Heusler compounds were predicted theoretically to be interesting candidates with large thermoelectric power factor. Here, we show that under doping conditions compatible with thermoelectric applications, these materials are prone to an unexpected magnetic instability detrimental to their thermoelectric performance. We rationalize the physics at the origin of this instability, provide guidelines for avoiding it and discuss its impact on the thermoelectric power factor. Doing so, we also point out the shortcomings of the rigid band approximation commonly used in high-throughput theoretical searches of new thermoelectrics., Comment: 16 pages, 20 figures; the Supplementary Material is included

Stable structures of layered SnS and SnSe and their associated electronic and vibrational spectra are predicted using first-principles DFT calculations. The calculations show that both materials undergo a phase transformation upon thinning whereby the in-plane lattice parameters converge to a pseudo-cubic phase, similar to the high-temperature behaviour observed for their bulk counterparts. The electronic properties of layered SnS and SnSe evolve to an almost symmetric dispersion whilst the gap changes from indirect to direct. Characteristic signatures in the phonon dispersion curves and surface phonon states where only atoms belonging to surface layers vibrate can also be observed for these materials.

The high pressure lattice dynamics of rhombohedral antimony have been studied by a combination of diffuse scattering and inelastic x-ray scattering. The evolution of the phonon behavior as function of pressure was analyzed by means of two theoretical approaches: density functional perturbation theory and symmetry-based phenomenological phase transition analysis. This paper focuses on the first structural phase transition, SbI-SbIV, and the role of vibrations in leading the transition. The phonon dispersion exhibits complex behaviour as one approaches the structural transition, with the branches, corresponding to the two transitions happening at high pressure in the Va elements (A7-to-BCC and A7-to-PC) both showing softening., Comment: 8 pages and 8 figures

First-principles studies of two-dimensional transition metal dichalcogenides have contributed considerably to the understanding of their dielectric, optical, elastic, and vibrational properties. The majority of works to date focus on a single material or physical property. Here we use a single first-principles methodology on the whole family of systems, to investigate in depth the relationships between different physical properties, the underlying symmetry and the composition of these materials, and observe trends. We compare to bulk counterparts to show strong interlayer effects in triclinic compounds. Previously unobserved relationships between these monolayer compounds become apparent. These trends can then be exploited by the materials science, nanoscience, and chemistry communities to better design devices and heterostructures for specific functionalities., Comment: 4 figures, 11 pages

London-van der Waals dispersion forces are a fundamental component of condensed matter systems, biological processes, and self-assembly. In this letter we propose a method to calculate the C6 coefficients that characterize dispersion forces in the non-retarded regime for molecules in a coherent superposition of excited states. Several ultrafast femtosecond pump-probe schemes are investigated. We apply the method to LiH molecules and show that their London-van der Waals interaction can change dramatically after the interaction with the pump pulse. The pulse modulates the C6 coefficients, and the interplay between polarization, orientation of the molecules, and the dipole fields gives rise to a rich variety of combinations.

We employ an ab-initio structure search algorithm to explore the configurational space of Bi in quasi two dimensions. A confinement potential restricts the movement of atoms within a pre-defined thickness during structure search calculations within the minima hopping method to find the stable and metastable forms of bilayer Bi. In addition to recovering the two known low-energy structures (puckered monoclinic and buckled hexagonal), our calculations predict three new structures of bilayer Bi. We call these structures the $\alpha$, $\beta$, and $\gamma$ phases of bilayer Bi, which are, respectively, 63, 72, and 83 meV/atom higher in energy than that of the monoclinic ground state, and thus potentially synthesizable using appropriate substrates. We also compare the structural, electronic, and vibrational properties of the different phases. The puckered monoclinic, buckled hexagonal, and $\beta$ phases exhibit a semiconducting energy gap, whereas $\alpha$ and $\gamma$ phases are metallic. We notice an unusual Mexican-hat type band dispersion leading to a van Hove singularity in the buckled hexagonal bilayer Bi. Notably, we find symmetry-protected topological Dirac points in the electronic spectrum of the $\gamma$ phase. The new structures suggest that bilayer Bi provides a novel playground to study distortion-mediated metal-insulator phase transitions.

We extend the density functional perturbation theory formalism to the case of non-collinear magnetism. The main problem comes with the exchange-correlation (XC) potential derivatives, which are the only ones that are affected by the non-collinearity of the system. Most of the present XC functionals are constructed at the collinear level, such that the off-diagonal (containing magnetization densities along $x$ and $y$ directions) derivatives cannot be calculated simply in the non-collinear framework. To solve this problem, we consider here possibilities to transform the non-collinear XC derivatives to a local collinear basis, where the $z$ axis is aligned with the local magnetization at each point. The two methods we explore are i) expanding the spin rotation matrix as a Taylor series, ii) evaluating explicitly the XC for the local density approximation through an analytical expression of the expansion terms. We compare the two methods and describe their practical implementation. We show their application for atomic displacement and electric field perturbations at the second order, within the norm-conserving pseudopotential methods.

Interest in the bulk transition metal dichalcogenides for their electronic, photovoltaic, and optical properties has grown and led to their use in many technological applications. We present a systematic investigation of their interlinked vibrational and dielectric properties, using density functional theory and density functional perturbation theory, studying the effects of the spin-orbit interaction and of the long-range e$^-$- e$^-$ correlation as part of our investigation. This study confirms that the spin-orbit interaction plays a small role in these physical properties, while the direct contribution of dispersion corrections is of crucial importance in the description of the interatomic force constants. Here, our analysis of the structural and vibrational properties, including the Raman spectra, compare well to experimental measurement. Three materials with different point groups are showcased and data trends on the full set of fifteen existing hexagonal, trigonal, and triclinic materials are demonstrated. This overall picture will enable the modeling of devices composed of these materials for novel applications., Comment: 11 pages, 6 figures

The direct manipulation of individual atoms in materials using scanning probe microscopy has been a seminal achievement of nanotechnology. Recent advances in imaging resolution and sample stability have made scanning transmission electron microscopy a promising alternative for single-atom manipulation of covalently bound materials. Pioneering experiments using an atomically focused electron beam have demonstrated the directed movement of silicon atoms over a handful of sites within the graphene lattice. Here, we achieve a much greater degree of control, allowing us to precisely move silicon impurities along an extended path, circulating a single hexagon, or back and forth between the two graphene sublattices. Even with manual operation, our manipulation rate is already comparable to the state-of-the-art in any atomically precise technique. We further explore the influence of electron energy on the manipulation rate, supported by improved theoretical modeling taking into account the vibrations of atoms near the impurities, and implement feedback to detect manipulation events in real time. In addition to atomic-level engineering of its structure and properties, graphene also provides an excellent platform for refining the accuracy of quantitative models and for the development of automated manipulation., Comment: 5 figures, 4 supporting figures

BoltzTraP2 is a software package for calculating a smoothed Fourier expression of periodic functions and the Onsager transport coefficients for extended systems using the linearized Boltzmann transport equation. It uses only the band and $k$-dependent quasi-particle energies, as well as the intra-band optical matrix elements and scattering rates, as input. The code can be used via a command-line interface and/or as a Python module. It is tested and illustrated on a simple parabolic band example as well as silicon. The positive Seebeck coefficient of lithium is reproduced in an example of going beyond the constant relaxation time approximation.

We report on the structural, electronic and magnetic properties of a monolayer of GdAg2, forming a moir\'e pattern on Ag(111). Combining scanning tunneling microscopy and ab-initio spin-polarized calculations, we show that the electronic band structure can be shifted linearly via thermal dependent strain of the intra-layer atomic distance in a range between 1-7%, leading to lateral hetero-structuring. Furthermore, the coupling of the incommensurable GdAg2 alloy layer to the Ag(111) substrate leads to spatially varying atomic relaxation causing subsurface layer buckling, texturing of the electronic and spin properties, and inhomogeneity of the magnetic anisotropy energy across the layer. These results provide perspectives for a control of electronic properties and magnetic ordering in atomically-thin layers.

We present time-resolved Kerr rotation measurements, showing spin lifetimes of over 100 ns at room temperature in monolayer MoSe$_2$. These long lifetimes are accompanied by an intriguing temperature dependence of the Kerr amplitude, which increases with temperature up to 50 K and then abruptly switches sign. Using ab initio simulations we explain the latter behavior in terms of the intrinsic electron-phonon coupling and the activation of transitions to secondary valleys. The phonon-assisted scattering of the photo-excited electron-hole pairs prepares a valley spin polarization within the first few ps after laser excitation. The sign of the total valley magnetization, and thus the Kerr amplitude, switches as a function of temperature, as conduction and valence band states exhibit different phonon-mediated inter-valley scattering rates. However, the electron-phonon scattering on the ps time scale does not provide an explanation for the long spin lifetimes. Hence, we deduce that the initial spin polarization must be transferred into spin states which are protected from the intrinsic electron-phonon coupling, and are most likely resident charge carriers which are not part of the itinerant valence or conduction band states., Comment: 18 pages, 17 figures

We investigate the chemical bonding characteristics of the transition metal dichalcogenides based on their static and dynamical atomic charges within Density Functional Theory. The dynamical charges of the trigonal transition metal dichalcogenides are anomalously large, while in their hexagonal counterparts, their sign is even counterintuitive i.e. the transition metal takes the negative charge. This phenomenon cannot be understood simply in terms of a change in the static atomic charge as it results from a local change of polarization. We present our theoretical understanding of these phenomena based on the perturbative response of the system to a static electric field and by investigating the hybridization of the molecular orbitals near the Fermi level. Furthermore, we establish a link between the sign of the Born effective charge and the $\pi$-backbonding in organic chemistry and propose an experimental procedure to verify the calculated sign of the dynamical charge in the transition metal dichalcogenides., Comment: 3 figures, 6 pages

Cubic SrTiO$_{3}$ becomes tetragonal below 105 K. The antiferrodistortive (AFD) distortion leads to clockwise and counter-clockwise rotation of adjacent TiO$_{6}$ octahedra. This insulator becomes a metal upon the introduction of extremely low concentration of n-type dopants. However, signatures of the structural phase transition in charge conduction have remained elusive. Employing the Montgomery technique, we succeed in resolving the anisotropy of charge conductivity induced by the AFD transition, in the presence of different types of dopants. We find that the slight lattice distortion ($<6 \times 10^{-4}$) gives rise to a twenty percent anisotropy in charge conductivity, in agreement with the expectations of band calculations. Application of uniaxial strain amplifies the detectable anisotropy by disfavoring one of the three possible tetragonal domains. In contrast with all other known anisotropic Fermi liquids, the anisotropy has opposite signs for elastic and inelastic scattering. Increasing the concentration of dopants leads to a drastic shift in the temperature of the AFD transition either upward or downward. The latter result puts strong constraints on any hypothetical role played by the AFD soft mode in the formation of Cooper pairs and the emergence of superconductivity in SrTiO$_3$., Comment: 6 pages with 5 figures

We present a theoretical investigation of the anomalous ferroelectricity of mixed-stack charge transfer molecular crystals, based on the Peierls-Hubbard model, and first principles calculations for its parameterization. This approach is first validated by reproducing the temperature-induced transition and the electronic polarization of TTF-CA, and then applied to a novel series of hydrogen-bonded crystals, for which room temperature ferroelectricity has recently been claimed. Our analysis shows that the hydrogen-bonded systems present a very low degree of charge transfer and hence support a very small polarization. A critical re-examination of experimental data supports our findings, shedding doubts on the ferroelectricity of these systems. More generally, our modelling allows the rationalization of general features of the ferroelectric transition in charge transfer crystals, and suggests design principles for materials optimization.

Lithium is one of the simplest metals, with negative charge carriers and a close reproduction of free electron dispersion. Experimentally, however, Li is one of a handful of elemental solids (along with Cu, Ag, Au etc.) where the sign of the Seebeck coefficient (S) is opposite to that of the carrier. This counterintuitive behavior still lacks a satisfactory interpretation. We calculate S fully from first-principles, within the framework of P.B. Allen's formulation of Boltzmann transport theory. Here it is crucial to avoid the constant relaxation time approximation, which gives a sign for S which is necessarily that of the carriers. Our calculated S are in excellent agreement with experimental data, up to the melting point. In comparison with another alkali metal Na, we demonstrate that within the simplest non-trivial model for the energy dependency of the electron lifetimes, the rapidly increasing density of states (DOS) across the Fermi energy is related to the sign of S in Li. The exceptional energy dependence of the DOS is beyond the free-electron model, as the dispersion is distorted by the Brillouin Zone edge, a stronger effect in Li than other Aliki metals. The electron lifetime dependency on energy is central, but the details of the electron-phonon interaction are found to be less important, contrary to what has been believed for several decades. Band engineering combined with the mechanism exposed here may open the door to new "ambipolar" thermoelectric materials, with a tunable sign for the thermopower even if either n- or p-type doping is impossible., Comment: 5 pages 4 figures