Probing the thermal resistance of solid-liquid interfaces in nanofluids with molecular dynamics.

The significance of interfacial thermal resistance in the thermal conductivity of nanofluids is not well understood, in part because of the absence of measurements of this quantity. Here, we study the interfacial thermal resistance for metal-oil nanofluids of interest as heat transfer fluids for concentrating solar power, using density functional theory and molecular dynamics simulations. Insights on the role of chemical interactions in determining the interfacial thermal resistance are revealed. The results presented here showcase a general picture in which the stronger the chemical interactions between species at the interface, the lower the associated interfacial thermal resistance. The implications toward nanofluid design are discussed. We show that, for this important family of metal-oil nanofluids, the interfacial thermal resistance values are low enough so that it is possible to afford a reduction in particle size, minimizing stability and rheological issues while still offering enhancement in the effective thermal conductivity with respect to the base fluid.

Thermodynamic analysis of the interaction between metal vacancies and hydrogen in bulk Cu.

Using grand canonical thermodynamic analysis with inputs from DFT calculations we calculated equilibrium molar fractions of copper vacancies (VCu), H interstitials (Hi) and their complexes in bulk Cu in a wide range of temperature and hydrogen pressure values. The results show that the equilibrium molar fractions of both VCu and Hi are low in most conditions of interest, in good agreement with available experimental data. Although Hi-VCu complexes have significantly lower formation energies than the isolated defects, the low molar fraction of H is predicted to have little impact on the rise in vacancy molar fraction for external hydrogen pressures below 100 bar. Only at relatively high hydrogen pressures exceeding 10 kbar in the presence of Cu vacancies, the H molar fraction was found to reach the same order of magnitude as the molar fraction of vacancies. These results put thermodynamic limits on the hydrogen-induced vacancy clustering and void formation in bulk Cu.

Linking in situ charge accumulation to electronic structure in doped SrTiO3 reveals design principles for hydrogen-evolving photocatalysts.

Recently, high solar-to-hydrogen efficiencies were demonstrated using La and Rh co-doped SrTiO3 (La,Rh:SrTiO3) incorporated into a low-cost and scalable Z-scheme device, known as a photocatalyst sheet. However, the unique properties that enable La,Rh:SrTiO3 to support this impressive performance are not fully understood. Combining in situ spectroelectrochemical measurements with density functional theory and photoelectron spectroscopy produces a depletion model of Rh:SrTiO3 and La,Rh:SrTiO3 photocatalyst sheets. This reveals remarkable properties, such as deep flatband potentials (+2 V versus the reversible hydrogen electrode) and a Rh oxidation state dependent reorganization of the electronic structure, involving the loss of a vacant Rh 4d mid-gap state. This reorganization enables Rh:SrTiO3 to be reduced by co-doping without compromising the p-type character. In situ time-resolved spectroscopies show that the electronic structure reorganization induced by Rh reduction controls the electron lifetime in photocatalyst sheets. In Rh:SrTiO3, enhanced lifetimes can only be obtained at negative applied potentials, where the complete Z-scheme operates inefficiently. La co-doping fixes Rh in the 3+ state, which results in long-lived photogenerated electrons even at very positive potentials (+1 V versus the reversible hydrogen electrode), in which both components of the complete device operate effectively. This understanding of the role of co-dopants provides a new insight into the design principles for water-splitting devices based on bandgap-engineered metal oxides.

Adsorption of Aspartic Acid on Ni{100}: A Combined Experimental and Theoretical Study.

Understanding the interaction of amino acids with metal surfaces is essential for the rational design of chiral modifiers able to confer enantioselectivity to metal catalysts. Here, we present an investigation of the adsorption of aspartic acid (Asp) on the Ni{100} surface, using a combination of synchrotron X-ray photoelectron spectroscopy (XPS), near-edge X-ray absorption fine structure, and density functional theory simulations. Based on the combined analysis of the experimental and simulated data, we can identify the dominant mode of adsorption as a pentadentate configuration with three O atoms at the bridge sites of the surfaces, and the remaining oxygen atom and the amino nitrogen are located on atop sites. From temperature-programmed XPS measurements, it was found that Asp starts decomposing above 400 K, which is significantly higher than typical decomposition temperatures of smaller organic molecules on Ni surfaces. Our results offer valuable insights into understanding the role of Asp as a chiral modifier of nickel catalyst surfaces in enantioselective hydrogenation reactions.

