Relationship between phonons and thermal expansion in Zn(CN)2 and Ni(CN)2 from inelastic neutron scattering and ab initio calculations

Zn(CN)2 and Ni(CN)2 are known for exhibiting anomalous thermal expansion over a wide temperature range. The volume thermal expansion coefficient for the cubic, three dimensionally connected material, Zn(CN)2, is negative (alpha(V) = −51  10(-6) K-1) while for Ni(CN)2, a tetragonal material, the thermal expansion coefficient is negative in the two dimensionally connected sheets (alpha(a) = −7  10(-6) K-1), but the overall thermal expansion coefficient is positive (alpha(V) = 48  10(-6) K-1). We have measured the temperature dependence of phonon spectra in these compounds and analyzed them using ab initio calculations. The spectra of the two compounds show large differences that cannot be explained by simple mass renormalization of the modes involving Zn (65.38 amu) and Ni (58.69 amu) atoms. This reflects the fact that the structure and bonding are quite different in the two compounds. The calculated pressure dependence of the phonon modes and of the thermal expansion coefficient, alpha(V), are used to understand the anomalous behavior in these compounds. Our ab initio calculations indicate that phonon modes of energy approx. 2 meV are major contributors to negative thermal expansion (NTE) in both the compounds. The low-energy modes of approx.8 and 13 meV in Zn(CN)2 also contribute significantly to the NTE in Zn(CN)2 and Ni(CN)2, respectively. The measured temperature dependence of the phonon spectra has been used to estimate the total anharmonicity of both compounds. For Zn(CN)2, the temperature-dependent measurements (total anharmonicity), along with our previously reported pressure dependence of the phonon spectra (quasiharmonic), is used to separate the explicit temperature effect at constant volume (intrinsic anharmonicity).

Symmetry-break in Voronoi tessellations

We analyse in a common framework the properties of the Voronoi tessellations resulting from regular 2D and 3D crystals and those of tessellations generated by Poisson distributions of points, thus joining on symmetry breaking processes and the approach to uniform random distributions of seeds. We perturb crystalline structures in 2D and 3D with a spatial Gaussian noise whose adimensional strength is α and analyse the statistical properties of the cells of the resulting Voronoi tessellations using an ensemble approach. In 2D we consider triangular, square and hexagonal regular lattices, resulting into hexagonal, square and triangular tessellations, respectively. In 3D we consider the simple cubic (SC), body-centred cubic (BCC), and face-centred cubic (FCC) crystals, whose corresponding Voronoi cells are the cube, the truncated octahedron, and the rhombic dodecahedron, respectively. In 2D, for all values α>0, hexagons constitute the most common class of cells. Noise destroys the triangular and square tessellations, which are structurally unstable, as their topological properties are discontinuous in α=0. On the contrary, the honeycomb hexagonal tessellation is topologically stable and, experimentally, all Voronoi cells are hexagonal for small but finite noise with α<0.12. Basically, the same happens in the 3D case, where only the tessellation of the BCC crystal is topologically stable even against noise of small but finite intensity. In both 2D and 3D cases, already for a moderate amount of Gaussian noise (α>0.5), memory of the specific initial unperturbed state is lost, because the statistical properties of the three perturbed regular tessellations are indistinguishable. When α>2, results converge to those of Poisson-Voronoi tessellations. In 2D, while the isoperimetric ratio increases with noise for the perturbed hexagonal tessellation, for the perturbed triangular and square tessellations it is optimised for specific value of noise intensity. The same applies in 3D, where noise degrades the isoperimetric ratio for perturbed FCC and BCC lattices, whereas the opposite holds for perturbed SCC lattices. This allows for formulating a weaker form of the Kelvin conjecture. By analysing jointly the statistical properties of the area and of the volume of the cells, we discover that also the cells shape heavily fluctuates when noise is introduced in the system. In 2D, the geometrical properties of n-sided cells change with α until the Poisson-Voronoi limit is reached for α>2; in this limit the Desch law for perimeters is shown to be not valid and a square root dependence on n is established, which agrees with exact asymptotic results. Anomalous scaling relations are observed between the perimeter and the area in the 2D and between the areas and the volumes of the cells in 3D: except for the hexagonal (2D) and FCC structure (3D), this applies also for infinitesimal noise. In the Poisson-Voronoi limit, the anomalous exponent is about 0.17 in both the 2D and 3D case. A positive anomaly in the scaling indicates that large cells preferentially feature large isoperimetric quotients. As the number of faces is strongly correlated with the sphericity (cells with more faces are bulkier), in 3D it is shown that the anomalous scaling is heavily reduced when we perform power law fits separately on cells with a specific number of faces.

