Source author record

Bartomeu Monserrat

Bartomeu Monserrat appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

32 published item(s)

preprint2025arXiv

Impact of electronic correlations on the superconductivity of high-pressure CeH$_9$

Rare-earth superhydrides have attracted considerable attention because of their high critical superconducting temperature under extreme pressures. They are known to have localized valence electrons, implying strong electronic correlations. However, such many-body effects are rarely included in first-principles studies of rare-earth superhydrides because of the complexity of their high-pressure phases. In this work, we use a combined density functional theory and dynamical mean-field theory approach to study both electrons and phonons in the prototypical rare-earth superhydride CeH$_9$, shedding light on the impact of electronic correlations on its critical temperature for phonon-mediated superconductivity. Our findings indicate that electronic correlations result in a larger electronic density at the Fermi level, a bigger superconducting gap, and softer vibrational modes associated with hydrogen atoms. Together, the inclusion of these correlation signatures increases the Migdal-Eliashberg superconducting critical temperature from 47 K to 96 K, close to the measured 95 K. Our results reconcile experimental observations and theoretical predictions for CeH$_9$ and herald a path towards the quantitative modeling of phonon-mediated superconductivity for interacting electron systems.

preprint2023arXiv

Temperature effects in topological insulators of transition metal dichalcogenide monolayers

We investigate the role of temperature on the topological insulating state of metal dichalcogenide monolayers, 1T$^\prime$-MX$_2$ (M=W, Mo and X=S, Se). Using first principles calculations based on density functional theory, we consider three temperature-related contributions to the topological band gap: electrons coupling with short-wavelength phonons, with long-wavelength phonons \textit{via} Fröhlich coupling, and thermal expansion. We find that electron-phonon coupling promotes the topology of the electronic structures of all 1T$^\prime$-MX$_2$ monolayers, while thermal expansion acts as a counteracting effect. Additionally, we derive the band renormalization from Fröhlich coupling in the two-dimensional context and observe its relatively modest contribution to 1T$^\prime$-MX$_2$ monolayers. Finally, we present a simplified physical picture to understand the "inverse Varshni" effect driven by band inversion in topological insulators. Our work reveals that, among the four 1T$^\prime$-MX$_2$ studied monolayers, MoSe$_2$ is a promising candidate for room temperature applications as its band gap displays remarkable resilience against thermal expansion, while the topological order of WS$_2$ can be tuned under the combined influence of strain and temperature. Both materials represent novel examples of temperature promoted topological insulators.

preprint2022arXiv

Finite Temperature Effects on the X-ray Absorption Spectra of Crystalline Aluminas from First Principles

By including phonon-assisted transitions within plane-wave DFT methods for calculating the X-ray absorption spectrum (XAS) we obtain the Al K-edge XAS at 300 K for two Al$_2$O$_3$ phases. The 300 K XAS reproduces the pre-edge peak for $α$-Al$_2$O$_3$, which is not visible at the static-lattice level of approximation. The 300 K XAS for $γ$-Al$_2$O$_3$ correctly describes two out of the three experimental peaks. We show that the second peak arises from 1s to mixed $s$-$p$ transitions and is absent in the 0 K XAS. This letter serves as a basis for future applications, as the method is generalizable to any atom and edge.

preprint2022arXiv

Nonuniform grids for Brillouin zone integration and interpolation

We present two developments for the numerical integration of a function over the Brillouin zone. First, we introduce a nonuniform grid, which we refer to as the Farey grid, that generalizes regular grids. Second, we introduce symmetry-adapted Voronoi tessellation, a general technique to assign weights to the points in an arbitrary grid. Combining these two developments, we propose a strategy to perform Brillouin zone integration and interpolation that provides a significant computational advantage compared to the usual approach based on regular uniform grids. We demonstrate our methodology in the context of first principles calculations with the study of Kohn anomalies in the phonon dispersions of graphene and MgB2, and in the evaluation of the electron-phonon driven renormalization of the band gaps of diamond and bismuthene. In the phonon calculations, we find speedups by a factor of 3 to 4 when using density functional perturbation theory, and by a factor of 6 to 7 when using finite differences in conjunction with supercells. As a result, the computational expense between density functional perturbation theory and finite differences becomes comparable. For electron-phonon coupling calculations we find even larger speedups. Finally, we also demonstrate that the Farey grid can be expressed as a combination of the widely used regular grids, which should facilitate the adoption of this methodology.

