The Blackett Memorial Lecture, 1991: Chemical Insights into High-Temperature Superconductors

The high-temperature superconductors are complex oxides, generally containing two-dimensional CuO2 sheets. Various families of the cuprate superconductors are described, paying special attention to aspects related to oxygen stoichiometry, phase stability, synthesis and chemical manipulation of charge carriers. Other aspects discussed are chemical applications of cuprates, possibly as gas sensors and copper-free oxide superconductors. All but the substituted Nd and Pr cuprates are hole-superconductors. Several families of cuprates show a nearly constant n(h) at maximum T(c). Besides this universality, the cuprates exhibit a number of striking common features. Based on Cu(2p) photoemission studies, it is found that the Cu-O charge-transfer energy, DELTA, and the Cu(3d)-O(2p) hybridization strength, t(pd), are key factors in the superconductivity of cuprates. The relative intensity of the satellite in the Cu(2p) core-level spectra, the polarizability of the CuO2 sheets as well as the hole concentration are related to DELTA/t(pd). These chemical bonding factors have to be explicitly taken into account in any model for superconductivity of the cuprates.

A simple method for the preparation of plasma-sprayable powders based on ZrO2

Plasma-sprayable powders of calcia, magnesia and yttria-stabilized zirconia have been prepared by using polyvinyl alcohol binders. The powders have been characterized for sprayability by spray coating on steer plates previously coated with an NiAl bond coat. The suitability of these coatings for thermal barrier applications have been examined. Thermal barrier and related properties, along with phase stability and mechanical properties, have been found to be good. Failure of the thermal barrier coating has been observed to occur at the interface between the bond coat and the substrate, due to the formation of a pile-up layer consisting of Fe-Zr-Al-O compound.

Thermodynamic modeling to analyse composition of carbonaceous coatings of MnO and other oxides of manganese grown by MOCVD

Equilibrium thermodynamic analysis has been applied to the low-pressure MOCVD process using manganese acetylacetonate as the precursor. ``CVD phase stability diagrams'' have been constructed separately for the processes carried out in argon and oxygen ambient, depicting the compositions of the resulting films as functions of CVD parameters. For the process conduced in argon ambient, the analysis predicts the simultaneous deposition of MnO and elemental carbon in 1: 3 molar proportion, over a range of temperatures. The analysis predicts also that, if CVD is carried out in oxygen ambient, even a very low flow of oxygen leads to the complete absence of carbon in the film deposited oxygen, with greater oxygen flow resulting in the simultaneous deposition of two different manganese oxides under certain conditions. The results of thermodynamic modeling have been verified quantitatively for low-pressure CVD conducted in argon ambient. Indeed, the large excess of carbon in the deposit is found to constitute a MnO/C nanocomposite, the associated cauliflower-like morphology making it a promising candidate for electrode material in supercapacitors. CVD carried out in oxygen flow, under specific conditions, leads to the deposition of more than one manganese oxide, as expected from thermodynamic analysis ( and forming an oxide-oxide nanocomposite). These results together demonstrate that thermodynamic analysis of the MOCVD process can be employed to synthesize thin films in a predictive manner, thus avoiding the inefficient trial-and-error method usually associated with MOCVD process development. The prospect of developing thin films of novel compositions and characteristics in a predictive manner, through the appropriate choice of CVD precursors and process conditions, emerges from the present work.

Process-kinetics in facing targets sputtering of multi-component oxide thin films

The kinetics of the processes in facing targets sputtering of multicomponent oxide films is presented. The novel configuration of the process exhibits an enhanced ionization efficiency. Discharge diagnostics performed using optical emission spectroscopy revealed strong dependence of plasma parameters on process conditions. Numerical simulation based on thermalization and diffusion of sputtered atoms has been performed to estimate the transport efficiency in off-axis mode. Composition, structure and epitaxial quality of YBa2Cu3O7-x films prepared was found to be strongly dependent on atomic flux ratios (of Cu/Y and Ba/Y) arriving at the substrate, resputtering effect and phase stability of YBa2Cu3O7-x These studies have been shown to be useful in understanding the complex processes that occur in sputtering of multicomponent films. (C) 1999 Elsevier Science S.A. All rights reserved.

Nano-size and miscibility gap

