Panchapakesan Ganesh

Origin of morphotropic phase boundaries in ferroelectrics.

A piezoelectric material is one that generates a voltage in response to a mechanical strain (and vice versa). The most useful piezoelectric materials display a transition region in their composition phase diagrams, known as a morphotropic phase boundary, where the crystal structure changes abruptly and the electromechanical properties are maximal. As a result, modern piezoelectric materials for technological applications are usually complex, engineered, solid solutions, which complicates their manufacture as well as introducing complexity in the study of the microscopic origins of their properties. Here we show that even a pure compound, in this case lead titanate, can display a morphotropic phase boundary under pressure. The results are consistent with first-principles theoretical predictions, but show a richer phase diagram than anticipated; moreover, the predicted electromechanical coupling at the transition is larger than any known. Our results show that the high electromechanical coupling in solid solutions with lead titanate is due to tuning of the high-pressure morphotropic phase boundary in pure lead titanate to ambient pressure. We also find that complex microstructures or compositions are not necessary to obtain strong piezoelectricity. This opens the door to the possible discovery of high-performance, pure-compound electromechanical materials, which could greatly decrease costs and expand the utility of piezoelectric materials.

Intrinsic thermodynamic and kinetic properties of Sb electrodes for Li-ion and Na-ion batteries: experiment and theory

A detailed comparative study between the electrochemical lithiation and sodiation of pure antimony (Sb), relating changes in structural, thermodynamic, kinetic and electrochemical properties has been carried out. For this purpose, a wide range of measurements using electrochemical (galvanostatic cycling, GITT, PITT), X-ray diffraction (XRD) and X-ray photoelectron spectroscopy (XPS) methods as well as density functional theory (DFT) based investigations have been undertaken. Assessment of the thermodynamics reveals that the reaction proceeds identically during the first and second cycles for Li whereas it differs between the first and subsequent cycles for Na as the reaction with Na proceeds through a different pathway associated with the formation of amorphous NaxSb phases. For the first time we rationalize the amorphization of NaxSb phases by the long ranged strain propagation due to Na-vacancy compared to Li–Sb. At full discharge, our XRD results show for the first time that a minor fraction of hexagonal Li3Sb forms concomitantly with cubic Li3Sb. The XRD results confirm that Sb crystallizes into hexagonal Na3Sb at full sodiation. The kinetics of the reaction is assessed by rate performance tests which highlight that both Li and Na can diffuse rapidly throughout micron thick films at room temperature. However, it is found that the (de)insertion of Li provides lower overpotentials and larger storage capacities compared to Na. The difference in rate performance is complemented by diffusion coefficient determinations near the 0 V region where both materials are crystallized into M3Sb (M = Li, Na). Interestingly, calculations show that the energy barrier for near-neighbor vacancy migration, predominant in these close-packed phases, is about twice for Na than for Li. Our analysis tries to relate the lower intrinsic diffusivity of Na compared to Li with the long-range strain propagation induced by the former, thereby leading to an intrinsic origin of differences in rates, mechanical properties and amorphization. Finally, the surface chemistry of Sb electrodes cycled in NaClO4 dissolved in pure PC with(out) the addition of 5 wt% EC or FEC shows presence of ethers and NaF for the EC- and FEC-based electrolytes, respectively, and SEI films rich in Na-based carbonates.

Signature of nearly icosahedral structures in liquid and supercooled liquid copper

A growing body of experiments display indirect evidence of icosahedral structures in supercooled liquid metals. Computer simulations provide more direct evidence but generally rely on approximate interatomic potentials of unproven accuracy. We use first-principles molecular dynamics simulations to generate realistic atomic configurations, providing structural detail not directly available from experiment, based on interatomic forces that are more reliable than conventional simulations. We analyze liquid copper, for which recent experimental results are available for comparison, to quantify the degree of local icosahedral and polytetrahedral order.

Accurate static and dynamic properties of liquid electrolytes for Li-ion batteries from ab initio molecular dynamics.

Lithium-ion batteries have the potential to revolutionize the transportation industry, as they did for wireless communication. A judicious choice of the liquid electrolytes used in these systems is required to achieve a good balance among high-energy storage, long cycle life and stability, and fast charging. Ethylene-carbonate (EC) and propylene-carbonate (PC) are popular electrolytes. However, to date, almost all molecular-dynamics simulations of these fluids rely on classical force fields, while a complete description of the functionality of Li-ion batteries will eventually require quantum mechanics. We perform accurate ab initio molecular-dynamics simulations of ethylene- and propylene-carbonate with LiPF(6) at experimental concentrations to build solvation models which explain available neutron scattering and nuclear magnetic resonance (NMR) results and to compute Li-ion solvation energies and diffusion constants. Our results suggest some similarities between the two liquids as well as some important differences. Simulations also provide useful insights into formation of solid-electrolyte interphases in the presence of electrodes in conventional Li-ion batteries.

Liquid-liquid transition in supercooled silicon determined by first-principles simulation.

First-principles molecular dynamics simulations reveal a liquid-liquid phase transition in supercooled elemental silicon. Two phases coexist below Tc approximately 1232 K and above pc approximately -12 kB. The low-density phase is nearly tetracoordinated, with a pseudogap at the Fermi surface, while the high-density phase is more highly coordinated and metallic in nature. The transition is observed through the formation of van der Waals loops in pressure-volume isotherms below Tc.

Aqueous proton transfer across single-layer graphene