Ensemble-Based Modeling of the NMR Spectra of Solid Solutions: Cation Disorder in Y2(Sn,Ti)2O7.

The sensitivity of NMR to the local environment, without the need for any long-range order, makes it an ideal tool for the characterization of disordered materials. Computational prediction of NMR parameters can be of considerable help in the interpretation and assignment of NMR spectra of solids, but the statistical representation of all possible chemical environments for a solid solution is challenging. Here, we illustrate the use of a symmetry-adapted configurational ensemble in the simulation of NMR spectra, in combination with solid-state NMR experiments. We show that for interpretation of the complex and overlapped lineshapes that are typically observed, it is important to go beyond a single-configuration representation or a simple enumeration of local environments. The ensemble method leads to excellent agreement between simulated and experimental spectra for Y2(Sn,Ti)2O7 pyrochlore ceramics, where the overlap of signals from different local environments prevents a simple decomposition of the experimental spectral lineshapes. The inclusion of a Boltzmann weighting confirms that the best agreement with experiment is obtained at higher temperatures, in the limit of full disorder. We also show that to improve agreement with experiment, in particular at low dopant concentrations, larger supercells are needed, which might require alternative simulation approaches as the complexity of the system increases. It is clear that ensemble-based modeling approaches in conjunction with NMR spectroscopy offer great potential for understanding configurational disorder, ultimately aiding the future design of functional materials.

Double-well potential energy surface in the interaction between h-BN and Ni(111).

Density functional theory calculations with non-local correlation functionals, properly accounting for dispersion forces, predict the presence of two minima in the interaction energy between h-BN and Ni(111). These can be described as a physisorbed state with no corrugation of the h-BN structure, and a chemisorbed state exhibiting noticeable corrugation and a shorter distance of h-BN to the metallic support. The latter corresponds indeed to the one reported in most experiments. The relative stability of the two minima depends on the specific density functional employed: of those investigated here only optB86b-vdW yields the correct order of stability. We also demonstrate that the effect of the metal support on the Raman frequency of the chemisorbed boron nitride monolayer cannot be reduced to the associated strain. This is important because the Raman frequency has been proposed as a signature to identify h-BN monolayers from multilayered samples. Our analysis shows that such signatures would be strongly dependent on the nature of the interaction between the support and h-BN.

Understanding the origin of disorder in kesterite-type chalcogenides A2ZnBQ4 (A = Cu, Ag; B = Sn, Ge; Q = S, Se): the influence of inter-layer interactions.

Semiconducting quaternary chalcogenides with A2ZnBQ4 stoichiometry, where A and B are monovalent and tetravalent metal ions and Q is a chalcogen (e.g. Cu2ZnSnS4 or CZTS) have recently attracted attention as potential solar-cell absorbers made from abundant and non-toxic elements. Unfortunately, they exhibit relatively poor sunlight conversion efficiencies, which has been linked to site disorder within the tetrahedral cation sub-lattice. In order to gain a better understanding of the factors controlling cation disorder in these chalcogenides, we have used powder neutron diffraction, coupled with Density Functional Theory (DFT) simulations, to investigate the detailed structure of A2ZnBQ4 phases, with A = Cu, Ag; B = Sn, Ge; and Q = S, Se. Both DFT calculations and powder neutron diffraction data demonstrate that the kesterite structure (space group: I4[combining macron]) is adopted in preference to the higher-energy stannite structure (space group: I4[combining macron]2m). The contrast between the constituent cations afforded by neutron diffraction reveals that copper and zinc cations are only partially ordered in the kesterites Cu2ZnBQ4 (B = Sn, Ge), whereas the silver-containing phases are fully ordered. The degree of cation order in the copper-containing phases shows a greater sensitivity to the identity of the B-cation than to the chalcogenide anion. DFT indicates that cation ordering minimises inter-planar Zn2+Zn2+ electrostatic interactions, while there is an additional intra-planar energy contribution associated with size mismatch. The complete Ag/Zn order in Ag2ZnBQ4 (B = Sn, Ge) phases can thus be related to the anisotropic expansion of the unit cell on replacing Cu with Ag.