Reports on the alloys formed from immiscible atoms when they are contained in a nano-sized system have initiated several research activities in the recent years. Bridging of the miscibility gap at nanoscale is significant as it has the potential to produce novel alloy materials with useful technological applications. Although the literature contains noticeable number of reports on the formation of solid solution between bulk immiscible atoms, several issues related to phase stability and microstructure remain unaddressed. This article discusses some of these issues using examples from the work done by the author's research group on isolated nanoparticles of bulk immiscible binary systems such as Ag-Ni, Ag-Fe and Ag-Co.

Computational design of model Re/Ru bearing Ni-base superalloys

It is well established that Re and Ru additions to Ni-base superalloys result in improved creep performance and phase stability. However, the role of Re and Ru and their synergetic effects are not well understood, and the first step in understanding these effects is to design alloys with controlled microstructural parameters. A computational approach was undertaken in the present work for designing model alloys with varying levels of Re and Ru. Thermodynamic and first principles calculations were employed complimentarily to design a set of alloys with varying Re and Ru levels, but which were constrained by constant microstructural parameters, i.e., phase fractions and lattice misfit across the alloys. Three ternary/quaternary alloys of type Ni-Al-xRe-yRu were thus designed. These compositions were subsequently cast, homogenized and aged. Experimental results suggest that while the measured volume fraction matches the predicted value in the Ru containing alloy, volume fraction is significantly higher than the designed value in the Re containing alloys. This is possibly due to errors in the thermodynamic database used to predict phase fraction and composition. These errors are also reflected in the mismatch between predicted and measured values of misfit.

Growth stress induced tunability of dielectric permittivity in thin films

Stress is inevitable during thin film growth. It is demonstrated here that the growth stress has a significant effect on the dielectric constant of high-k thin films. ZrO2 thin films were deposited on Ge by reactive direct current sputtering. Stress in these films was measured using in-situ curvature measurement tool. The growth stress was tuned from -2.8 to 0.1 GPa by controlling deposition rate. Dielectric permittivity of ZrO2 depends on temperature, phase, and stress. The correct combination of parameters-phase, texture, and stress-is shown to yield films with an equivalent oxide thickness of 8 angstrom. Growth stresses are shown to affect the dielectric constant both directly by affecting lattice parameter and indirectly through the effect on phase stability of ZrO2. (c) 2016 AIP Publishing LLC.

First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

Polaron hopping in olivine phosphates studied by nuclear resonant scattering

Valence fluctuations of Fe²⁺ and Fe³⁺ were studied in a solid solution of Li_xFePO₄ by nuclear resonant forward scattering of synchrotron x rays while the sample was heated in a diamond-anvil pressure cell. The spectra acquired at different temperatures and pressures were analyzed for the frequencies of valence changes using the Blume-Tjon model of a system with a fluctuating Hamil- tonian. These frequencies were analyzed to obtain activation energies and an activation volume for polaron hopping. There was a large suppression of hopping frequency with pressure, giving an anomalously large activation volume. This large, positive value is typical of ion diffusion, which indicates correlated motions of polarons, and Li⁺ ions that alter the dynamics of both.

In a parallel study of Na_xFePO₄, the interplay between sodium ordering and electron mobility was investigated using a combination of synchrotron x-ray diffraction and nuclear resonant scattering. Conventional Mossbauer spectra were collected while the sample was heated in a resistive furnace. An analysis of the temperature evolution of the spectral shapes was used to identify the onset of fast electron hopping and determine the polaron hopping rate. Synchrotron x-ray diffraction measurements were carried out in the same temperature range. Reitveld analysis of the diffraction patterns was used to determine the temperature of sodium redistribution on the lattice. The diffraction analysis also provides new information about the phase stability of the system. The temperature evolution of the iron site occupancies from the Mossbauer measurements, combined with the synchrotron diffraction results give strong evidence for a relationship between the onset of fast electron dynamics and the redistribution of sodium in the lattice.

Measurements of activation barriers for polaron hopping gave fundamental insights about the correlation between electronic carriers and mobile ions. This work established that polaron-ion interactions can alter the local dynamics of electron and ion transport. These types of coupled processes may be common in many materials used for battery electrodes, and new details concerning the influence of polaron-ion interactions on the charge dynamics are relevant to optimizing their electrochemical performance.

Análise global da estabilidade termodinâmica de misturas: um estudo com o método do conjunto gerador