preprint2022arXiv

Temperature Effect on Charge-state Transition Levels of Defects in Semiconductors

Defects are crucial in determining the overall physical properties of semiconductors. Generally, the charge-state transition level (TEL), one of the key physical quantities that determines the dopability of defects in semiconductors, is temperature dependent. However, little is known about the temperature dependence of TEL, and, as a result, almost all existing defect theories in semiconductors are built on a temperature-independent approximation. In this article, by deriving the basic formulas for temperature-dependent TEL, we have established two fundamental rules for the temperature dependence of TEL in semiconductors. Based on these rules, surprisingly, it is found that the temperature dependences of TEL for different defects are rather diverse: it can become shallower, deeper, or stay unchanged. This defect-specific behavior is mainly determined by the synergistic or opposing effects between free energy corrections (determined by the local volume change around the defect during a charge-state transition) and band edge changes (which differ for different semiconductors). These basic formulas and rules, confirmed by a large number of state-of-the-art temperature-dependent defect calculations in GaN, may potentially be widely adopted as guidelines for understanding or optimizing doping behaviors in semiconductors at finite temperatures.

preprint2022arXiv

Tunable photostriction of halide perovskites through energy dependent photoexcitation

Halide perovskites exhibit giant photostriction, that is, volume or shape changes upon illumination. However, the microscopic origin of this phenomenon remains unclear and there are experimental reports of both light-induced lattice expansion and contraction. In this paper we establish a general method, based on first-principles calculations and molecular orbital theory, which provides a microscopic picture of photostriction in insulators based on the orbital characters of their electronic bands near the Fermi level. For lead-halide perovskites, we find that different valence states have different bonding characters, leading to opposing strengthening or weakening of bonds depending on the photoexcitation energy. The overall trend is that light induces lattice contraction at low excitation energies, while giant lattice expansion occurs at high excitation energies, rationalizing experimental reports.

preprint2021arXiv

Degenerate topological line surface phonons in quasi-1D double helix crystal SnIP

Degenerate points/lines in the bulk band structures of crystals have become a staple of the growing number of topological materials. The bulk-boundary correspondence provides a relation between bulk topology and surface states. While line degeneracies of bulk excitations have been extensively characterized, line degeneracies of surface states are not well understood. We show that SnIP, a quasi-one-dimensional van der Waals material with a double helix crystal structure, exhibits topological nodal rings/lines in both the bulk phonon modes and their corresponding surface states. Using a combination of first-principles calculations, symmetry-based indicator theories and Zak phase analysis, we find that two neighbouring bulk nodal rings form doubly degenerate lines in their drumhead-like surface states, which are protected by the combination of time-reversal and glide mirror symmetries $\mathcal{T}\bar{M}_y$. Our results indicate that surface degeneracies can be generically protected by symmetries such as $\mathcal{T}\bar{M}_y$, and phonons provide an ideal platform to explore such degeneracies.

preprint2021arXiv

Manipulation and braiding of Weyl nodes using symmetry-constrained phase transitions

Weyl semimetals are arguably the most paradigmatic form of a gapless topological phase. While the stability of Weyl nodes, as quantified by their topological charge, has been extensively investigated, recent interest has shifted to the manipulation of the location of these Weyl nodes for non-Abelian braiding. To accomplish this braiding it is necessary to drive significant Weyl node motion using realistic experimental parameter changes. We show that a family of phase transitions characterized by certain symmetry constraints impose that the Weyl nodes have to reorganise by a large amount, shifting from one high symmetry plane to another. Additionally, for a subset of pairs of nodes with nontrivial Euler class topology, this reorganization can only occur through a braiding process with adjacent nodes. As a result, the Weyl nodes are forced to move a large distance across the Brillouin zone and to braid, all driven by small temperature changes, a process we illustrate with Cd$_2$Re$_2$O$_7$. Our work opens up routes to readily manipulate Weyl nodes using only slight external parameter changes, paving the way for the practical realization of reciprocal space braiding.

preprint2021arXiv

Multi-gap topology and non-Abelian braiding of phonons from first principles