Screening heteroatom distributions in zeotype materials using an effective Hamiltonian approach: the case of aluminogermanate PKU-9.

We introduce a method to allow the screening of large configurational spaces of heteroatom distributions in zeotype materials. Based on interatomic potential calculations of configurations containing up to two heteroatoms per cell, we parameterize an atomistic effective Hamiltonian to describe the energy of multiple substitutions, with consideration of both short- and long-range interactions. Then, the effective Hamiltonian is used to explore the full configurational space at other compositions, allowing the identification of the most stable structures for further analysis. We illustrate our approach with the aluminogermanate PKU-9, where we show that increasing the aluminium concentration changes the likely siting of Al, in agreement with experiment.

Critical Role of Dynamic Flexibility in Ge-Containing Zeolites: Impact on Diffusion.

Incorporation of germanium in zeolites is well known to confer static flexibility to their framework, by stabilizing the formation of small rings. In this work, we show that the flexibility associated to Ge atoms in zeolites goes beyond this static effect, manifesting also a clear dynamic nature, in the sense that it leads to enhanced molecular diffusion. Our study combines experimental and theoretical methods providing evidence for this effect, which has not been described previously, as well as a rationalization for it, based on atomistic grounds. We have used both pure-silica and silico-germanate ITQ-29 (LTA topology) zeolites as a case study. Based on our simulations, we identify the flexibility associated to the pore breathing-like behavior induced by the Ge atoms, as the key factor leading to the enhanced diffusion observed experimentally in Ge-containing zeolites.

Modelling a Linker Mix-and-Match Approach for Controlling the Optical Excitation Gaps and Band Alignment of Zeolitic Imidazolate Frameworks.

Tuning the electronic structure of metal-organic frameworks is the key to extending their functionality to the photocatalytic conversion of absorbed gases. Herein we discuss how the band edge positions in zeolitic imidazolate frameworks (ZIFs) can be tuned by mixing different imidazole-based linkers within the same structure. We present the band alignment for a number of known and hypothetical Zn-based ZIFs with respect to the vacuum level. Structures with a single type of linker exhibit relatively wide band gaps; however, by mixing linkers of a low-lying conduction edge with linkers of a high-lying valence edge, we can predict materials with ideal band positions for visible-light water splitting and CO2 reduction photocatalysis. By introducing copper in the tetrahedral position of the mixed-linker ZIFs, it would be possible to increase both photo-absorption and the electron-hole recombination times.

Engineering the electronic bandgaps and band edge positions in carbon-substituted 2D boron nitride: a first-principles investigation.

Modification of graphene to open a robust gap in its electronic spectrum is essential for its use in field effect transistors and photochemistry applications. Inspired by recent experimental success in the preparation of homogeneous alloys of graphene and boron nitride (BN), we consider here engineering the electronic structure and bandgap of C2xB1-xN1-x alloys via both compositional and configurational modification. We start from the BN end-member, which already has a large bandgap, and then show that (a) the bandgap can in principle be reduced to about 2 eV with moderate substitution of C (x < 0.25); and (b) the electronic structure of C2xB1-xN1-x can be further tuned not only with composition x, but also with the configuration adopted by C substituents in the BN matrix. Our analysis, based on accurate screened hybrid functional calculations, provides a clear understanding of the correlation found between the bandgap and the level of aggregation of C atoms: the bandgap decreases most when the C atoms are maximally isolated, and increases with aggregation of C atoms due to the formation of bonding and anti-bonding bands associated with hybridization of occupied and empty defect states. We determine the location of valence and conduction band edges relative to vacuum and discuss the implications on the potential use of 2D C2xB1-xN1-x alloys in photocatalytic applications. Finally, we assess the thermodynamic limitations on the formation of these alloys using a cluster expansion model derived from first-principles.

