Jesús Carrete

Thermal conductivity and phonon linewidths of monolayer MoS2 from first principles

Using ab initio calculations, we have investigated the phonon linewidths and the thermal conductivity (κ) of monolayer MoS2. κ for a typical sample size of 1 μm is 83 W/m K at room temperature in the completely rough edge limit, suggesting κ is not a limiting factor for the electronic application of monolayer MoS2. κ can be further increased by 30% in 10 μm sized samples. Due to strong anharmonicity, isotope enhancement of room temperature κ is only 10% for 1 μm sized samples. However, linewidths can be significantly reduced, for instance, for Raman active modes A1g and E2g1, in isotopically pure samples.

Molecular Dynamics Simulation of the Structure and Dynamics of Water–1-Alkyl-3-methylimidazolium Ionic Liquid Mixtures

We have performed extensive molecular dynamic simulations to analyze the influence of cation and anion natures, and of water concentration, on the structure and dynamics of water-1-alkyl-3-methylimidazolium ionic liquid mixtures. The dependence on water concentration of the radial distribution functions, coordination numbers, and hydrogen bonding degree between the different species has been systematically analyzed for different lengths of the cation alkyl chain (alkyl = ethyl, butyl, hexyl, and octyl) and several counterions. These include two halogens of different sizes and positions in Hoffmeister series, Cl(-) and Br(-), and the highly hydrophobic inorganic anion PF(6)(-) throughout its whole solubility regime. The formation of water clusters in the mixture has been verified, and the influences of both anion hydrophobicity and cation chain length on the structure and size of these clusters have been analyzed. The water cluster size is shown to be relatively independent of the cation chain length, but strongly dependent on the hydrophobicity of the anion, which also determines critically the network formation of water and therefore the miscibility of the ionic liquid. The greater influence of the anion relative to the cation one is seen to be reflected in all the analyzed physical properties. Finally, single-particle dynamics in IL-water mixtures is considered, obtaining the self-diffusion coefficients and the velocity autocorrelation functions of water molecules in the mixture, and analyzing the effect of cation, anion, and water concentration on the duration of the ballistic regime and on the time of transition to the diffusive regime. Complex non-Markovian behavior was detected at intermediate times within an interval progressively shorter as water concentration increases.

Finding Unprecedentedly Low-Thermal-Conductivity Half-Heusler Semiconductors via High-Throughput Materials Modeling

The lattice thermal conductivity ({\kappa}{\omega}) is a key property for many potential applications of compounds. Discovery of materials with very low or high {\kappa}{\omega} remains an experimental challenge due to high costs and time-consuming synthesis procedures. High-throughput computational pre-screening is a valuable approach for significantly reducing the set of candidate compounds. In this article, we introduce efficient methods for reliably estimating the bulk {\kappa}{\omega} for a large number of compounds. The algorithms are based on a combination of machine-learning algorithms, physical insights, and automatic ab-initio calculations. We scanned approximately 79,000 half-Heusler entries in the AFLOWLIB.org database. Among the 450 mechanically stable ordered semiconductors identified, we find that {\kappa}{\omega} spans more than two orders of magnitude- a much larger range than that previously thought. {\kappa}{\omega} is lowest for compounds whose elements in equivalent positions have large atomic radii. We then perform a thorough screening of thermodynamical stability that allows to reduce the list to 77 systems. We can then provide a quantitative estimate of {\kappa}{\omega} for this selected range of systems. Three semiconductors having {\kappa}{\omega} < 5 W /(m K) are proposed for further experimental study.

MD simulations of the formation of stable clusters in mixtures of alkaline salts and imidazolium-based ionic liquids.

Structural and dynamical properties of room-temperature ionic liquids containing the cation 1-butyl-3-methylimidazolium ([BMIM](+)) and three different anions (hexafluorophosphate, [PF6](-), tetrafluoroborate, [BF4](-), and bis(trifluoromethylsulfonyl)imide, [NTf2](-)) doped with several molar fractions of lithium salts with a common anion at 298.15 K and 1 atm were investigated by means of molecular dynamics simulations. The effect of the size of the salt cation was also analyzed by comparing these results with those for mixtures of [BMIM][PF6] with NaPF6. Lithium/sodium solvation and ionic mobilities were analyzed via the study of radial distribution functions, coordination numbers, cage autocorrelation functions, mean-square displacements (including the analysis of both ballistic and diffusive regimes), self-diffusion coefficients of all the ionic species, velocity and current autocorrelation functions, and ionic conductivity in all the ionic liquid/salt systems. We found that lithium and sodium cations are strongly coordinated in two different positions with the anion present in the mixture. Moreover, [Li](+) and [Na](+) cations were found to form bonded-like, long-lived aggregates with the anions in their first solvation shell, which act as very stable kinetic entities within which a marked rattling motion of salt ions takes place. With very long MD simulation runs, this phenomenon is proved to be on the basis of the decrease of self-diffusion coefficients and ionic conductivities previously reported in experimental and computational results.

Unexpected High-Temperature Stability of β-Zn4Sb3 Opens the Door to Enhanced Thermoelectric Performance

β-Zn4Sb3 has one of the highest ZT reported for binary compounds, but its practical applications have been hindered by a reported poor stability. Here we report the fabrication of nearly dense single-phase β-Zn4Sb3 and a study of its thermoelectric transport coefficients across a wide temperature range. Around 425 K we find an abrupt decrease of its thermal conductivity. Past this point, Zn atoms can migrate from crystalline sites to interstitial positions; β-Zn4Sb3 becomes metastable and gradually decomposes into Zn(hcp) and ZnSb. However, above 565 K it recovers its stability; in fact, the damage caused by decomposition can be repaired completely. This is key to its excellent thermoelectric performance at high temperature: the maximum ZT reaches 1.4. Molecular dynamics simulations are used to shed light on the microscopic behavior of the material.

