Search

This work introduces a completely rewritten version of the program RMCProfile (version 7), big-box, reverse Monte Carlo modelling software for analysis of total scattering data. The major new feature of RMCProfile7 is the ability to refine multiple phases simultaneously, which is relevant for many current research areas such as energy materials, catalysis and engineering. Other new features include improved support for molecular potentials and rigid-body refinements, as well as multiple different data sets. An empirical resolution correction and calculation of the pair distribution function as a back-Fourier transform are now also available. RMCProfile7 is freely available for download at https://rmcprofile.ornl.gov/., (© Wojciech A. Sławiński et al. 2024.)

Metal-organic framework (MOF) glasses are an emerging class of glasses which complement traditional inorganic, organic and metallic counterparts due to their hybrid nature. Although a few zeolitic imidazolate frameworks have been made into glasses, how to melt and quench the largest subclass of MOFs, metal carboxylate frameworks, into glasses remains challenging. Here, we develop a strategy by grafting the zwitterions on the carboxylate ligands and incorporating organic acids in the framework channels to enable the glass formation. The charge delocalization of zwitterion-acid subsystem and the densely filled channels facilitate the coordination bonding mismatch and thus reduce the melting temperature. Following melt-quenching realizes the glass formation of a family of carboxylate MOFs (UiO-67, UiO-68 and DUT-5), which are usually believed to be un-meltable. Our work opens up an avenue for melt-quenching porous molecular solids into glasses., (© 2024. The Author(s).)

Samarium hexaboride, SmB 6 , is a negative thermal expansion (NTE) material whose structure is similar to other known NTE materials such as the family of Prussian blues. In the Prussian blues, NTE is due to a phonon mechanism, but we recently showed from DFT calculations that this is unlikely in SmB 6 (Li et al. , Phys. Chem. Chem. Phys . 2023, 25 , 10749). We now report experimental X-ray diffraction and pair distribution function analysis of this material in the temperature range 20-300 K. The interatomic distances shown by both methods are consistent with the NTE instead arising from an electronic effect, by which the samarium atoms lose electrons and thus have a smaller ionic radius as the temperature increases.

In this paper, we use weak measurement to study the mixed magneto-optical Kerr effects (MOKEs) in a magnetic single layer of Co and in organic/Co bilayer films. The relationship between the amplified shift and the MOKE parameters is theoretically established and then experimentally observed as a function of both the magnetization intensity and the magnetization direction in magnetic thin films with an arbitrary magnetization. Furthermore, we experimentally observe a magnetic coupling at the organic/ferromagnetic interface. Due to this being a robust coupling, there is a significant dependence of the amplified shift with the thickness of the organic layer. This work has value for its application in the magnetic measurement technique for the more general case based on weak measurement.

Negative thermal expansion (NTE) in ZrW 2 O 8 was investigated using a flexibility analysis of ab initio phonons. It was shown that no previously proposed mechanism adequately describes the atomic-scale origin of NTE in this material. Instead it was found that the NTE in ZrW 2 O 8 is driven, not by a single mechanism, but by wide bands of phonons that resemble vibrations of near-rigid WO 4 units and Zr-O bonds at low frequency, with deformation of O-W-O and O-Zr-O bond angles steadily increasing with increasing NTE-phonon frequency. It is asserted that this phenomenon is likely to provide a more accurate explanation for NTE in many complex systems not yet studied.

We have performed first-principles calculations of the structure and lattice dynamics in the metal hexaborides SmB 6 , CaB 6 , SrB 6 and BaB 6 using Density Functional Theory in an attempt to understand the negative thermal expansion in the first of these materials. The focus is on the role of Rigid Unit Modes involving rotations of the B 6 octahedra similar to the rotations of structural polyhedra connected by bonds in Zn(CN) 2 , Prussian Blue and Si(NCN) 2 . However, it was found that there is very low flexibility of the network of connected B 6 octahedra, and the lattice dynamics do not support negative thermal expansion except possibly at very low temperature. Thus the negative thermal expansion observed in SmB 6 probably has an electronic origin.