Non-Abelian states of matter, in which the final state depends on the order of the interchanges of two quasiparticles, can encode information immune from environmental noise with the potential to provide a robust platform for topological quantum computation. We demonstrate that phonons can carry non-Abelian frame charges at the band crossing points of their frequency spectrum, and that external stimuli can drive their braiding. We present a general framework to understand the topological configurations of phonons from first principles calculations using a topological invariant called Euler class, and provide a complete analysis of phonon braiding by combining different topological configurations. Taking a well-known dielectric material, Al$_2$O$_3$, as a representative example, we demonstrate that electrostatic doping gives rise to phonon band inversions that can induce redistribution of the frame charges, leading to non-Abelian braiding of phonons. Our work provides a new quasiparticle platform for realizable non-Abelian braiding in reciprocal space, and expands the toolset for studying braiding processes.

preprint2021arXiv

Phonons as a platform for non-Abelian braiding and its manifestation in layered silicates

Topological phases of matter have revolutionised the fundamental understanding of band theory and hold great promise for next-generation technologies such as low-power electronics or quantum computers. Single-gap topologies have been extensively explored, and a large number of materials have been theoretically proposed and experimentally observed. These ideas have recently been extended to multi-gap topologies with band nodes that carry non-Abelian charges, characterised by invariants that arise by the momentum space braiding of such nodes. However, the constraints placed by the Fermi-Dirac distribution to electronic systems have so far prevented the experimental observation of multi-gap topologies in real materials. Here, we show that multi-gap topologies and the accompanying phase transitions driven by braiding processes can be readily observed in the bosonic phonon spectra of known monolayer silicates. The associated braiding process can be controlled by means of an electric field and epitaxial strain, and involves, for the first time, more than three bands. Finally, we propose that the band inversion processes at the $Γ$ point can be tracked by following the evolution of the Raman spectrum, providing a clear signature for the experimental verification of the band inversion accompanied by the braiding process.

preprint2021arXiv

Topological continuum charges of acoustic phonons in 2D

We analyze the band topology of acoustic phonons in 2D materials by considering the interplay of spatial and internal symmetries with additional constraints that arise from the physical context. These supplemental constraints trace back to the Nambu-Goldstone theorem and the requirements of structural stability. We show that this interplay can give rise to previously unaddressed non-trivial nodal charges that are associated with the crossing of the acoustic phonon branches at the center ($Γ$-point) of the phononic Brillouin zone. We moreover apply our perspective to the concrete context of graphene, where we demonstrate that the phonon spectrum harbors these kinds of non-trivial nodal charges. Apart from its fundamental appeal, this analysis is physically consequential and dictates how the phonon dispersion is affected when graphene is grown on a substrate. Given the generality of our framework, we anticipate that our strategy that thrives on combining physical context with insights from topology should be widely applicable in characterizing systems beyond electronic band theory.

preprint2020arXiv

Boron phosphide as a \emph{p}-type transparent conductor: optical absorption and transport through electron-phonon coupling

Boron phosphide has recently been identified as a potential high hole mobility transparent conducting material. This promise arises from its low hole effective masses. However, BP has a relatively small 2 eV indirect band gap which will affect its transparency. In this work, we computationally study both optical absorption across the indirect gap and phonon-limited electronic transport to quantify the potential of boron phosphide as a \emph{p}-type transparent conductor. We find that phonon-mediated indirect optical absorption is weak in the visible spectrum and that the phonon-limited hole mobility is very high (around 900 cm$^2$/Vs) at room temperature. This exceptional mobility comes from a combination of low hole effective mass and very weak scattering by polar phonon modes. We rationalize the weak scattering by the less ionic bonding in boron phosphide compared to oxides. We suggest this could be a general advantage of non-oxides for \emph{p}-type transparent conducting applications. Using our computed properties, we assess the transparent conductor figure of merit of boron phosphide and shows that it exceeds by one order of magnitude that of established \emph{p}-type transparent conductors, confirming the potential of this material.

preprint2020arXiv

Extracting ice phases from liquid water: why a machine-learning water model generalizes so well

We investigate the structural similarities between liquid water and 53 ices, including 20 knowncrystalline phases. We base such similarity comparison on the local environments that consist of atoms within a certain cutoff radius of a central atom. We reveal that liquid water explores the localenvironments of the diverse ice phases, by directly comparing the environments in these phases using general atomic descriptors, and also by demonstrating that a machine-learning potential trained on liquid water alone can predict the densities, the lattice energies, and vibrational properties of theices. The finding that the local environments characterising the different ice phases are found in water sheds light on water phase behaviors, and rationalizes the transferability of water models between different phases.