Solvation of lithium salts in protic ionic liquids: a molecular dynamics study.

The structure of solutions of lithium nitrate in a protic ionic liquid with a common anion, ethylammonium nitrate, at room temperature is investigated by means of molecular dynamics simulations. Several structural properties, such as density, radial distribution functions, hydrogen bonds, spatial distribution functions, and coordination numbers, are analyzed in order to get a picture of the solvation of lithium cations in this hydrogen-bonded, amphiphilically nanostructured environment. The results reveal that the ionic liquid mainly retains its structure upon salt addition, the interaction between the ammonium group of the cation and the nitrate anion being only slightly perturbed by the addition of the salt. Lithium cations are solvated by embedding them in the polar nanodomains of the solution formed by the anions, where they coordinate with the latter in a solid-like fashion reminiscent of a pseudolattice structure. Furthermore, it is shown that the average coordination number of [Li](+) with the anions is 4, nitrate coordinating [Li](+) in both monodentate and bidentate ways, and that in the second coordination layer both ethylammonium cations and other lithiums are also found. Additionally, the rattling motion of lithium ions inside the cages formed by their neighboring anions, indicative of the so-called caging effect, is confirmed by the analysis of the [Li](+) velocity autocorrelation functions. The overall picture indicates that the solvation of [Li](+) cations in this amphiphilically nanostructured environment takes place by means of a sort of inhomogeneous nanostructural solvation, which we could refer to as nanostructured solvation, and which could be a universal solvation mechanism in ionic liquids.

Surface roughness and thermal conductivity of semiconductor nanowires: Going below the Casimir limit

By explicitly considering surface roughness at the atomic level, we quantitatively show that the thermal conductivity of Si nanowires can be lower than Casimirs classical limit. However, this violation only occurs for deep surface degradation. For shallow surface roughness, the Casimir formula is shown to yield a good approximation to the phonon mean free paths and conductivity, even for nanowire diameters as thin as 2.22 nm. Our exact treatment of roughness scattering is in stark contrast with a previously proposed perturbative approach, which is found to overpredict scattering rates by an order of magnitude. The obtained results suggest that a complete theoretical understanding of some previously published experimental results is still lacking.

Physically founded phonon dispersions of few-layer materials and the case of borophene

ABSTRACT By building physically sound interatomic force constants, we offer evidence of the universal presence of a quadratic phonon branch in all unstrained 2D materials, thus contradicting much of the existing literature. Through a reformulation of the interatomic force constants (IFCs) in terms of internal coordinates, we find that a delicate balance between the IFCs is responsible for this quadraticity. We use this approach to predict the thermal conductivity of Pmmn borophene, which is comparable to that of , and displays a remarkable in-plane anisotropy. These qualities may enable the efficient heat management of borophene devices in potential nanoelectronic applications. IMPACT STATEMENT The newly found universality of quadratic dispersion will change the way 2D-material phonons are calculated. Predicted results for borophene shall become a fundamental reference for future research on this material. GRAPHICAL ABSTRACT

Mixtures of protic ionic liquids and molecular cosolvents: a molecular dynamics simulation.

In this work, the effect of molecular cosolvents (water, ethanol, and methanol) on the structure of mixtures of these compounds with a protic ionic liquid (ethylammonium nitrate) is analyzed by means of classical molecular dynamics simulations. Included are as-yet-unreported measurements of the densities of these mixtures, used to test our parameterized potential. The evolution of the structure of the mixtures throughout the concentration range is reported by means of the calculation of coordination numbers and the fraction of hydrogen bonds in the system, together with radial and spatial distribution functions for the various molecular species and molecular ions in the mixture. The overall picture indicates a homogeneous mixing process of added cosolvent molecules, which progressively accommodate themselves in the network of hydrogen bonds of the protic ionic liquid, contrarily to what has been reported for their aprotic counterparts. Moreover, no water clustering similar to that in aprotic mixtures is detected in protic aqueous mixtures, but a somehow abrupt replacing of [NO3](-) anions in the first hydration shell of the polar heads of the ionic liquid cations is registered around 60% water molar concentration. The spatial distribution functions of water and alcohols differ in the coordination type, since water coordinates with [NO3](-) in a bidentate fashion in the equatorial plane of the anion, while alcohols do it in a monodentate fashion, competing for the oxygen atoms of the anion. Finally, the collision times of the different cosolvent molecules are also reported by calculating their velocity autocorrelation functions, and a caging effect is observed for water molecules but not in alcohol mixtures.

Twisting phonons in complex crystals with quasi-one-dimensional substructures

A variety of crystals contain quasi-one-dimensional substructures, which yield distinctive electronic, spintronic, optical and thermoelectric properties. There is a lack of understanding of the lattice dynamics that influences the properties of such complex crystals. Here we employ inelastic neutron scatting measurements and density functional theory calculations to show that numerous low-energy optical vibrational modes exist in higher manganese silicides, an example of such crystals. These optical modes, including unusually low-frequency twisting motions of the Si ladders inside the Mn chimneys, provide a large phase space for scattering acoustic phonons. A hybrid phonon and diffuson model is proposed to explain the low and anisotropic thermal conductivity of higher manganese silicides and to evaluate nanostructuring as an approach to further suppress the thermal conductivity and enhance the thermoelectric energy conversion efficiency. This discovery offers new insights into the structure-property relationships of a broad class of materials with quasi-one-dimensional substructures for various applications.