We have studied the two phases of the molecular crystal acetylene, C 2 H 2 , using calculations of the lattice dynamics by Density Functional Theory methods. together with the use of classical molecular dynamics (MD) simulation methods. The two phases share the same simple face-centred cubic lattice arrangement of the molecular centres of mass, but with different molecular orientations. We show that the higher-temperature phase has lower phonon frequencies and hence higher entropy, giving thermodynamic stability at higher temperature. The calculated lattice dynamics of this phase show instabilities associated with phonons involving pure rotations of the molecules. The MD results show large amplitudes of librational motion in this phase. The MD simulations also showed a potential phase transition to a structure of tetragonal symmetry. The picture that emerges is that the phase transition in acetylene is a very rare example of one that encompasses elements of three types of transition: displacive, order-disorder and reconstructive.

We review a set of ideas concerning the flexibility of network materials, broadly defined as structures in which atoms form small polyhedral units that are connected at corners. One clear example is represented by the family of silica polymorphs, with structures composed of corner-lined SiO 4 tetrahedra. The Rigid Unit Mode (RUM) is defined as any normal mode in which the structural polyhedra can translate and/or rotate without distortion, and since forces associated with changing the size and shape of the polyhedra are much stronger than those associated with rotations of two polyhedra around a shared vertex, the RUMs might be expected to have low frequencies compared to all other phonon modes. In this paper we discuss the flexibility of network structures, and how RUMs can arise in such structures, both in principle and in a number of specific examples of real systems. We also discuss applications of the RUM model, particularly for our understanding of phenomena such as displacive phase transitions and negative thermal expansion in network materials.&#xD., (© 2023 IOP Publishing Ltd.)

The orientational disorder in crystalline sulfur hexafluoride, SF 6 , has been studied using a combination of neutron total scattering and the reverse Monte Carlo method. Analysis of the atomic configurations has shown the extent of the disorder through the evaluation of the S-F bond orientational distribution function, consistent with, but improving upon, the results of earlier neutron powder diffraction data. The correlations between orientations of neighbouring molecules have been studied through analysis of the distributions of F-F distances, showing that nearest-neighbour F-F close contacts are avoided, consistent with previous molecular dynamics simulation results. The results present a new case study of the application of neutron total scattering and the reverse Monte Carlo methods for the study of orientational disorder, where in this instance the disorder arises from orientational frustration rather than from a mismatch of molecular and site symmetries., (Creative Commons Attribution license.)

We report results from a study of the crystal and magnetic structures of strontium-doped BiFeO 3 using neutron powder diffraction and the Rietveld method. Measurements were obtained over a wide range of temperatures from 300-800 K for compositions between 10%-16% replacement of bismuth by strontium. The results show a clear variation of the two main structural deformations-symmetry-breaking rotations of the FeO 6 octahedra and polar ionic displacements that give ferroelectricity-with chemical composition, but relatively little variation with temperature. On the other hand, the antiferromagnetic order shows a variation with temperature and a second-order phase transition consistent with the classical Heisenberg model. There is, however, very little variation in the behaviour of the antiferromagnetism with chemical composition, and hence with the degree of the structural symmetry-breaking distortions. We therefore conclude that there is no significant coupling between antiferromagnetism and ferroelectricity in Sr-doped BiFeO 3 and, by extension, in pure BiFeO 3 ., (Creative Commons Attribution license.)

We report the results of a neutron powder diffraction study of the phase transitions in deuterated methylammonium lead iodide, with a focus on the system of orientational distortions of the framework of PbI 6 octahedra. The results are analysed in terms of symmetry-adapted lattice strains and normal mode distortions. The higher-temperature cubic-tetragonal phase transition at 327 K is weakly discontinuous and nearly tricritical. The variations of rotation angles and spontaneous strains with temperature are consistent with a standard Landau theory treatment. The lower-temperature transition to the orthorhombic phase at 165 K is discontinuous, with two systems of octahedral rotations and internal distortions that together can be described by 5 order parameters of different symmetry. In this paper we quantify the various symmetry-breaking distortions and their variation with temperature, together with their relationship to the spontaneous strains, within the formalism of Landau theory. A number of curious results in the low-temperature phase are identified, particularly regarding distortion amplitudes that decrease rather than increase with lowering temperature., (Creative Commons Attribution license.)