preprint2020arXiv

First principles modeling of exciton-polaritons in polydiacetylene chains

Exciton-polaritons in organic materials are hybrid states that result from the strong interaction of photons and the bound excitons that these materials host. Organic polaritons hold great interest for optoelectronic applications, however progress towards this end has been impeded by the lack of a first principles approach that quantifies light-matter interactions in these systems, and which would allow the formulation of molecular design rules. Here we develop such a first principles approach, quantifying light-matter interactions. We exemplify our approach by studying variants of the conjugated polymer polydiacetylene, and we show that a large polymer conjugation length is critical towards strong exciton-photon coupling, hence underlying the importance of pure structures without static disorder. By comparing to our experimental reflectivity measurements, we show that the coupling of excitons to vibrations, manifested by phonon side bands in the absorption, has a strong impact on the magnitude of light-matter coupling over a range of frequencies. Our approach opens the way towards a deeper understanding of polaritons in organic materials, and we highlight that a quantitatively accurate calculation of the exciton-photon interaction would require accounting for all sources of disorder self-consistently.

preprint2020arXiv

Sub-picosecond photo-induced displacive phase transition in two-dimensional MoTe$_2$

Photo-induced phase transitions (PIPTs) provide an ultrafast, energy-efficient way for precisely manipulating the topological properties of transition-metal ditellurides, and can be used to stabilize a topological phase in an otherwise semiconducting material. Using first-principles calculations, we demonstrate that the PIPT in monolayer MoTe$_2$ from the semiconducting 2H phase to the topological 1T$'$ phase can be triggered purely by electronic excitations that soften multiple lattice vibrational modes. These softenings, driven by a Peierls-like mechanism within the conduction bands, lead to structural symmetry breaking within sub-picosecond timescales, which is shorter than the timescale of a thermally driven phase transition. The transition is predicted to be triggered by photons with energies over $1.96$\,eV, with an associated excited carrier density of $3.4\times10^{14}$\,cm$^{-2}$, which enables a controllable phase transformation by varying the laser wavelength. Our results provide insight into the underlying physics of the phase transition in 2D transition-metal ditellurides, and show an ultrafast phase transition mechanism for manipulation of the topological properties of 2D systems.

preprint2020arXiv

The impact of exciton delocalization on exciton-vibration interactions in organic semiconductors

Organic semiconductors exhibit properties of individual molecules and extended crystals simultaneously. The strongly bound excitons they host are typically described in the molecular limit, but excitons can delocalize over many molecules, raising the question of how important the extended crystalline nature is. Using accurate Green's function based methods for the electronic structure and non-perturbative finite difference methods for exciton-vibration coupling, we describe exciton interactions with molecular and crystal degrees of freedom concurrently. We find that the degree of exciton delocalization controls these interactions, with thermally activated crystal phonons predominantly coupling to delocalized states, and molecular quantum fluctuations predominantly coupling to localized states. Based on this picture, we quantitatively predict and interpret the temperature and pressure dependence of excitonic peaks in the acene series of organic semiconductors, which we confirm experimentally, and we develop a simple experimental protocol for probing exciton delocalization. Overall, we provide a unified picture of exciton delocalization and vibrational effects in organic semiconductors, reconciling the complementary views of finite molecular clusters and periodic molecular solids.

preprint2016arXiv

Correlation effects on electron-phonon coupling in semiconductors: many-body theory along thermal lines

A method is proposed for the inclusion of electron correlation in the calculation of the temperature dependence of band structures arising from electron-phonon coupling. It relies on an efficient exploration of the vibrational phase space along the recently introduced thermal lines. Using the $G_0W_0$ approximation, the temperature dependence of the direct gaps of diamond, silicon, lithium fluoride, magnesium oxide, and titanium dioxide is calculated. Within the proposed formalism, a single calculation at each temperature of interest is sufficient to obtain results of the same accuracy as in alternative, more expensive methods. It is shown that many-body contributions beyond semi-local density functional theory modify the electron-phonon coupling strength by almost $50$% in diamond, silicon, and titanium dioxide, but by less than $5$% in lithium flouride and magnesium oxide. The results reveal a complex picture regarding the validity of semi-local functionals for the description of electron-phonon coupling.