Optimizing the methods of synthesis for barium hexagonal ferrite: an experimental and theoretical study

M-type barium hexaferrite (BaM) is a hard ferrite, crystallizing in space group P6(3)/mmc possessing a hexagonal magneto-plumbite structure, which consists of alternate hexagonal and spinel blocks. The structure of BaM is thus related to those of garnet and spinel ferrite. However the material has proved difficult to synthesize. By taking into account the presence of the spinel block in barium hexagonal ferrite, highly efficient new synthetic methods were devised with routes significantly different from existing ones. These successful variations in synthetic methods have been derived by taking into account a detailed investigation of the structural features of barium hexagonal ferrite and the least change principle whereby configuration changes are kept to a minimum. Thus considering the relevant mechanisms has helped to improve the synthesis efficiencies for both hydrothermal and co-precipitation methods by choosing conditions that invoke the formation of the cubic block or the less stable Fe3O4. The role played by BaFe2O4 in the synthesis is also discussed. The distribution of iron from reactants or intermediates among different sites was also successfully explained. The proposed mechanisms are based on the principle that the cubic block must be self-assembled to form the final product. Thus, it is believed that these formulated mechanisms should be helpful in designing experiments to obtain a deeper understanding of the synthesis process and to investigate the substitution of magnetic ions with doping ions.

Determination of tube theory parameters using a simple grid model as an example

Although the tube theory is successful in describing entangled polymers qualitatively, a more quantitative description requires precise and consistent definitions of its parameters. Here we investigate the simplest model of entangled polymers, namely a single Rouse chain in a cubic lattice of line obstacles, and illustrate the typical problems and uncertainties of the tube theory. In particular we show that in general one needs 3 entanglement related parameters, but only 2 combinations of them are relevant for the long-time dynamics. Conversely, the plateau modulus can not be determined from these two parameters and requires a more detailed model of entanglements with explicit entanglement forces, such as the slipsprings model. It is shown that for the grid model the Rouse time within the tube is larger than the Rouse time of the free chain, in contrast to what the standard tube theory assumes.

Cyclic prefixed OQAM-OFDM and its application to single-carrier FDMA

Single-carrier frequency division multiple access (SC-FDMA) has appeared to be a promising technique for high data rate uplink communications. Aimed at SC-FDMA applications, a cyclic prefixed version of the offset quadrature amplitude modulation based OFDM (OQAM-OFDM) is first proposed in this paper. We show that cyclic prefixed OQAMOFDM CP-OQAM-OFDM) can be realized within the framework of the standard OFDM system, and perfect recovery condition in the ideal channel is derived. We then apply CP-OQAMOFDM to SC-FDMA transmission in frequency selective fading channels. Signal model and joint widely linear minimum mean square error (WLMMSE) equalization using a prior information with low complexity are developed. Compared with the existing DFTS-OFDM based SC-FDMA, the proposed SC-FDMA can significantly reduce envelope fluctuation (EF) of the transmitted signal while maintaining the bandwidth efficiency. The inherent structure of CP-OQAM-OFDM enables low-complexity joint equalization in the frequency domain to combat both the multiple access interference and the intersymbol interference. The joint WLMMSE equalization using a prior information guarantees optimal MMSE performance and supports Turbo receiver for improved bit error rate (BER) performance. Simulation resultsconfirm the effectiveness of the proposed SC-FDMA in termsof EF (including peak-to-average power ratio, instantaneous-toaverage power ratio and cubic metric) and BER performances.

Approximation by harmonic polynomials in star-shaped domains and exponential convergence of Trefftz hp-dGFEM

We study the approximation of harmonic functions by means of harmonic polynomials in two-dimensional, bounded, star-shaped domains. Assuming that the functions possess analytic extensions to a delta-neighbourhood of the domain, we prove exponential convergence of the approximation error with respect to the degree of the approximating harmonic polynomial. All the constants appearing in the bounds are explicit and depend only on the shape-regularity of the domain and on delta. We apply the obtained estimates to show exponential convergence with rate O(exp(−b square root N)), N being the number of degrees of freedom and b>0, of a hp-dGFEM discretisation of the Laplace equation based on piecewise harmonic polynomials. This result is an improvement over the classical rate O(exp(−b cubic root N )), and is due to the use of harmonic polynomial spaces, as opposed to complete polynomial spaces.

Thiocyanate complexes of uranium in multiple oxidation states: a combined structural, magnetic, spectroscopic, spectroelectrochemical, and theoretical study