A DFT study of the structures, stabilities and redox behaviour of the major surfaces of magnetite Fe3O4.

The renewed interest in magnetite (Fe3O4) as a major phase in different types of catalysts has led us to study the oxidation-reduction behaviour of its most prominent surfaces. We have employed computer modelling techniques based on the density functional theory to calculate the geometries and surface free energies of a number of surfaces at different compositions, including the stoichiometric plane, and those with a deficiency or excess of oxygen atoms. The most stable surfaces are the (001) and (111), leading to a cubic Fe3O4 crystal morphology with truncated corners under equilibrium conditions. The scanning tunnelling microscopy images of the different terminations of the (001) and (111) stoichiometric surfaces were calculated and compared with previous reports. Under reducing conditions, the creation of oxygen vacancies in the surface leads to the formation of reduced Fe species in the surface in the vicinity of the vacant oxygen. The (001) surface is slightly more prone to reduction than the (111), due to the higher stabilisation upon relaxation of the atoms around the oxygen vacancy, but molecular oxygen adsorbs preferentially at the (111) surface. In both oxidized surfaces, the oxygen atoms are located on bridge positions between two surface iron atoms, from which they attract electron density. The oxidised state is thermodynamically favourable with respect to the stoichiometric surfaces under ambient conditions, although not under the conditions when bulk Fe3O4 is thermodynamically stable with respect to Fe2O3. This finding is important in the interpretation of the catalytic properties of Fe3O4 due to the presence of oxidised species under experimental conditions.

Theoretical Investigation of the Lattice Thermal Conductivities of II-IV-V2 Pnictide Semiconductors.

Ternary pnictide semiconductors with II-IV-V2 stoichiometry hold potential as cost-effective thermoelectric materials with suitable electronic transport properties, but their lattice thermal conductivities (κ) are typically too high. Insights into their vibrational properties are therefore crucial to finding strategies to reduce κ and achieve improved thermoelectric performance. We present a theoretical exploration of the lattice thermal conductivities for a set of pnictide semiconductors with ABX2 composition (A = Zn, Cd; B = Si, Ge, Sn; and X = P, As) using machine-learning-based regression algorithms to extract force constants from a reduced number of density functional theory simulations and then solving the Boltzmann transport equation for phonons. Our results align well with available experimental data, decreasing the mean absolute error by ∼3 W m-1 K-1 with respect to the best previous set of theoretical predictions. Zn-based ternary pnictides have, on average, more than double the thermal conductivity of the Cd-based compounds. Anisotropic behavior increases with the mass difference between A and B cations, but while the nature of the anion does not affect the structural anisotropy, the thermal conductivity anisotropy is typically higher for arsenides than for phosphides. We identify compounds such as CdGeAs2, for which nanostructuring to an affordable range of particle sizes could lead to κ values low enough for thermoelectric applications.

Anion-Dependent Strength Scale of Interactions in Ionic Liquids from X-ray Photoelectron Spectroscopy, Ab Initio Molecular Dynamics, and Density Functional Theory.