Hybrid organic-inorganic perovskite (HOIP) ferroelectrics are attracting considerable interest because of their high performance, ease of synthesis, and lightweight. However, the intrinsic thermodynamic origins of their ferroelectric transitions remain insufficiently understood. Here, we identify the nature of the ferroelectric phase transitions in displacive [(CH 3 ) 2 NH 2 ][Mn(N 3 ) 3 ] and order-disorder type [(CH 3 ) 2 NH 2 ][Mn(HCOO) 3 ] via spatially resolved structural analysis and ab initio lattice dynamics calculations. Our results demonstrate that the vibrational entropy change of the extended perovskite lattice drives the ferroelectric transition in the former and also contributes importantly to that of the latter along with the rotational entropy change of the A-site. This finding not only reveals the delicate atomic dynamics in ferroelectric HOIPs but also highlights that both the local and extended fluctuation of the hybrid perovskite lattice can be manipulated for creating ferroelectricity by taking advantages of their abundant atomic, electronic, and phononic degrees of freedom.

We report calculations of the crystal structures and lattice dynamics of the tetragonal and orthorhombic phases of the molecular crystal hydrogen cyanide, HCN, using density functional theory methods. By treating negative pressure as a proxy for raising temperature we show that the ferroelastic phase transition involves softening of a transverse acoustic mode, and confirm that the phase transition is discontinuous. Analysis of the complete phonon spectrum shows that the acoustic modes are responsible both for the very large thermal expansion seen in HCN and also for the thermodynamic driving force for the phase transition., (© 2021 IOP Publishing Ltd.)

We report a study of the orientational order and phase transitions in crystalline deuterated methane, carried out using neutron total scattering and the reverse Monte Carlo method. The resultant atomic configurations are consistent with the average structures obtained from Rietveld refinement of the powder diffraction data, but additionally enable us to determine the C-D bond orientational distribution functions (ODF) for the disordered molecules in the high-temperature phase, and for both ordered and disordered molecules in the intermediate-temperature phase. We show that this approach gives more accurate information than can been obtained from fitting a bond ODF to diffraction data. Given the resurgence of interest in orientationally-disordered crystals, we argue that the approach of total scattering with the RMC method provides a unique quantification of orientational order and disorder., (© 2021 IOP Publishing Ltd.)

Magic-size clusters are ultra-small colloidal semiconductor systems that are intensively studied due to their monodisperse nature and sharp UV-vis absorption peak compared with regular quantum dots. However, the small size of such clusters (<2 nm), and the large surface-to-bulk ratio significantly limit characterisation techniques that can be utilised. Here we demonstrate how a combination of EXAFS and XANES analyses can be used to obtain information about sample stoichiometry and cluster symmetry. Investigating two types of clusters that show sharp UV-vis absorption peaks at 311 nm and 322 nm, we found that both samples possess approximately 2 : 1 Cd : S ratio and have similar nearest-neighbour structural arrangements. However, both samples demonstrate a significant departure from the tetrahedral structural arrangement, with an average bond angle determined to be around 106.1° showing a bi-fold bond angle distribution. Our results suggest that both samples are quasi-isomers - their core structures have identical chemical compositions, but different atomic arrangements with distinct bond angle distributions.

We report an ab initio calculation of crystal structure and lattice dynamics of cubic silicon dicarbodiimide, Si(NCN) 2 , using density functional theory methods. The calculations reveal a low-energy spectrum of rigid unit modes that are shown to be associated with negative thermal expansion. Comparisons are drawn with the closely-related materials Zn(CN) 2 and the cubic-cristobalite phase of SiO 2 . Instabilities in the spectrum of rigid unit modes point to the existence of disorder of the positions and orientations of the dicarbodiimide molecular anions., (© 2020 IOP Publishing Ltd.)