A comprehensive study of the complexes A4[U(NCS)8] (A = Cs, Et4N, nBu4N) and A3[UO2(NCS)5] (A = Cs, Et4N) is described, with the crystal structures of [nBu4N]4[U(NCS)8]·2MeCN and Cs3[UO2(NCS)5]·O0.5 reported. The magnetic properties of square antiprismatic Cs4[U(NCS)8] and cubic [Et4N]4[U(NCS)8] have been probed by SQUID magnetometry. The geometry has an important impact on the low-temperature magnetic moments: at 2 K, μeff = 1.21 μB and 0.53 μB, respectively. Electronic absorption and photoluminescence spectra of the uranium(IV) compounds have been measured. The redox chemistry of [Et4N]4[U(NCS)8] has been explored using IR and UV–vis spectroelectrochemical methods. Reversible 1-electron oxidation of one of the coordinated thiocyanate ligands occurs at +0.22 V vs Fc/Fc+, followed by an irreversible oxidation to form dithiocyanogen (NCS)2 which upon back reduction regenerates thiocyanate anions coordinating to UO22+. NBO calculations agree with the experimental spectra, suggesting that the initial electron loss of [U(NCS)8]4– is delocalized over all NCS– ligands. Reduction of the uranyl(VI) complex [Et4N]3[UO2(NCS)5] to uranyl(V) is accompanied by immediate disproportionation and has only been studied by DFT methods. The bonding in [An(NCS)8]4– (An = Th, U) and [UO2(NCS)5]3– has been explored by a combination of DFT and QTAIM analysis, and the U–N bonds are predominantly ionic, with the uranyl(V) species more ionic that the uranyl(VI) ion. Additionally, the U(IV)–NCS ion is more ionic than what was found for U(IV)–Cl complexes.

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.

Dynamic simulation studies of structural formation and transition in electro-magneto-rheological fluids

The computer simulation method has been used to study the structural formation and transition of electro-magneto-rheological (EMR) fluids under compatible electric and magnetic fields. When the fields are applied simultaneously and perpendicularly to each other, the particles rapidly arrange into two-dimensional close-packed layer structures parallel to both fields. The layers then combine together to form thicker sheet-like structures, which finally relax into three-dimensional close-packed structures with the help of the thermal fluctuations. On the other hand, if the electric field is applied firstly to induce the body-centered tetragonal (BCT) columns in the system, and then the magnetic field is applied in the perpendicular direction. the BCT to face-centered cubic (FCC) structure transition is observed in very short time. Following that. the structure keeps on evolving due to the demagnetization effect and finally form the three-dimensional close-packed structures.

Simulation of field-induced structural formation and transition in electromagnetorheological suspensions

A computer simulation method has been used to study the three-dimensional structural formation and transition of eleetromagnetorheological (EMR) suspensions under compatible electric and magnetic fields. When the fields are applied simultaneously and perpendicularly to each other, the particles rapidly arrange into single layer structures parallel to both fields. In each layer, there is a two-dimensional hexagonal lattice. The single layers then combine together to form thicker sheetlike structures. With the help of the thermal fluctuations, the thicker structures relax into three-dimensional close-packed structures, which may be face-centered cubic (fcc), hexagonal close-packed (hup) lattices, or, more probably, the mixture of them, depending on the initial configurations and the thermal fluctuations. On the other hand, if the electric field is applied first to induce the body-centered tetragonal (bct) columns in the system, and then the magnetic field is applied in the perpendicular direction, the bet to fee structure transition is observed in a very short time. Following that, the structure keeps on evolving due to the demagnetization effect and finally forms close-packed structures with fee and hcp lattice character. The simulation results are in agreement with the theoretical and experimental results.

High-temperature order–disorder transitions in the skutterudites CoGe1.5Q1.5 (Q=S, Te)

The temperature dependence of anion ordering in the skutterudites CoGe1.5Q1.5 (Q=S, Te) has been investigated by powder neutron diffraction. Both materials adopt a rhombohedral structure at room temperature (space group R-3 ) in which the anions are ordered trans to each other within Ge2Q2 rings. In CoGe1.5S1.5, anion ordering is preserved up to the melting point of 950 °C. However, rhombohedral CoGe1.5Te1.5 undergoes a phase transition at 610 °C involving a change to cubic symmetry (space group Im-3). In the high-temperature modification, there is a statistical distribution of anions over the available sites within the Ge2Te2 rings. The structural transition involves a reduction in the degree of distortion of the Ge2Te2 rings which progressively transform from a rhombus to a rectangular shape. The effect of this transition on the thermoelectric properties has been investigated.

Simulation of long-term thermal characteristics of three Estonian lakes