Proton transfer across single-layer graphene proceeds with large computed energy barriers and is therefore thought to be unfavourable at room temperature unless nanoscale holes or dopants are introduced, or a potential bias is applied. Here we subject single-layer graphene supported on fused silica to cycles of high and low pH, and show that protons transfer reversibly from the aqueous phase through the graphene to the other side where they undergo acid–base chemistry with the silica hydroxyl groups. After ruling out diffusion through macroscopic pinholes, the protons are found to transfer through rare, naturally occurring atomic defects. Computer simulations reveal low energy barriers of 0.61–0.75 eV for aqueous proton transfer across hydroxyl-terminated atomic defects that participate in a Grotthuss-type relay, while pyrylium-like ether terminations shut down proton exchange. Unfavourable energy barriers to helium and hydrogen transfer indicate the process is selective for aqueous protons.

Development of a ReaxFF potential for carbon condensed phases and its application to the thermal fragmentation of a large fullerene

In this article, we report the development of a ReaxFF reactive potential that can accurately describe the chemistry and dynamics of carbon condensed phases. Density functional theory (DFT)-based calculations were performed to obtain the equation of state for graphite and diamond and the formation energies of defects in graphene and amorphous phases from fullerenes. The DFT data were used to reparametrize ReaxFFCHO, resulting in a new potential called ReaxFFC-2013. ReaxFFC-2013 accurately predicts the atomization energy of graphite and closely reproduces the DFT-based energy difference between graphite and diamond, and the barrier for transition from graphite to diamond. ReaxFFC-2013 also accurately predicts the DFT-based energy barrier for Stone-Wales transformation in a C60(Ih) fullerene through the concerted rotation of a C2 unit. Later, MD simulations of a C180 fullerene using ReaxFFC-2013 suggested that the thermal fragmentation of these giant fullerenes is an exponential function of time. An Arrhenius-type equation was fit to the decay rate, giving an activation energy of 7.66 eV for the loss of carbon atoms from the fullerene. Although the decay of the molecule occurs primarily via the loss of C2 units, we observed that, with an increase in temperature, the probability of loss of larger fragments increases. The ReaxFFC-2013 potential developed in this work, and the results obtained on fullerene fragmentation, provide an important step toward the full computational chemical modeling of coal pyrolysis, soot incandescence, high temperature erosion of graphitic rocket nozzles, and ablation of carbon-based spacecraft materials during atmospheric reentry.

Understanding controls on interfacial wetting at epitaxial graphene: Experiment and Theory

The interaction of interfacial water with graphitic carbon at the atomic scale is studied as a function of the hydrophobicity of epitaxial graphene. High resolution x-ray reflectivity shows that the graphene-water contact angle is controlled by the average graphene thickness, due to the fraction of the film surface expressed as the epitaxial buffer layer whose contact angle (contact angle {Theta}{sub c} = 73{sup o}) is substantially smaller than that of multilayer graphene ({Theta}{sub c} = 93{sup o}). Classical and ab initio molecular dynamics simulations show that the reduced contact angle of the buffer layer is due to both its epitaxy with the SiC substrate and the presence of interfacial defects. This insight clarifies the relationship between interfacial water structure and hydrophobicity, in general, and suggests new routes to control interface properties of epitaxial graphene.

Understanding the origin of high-rate intercalation pseudocapacitance in Nb2O5 crystals

Pseudocapacitors aim to maintain the high power density of supercapacitors while increasing the energy density towards those of energy dense storage systems such as lithium ion batteries. Recently discovered intercalation pseudocapacitors (e.g. Nb2O5) are particularly interesting because their performance is seemingly not limited by surface reactions or structures, but instead determined by the bulk crystalline structure of the material. We study ordered polymorphs of Nb2O5 and detail the mechanism for the intrinsic high rates and energy density observed for this class of materials. We find that the intercalating atom (lithium) forms a solid solution adsorbing at specific sites in a network of quasi-2D NbOx faces with x = {1.3, 1.67, or 2}, donating electrons locally to its neighboring atoms, reducing niobium. Open channels in the structure have low diffusion barriers for ions to migrate between these sites (Eb ∼ 0.28–0.44 eV) comparable to high-performance solid electrolytes. Using a combination of complementary theoretical methods we rationalize this effect in LixNb2O5 for a wide range of compositions (x) and at finite temperatures. Multiple adsorption sites per unit-cell with similar adsorption energies and local charge transfer result in high capacity and energy density, while the interconnected open channels lead to low cost diffusion pathways between these sites, resulting in high power density. The nano-porous structure exhibiting local chemistry in a crystalline framework is the origin of high-rate pseudocapacitance in this new class of intercalation pseudocapacitor materials. This new insight provides guidance for improving the performance of this family of materials.

Binding and Diffusion of Lithium in Graphite: Quantum Monte Carlo Benchmarks and Validation of van der Waals Density Functional Methods

Highly accurate diffusion quantum Monte Carlo (QMC) studies of the adsorption and diffusion of atomic lithium in AA-stacked graphite are compared with van der Waals-including density functional theory (DFT) calculations. Predicted QMC lattice constants for pure AA graphite agree with experiment. Pure AA-stacked graphite is shown to challenge many van der Waals methods even when they are accurate for conventional AB graphite. Highest overall DFT accuracy, considering pure AA-stacked graphite as well as lithium binding and diffusion, is obtained by the self-consistent van der Waals functional vdW-DF2, although errors in binding energies remain. Empirical approaches based on point charges such as DFT-D are inaccurate unless the local charge transfer is assessed. The results demonstrate that the lithium-carbon system requires a simultaneous highly accurate description of both charge transfer and van der Waals interactions, favoring self-consistent approaches.