preprint2016arXiv

Hexagonal structure of phase III of solid hydrogen

A hexagonal structure of solid molecular hydrogen with $P6_122$ symmetry is calculated to be more stable below about 200 GPa than the monoclinic $C2/c$ structure identified previously as the best candidate for phase III. We find that the effects of nuclear quantum and thermal vibrations play a central role in the stabilization of $P6_122$. The $P6_122$ and $C2/c$ structures are very similar and their Raman and infra-red data are in good agreement with experiment. However, our calculations show that the hexagonal $P6_122$ structure provides better agreement with the available x-ray diffraction data than the $C2/c$ structure at pressures below about 200 GPa. We suggest that two phase-III-like structures may be formed at high pressures, hexagonal $P6_122$ below about 200 GPa and monoclinic $C2/c$ at higher pressures.

preprint2016arXiv

Temperature Dependence of the Energy Levels of Methylammonium Lead Iodide Perovskite from First Principles

Environmental effects and intrinsic energy-loss processes lead to fluctuations in the operational temperature of solar cells, which can profoundly influence their power conversion efficiency. Here we determine from first principles the effects of temperature on the band gap and band edges of the hybrid pervoskite CH$_3$NH$_3$PbI$_3$ by accounting for electron-phonon coupling and thermal expansion. From $290$ to $380$ K, the computed band gap change of $40$ meV coincides with the experimental change of $30$-$40$ meV. The calculation of electron-phonon coupling in CH$_3$NH$_3$PbI$_3$ is particularly intricate, as the commonly used Allen-Heine-Cardona theory overestimates the band gap change with temperature, and excellent agreement with experiment is only obtained when including high-order terms in the electron-phonon interaction. We also find that spin-orbit coupling enhances the electron-phonon coupling strength, but that the inclusion of nonlocal correlations using hybrid functionals has little effect. We reach similar conclusions in the metal-halide perovskite CsPbI$_3$. Our results unambiguously confirm for the first time the importance of high-order terms in the electron-phonon coupling by direct comparison with experiment.

preprint2016arXiv

Temperature effects in the band structure of topological insulators

We study the effects of temperature on the band structure of the Bi$_2$Se$_3$ family of topological insulators using first-principles methods. Increasing temperature drives these materials towards the normal state, with similar contributions from thermal expansion and from electron-phonon coupling. The band gap changes with temperature reach $0.3$ eV at $600$ K, of similar size to the changes caused by electron correlation. Our results suggest that temperature-induced topological phase transitions should be observable near critical points of other external parameters.

preprint2016arXiv

Vibrational renormalisation of the electronic band gap in hexagonal and cubic ice

Electron-phonon coupling in hexagonal and cubic water ice is studied using first-principles quantum mechanical methods. We consider 29 distinct hexagonal and cubic ice proton-orderings with up to 192 molecules in the simulation cell to account for proton-disorder. We find quantum zero-point vibrational corrections to the minimum electronic band gaps ranging from -1.5 to -1.7 eV, which leads to improved agreement between calculated and experimental band gaps. Anharmonic nuclear vibrations play a negligible role in determining the gaps. Deuterated ice has a smaller band-gap correction at zero-temperature of -1.2 to -1.4eV. Vibrations reduce the differences between the electronic band gaps of different proton-orderings from around 0.17 eV to less than 0.05 eV, so that the electronic band gaps of hexagonal and cubic ice are almost independent of the proton-ordering when quantum nuclear vibrations are taken into account. The comparatively small reduction in the band gap over the temperature range 0-240 K of around 0.1 eV does not depend on the proton ordering, or whether the ice is protiated or deuterated, or hexagonal or cubic. We explain this in terms of the atomistic origin of the strong electron-phonon coupling in ice.

preprint2015arXiv

Anharmonic nuclear motion and the relative stability of hexagonal and cubic ice

We use extensive first-principles quantum mechanical calculations to show that, although the static lattice and harmonic vibrational energies are almost identical, the anharmonic vibrational energy of hexagonal ice is significantly lower than that of cubic ice. This difference in anharmonicity is crucial, stabilising hexagonal ice compared with cubic ice by at least 1.4 meV/H2O, in agreement with experimental estimates. The difference in anharmonicity arises predominantly from molecular O-H bond stretching vibrational modes and is related to the different stacking of atomic layers.