A one-dimensional surface energy-balance lake model, coupled to a thermodynamic model of lake ice, is used to simulate variations in the temperature of and evaporation from three Estonian lakes: Karujärv, Viljandi and Kirjaku. The model is driven by daily climate data, derived by cubic-spline interpolation from monthly mean data, and was run for periods of 8 years (Kirjaku) up to 30 years (Viljandi). Simulated surface water temperature is in good agreement with observations: mean differences between simulated and observed temperatures are from −0.8°C to +0.1°C. The simulated duration of snow and ice cover is comparable with observed. However, the model generally underpredicts ice thickness and overpredicts snow depth. Sensitivity analyses suggest that the model results are robust across a wide range (0.1–2.0 m−1) of lake extinction coefficient: surface temperature differs by less than 0.5°C between extreme values of the extinction coefficient. The model results are more sensitive to snow and ice albedos. However, changing the snow (0.2–0.9) and ice (0.15–0.55) albedos within realistic ranges does not improve the simulations of snow depth and ice thickness. The underestimation of ice thickness is correlated with the overestimation of snow cover, since a thick snow layer insulates the ice and limits ice formation. The overestimation of snow cover results from the assumption that all the simulated winter precipitation occurs as snow, a direct consequence of using daily climate data derived by interpolation from mean monthly data.

Synthesis, characterization and physical properties of the skutterudites YbxFe2Ni2Sb12 (0≤x≤0.4)

The skutterudites YbxFe2Ni2Sb12 (0≤x≤0.4) have been prepared by solid-state reaction and characterised by powder X-ray diffraction. The compounds crystallise in the cubic space group Im View the MathML source3¯ (a≈9.1 Å) with Yb atoms partially filling the voids in the skutterudite framework. A neutron time-of-flight diffraction experiment for Fe2Ni2Sb12 confirms the disorder of Fe and Ni atoms on the transition-metal site. Electrical resistivity, Seebeck coefficient and thermal conductivity measurements indicate that the thermoelectric performance of the skutterudites shows a marked dependence on the Yb content. Magnetic measurements over the temperature range 2≤T/K≤300 show paramagnetic behaviour for all compounds. Decomposition studies under an oxidising atmosphere at elevated temperatures have also been carried out by thermogravimetric analysis.

Na5(Ga4S)(GaS4)3·6H2O: a three-dimensional thiogallate containing a novel octahedral building block

The new thiogallate Na5(Ga4S)(GaS4)3·6H2O has been prepared solvothermally, using 3,5-dimethyl pyridine as a solvent, and characterised by powder and single crystal X-ray diffraction. This material, which exhibits a three-dimensional crystal structure, crystallises in the cubic space group View the MathML sourceF4¯3c (a = 17.557(4) Å). The crystal structure contains octahedral building blocks [Ga4S (GaS4)6]20−, linked into a three-dimensional network with a perovskite-type topology, and sodium hydrate clusters, [Na5(H2O)6]5+, filling the cavities in the [Ga4S(GaS4)6/2]5− framework. UV–Vis diffuse reflectance measurements indicate that this material is a wide band gap semiconductor, with a band gap of ca. 4.4 eV.

The atomic structure of low-index surfaces of the intermetallic compound InPd

The intermetallic compound InPd (CsCl type of crystal structure with a broad compositional range) is considered as a candidate catalyst for the steam reforming of methanol. Single crystals of this phase have been grown to study the structure of its three low-index surfaces under ultra-high vacuum conditions, using low energy electron diffraction (LEED), X-ray photoemission spectroscopy (XPS), and scanning tunneling microscopy (STM). During surface preparation, preferential sputtering leads to a depletion of In within the top few layers for all three surfaces. The near-surface regions remain slightly Pd-rich until annealing to ∼580 K. A transition occurs between 580 and 660 K where In segregates towards the surface and the near-surface regions become slightly In-rich above ∼660 K. This transition is accompanied by a sharpening of LEED patterns and formation of flat step-terrace morphology, as observed by STM. Several superstructures have been identified for the different surfaces associated with this process. Annealing to higher temperatures (≥750 K) leads to faceting via thermal etching as shown for the (110) surface, with a bulk In composition close to the In-rich limit of the existence domain of the cubic phase. The Pd-rich InPd(111) is found to be consistent with a Pd-terminated bulk truncation model as shown by dynamical LEED analysis while, after annealing at higher temperature, the In-rich InPd(111) is consistent with an In-terminated bulk truncation, in agreement with density functional theory (DFT) calculations of the relative surface energies. More complex surface structures are observed for the (100) surface. Additionally, individual grains of a polycrystalline sample are characterized by micro-spot XPS and LEED as well as low-energy electron microscopy. Results from both individual grains and “global” measurements are interpreted based on comparison to our single crystals findings, DFT calculations and previous literature.