The counterintuitive phenomenon of pressure-induced softening in materials is likely to be caused by the same dynamical behavior that produces negative thermal expansion. Through a combination of molecular dynamics simulation on an idealized model and neutron diffraction at variable temperature and pressure, we show the existence of extraordinary and unprecedented pressure-induced softening in the negative thermal expansion material scandium fluoride ScF&#95;{3}. The pressure derivative of the bulk modulus B, B^{'}=(∂B/∂P)&#95;{P=0}, reaches values as low as -220±30 at 50 K, and is constant at -50 between 150 and 250 K.

Orientational disorder of the molecularCO32-anions in BaCO 3 , which occurs naturally as the mineral witherite, has been studied using a combination of neutron total scattering analysed by the reverse Monte Carlo method and molecular dynamics simulations. The primary focus is on the phase transition to the cubic phase, which assumes a rocksalt structure (Strukturbericht type B1) with highly disordered orientations consistent with the mismatch between the site (m3¯m) and molecular (3/ m ) symmetries. Both experiment and simulation show a high degree of disorder, with the C-O bond orientation distribution never exceeding 25% variation from that of a completely uniform distribution, although there are differences between the two methods regarding the nature of these variations. Molecular dynamics simulations are also reported for the analogous phase transitions in the very important mineral calcite, CaCO 3 . The combination of the simulations and comparison with BaCO 3 shows that the properties of calcite at all temperatures within its stability field are affected mostly by the onset of orientational disorder associated with the high-temperature cubic phase, even though this lies outside the stability field of calcite. This is a new understanding of calcite, which previously had been interpreted purely in terms of the phase transition to an intermediate partially-disordered phase. Finally, we also found that witherite itself appears to support the development of orientational disorder on heating, with the simulations showing a sequence of phase transitions that explain the much larger thermal expansion of one axis., (Creative Commons Attribution license.)

Zero thermal expansion (ZTE) and zero linear compressibility (ZLC) are unique and rare properties. Materials combining ZTE and ZLC will have promising prospects. A novel route is proposed in this work to design the coexistence of uniaxial-ZTE and ZLC based on layered hybrid semiconductors [ZnTe(L)0.5] [L = N2H4, ethylenediamine (en), propyldiamine (pda)]. In the framework of [ZnTe(L)0.5], the organic chain contains the attractive and repulsive interactions that arise from the different organic components. It is demonstrated that changing the length of the organic chain can effectively regulate the interaction between different organic components and then achieve the uniaxial-ZTE and ZLC or the desired thermal expansion and compression behaviors. The origin of the coexistence of abnormal axial responses has been traced from thermodynamic formalisms, model Grüneisen parameters and specific vibration modes. It is found that low-energy phonons play an important internal role in realizing the multi-peculiar properties.

The many-body polarization energy is the major source of nonadditivity in strongly polar systems such as water. This nonadditivity is often considerable and must be included, if only in an average manner, to correctly describe the physical properties of the system. Models for the polarization energy are usually parametrized using experimental data, or theoretical estimates of the many-body effects. Here we show how many-body polarization models can be developed for water complexes using data for the monomer and dimer only using ideas recently developed in the field of intermolecular perturbation theory and state-of-the-art approaches for calculating distributed molecular properties based on the iterated stockholder atoms (ISA) algorithm. We show how these models can be calculated, and we validate their accuracy in describing the many-body nonadditive energies of a range of water clusters. We further investigate their sensitivity to the details of the polarization damping models used. We show how our very best polarization models yield many-body energies that agree with those computed with coupled-cluster methods, but at a fraction of the computational cost.

Four types of magic-size CdS clusters and three different CdS quantum dots have been studied using the technique of X-ray total scattering and pair distribution function analysis. We found that the CdS quantum dots could be modelled as a mixed phase of atomic structures based on the two bulk crystalline phases, which is interpreted as representing the effects of random stacking of layers. However, the results for the magic-size clusters are significantly different. On one hand, the short-range features in the pair distribution function reflect the bulk, indicating that these structures are based on the same tetrahedral coordination found in the bulk phases (and therefore excluding new types of structures such as cage-like arrangements of atoms). But on the other hand, the longer-range atomic structure clearly does not reflect the layer structures found in the bulk and the quantum dots. We compare the effect of two ligands, phenylacetic acid and oleic acid, showing that in one case the ligand has little effect on the atomic structure of the magic-size nanocluster, and in another it has a significant effect.