Using a combination of experiments and calculations, we have gained new insights into the nature of anion-cation interactions in ionic liquids (ILs). An X-ray photoelectron spectroscopy (XPS)-derived anion-dependent electrostatic interaction strength scale, determined using XPS core-level binding energies for IL cations, is presented here for 39 different anions, with at least 18 new anions included. Linear correlations of experimental XPS core-level binding energies for IL cations with (a) calculated core binding energies (ab initio molecular dynamics (AIMD) simulations were used to generate high-quality model IL structures followed by single-point density functional theory (DFT) to obtain calculated core binding energies), (b) experimental XPS core-level binding energies for IL anions, and (c) other anion-dependent interaction strength scales led to three main conclusions. First, the effect of different anions on the cation can be related to ground-state interactions. Second, the variations of anion-dependent interactions with the identity of the anion are best rationalized in terms of electrostatic interactions and not occupied valence state/unoccupied valence state interactions or polarizability-driven interactions. Therefore, the XPS-derived anion-dependent interaction strength scale can be explained using a simple electrostatic model based on electrostatic site potentials. Third, anion-probe interactions, irrespective of the identity of the probe, are primarily electrostatic, meaning that our electrostatic interaction strength scale captures some inherent, intrinsic property of anions independent of the probe used to measure the interaction strength scale.

A computational study of magnesium incorporation in the bulk and surfaces of hydroxyapatite.

We have used a combination of static lattice energy minimization and molecular dynamics simulations to investigate the thermodynamics of Mg incorporation into the bulk and hydrated surfaces of hydroxyapatite (HA). In agreement with recent experimental and theoretical work, our simulations show that the incorporation of low levels of Mg in the Ca (II) site is preferred with respect to incorporation in Ca (I) sites. However, we predict that Mg in the HA bulk material is metastable both with respect to the Mg/Ca exchange with aqueous solution and with respect to separation into bulk phases of magnesium phosphate [Mg3(PO4)2] and magnesium hydroxide [Mg(OH)2]. This finding suggests that Mg siting in the HA bulk is at least partially controlled by kinetics rather than by thermodynamics during crystal growth, which can explain the discrepancies found in the literature about the preferential substitution site. Finally, we found that Mg incorporation from solution into the hydrated surfaces, rather than the bulk material, is energetically favorable, in particular in the (011̅0) plane where cation sites are exposed to solution, thereby enabling the favorable interaction of Mg with water.

Hydrogen Bonds in Lead Halide Perovskites: Insights from Ab Initio Molecular Dynamics.

Hydrogen bonds (HBs) play an important role in the rotational dynamics of organic cations in hybrid organic/inorganic halide perovskites, thus affecting the structural and electronic properties of the perovskites. However, the properties and even the existence of HBs in these perovskites are not well established. In this study, we investigate HBs in perovskites MAPbBr3 (MA+ = CH3NH3+), FAPbI3 (FA+ = CH(NH2)2+), and their solid solution with composition (FAPbI3)7/8(MAPbBr3)1/8, using ab initio molecular dynamics and electronic structure calculations. We consider HBs donated by X-H fragments (X = N and C) of the organic cations and accepted by the halides (Y = Br and I) and characterize their properties based on pair distribution functions and on a combined distribution function of the hydrogen-acceptor distance with the donor-hydrogen-acceptor angle. By analyzing these functions, we establish geometrical criteria for HB existence based on the hydrogen-acceptor (H-Y) distance and donor-hydrogen-acceptor angle (X-H-Y). The distance condition is defined as d(H - Y) < 3 Å for N-H-donated HBs and d(H - Y) < 4 Å for C-H-donated HBs. The angular condition is 135° < (X - H - Y) < 180° for both types of HBs. A HB is considered to be formed when both angular and distance conditions are simultaneously satisfied. At the simulated temperature (350 K), the HBs dynamically break and form. We compute the time correlation functions of HB existence and HB lifetimes, which range between 0.1 and 0.3 ps at that temperature. The analysis of HB lifetimes indicates that N-H-Br bonds are relatively stronger than N-H-I bonds, while C-H-Y bonds are weaker, with a minimal influence from the halide and cation. To evaluate the impact of HBs on the vibrational spectra, we present the power spectrum in the region of N-H and C-H stretching modes, comparing them with the normal mode frequencies of isolated cations. We show that the peaks associated with N-H stretching modes in perovskites are redshifted and asymmetrically deformed, while the C-H peaks do not exhibit these effects.