O cálculo do equilíbrio de fases é um problema de grande importância em processos da engenharia, como, por exemplo, na separação por destilação, em processos de extração e simulação da recuperação terciária de petróleo, entre outros. Mas para resolvê-lo é aconselhável que se estude a priori a estabilidade termodinâmica do sistema, a qual consiste em determinar se uma dada mistura se apresenta em uma ou mais fases. Tal problema pode ser abordado como um problema de otimização, conhecido como a minimização da função distância do plano tangente à energia livre de Gibbs molar, onde modelos termodinâmicos, de natureza não convexa e não linear, são utilizados para descrevê-lo. Esse fato tem motivado um grande interesse em técnicas de otimização robustas e eficientes para a resolução de problemas relacionados com a termodinâmica do equilíbrio de fases. Como tem sido ressaltado na literatura, para proporcionar uma completa predição do equilíbrio de fases, faz-se necessário não apenas a determinação do minimizador global da função objetivo do teste de estabilidade, mas também a obtenção de todos os seus pontos estacionários. Assim, o desenvolvimento de metodologias para essa tarefa desafiadora tem se tornado uma nova área de pesquisa da otimização global aplicada à termodinâmica do equilíbrio, com interesses comuns na engenharia química e na engenharia do petróleo. O foco do presente trabalho é uma nova metodologia para resolver o problema do teste de estabilidade. Para isso, usa-se o chamado método do conjunto gerador para realizar buscas do tipo local em uma rede de pontos previamente gerada por buscas globais efetuadas com uma metaheurística populacional, no caso o método do enxame de partículas.Para se obter mais de um ponto estacionário, minimizam-se funções de mérito polarizadas, cujos pólos são os pontos previamente encontrados. A metodologia proposta foi testada na análise de quatorze misturas polares previamente consideradas na literatura. Os resultados mostraram que o método proposto é robusto e eficiente a ponto de encontrar, além do minimizador global, todos os pontos estacionários apontados previamente na literatura, sendo também capaz de detectar, em duas misturas ternárias estudadas, pontos estacionários não obtidos pelo chamado método de análise intervalar, uma técnica confiável e muito difundida na literatura. A análise do teste de estabilidade pela simples utilização do método do enxame de partículas associado à técnica de polarização mencionada acima, para a obtenção de mais de um ponto estacionário (sem a busca local feita pelo método do conjunto gerador em uma dada rede de pontos), constitui outra metodologia para a resolução do problema de interesse. Essa utilização é uma novidade secundária deste trabalho. Tal metodologia simplificada exibiu também uma grande robustez, sendo capaz de encontrar todos os pontos estacionários pesquisados. No entanto, quando comparada com a abordagem mais geral proposta aqui, observou-se que tal simplificação pode, em alguns casos onde a função de mérito apresenta uma geometria mais complexa, consumir um tempo de máquina relativamente grande, dessa forma é menos eficiente.

Formulação teórica dos fundamentos da otimização global topográfica com análise de desempenho e aplicações à estabilidade de fases de misturas termodinâmicas

Métodos de otimização que utilizam condições de otimalidade de primeira e/ou segunda ordem são conhecidos por serem eficientes. Comumente, esses métodos iterativos são desenvolvidos e analisados à luz da análise matemática do espaço euclidiano n-dimensional, cuja natureza é de caráter local. Consequentemente, esses métodos levam a algoritmos iterativos que executam apenas as buscas locais. Assim, a aplicação de tais algoritmos para o cálculo de minimizadores globais de uma função não linear,especialmente não-convexas e multimodais, depende fortemente da localização dos pontos de partida. O método de Otimização Global Topográfico é um algoritmo de agrupamento, que utiliza uma abordagem baseada em conceitos elementares da teoria dos grafos, a fim de gerar bons pontos de partida para os métodos de busca local, a partir de pontos distribuídos de modo uniforme no interior da região viável. Este trabalho tem dois objetivos. O primeiro é realizar uma nova abordagem sobre método de Otimização Global Topográfica, onde, pela primeira vez, seus fundamentos são formalmente descritos e suas propriedades básicas são matematicamente comprovadas. Neste contexto, propõe-se uma fórmula semi-empírica para calcular o parâmetro chave deste algoritmo de agrupamento, e, usando um método robusto e eficiente de direções viáveis por pontos-interiores, estendemos o uso do método de Otimização Global Topográfica a problemas com restrições de desigualdade. O segundo objetivo é a aplicação deste método para a análise de estabilidade de fase em misturas termodinâmicas,o qual consiste em determinar se uma dada mistura se apresenta em uma ou mais fases. A solução deste problema de otimização global é necessária para o cálculo do equilíbrio de fases, que é um problema de grande importância em processos da engenharia, como, por exemplo, na separação por destilação, em processos de extração e simulação da recuperação terciária de petróleo, entre outros. Além disso, afim de ter uma avaliação inicial do potencial dessa técnica, primeiro vamos resolver 70 problemas testes, e então comparar o desempenho do método proposto aqui com o solver MIDACO, um poderoso software recentemente introduzido no campo da otimização global.