preprint2015arXiv

Giant electron-phonon interactions in molecular crystals and the importance of non-quadratic coupling

We investigate electron-phonon coupling in the molecular crystals CH$_4$, NH$_3$, H$_2$O, and HF, using first-principles quantum mechanical calculations. We find vibrational corrections to the electronic band gaps at zero temperature of -1.97 eV, -1.01 eV, -1.52 eV, and -1.62 eV, respectively, which are comparable in magnitude to those from electron-electron correlation effects. Microscopically, the strong electron-phonon coupling arises in roughly equal measure from the almost dispersionless high-frequency molecular modes and from the lower frequency lattice modes. We also highlight the limitations of the widely used Allen-Heine-Cardona theory, which gives significant discrepancies compared to our more accurate treatment.

preprint2015arXiv

Lattice dynamics and electron-phonon coupling calculations using non-diagonal supercells

We study the direct calculation of total energy derivatives for lattice dynamics and electron-phonon coupling calculations using supercell matrices with non-zero off-diagonal elements. We show that it is possible to determine the response of a periodic system to a perturbation characterized by a wave vector with reduced fractional coordinates $(m_1/n_1,m_2/n_2,m_3/n_3)$ using a supercell containing a number of primitive cells equal to the least common multiple of $n_1$, $n_2$, and $n_3$. If only diagonal supercell matrices are used, a supercell containing $n_1n_2n_3$ primitive cells is required. We demonstrate that the use of non-diagonal supercells significantly reduces the computational cost of obtaining converged zero-point energies and phonon dispersions for diamond and graphite. We also perform electron-phonon coupling calculations using the direct method to sample the vibrational Brillouin zone with grids of unprecedented size, which enables us to investigate the convergence of the zero-point renormalization to the thermal and optical band gaps of diamond.

preprint2015arXiv

Quantum Monte Carlo study of the phase diagram of solid molecular hydrogen at extreme pressures

Establishing the phase diagram of hydrogen is a major challenge for experimental and theoretical physics. Experiment alone cannot establish the atomic structure of solid hydrogen at high pressure, because hydrogen scatters X-rays only weakly. Instead our understanding of the atomic structure is largely based on density functional theory (DFT). By comparing Raman spectra for low-energy structures found in DFT searches with experimental spectra, candidate atomic structures have been identified for each experimentally observed phase. Unfortunately, DFT predicts a metallic structure to be energetically favoured at a broad range of pressures up to 400 GPa, where it is known experimentally that hydrogen is nonmetallic. Here we show that more advanced theoretical methods (diffusion quantum Monte Carlo calculations) find the metallic structure to be uncompetitive, and predict a phase diagram in reasonable agreement with experiment. This greatly strengthens the claim that the candidate atomic structures accurately model the experimentally observed phases.

preprint2015arXiv

Vibrational averages along thermal lines

A method is proposed for the calculation of vibrational quantum and thermal expectation values of physical properties from first principles. Thermal lines are introduced: these are lines in configuration space parametrized by temperature, such that the value of any physical property along them is approximately equal to the vibrational average of that property. The number of sampling points needed to explore the vibrational phase space is reduced by up to an order of magnitude when the full vibrational density is replaced by thermal lines. Calculations of the vibrational averages of several properties and systems are reported, namely the internal energy and the electronic band gap of diamond and silicon, and the chemical shielding tensor of L-alanine. Thermal lines pave the way for complex calculations of vibrational averages, including large systems and methods beyond semi-local density functional theory.

preprint2014arXiv

Comparing electron-phonon coupling strength in diamond, silicon and silicon carbide: First-principles study

Renormalization of the electronic band gap due to electron-phonon coupling in the tetrahedral semiconductors diamond, silicon and cubic silicon carbide is studied from first principles. There is a marked difference between the coupling of the vibrational state to the valence band maximum and to the conduction band minimum. The strength of phonon coupling to the valence band maximum is similar between the three systems and is dominated by vibrations that change the bond length. The coupling strength to the conduction band minimum differs significantly in diamond, silicon carbide and silicon. In diamond, the coupling is dominated by six small pockets of vibrational states in the phonon Brillouin zone, that are ultimately responsible for the stronger electron-phonon coupling in this material. Our results represent a first step towards the development of an a priori understanding of electron-phonon coupling in semiconductors and insulators, that should aid the design of materials with tailored electron-phonon coupling properties.