Mg/Ca partitioning between aqueous solution and aragonite mineral: a molecular dynamics study.

We have calculated the concentrations of Mg in the bulk and surfaces of aragonite CaCO(3) in equilibrium with aqueous solution, based on molecular dynamics simulations and grand-canonical statistical mechanics. Mg is incorporated in the surfaces, in particular in the (001) terraces, rather than in the bulk of aragonite particles. However, the total Mg content in the bulk and surface of aragonite particles was found to be too small to account for the measured Mg/Ca ratios in corals. We therefore argue that most Mg in corals is either highly metastable in the aragonite lattice, or is located outside the aragonite phase of the coral skeleton, and we discuss the implications of this finding for Mg/Ca paleothermometry.

Density functional theory study of rutile VO2 surfaces.

We present the results of a density functional theory investigation of the surfaces of rutile-like vanadium dioxide, VO(2)(R). We calculate the surface energies of low Miller index planes and find that the most stable surface orientation is the (110). The equilibrium morphology of a VO(2)(R) particle has an acicular shape, laterally confined by (110) planes and topped by (011) planes. The redox properties of the (110) surface are investigated by calculating the relative surface free energies of the non-stoichiometric compositions as a function of oxygen chemical potential. It is found that the VO(2)(110) surface is oxidized with respect to the stoichiometric composition, not only at ambient conditions but also at the more reducing conditions under which bulk VO(2) is stable in comparison with bulk V(2)O(5). The adsorbed oxygen forms surface vanadyl species much more favorably than surface peroxo species.

Tin-Substituted Chalcopyrite: An n-Type Sulfide with Enhanced Thermoelectric Performance.

The dearth of n-type sulfides with thermoelectric performance comparable to that of their p-type analogues presents a problem in the fabrication of all-sulfide devices. Chalcopyrite (CuFeS2) offers a rare example of an n-type sulfide. Chemical substitution has been used to enhance the thermoelectric performance of chalcopyrite through preparation of Cu1-x Sn x FeS2 (0 ≤ x ≤ 0.1). Substitution induces a high level of mass and strain field fluctuation, leading to lattice softening and enhanced point-defect scattering. Together with dislocations and twinning identified by transmission electron microscopy, this provides a mechanism for scattering phonons with a wide range of mean free paths. Substituted materials retain a large density-of-states effective mass and, hence, a high Seebeck coefficient. Combined with a high charge-carrier mobility and, thus, high electrical conductivity, a 3-fold improvement in power factor is achieved. Density functional theory (DFT) calculations reveal that substitution leads to the creation of small polarons, involving localized Fe2+ states, as confirmed by X-ray photoelectron spectroscopy. Small polaron formation limits the increase in carrier concentration to values that are lower than expected on electron-counting grounds. An improved power factor, coupled with substantial reductions (up to 40%) in lattice thermal conductivity, increases the maximum figure-of-merit by 300%, to zT ≈ 0.3 at 673 K for Cu0.96Sn0.04FeS2.

Understanding X-ray Photoelectron Spectra of Ionic Liquids: Experiments and Simulations of 1-Butyl-3-methylimidazolium Thiocyanate.

We demonstrate a combined experimental and computational approach to probe the electronic structure and atomic environment of an ionic liquid, based on core level binding energies. The 1-butyl-3-methylimidazolium thiocyanate [C4C1Im][SCN] ionic liquid was studied using ab initio molecular dynamics, and results were compared against previously published and new experimental X-ray photoelectron spectroscopy (XPS) data. The long-held assumption that initial-state effects in XPS dominate the measured binding energies is proven correct, which validates the established premise that the ground state electronic structure of the ionic liquid can be inferred directly from XPS measurements. A regression model based upon site electrostatic potentials and intramolecular bond lengths is shown to account accurately for variations in core-level binding energies within the ionic liquid, demonstrating the important effect of long-range interactions on the core levels and throwing into question the validity of traditional single ion pair ionic liquid calculations for interpreting XPS data.