The Rigid Unit Mode model was initially developed to understand the origin of displacive phase transitions in network silicate structures. Here, we review the successes of the model, and consider how it might apply to a wider range of network structures. This article is part of the theme issue 'Mineralomimesis: natural and synthetic frameworks in science and technology'.

The hybrid perovskites are coordination frameworks with the same topology as the inorganic perovskites, but with properties driven by different chemistry, including host-framework hydrogen bonding. Like the inorganic perovskites, these materials exhibit many different phases, including structures with potentially exploitable functionality. However, their phase transformations under pressure are more complex and less well understood. We have studied the structures of manganese and cobalt guanidinium formate under pressure using single-crystal X-ray and powder neutron diffraction. Under pressure, these materials transform to a rhombohedral phase isostructural to cadmium guanidinium formate. This transformation accommodates the reduced cell volume while preserving the perovskite topology of the framework. Using density-functional theory calculations, we show that this behaviour is a consequence of the hydrogen-bonded network of guanidinium ions, which act as struts protecting the metal formate framework against compression within their plane. Our results demonstrate more generally that identifying suitable host-guest hydrogen-bonding geometries may provide a route to engineering hybrid perovskite phases with desirable crystal structures. This article is part of the theme issue 'Mineralomimesis: natural and synthetic frameworks in science and technology'.

The crystal structure of the low-temperature [Formula: see text] phase of crystalline malononitrile, CH 2 (CN) 2 (stable phase below 260 K), has been determined using Rietveld refinement on neutron powder diffraction data. The [Formula: see text] phase has a slightly lower density than the other three low-pressure phases, and unlike those phases it has a polar structure. The transition from the [Formula: see text] to [Formula: see text] phase involves a major reconstruction of the structure, including establishing a network of hydrogen bonds. DFT simulations of the structure and phonon dispersion curves of both [Formula: see text] and [Formula: see text] phases give free energy curves consistent with the phase transitions. It is noted that the transition from the [Formula: see text] to [Formula: see text] phase at 260 K is facilitated by the entropy arising from the low-frequency phonons associated with the soft mode for the sequence of [Formula: see text]-[Formula: see text]-[Formula: see text] phase transitions.

A hallmark of a thermodynamic phase transition is the qualitative change of system thermodynamic properties such as energy and heat capacity. On the other hand, no phase transition is thought to operate in the supercritical state of matter and, for this reason, it was believed that supercritical thermodynamic properties vary smoothly and without any qualitative changes. Here, we perform extensive molecular dynamics simulations in a wide temperature range and find that a deeply supercritical state is thermodynamically heterogeneous, as witnessed by different temperature dependence of energy, heat capacity and its derivatives at low and high temperature. The evidence comes from three different methods of analysis, two of which are model-independent. We propose a new definition of the relative width of the thermodynamic crossover and calculate it to be in the fairly narrow relative range of 13%-20%. On the basis of our results, we relate the crossover to the supercritical Frenkel line.

Physical properties of an interacting system are governed by collective excitations, but their nature at extreme supercritical conditions is unknown. Here, we present direct evidence for propagating solid-like longitudinal phonon-like excitations with wavelengths extending to interatomic separations deep in the supercritical state at temperatures up to 3,300 times the critical temperature. We observe that the crossover of dispersion curves develops at k points reducing with temperature. We interpret this effect as the crossover from the collective phonon to the collisional mean-free path regime of particle dynamics and find that the crossover points are close to both the inverse of the shortest available wavelength in the system and to the particle mean free path inferred from experiments and theory. Notably, both the shortest wavelength and mean free path scale with temperature with the same power law, lending further support to our findings.