Double-ceramic-layer thermal barrier coatings of La2Zr2O7/YSZ deposited by electron beam-physical vapor deposition

Double-ceramic-layer(DCL) thermal barrier coatings (TBCs) of La2Zr2O7 (LZ) and yttria stabilized zirconia (YSZ) were deposited by electron beam-physical vapor deposition (EB-PVD). The composition, crystal structure, surface and cross-sectional morphologies and cyclic oxidation behavior of the DCL coating were studied. Both the X-ray diffraction (XRD) and thermogravimetric-differential thermal analysis (TG-DTA) prove that LZ and YSZ have good chemical applicability to form a DCL coating. The thermal cycling test at 1373 K in an air furnace indicates the DCL coating has a much longer lifetime than the single layer LZ coating. and even longer than that of the single layer YSZ coating. The failure of the DCL coating is a result of both the bond coat oxidation and the thermal strain between bond coat and ceramic layer generated by the thermal expansion mismatch.

Preparation and characterization of La2Zr2O7 coating with the addition of Y2O3 by EB-PVD

Thermal barrier coatings (TBCs) of La2Zr2O7 (LZ) with the addition of 3 wt.% Y2O3 (LZ3Y) were deposited by electron beam-physical vapor deposition (EB-PVD). The phase stabilities, thermophysical and mechanical properties, and chemical compositions of these ceramics and coatings were studied in detail. The phase stability and thermal expansion behavior of LZ3Y bulk material are identical to those of LZ bulk material, but the mechanical properties of the former are superior to those of the latter. Elemental analysis and X-ray diffraction indicate that compositional deviation of LZ coating can be optimized after doping by 3 wt.% Y2O3, Y2O3 acts as a dopant as well as a process regulator. The optimal composition of LZ3Y coating could be effectively achieved by the addition of excess Y2O3 into the ingot and by properly controlling the current of electron beam (i.e. similar to 650 mA).

Thermal barrier coating of lanthanum-zirconium-cerium composite oxide made by electron beam-physical vapor deposition

Lanthanum-zirconium-cerium composite oxide (La-2(Zr0.7Ce0.3)(2)O-7, LZ7C3) as a candidate material for thermal barrier coatings (TBCs) was prepared by electron beam-physical vapor deposition (EB-PVD). The composition, crystal structure, thermophysical properties, surface and cross-sectional morphologies and cyclic oxidation behavior of the LZ7C3 coating were studied. The results indicated that LZ7C3 has a high phase stability between 298 K and 1573 K, and its linear thermal expansion coefficient (TEC) is similar to that of zirconia containing 8 wt% yttria (8YSZ). The thermal conductivity of LZ7C3 is 0.87 W m(-1) K-1 at 1273 K, which is almost 60% lower than that of 8YSZ. The deviation of coating composition from the ingot can be overcome by the addition of excess CeO2 and ZrO2 during ingot preparation or by adjusting the process parameters.

Studies on the doubly charged cluster ions of sodium and potassium nitrates by electrospray ionization tandem mass spectrometry

Doubly charged cluster ions, besides singly charged cluster ions, from sodium and potassium nitrates were produced evidently under normal source capillary temperature of 200 degrees C in both positive and negative ion electrospray ionization (ESI) ion trap mass spectrometry. The fragmentation pathways for doubly charged cluster ions were studied in detail using ESI tandem mass spectrometry and two pathways were observed depending on the cluster sizes of alkali metal nitrates. In addition, factors that affect the formation of cluster ions were also interrogated.