preprint2014arXiv

Dissociation of high-pressure solid molecular hydrogen: Quantum Monte Carlo and anharmonic vibrational study

A theoretical study is reported of the molecular-to-atomic transition in solid hydrogen at high pressure. We use the diffusion quantum Monte Carlo method to calculate the static lattice energies of the competing phases and a density-functional-theory-based vibrational self-consistent field method to calculate anharmonic vibrational properties. We find a small but significant contribution to the vibrational energy from anharmonicity. A transition from the molecular Cmca-12 direct to the atomic I4_1/amd phase is found at 374 GPa. The vibrational contribution lowers the transition pressure by 91 GPa. The dissociation pressure is not very sensitive to the isotopic composition. Our results suggest that quantum melting occurs at finite temperature.

preprint2014arXiv

Temperature effects in first-principles solid state calculations of the chemical shielding tensor made simple

We study the effects of atomic vibrations on the solid-state chemical shielding tensor using first principles density functional theory calculations. At the harmonic level, we use a Monte Carlo method and a perturbative expansion. The Monte Carlo method is accurate but computationally expensive, while the perturbative method is computationally more efficient, but approximate. We find excellent agreement between the two methods for both the isotropic shift and the shielding anisotropy. The effects of zero-point quantum mechanical nuclear motion are important up to relatively high temperatures: at 500 K they still represent about half of the overall vibrational contribution. We also investigate the effects of anharmonic vibrations, finding that their contribution to the zero-point correction to the chemical shielding tensor is small. We exemplify these ideas using magnesium oxide and the molecular crystals L-alanine and beta-aspartyl-L-alanine. We therefore propose as the method of choice to incorporate the effects of temperature in solid state chemical shielding tensor calculations the perturbative expansion within the harmonic approximation. This approach is accurate and requires a computational effort that is about an order of magnitude smaller than that of dynamical or Monte Carlo approaches, so these effects might be routinely accounted for.

preprint2013arXiv

Anharmonic vibrational properties in periodic systems: energy, electron-phonon coupling, and stress

A unified approach is used to study vibrational properties of periodic systems with first-principles methods and including anharmonic effects. Our approach provides a theoretical basis for the determination of phonon-dependent quantities at finite temperatures. The low-energy portion of the Born-Oppenheimer energy surface is mapped and used to calculate the total vibrational energy including anharmonic effects, electron-phonon coupling, and the vibrational contribution to the stress tensor. We report results for the temperature dependence of the electronic band gap and the linear coefficient of thermal expansion of diamond, lithium hydride, and lithium deuteride.

preprint2013arXiv

Electron-Phonon Coupling and the Metalization of Solid Helium at Terapascal Pressures

Solid He is studied in the pressure and temperature ranges 1-40 TPa and 0-10,000 K using first-principles methods. Anharmonic vibrational properties are calculated within a self-consistent field framework, including the internal and free energies, density-pressure relation, stress tensor, thermal expansion, and the electron-phonon coupling renormalization of the electronic band gap. We find that an accurate description of electron-phonon coupling requires us to use a non-perturbative approach. The metalization pressure of 32.9 TPa at 0 K is larger than found previously. The vibrational effects are large; for example at P=30 TPa the band gap is increased by 2.8 eV by electron-phonon coupling and a further 0.1 eV by thermal expansion compared to the static value. The implications of the calculated metalization pressure for the cooling of white dwarfs are discussed.

preprint2013arXiv

Theory of temperature dependent phonon-renormalized properties

We present a general harmonic theory for the temperature dependence of phonon-renormalized properties of solids. Firstly, we formulate a perturbation theory in phonon-phonon interactions to calculate the phonon renormalization of physical quantities. Secondly, we propose two new schemes for extrapolating phonon zero-point corrections from temperature dependent data that improve the accuracy by an order of magnitude compared to previous approaches. Finally, we consider the low-temperature limit of the class of observables that includes the electronic band gap, obtaining a $T^4$ dependence in three dimensions, $T^2$ in two dimensions, and $T^{3/2}$ in one dimension.