Structural isomerism of colloidal semiconductor nanocrystals has been largely unexplored. Here, we report one pair of structural isomers identified for colloidal nanocrystals which exhibit thermally-induced reversible transformations behaving like molecular isomerization. The two isomers are CdS magic-size clusters with sharp absorption peaks at 311 and 322 nm. They have identical cluster masses, but slightly different structures. Furthermore, their interconversions follow first-order unimolecular reaction kinetics. We anticipate that such isomeric kinetics are applicable to a variety of small-size functional nanomaterials, and that the methodology developed for our kinetic study will be helpful to investigate and exploit solid-solid transformations in other semiconductor nanocrystals. The findings on structural isomerism should stimulate attention toward advanced design and synthesis of functional nanomaterials enabled by structural transformations.

Background: Sporadic fatal insomnia (sFI) is a rapid progressive neurodegenerative disease characterised by gradual to perpetual insomnia, followed by dysautonomia, coma and death. 1 The cause of sFI was recently mapped to a mutation in a protein, the prion, found in the human brain. It is the unfolding of the prion that leads to the generation of toxic oligomers that destroy brain tissue and function. Recent studies have confirmed that a methionine mutation at codon 129 of the human Prion is characteristic of sFI. Current treatment slows down the progression of the disease, but no cure has been found, yet., Methods: We used Molecular Docking and Molecular Dynamics simulation methods, to study the toxic Fatal-Insomnia-prion conformations at local unfolding. The idea was to determine these sites and to stabilise these regions against unfolding and miss-folding, using a small ligand, based on a phenothiazine "moiety"., Conclusion: As a result we here discuss current fatal insomnia therapy and present seven novel possible compounds for in vitro and in vivo screening.

We study the thermomechanical properties of matter under extreme conditions deep in the supercritical state, at temperatures exceeding the critical one by up to four orders of magnitude. We calculate the Grüneisen parameter γ and find that on isochores it decreases with temperature from 3 to 1, depending on the density. Our results indicate that from the perspective of thermomechanical properties, the supercritical state is characterized by a wide range of γ's which includes solidlike values-an interesting finding in view of the common perception of the supercritical state as being an intermediate state between gases and liquids. We rationalize this result by considering the relative weights of oscillatory and diffusive components of the supercritical system below the Frenkel line. We also find that γ is nearly constant at the Frenkel line above the critical point and explain this universality in terms of the pressure and temperature scaling of system properties along the lines where particle dynamics changes qualitatively.

Fundamental understanding of strongly interacting systems necessarily involves collective modes, but their nature and evolution is not generally understood in dynamically disordered and strongly interacting systems such as liquids and supercritical fluids. We report the results of extensive molecular dynamics simulations and provide direct evidence that liquids develop a gap in a solidlike transverse spectrum in the reciprocal space, with no propagating modes between zero and a threshold value. In addition to the liquid state, this result importantly applies to the supercritical state of matter. We show that the emerging gap increases with the inverse of liquid relaxation time and discuss how the gap affects properties of liquid and supercritical states.

The molecular crystals adamantane, C 10 H 16 , and adamantanecarboxylic acid, C 10 H 15 COOH, undergo order-disorder phase transitions at 208 and 250 K, respectively. Reverse Monte Carlo refinement of total neutron scattering data collected from deuterated samples immediately above these phase transitions shows that the high-temperature phases are well described by models in which the adamantyl groups are disordered over two sites. No correlation between the orientations of neighbouring molecules is observed. These results demonstrate that the intermolecular potential energy of these materials depends strongly on the orientation of the reference molecule but only very weakly on the orientations of its neighbours., (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

We develop an approach to liquid thermodynamics based on collective modes. We perform extensive molecular-dynamics simulations of noble, molecular, and metallic liquids, and we provide direct evidence that liquid energy and specific heat are well-described by the temperature dependence of the Frenkel (hopping) frequency. The agreement between predicted and calculated thermodynamic properties is seen in the notably wide range of temperature spanning tens of thousands of Kelvin. The range includes both subcritical liquids and supercritical fluids. We discuss the structural crossover and interrelationships between the structure, dynamics, and thermodynamics of liquids and supercritical fluids.

Carbon capture and storage (CCS) offers a possible solution to curb the CO 2 emissions from stationary sources in the coming decades, considering the delays in shifting energy generation to carbon neutral sources such as wind, solar and biomass. The most mature technology for post-combustion capture uses a liquid sorbent, amine scrubbing. However, with the existing technology, a large amount of heat is required for the regeneration of the liquid sorbent, which introduces a substantial energy penalty. The use of alternative sorbents for CO 2 capture, such as the CaO-CaCO 3 system, has been investigated extensively in recent years. However there are significant problems associated with the use of CaO based sorbents, the most challenging one being the deactivation of the sorbent material. When sorbents such as natural limestone are used, the capture capacity of the solid sorbent can fall by as much as 90 mol% after the first 20 carbonation-regeneration cycles. In this study a variety of techniques were employed to understand better the cause of this deterioration from both a structural and morphological standpoint. X-ray and neutron PDF studies were employed to understand better the local surface and interfacial structures formed upon reaction, finding that after carbonation the surface roughness is decreased for CaO. In situ synchrotron X-ray diffraction studies showed that carbonation with added steam leads to a faster and more complete conversion of CaO than under conditions without steam, as evidenced by the phases seen at different depths within the sample. Finally, in situ X-ray tomography experiments were employed to track the morphological changes in the sorbents during carbonation, observing directly the reduction in porosity and increase in tortuosity of the pore network over multiple calcination reactions.

Negative thermal expansion (NTE) is the phenomenon in which materials shrink rather than expand on heating. Although NTE had been previously observed in a few simple materials at low temperature, it was the realisation in 1996 that some materials have NTE over very wide ranges of temperature that kick-started current interest in this phenomenon. Now, nearly two decades later, a number of families of ceramic NTE materials have been identified. Increasingly quantitative studies focus on the mechanism of NTE, through techniques such as high-pressure diffraction, local structure probes, inelastic neutron scattering and atomistic simulation. In this paper we review our understanding of vibrational mechanisms of NTE for a range of materials. We identify a number of different cases, some of which involve a small number of phonons that can be described as involving rotations of rigid polyhedral groups of atoms, others where there are large bands of phonons involved, and some where the transverse acoustic modes provide the main contribution to NTE. In a few cases the elasticity of NTE materials has been studied under pressure, identifying an elastic softening under pressure. We propose that this property, called pressure-induced softening, is closely linked to NTE, which we can demonstrate using a simple model to describe NTE materials. There has also been recent interest in the role of intrinsic anharmonic interactions on NTE, particularly guided by calculations of the potential energy wells for relevant phonons. We review these effects, and show how anhamonicity affects the response of the properties of NTE materials to pressure.

We report total neutron scattering measurements on the metal-organic perovskite analogue dimethylammonium manganese(ii) formate, (CD3)2ND2[Mn(DCO2)3]. Reverse Monte Carlo modelling shows that, in both the disordered high-temperature and ordered low-temperature phases, the ammonium moiety forms substantially shorter hydrogen bonds (N...O = 2.4 Å and 2.6 Å) than are visible in the average crystal structures. These bonds result from a pincer-like motion of two adjacent formate ions about the dimethylammonium ion in such a way that the framework can adjust independently to the positions of nearest-neighbour dimethylammonium ions. At low temperatures the shortest hydrogen bond is less favourable, apparently because it involves a greater distortion of the framework. Furthermore, in the high-temperature phase, in addition to the three disordered nitrogen positions expected from the average crystal structure, there appear to be also smaller probability maxima between these positions, corresponding to orientations in which the dimethylammonium is hydrogen-bonded to the two oxygen atoms of a single formate ion. The spontaneous strain across the phase transition reveals a contraction of the framework about the dimethylammonium cation, continuing as the material is cooled below the transition temperature. These results provide direct evidence of the local atomic structure of the guest-framework hydrogen bonding, and in particular the distortions of the framework responsible for the phase transition in this system.

We report the results of a series of molecular dynamics simulations on a number of zinc zeolitic imidazolate framework (ZIF) structures together with some lattice dynamics calculations on ZIF-4, providing information about the flexibilities of these structures. The simulations have used a force field we developed based on ab initio calculations of clusters of ligands and metal cations. We have shown that there are instabilities of the structures of some ZIF structures at low temperatures and high pressures. A rigidity analysis based on the Rigid Unit Mode model shows considerable degree of network flexibility, including a significant elastic flexibility.

A simulation study of negative thermal expansion in Y2W3O12 was carried out using calculations of phonon dispersion curves through the application of density functional perturbation theory. The mode eigenvectors were mapped onto flexibility models and results compared with calculations of the mode Grüneisen parameters. It was found that many lower-frequency phonons contribute to negative thermal expansion in Y2W3O12, all of which can be described in terms of rotations of effectively rigid WO4 tetrahedra and Y-O rods. The results are strikingly different from previous phonon studies of higher-symmetry materials that show negative thermal expansion.

A new dynamic line, the Frenkel line, has recently been proposed to separate the supercritical state into rigid-liquid and nonrigid gaslike fluid. The location of the Frenkel line on the phase diagram is unknown for real fluids. Here we map the Frenkel line for three important systems: CO(2), H(2)O, and CH(4). This provides an important demarcation on the phase diagram of these systems, the demarcation that separates two distinct physical states with liquidlike and gaslike properties. We find that the Frenkel line can have a similar trend as the melting line above the critical pressure. Moreover, we discuss the relationship between unexplained solubility maxima and Frenkel line, and we propose that the Frenkel line corresponds to the optimal conditions for solubility.

The mechanism of negative thermal expansion (NTE) in MOF-5 has been studied using a rigid unit mode analysis of the phonons responsible for NTE. In addition to confirming the role of the previously proposed optic 'trampoline mode' mechanism [Lock et al., Dalton Trans., 2013, 42, 1933-2320], we have identified a new acoustic mode mechanism that makes a major contribution to the macroscopic NTE of this framework.

By exploring a simple model of a negative thermal expansion (NTE) system, we introduce a phenomenological expression to describe the temperature dependence of the pressure-induced softening in NTE structures.

Electronic effects have been shown to be important in high-energy radiation damage processes where a high electronic temperature is expected, yet their effects are not currently understood. Here, we perform molecular dynamics simulations of high-energy collision cascades in α-iron using a coupled two-temperature molecular dynamics (2T-MD) model that incorporates both the effects of electronic stopping and electron-phonon interaction. We subsequently compare it with the model employing electronic stopping only, and find several interesting novel insights. The 2T-MD results in both decreased damage production in the thermal spike and faster relaxation of the damage at short times. Notably, the 2T-MD model gives a similar amount of final damage at longer times, which we interpret to be the result of two competing effects: a smaller amount of short-time damage and a shorter time available for damage recovery.

We measured the density of vibrational states (DOS) and the specific heat of various glassy and crystalline polymorphs of SiO2. The typical (ambient) glass shows a well-known excess of specific heat relative to the typical crystal (α-quartz). This, however, holds when comparing a lower-density glass to a higher-density crystal. For glassy and crystalline polymorphs with matched densities, the DOS of the glass appears as the smoothed counterpart of the DOS of the corresponding crystal; it reveals the same number of the excess states relative to the Debye model, the same number of all states in the low-energy region, and it provides the same specific heat. This shows that glasses have higher specific heat than crystals not due to disorder, but because the typical glass has lower density than the typical crystal.

Two solid phases of cyclohexane have been investigated over a temperature range spanning 13-266 K on a powdered, perdeuterated sample using neutron total scattering. Phase II has an ordered structure (C2/c) that forms below 186 K. Between 186 and 280 K it exists as a plastic solid-phase I (Fm3m), where the molecules are rotationally disordered about the lattice points of the face-centred cubic cell. Data-dependent atomistic configurations that represent the 'instantaneous' crystal structure have been generated from the total scattering data using reverse Monte Carlo refinement. Analysis of the local structure reveals that instantaneous distortions in phase I resemble the average structure of phase II.