^{1}Department of Materials and The Thomas Young Centre, Imperial College London, London, United Kingdom^{2}Department of Physics, Imperial College London, London, United Kingdom^{3}School of Physical Sciences, University of Kent, Canterbury, United Kingdom^{4}Department of Chemistry, University of Warwick, Coventry, United Kingdom

Uniaxial negative thermal expansion (NTE) is known to occur in low *n* members of the A_{n+1}B_{n}O_{3n+1} Ruddlesden–Popper (RP) layered perovskite series with a frozen rotation of BO_{6} octahedra about the layering axis. Previous work has shown that this NTE arises due to the combined effects of a close proximity to a transition to a competing phase, so called “symmetry trapping”, and highly anisotropic elastic compliance specific to the symmetry of the NTE phase. We extend this analysis to the broader RP family (*n* = 1, 2, 3, 4, …, ∞), demonstrating that by changing the fraction of layer interface in the structure (i.e., the value of 1/*n*) one may control the anisotropic compliance that is necessary for the pronounced uniaxial NTE observed in these systems. More detailed analysis of how the components of the compliance matrix develop with 1/*n* allows us to identify different regimes, linking enhancements in compliance between these regimes to the crystallographic degrees of freedom in the structure. We further discuss how the perovskite layer thickness affects the frequencies of soft zone boundary modes with large negative Grüneisen parameters, associated with the aforementioned phase transition, that constitute the thermodynamic driving force for NTE. This new insight complements our previous work—showing that chemical control may be used to switch from positive to negative thermal expansion in these systems—since it makes the layer thickness, *n*, an additional design parameter that may be used to engineer layered perovskites with tuneable thermal expansion. In these respects, we predict that, with appropriate chemical substitution, the *n* = 1 phase will be the system in which the most pronounced NTE could be achieved.

## 1. Introduction

Ruddlesden–Popper (RP) oxides are an intriguing class of ceramic materials. They have the basic formula A_{n+1}B_{n}O_{3n+1} and consist of a perovskite block of *n* corner sharing BO_{6} octahedra separated by an AO rock salt layer. Blocks of octahedra are stacked perpendicular to the long crystallographic axis making this layering axis structurally distinct from the two in-plane axes. Neighboring blocks are de-phased from each other by a lattice translation of (0.5, 0.5, 0.5), and the aristotypical symmetry for any value of *n* is *I*4/*mmm* (Figure 1). As with the perovskites, the A-site chemistry is dominated by larger alkali, alkali-earth, and rare earth metals, and the B-site by transition metals. In the limit *n* = ∞, the perovskite structure is recovered. While in practice most chemistries are found to predominantly exhibit the *n* = 1, 2 phases only (Palgrave et al., 2012), in principle any value of *n* between 1 and ∞ is possible; *n* = 3 structures have been synthesized by careful compositional control (Battle et al., 1998) and although *n* > 3 phases are often predicted to be unstable to decomposition (McCoy et al., 1997), epitaxial growth techniques have allowed the synthesis of *n* = 2−5 (Haeni et al., 2001), *n* = 6 (Yan et al., 2011), and *n* = 10 (Lee et al., 2013) structures.

**Figure 1**. In the A_{n+1}B_{n}O_{3n+1} Ruddlesden–Popper series, blocks consisting of *n* layers of ABO_{3} perovskite structure are separated by a single layer of AO rock salt structure, with BO_{6} octahedra in the next block displaced by a (0.5, 0.5, 0.5) lattice translation. In the *n* = ∞ limit the pure ABO_{3} perovskite structure is recovered.

One of the most explored systems is the *n* = 1 A_{2}CuO_{4} on account of its high-temperature superconductivity, where doping of divalent A = Ba and Sr with trivalent rare earth cations has been extensively investigated (Dwivedi and Cormack, 1991). Superconductivity in these systems is not limited to the cuprates, and there has been substantial interest in Sr_{2}RuO_{4} (Mackenzie and Maeno, 2003) on account of its superconducting phase transition below 0.8 K, and in Sr_{3}Ru_{2}O_{7} for its nematic orbitally-ordered phase (Borzi et al., 2007). The doped nickelates have also been much studied due to their believed proximity to a superconducting phase transition, and their charge ordering physics (Yoshizawa et al., 2000). More recently, the *n* = 2 member of the RP family has received much attention on the account of a new form of improper ferroelectricity predicated in Ca_{3}Mn_{2}O_{7} and Ca_{3}Ti_{2}O_{7}, termed hybrid improper ferroelectricity (Benedek and Fennie, 2011). This mechanism circumvents the so-called *d*^{0} criterion for ferroelectricity, as it does not require an off-centring of cations to drive the phase transition. Instead, this off-centring (*P*) may occur as a slave process driven by an octahedral tilt (*R*_{1}) and rotation mode (*R*_{2}) of the parent structure that are inherently unstable in some of these systems. This leads to a so-called trilinear term β*R*_{1}*R*_{2}*P* in the free energy expansion about the parent structure (Benedek et al., 2015) which, regardless of the sign of the coefficient β, as *R*_{1} and *R*_{2} are inherently unstable, leads to a non-zero value of the polarization *P* (either positive or negative depending on the sign of β).

Our contribution to this field of hybrid improper ferroelectricity was to provide experimental confirmation of this mechanism for the case of Ca_{3}Ti_{2}O_{7} (Senn et al., 2015). However, our high-resolution powder diffraction data for Ca_{3}Mn_{2}O_{7} revealed an added complexity. What was believed to be a single phase at room temperature, having the polar symmetry *A*2_{1}*am*, was in fact a mixture of this and a phase that we identified as having *Acaa* symmetry. Crucially, this phase only has a single octahedral rotation that is out-of-phase rather than in-phase, and the octahedra remain untilted along the *c*-axis. No hybrid improper ferroelectric mechanism is therefore possible. However, over the large phase coexistence region, which spans a temperature range of 120 K, we did observe pronounced uniaxial negative thermal expansion (NTE) along the *c*-axis in the *Acaa* phase. This had not been observed previously in the *n* = 2 system and, although reported in the literature in an analogous *n* = 1 system (Takahashi and Kamegashira, 1993), its significance had not been noted. We were able to explain this NTE phenomenon as being driven by a leftover degree of freedom, the octahedral tilt in the *Acaa* phase, which remains dynamic.

NTE is a rare property, that when it does occur is known to be caused by a diverse range of mechanisms in different materials. Even within inorganic perovskite-based systems, NTE has been found to originate due to coupling of the lattice parameters to: charge ordering (Azuma et al., 2011), ferroelectric ordering (Chen et al., 2013), and magnetic and orbital ordering via an invar-like mechanism (Yoshida et al., 2005, Qi et al., 2010). In framework structures, formed from connected strongly-bonded polyhedral units, NTE has been explained by transverse vibrations of these units, known as rigid-unit modes (RUMs) (Dove et al., 1995, Heine et al., 1999). We argued in *Acaa* Ca_{3}Mn_{2}O_{7} that certain vibrational modes with RUM character would have negative Grüneisen parameters and be soft on account of the proximity of the system to the symmetry-forbidden phase transition to *A*2_{1}*am* (Senn et al., 2015). Using this idea of “trapping” a soft mode in the *Acaa* phase of Ca_{3}Mn_{2}O_{7} to systematically control and tune the uniaxial thermal expansion properties of the solid solution Ca_{3−x}Sr_{x}Mn_{2}O_{7} (Senn et al., 2016), we were able to demonstrate that this is a property exclusively of the *Acaa* phase in these materials, and that NTE is enhanced as the system approaches the *A*2_{1}*am* phase boundary as a function of chemical composition *x*. Although other effects operate in related materials, in this study we restrict our discussion to NTE driven by the coupling of the cell parameters to soft lattice modes since it is the most appropriate mechanism to describe our system.

The presence of dynamic octahedral tilts in this Ca_{3−x}Sr_{x}Mn_{2}O_{7} system explained the thermodynamic driving force for NTE along the layering axis. However, the question remained open of why NTE was only observed in this *Acaa* phase with a frozen in-plane rotation and not in the high-symmetry *I*4/*mmm* phase or related ABO_{3} perovskite phases, where dynamic octahedral tilts would still operate. We were able to answer this question in a recent computational study using density functional theory (DFT) and working within the quasi-harmonic approximation (QHA) to reproduce experimentally measured uniaxial NTE in the *I*4_{1}/*acd* phase of *n* = 1 Ca_{2}MnO_{4} (Ablitt et al., 2017).

Equation (1) (Grüneisen and Goens, 1924) describes the thermal expansion, α_{η}(*T*), at temperature, *T*, of the three cell axes of a tetragonal crystal (η = 1, 2, 3 where α_{1} = α_{2} by symmetry). Equation (1) is explained in detail in Appendix 1 and the concept of a Φ vector driving bulk PTE being transformed by a highly anisotropic * s* into uniaxial NTE (Barron and Munn, 1967) is depicted pictorially in Figure 2. In this picture, the anisotropic thermal expansion is separated into a thermodynamic driving force vector, Φ—arising from the lattice dynamics—that is transformed by the anisotropic elastic compliance matrix,

*. By computing the compliance matrix for our NTE phase, we were able to extract the thermodynamic driving force vector from our QHA simulation and found that the effect from dynamic tilts alone would not predict NTE over the wide temperature range observed in experiment. It is only when this thermodynamic driving force is transformed by the highly anisotropic elastic compliance of the layered RP phase that our simulations demonstrated uniaxial NTE of a magnitude and over a temperature range comparable to experiment. Comparing the compliance matrices computed for different phases, we found that particularly high anisotropic compliance is unique to the NTE phase of the RP structure and we linked this anisotropy to combined in-plane (frozen rotations) and out-of-plane (the AO layer) symmetry breaking in the NTE phase.*

**s****Figure 2**. Equation (1) is illustrated on axes describing normal cell deformations using Voigt notation (e.g., where the vector [ϵ, ϵ, 0] corresponds to normal strains of ϵ of the *a* and *b* axes with the *c* axis unstrained). A 3 × 3 elastic compliance matrix * s* is shown here in its quadratic form as an ellipsoid projected onto the ($\stackrel{\u0304}{1},1,0$) plane.

*may transform a vector describing the thermodynamic driving force for thermal expansion,*

**s****Φ**, in the positive quadrant (corresponding to bulk positive thermal expansion), into a thermal expansion vector,

**α**, in the quadrant corresponding to uniaxial negative thermal expansion of the

*c*axis. The direction of

**α**is given in the figure by utilizing the radius normal property of the representation ellipsoid of the matrix

*.*

**s**Until now, our research has focused on understanding uniaxial NTE in the low *n* RP compounds. In the present paper, we focus on extending these concepts to predict how the material properties that we have linked to NTE develop as a function of *n*. Studying the Ca_{n+1}Ge_{n}O_{3n+1} system within the framework of DFT, we find that the magnitude of anisotropic elastic compliance is dependent upon the proportion of CaGeO_{3}:CaO interface in the structure—which may be conveniently expressed by the fraction 1/*n*. This high compliance is then maximized with the highest proportion of interface (*n* = 1). To provide an explanation for this key result, we analyze how the components of the compliance matrix vary with 1/*n*. We identify a series of regimes in which groups of structures display similar elastic behavior based on the atomic degrees of freedom allowed by symmetry, and propose mechanisms by which these internal degrees of freedom couple to cell strains. The most important of these is the atomic “corkscrew” mechanism that operates at this interface in the NTE phase. We go on to investigate how the frequencies of the softest (lowest frequency) octahedral tilt modes, which provide the thermodynamic driving force for NTE, vary as a function of *n* in the high-symmetry and NTE phases. We find that a higher proportion of interface causes these phonons to stiffen (increase in frequency) and we infer that within a given chemical composition, the layer thickness *n* provides a structural constraint on an approximate temperature window for which the NTE phase will be stable—where increasing *n* drives this window to higher temperatures. One may thus use this insight, combined with our previous discovery that chemical substitution within a given structure may be used to tune the thermodynamic driving force for NTE, to use layer thickness *n* and composition as design parameters to engineer RP phases with optimal thermal expansion properties.

The layout of the paper is as follows: section 2 gives details of simulation parameters used in this work; section 3 is then split into four subsections: section 3.1 presents the key result of the paper, showing how the magnitude of anisotropic compliance, linked with uniaxial NTE, varies with 1/*n*; for the more interested reader, sections 3.2 and 3.3 then take a step back and analyze the origin of the result in section 3.1 in more detail, in section 3.2 by analyzing how the different elements of the compliance matrix evolve with 1/*n* and in section 3.3 by identifying different compliance regimes linked to crystallographic degrees of freedom in the structure that couple to strain; finally in section 3.4 we consider how the thermodynamic driving force for NTE varies with 1/*n* by presenting lattice dynamics calculations investigating phonons corresponding to tilts of GeO_{6} octahedra in Ca_{n+1}Ge_{n}O_{3n+1}, and compare these results against experimental phase diagrams constructed with data taken from the literature of the analogous Ca_{n+1}Mn_{n}O_{3n+1} system. Additionally, in the associated Appendices file, Appendix 1 gives a brief overview of the mathematical concepts relevant for the study of thermal expansion in an anisotropic material and Appendix 2 outlines the symmetry of the phases simulated throughout this study.

## 2. Methods

Calculations were performed using CASTEP, a plane-wave DFT code, version 7.0.3 (Clark et al., 2005). A plane wave cut-off energy of 1,400 eV was employed for all calculations with electron density stored on a grid twice as dense. A 7 × 7 × 2 Monkhorst–Pack grid of kpoints shifted away from the Γ-point was used for calculations of the 14 atom *I*4/*mmm* phase of (*n* = 1) Ca_{2}GeO_{4}, with grids of equivalent reciprocal space density used for all other structures (high-symmetry and rotation phases for *n* = 1, 2, 3, 4, ∞—see Appendix 2). Norm-conserving pseudopotentials, generated on-the-fly using CASTEP version 16.0, were used for all calculations and the associated pseudopotential strings may be found Table S1 in the Supplementary Information. All calculations used the PBEsol exchange-correlation functional (Perdew et al., 2008). Absolute energies were converged to an accuracy of 0.5 meV/atom with respect to k-point grid density and plane wave cut-off energy. Geometric relaxations were performed with a force tolerance of 10^{−4} eV/Å and a stress tolerance of 10 MPa.

We expect our Ca_{n+1}Ge_{n}O_{3n+1} system to be well-described by conventional GGA density functionals. There are other members of the chemical space that might require more careful consideration in terms of the appropriate methodology, such as hybrid functionals, DFT+U, or potentially even DMFT in order to accurately describe the physics associated with localized d and f-electrons.

Elastic constants were computed by fitting second order polynomials to the energies of cells with applied strains of ±0.2, 0.4% from the fully relaxed cell, where the internal degrees of freedom (the atomic positions) were free to relax. The quadratic terms to these fits were used to construct terms within the elastic constant matrix, **c**, and this matrix inverted to compute the elastic compliance matrix, * s* (see Appendix 1 for the definition of

*studied).*

**s**Bulk moduli, *K*, were computed by allowing the cell and all internal degrees of freedom to relax in response to hydrostatic pressures in the range −2 to +2 GPa. The bulk modulus was then found by fitting the computed relaxed volume, *V*, as a function of the external pressure, *P*, to the equation: $K=-\frac{d(ln\left[V\right])}{dP}$. The bulk compressibility, β, is then given by β = *K*^{−1}.

Density functional perturbation theory (DFPT) was used within CASTEP (Refson et al., 2006) to perform phonon calculations. In the present study, only phonon frequencies computed at single, high symmetry q-points are reported. In *n* = 1, 3, ∞ phases, these are at the *X* (½, ½, 0) and *P* (½, ½, ½) points in the *I*4/*mmm* high-symmetry phase (labeled *M* and *R*, respectively, in *n* = ∞ $Pm\stackrel{\u0304}{3}m$ ABO_{3}) and for *n* = 2, 4 phases at the *X*-point in *I*4/*mmm*. Phonons in child rotation phases were always computed at the Γ-point. In every compound studied, the initial structure was the highest symmetry phase that was fully relaxed. All child phases were found by freezing unstable phonons (modes with imaginary frequencies) into the structure with small amplitudes and allowing this child structure to relax. The lattice parameters and cell energies relative to the high-symmetry parent of all relaxed structures may be found in Table S2.

## 3. Results and Discussion

### 3.1. Anisotropy of the Compliance Matrix

It was previously shown in first-principles calculations performed on the NTE phase of Ca_{2}GeO_{4} (i.e., the *I*4_{1}/*acd* rotation phase) that a highly anisotropic * s* is an essential ingredient for uniaxial NTE in this system (Ablitt et al., 2017). κ, as defined in Equation (2), is the ratio of the highest (

*s*

_{H}) and lowest (

*s*

_{L}) eigenvalues of

*and gives a measure of the anisotropy of*

**s***; where higher κ indicates that*

**s***is more conducive to uniaxial NTE. Figure 3 therefore shows how*

**s***s*

_{H},

*s*

_{L}, and κ evolve with varying

*n*for high-symmetry and rotation phases in the Ca

_{n+1}Ge

_{n}O

_{3n+1}series.

**Figure 3. (A)** Lowest, *s*_{L} and highest, *s*_{H}, eigenvalues to the compliance matrix, * s*, and

**(B)**the anisotropy ratio κ plotted for high-symmetry and rotation phases in the Ca

_{n+1}Ge

_{n}O

_{3n+1}series against the mole fraction of CaGeO

_{3}:CaO interface (1/

*n*). Interpolations between values in the CaGeO

_{3}(1/

*n*= 0) and Ca

_{2}GeO

_{4}(1/

*n*= 1) structures are plotted for both high-symmetry and rotation phases (blue dot-dashed and red dotted lines, respectively).

In the composite mechanics community, the elastic properties of laminates are typically described by the properties of the constituent phases, weighted by the relative fraction of that phase (Sarlosi and Bocko, 2016). However, in the functional oxides community, it is well known that the local structure of atomic interfaces within a crystal often dictate the physical properties of the entire material 2012. There are therefore two ways to conceive the layered RP structure shown in Figure 1: (i) as being comprised of constituent CaGeO_{3} and CaO phases or (ii) as being comprised of CaGeO_{3} and the CaGeO_{3}:CaO interface.

From Equation (3), it is clear that (i) the mole fraction of CaO in Ca_{n+1}Ge_{n}O_{3n+1} is given by the ratio 1/(*n* + 1) and (ii) the mole fraction of Ca_{2}GeO_{4}, which represents the fraction of CaGeO_{3}:CaO interface in the structure, is given by the ratio 1/*n*. Therefore, to reflect the importance of the interface, *s*_{H}, *s*_{L}, and κ are plotted as a function of 1/*n* in the subplots on Figure 3. Straight lines have also been plotted interpolating between values for the CaGeO_{3} (1/*n* = 0) and Ca_{2}GeO_{4} (1/*n* = 1) end members to show how well the structure may be considered as a mixture of these two constituents in the high-symmetry (blue dot-dashed) and rotation (red dotted) phase series.

The least compliant eigenvector, *s*_{L}, corresponds to isotropic expansion/contraction for all structures investigated (see Table S3) and thus is closely linked to the bulk compressibility, β. Figure 3 shows that *s*_{L} increases linearly with higher Ca_{2}GeO_{4} mole fraction but is invariant to changes in symmetry for a given *n*. All values for *s*_{L} lie on the line interpolating between CaGeO_{3} and Ca_{2}GeO_{4} high-symmetry end members regardless of phase symmetry implying that *s*_{L} is determined mainly by the composition.

*s*_{H} also increases in magnitude with Ca_{2}GeO_{4} content for all RP phases. However, unlike *s*_{L}, *s*_{H} is greatly enhanced in the phase with a frozen rotation compared to the high-symmetry phase, and the rate of increase in *s*_{H} for rotation phases with 1/*n* is also greater in the rotation phase than in the high-symmetry parent. For all the tetragonal phases studied, the eigenvector *s*_{H} lies in a strain direction corresponding to a cooperative increase in in-plane lattice parameters, *ab*, and decrease in the lattice parameter along the layering axis, *c*, or visa versa. In a previous work, we proposed an atomic mechanism to facilitate a large compliance eigenvector in RP phases with a frozen octahedral rotation that relies on combined in-plane and out-of-plane symmetry breaking at the CaGeO_{3}:CaO layer interface to closely couple the *ab* and *c* axes (Ablitt et al., 2017). Since this mechanism operates at the CaGeO_{3}:CaO interface, it is interesting to note that *s*_{H} in the rotation phase is linearly dependent upon the mole fraction of this interface in the structure, increasing as this interface fraction becomes greater, and thus *s*_{H} for intermediate values of 1/*n* may be easily predicted by interpolating between the *s*_{H} values for CaGeO_{3} (with no interface) and Ca_{2}GeO_{4} (maximum interface) rotation phases.

This steeper increase in *s*_{H} for rotation phases than high-symmetry phases with interface mole fraction (1/*n*) thus manifests as a large enhancement in κ between the child structure and parent, where the magnitude of this enhancement increases greatly with 1/*n*, reaching a maximum in the *n* = 1 structure. The key result of this analysis of the compliances is hence that this *n* = 1 structure is the best in the Ruddelsden-Popper series for facilitating uniaxial NTE.

### 3.2. Elastic Compliances

Figure 3 showed how the eigenvalues of * s* vary with the CaGeO

_{3}:CaO interface fraction (1/

*n*). In this section we take a step back and analyse how the individual components of the compliance matrix,

*s*

_{ij}, vary with 1/

*n*. In the second half of the section, we assess the quality of the two interpolations, (i) between CaGeO

_{3}and CaO constituents and (ii) between CaGeO

_{3}and CaGeO

_{3}:CaO interface constituents, to predict the compliance components of intermediate values of 1/

*n*.

Figure 4 shows the elastic compressibility, β, and components of the elastic compliance matrix, * s*, computed for fully relaxed high-symmetry and rotation phases in the Ca

_{n+1}Ge

_{n}O

_{3n+1}series. Since all phases are (pseudo-)tetragonal, only the four symmetrically distinct

*s*

_{ij}components identified in Equation (1) are plotted. The full

*matrix, associated eigenvalues and eigenvectors and β may be found for all structures in Table S3.*

**s****Figure 4. (A)** The bulk compressibility, β and **(B**–**E)** components *s*_{ij} of the elastic compliance matrix for a tetragonal material plotted for high-symmetry and rotation phases in the Ca_{n+1}Ge_{n}O_{3n+1} series against the mole fraction of CaGeO_{3}:CaO interface (1/*n*). Interpolations between values in the CaGeO_{3} (1/*n* = 0) and Ca_{2}GeO_{4} (1/*n* = 1) structures are plotted for both high-symmetry and rotation phases (blue dot-dashed and red dotted lines, respectively). A third curve shows the interpolation between values in the $Pm\stackrel{\u0304}{3}m$ CaGeO_{3} phase and CaO rock salt structure computed as a function of CaO mole fraction 1/(*n* + 1) (cyan dashed).

The bulk elastic compressibility, β, increases linearly with 1/*n* but only very slight enhancement in compressibility is seen between the high-symmetry and rotation phases for a given *n*. Differences in β can therefore not be used to explain why uniaxial NTE is common in low *n* RP rotation phases but not in parent *I*4/*mmm* phases.

The normal compliance components, *s*_{11} and *s*_{33}, also increase with 1/*n* but unlike in β there is a significant enhancement in the rotation phase compared to the high-symmetry parent, with both the magnitude and gradient with respect to 1/*n* greater in the rotation phase.

The sign of the off-diagonal compliance components, *s*_{12} and *s*_{13}, that couple normal stresses to normal strains between axes, are negative for all compounds. This indicates that all materials have all positive Poisson ratios, ν_{ij}, where ν_{ij} describes the normal strain of axis *j* in response to a strain of axis *i* (${\nu}_{ij}=-\frac{{\epsilon}_{j}}{{\epsilon}_{i}}$). Most materials have ν_{ij} > 0, so these NTE RP phases are not auxetic (ν_{ij} < 0), even though auxetic materials have been linked with materials that exhibit anisotropic NTE (Wang et al., 2017).

Despite the negative sign, the behavior of *s*_{13} is similar to that of *s*_{11} and *s*_{33}: compliance increases with 1/*n* and there is a large enhancement in both the magnitude and the gradient increase with 1/*n* in the NTE phase compared with the high-symmetry parent. *s*_{12}, on the other hand, displays the opposite trend since the magnitude of coupling decreases with 1/*n* and going from the high-symmetry to rotation phases.

As in Figure 3, straight lines have been plotted on the subplots in Figure 4 interpolating between values for the CaGeO_{3} (1/*n* = 0) and Ca_{2}GeO_{4} (1/*n* = 1) end members. A third dashed cyan line has been added to interpolate between the pure high-symmetry CaGeO_{3} and CaO rock salt constituent phases. Because the mole fraction of CaO is actually expressed as 1/(*n* + 1) (Equation 3) these lines appear curved when plotted against the 1/*n* x-axis.

The trend in β follows that that would be predicted by modeling the RP series as a laminate of CaGeO_{3} and CaO, suggesting that β is determined predominantly by the composition. Since bulk volume thermal expansion, α, is proportional to the bulk compressibility, β, this implies that the magnitude of α is heavily dependent on chemistry. This result echoes recent work showing that experimental measurements of many thermodynamic properties of RP structures may be predicted by interpolating between values of their chemical constituents (Glasser, 2017).

Whereas, β could be approximated well as a function of CaO content for RP phases, *s*_{11} of high-symmetry phases increases above that predicted by the cyan curve. This indicates that even in the high-symmetry phase, the CaGeO_{3} and CaO layers do not behave independently and are affected by the interface between them. The prediction for *s*_{33} based on the CaO content is quite good, which may be because *s*_{33} corresponds to deformations along the layering axis (with the *a* and *b* lattice parameters fixed) and therefore the different constituent layers are being stretched in series^{1}. For both *s*_{11} and *s*_{33}, the rotation phases follow a linear relationship with 1/*n* and therefore may be considered dependent upon the fraction of CaGeO_{3}:CaO interface in the structure (red dotted line). However, in the high-symmetry phase the *n* > 1 values for both normal compliance components increase slightly beyond that predicted by interpolating between the extreme CaGeO_{3} and Ca_{2}GeO_{4} values (blue dot-dashed line). This is surprising since it is not immediately obvious how the structure of higher *n* compounds is different to local regions of CaGeO_{3} and Ca_{2}GeO_{4} and therefore what additional compliance mechanisms could operate.

For both *s*_{12} and *s*_{13}, modeling the compliance according to the mole fraction of CaO is a poor approximation, so much so that this prediction actually gives the wrong sign of the change in *s*_{12} with 1/*n*.

### 3.3. Compliance Enhancement Mechanisms

In section 3.2 we showed that certain elastic properties, such as the bulk compressibility, β, are insensitive to small changes in crystal symmetry and may be accurately predicted by interpolating between the value of β in CaGeO_{3} and CaO end member structures based on the mole fraction of CaO. However, components of the anisotropic compliance matrix, *s*_{ij}, typically differ in magnitude between high and low symmetry phases and are generally more compliant than a CaGeO_{3}:CaO interpolation predicts. Figure 5 shows the same plot as in Figure 4E (*s*_{13} vs. 1/*n*) but with annotations decomposing the *s*_{13} behavior of different structures into regimes of increasingly enhanced compliance. By separating the compliance regimes in this way, in this section we discuss the atomic displacements allowed in each regime by the phase symmetry and thus propose atomic mechanisms that may explain these enhancements in the *s*_{13} axis coupling parameter. In many cases (although not discussed here) this analysis may be used to explain the different regimes of the *s*_{11}, *s*_{33}, and *s*_{12} components in Figure 4.

**Figure 5**. The *s*_{13} term of the elastic compliance matrix (corresponding to coupling between an in-plane axis, *a* or *b*, and the layering axis, *c*) for a tetragonal material plotted for high-symmetry and rotation phases in the Ca_{n+1}Ge_{n}O_{3n+1} series against the mole fraction of CaGeO_{3}:CaO interface (1/*n*). Numbered arrows show the enhancement in compliance between different regimes represented by lines.

Taking the value of *s*_{13} that would be predicted by interpolating between values in the CaGeO_{3} and CaO constituent structures as a base (the dashed cyan curve in Figure 5), arrow ① represents an increase in the coupling between in-plane (*a* and *b*) axes and the layering axis (*c*) in the *n* = 1 Ca_{2}GeO_{4} *I*4/*mmm* phase. In a pure cubic ABO_{3} perovskite, the A cations by symmetry have the same *z* position as the apical O anions. However, the inclusion of the AO layer in high-symmetry RP phases causes symmetry breaking along *c* at the ABO_{3}:AO interface such that the apical O and interfacial A ions are no longer restricted to the same *z* coordinate, leading to a so-called “rumpling” of the AO layer. We propose that this rumpling facilitates a mechanism for enhanced *s*_{13} coupling illustrated in Figure 6A: As the in-plane, *ab*, axes are strained, the interstitial void between BO_{6} octahedra below the interfacial A cation changes in size, but the rumpling adds a degree of freedom to the *z* coordinate of the A cation which may thus move further into/out of the void in response to the in-plane strain. This thus couples the in-plane, *ab*, axes to internal displacements along the layering axis *c* and therefore to the layering axis itself.

**Figure 6**. Strain coupling mechanisms that make structures particularly compliant to cooperative strains. Since these mechanisms couple the *a* (and/or *b*) and *c* axes, they are typified by large *s*_{13} off-diagonal components of the compliance matrix. **(A)** Mechanism for increased compliance in *I*4/*mmm* RP structures where the rumpling between the interfacial A and apical O ions means that A is free to displace independently of O in response to in-plane strains; **(B)** corkscrew mechanism in RP rotation phases—the structure is able to couple in-plane tensile strains to compressive strains along *c* via internal displacements assuming all nearest neighbor B–O bonds and the two shortest A–O bonds do not deform (shown as rigid struts labeled *r*_{1}–*r*_{4}) by changing the angle of in-plane and out-of-plane hinges (described by the angles θ and α as shown); **(C)** wine rack structure where tensile strains along *a* may couple to compressive strains along *c* in a system of rigid struts in a trellis structure by only changing the angle ξ.

It was commented in the preceding section that the enhancement in compliance from the interpolation between CaGeO_{3} and Ca_{2}GeO_{4} values to *n* > 1 high-symmetry RP structures (shown by arrow ② in Figure 5) is surprising since all structures consist of a CaGeO_{3}:CaO interface (with rumpling of the CaO *z* positions) and blocks of CaGeO_{3} (that one might expect to behave as bulk cubic CaGeO_{3}). Close inspection of the *n* = 1 and *n* = 2 *I*4/*mmm* phases in Figure 1 allows us to see that in A_{2}BO_{4} structures, the length of all apical B–O bonds are equal due to the mirror symmetry plane lying in each BO_{6} layer. Similarly the angle between epitaxial B–O and apical B–O bonds must be 90° by the same reasoning. However, in A_{3}B_{2}O_{7} *I*4/*mmm* phases, this restriction that the perovskite blocks must contain a mirror plane at the centre only means that the two outer apical B–O bond lengths and BO_{6} internal angles must be equal and the apical B–O bond lengths between the two inner BO_{6} layers must be equal. There is no restriction by symmetry that all apical B–O bond lengths must be equal or in fact that all BO_{6} internal angles must be 90°. These weaker restrictions create internal degrees of structural freedom that may facilitate greater compliance, since there is greater freedom for the atoms to move in response to external strains. In structures relaxed using DFT, we find that there are slight differences in these two bondlengths: 1.87 and 1.90 Å for the outer and inner apical B–O bonds, respectively, and the angle between outer apical B–O and epitaxial B–O bonds is 91.2°. This same argument may be applied to all *n* > 1 *I*4/*mmm* phases. However, in the special case of the *n* = ∞ series end member, there is no AO layer to break any symmetry along *c* and thus all B–O bond lengths are equal and all O–B–O angles are 90°. That this difference between *s*_{13} of *n* = 2, 3, 4 high-symmetry phases and the CaGeO_{3}:Ca_{2}GeO_{4} interpolation is not also seen in the significantly more compliant rotation phases may be because in these rotation phases, additional compliance mechanisms operate that dwarf the effect described by arrow ②.

Arrow ③ in Figure 5 represents the increase in coupling between the *a*/*b* and *c* axes from the cubic CaGeO_{3} phase to an *I*4/*mcm* phase with anti-phase octahedral rotation about the *c* axis. From a structural symmetry point of view, this may again come as a surprise since the *I*4/*mcm* phase only has an internal degree of freedom in the *ab* plane and not along the *c* axis that this plane couples to. Therefore, one might expect large changes in in-plane strain in response to biaxial stress from rigid BO_{6} octahedra rotating, but not also large changes in *c*. However, DFT studies on LaAlO_{3} (Hatt and Spaldin, 2010) and LaNiO_{3} (Weber et al., 2016) show that while the application of compressive biaxial strain to *I*4/*mcm* phases causes the rotation angle to increase, it also leads to a large tetragonal distortion of the BO_{6} units with compression of the epitaxial B–O bonds and extension of the apical B–O bonds.

Finally, arrow ④ in Figure 5 represents the increase in *s*_{13} coupling between rotation phases of CaGeO_{3} perovskite and RP phases with a frozen octahedral rotation. In these phases, the rumpling of the interfacial AO layer identified in the high-symmetry structure is still present. However, whereas in the high-symmetry structure, in-plane strains necessarily involved deformations of stiff epitaxial B–O bonds, in this lower-symmetry phase, the frozen in-plane octahedral rotation adds an internal degree of freedom in the in-plane epitaxial O positions. There are thus internal degrees of freedom in both in-plane and layering axes in RP rotation phases. In our previous paper (Ablitt et al., 2017) we proposed a “corkscrew” mechanism to explain high *s*_{13} coupling in RP1 rotation phases but not in *I*4/*mcm* ABO_{3} phases. This coupling mechanism in theory allows *a* and *c* to deform cooperatively without extending the four most stiff cation-anion bonds identified in the *I*4_{1}/*acd* Ca_{2}GeO_{4} structure, by only changing two bond angles, labeled θ and α in Figure 6B. We call this mechanism corkscrew since an in-plane rotation leads to an extension along the rotation axis: Figure 6B shows how the in-plane rotation angle θ may decrease in response to a biaxial expansion, such that this in turn pulls the stiff O-A bond (shown as a rigid rod), decreasing the angle α and thus increasing the rumpling of the rock salt layer and forcing contraction along the *c*-axis.

The net result of these interfacial strain coupling mechanisms in RP phases, yielding an enhanced off-diagonal *s*_{13} compliance term, is rather reminiscent of the “wine-rack” mechanism such as that which operates in methanol monohydrate (Fortes et al., 2011). In the wine-rack trellis structure, shown schematically in Figure 6C with rigid struts but flexible hinges, *a* couples strongly to *c*, mediated by the hinge angle, ξ, in the *ac* plane. Following the method used to analyse the wine-rack (Grima and Evans, 2000) and other idealized geometries (Smith et al., 2000, Grima et al., 2012), we were able to derive the mechanical properties that our pure corkscrew mechanism would exhibit under the conditions that the four *stiff* bonds shown in Figure 6B indeed remained rigid and all resistance to strain came from a harmonic potential in the θ and α-hinges (Ablitt et al., 2018). Under these restrictions, a corkscrew model would have an *s*_{13} compliance parameter as a function of RP layer thickness, *n*, given by

where *k*_{θ} and *k*_{α} represent the harmonic stiffness of the θ and α-hinges, respectively, and f(θ, α) is an expression of trigonometric functions of θ and α

Under the constraints of this model, α is explicitly dependent on θ by the equation

and therefore given the bond lengths *r*_{1}-*r*_{4} (as defined in Figure 6B), the value of θ fully determines the structure.

In the limit that *k*_{α} >> *k*_{θ}, the AO interface is stiffer than the ABO_{3} perovskite blocks and *s*_{13} loses dependence on *n*. However, in the opposing limit that *k*_{θ} >> *k*_{α} and changing the in-plane rotation angle is the main obstacle to strain, *s*_{13}∝1/*n*. The real system is closer to the *k*_{θ} >> *k*_{α} limit as we observe linear behavior of *s*_{13} with 1/*n*. Although this model has been derived assuming that all equatorial B–O bond lengths (*r*_{1}) and all apical B–O bond lengths (*r*_{2}) are equal, the bond lengths *r*_{2} and *r*_{3} do not feature in Equation (4) and *r*_{1} and θ only need refer to the in-plane bond lengths and rotation angle in the outer layer of the perovskite block. Therefore, this mechanism is compatible with a distribution of possible bond lengths and rotation angles in different layers of a perovskite block if *n* > 1.

We accept the limitations of such simple models but developing the study of how symmetry-allowed local distortions can give rise to new compliance mechanisms in different crystallographic phases, such as those identified using arrows ①–④ in Figure 5, may allow for the prediction of phases with high cross compliances by symmetry alone, before needing to explicitly compute the elastic constants. Calculations of elastic constants are, in turn, frequently less expensive than performing lattice dynamics calculations across the full Brillouin zone to compute the thermodynamic driving force for anisotropic thermal expansion, Φ.

Therefore, looking for other materials with such high cross compliances, using symmetry analysis as a guide to narrow the pool of structures, may prove a more general method for searching for novel NTE materials. Indeed, by considering this analysis and our thermodynamic criteria requiring a proximity to a competing phase transition to provide Φ, we have already been able to identify (Ablitt et al., 2017) the layered double perovskites Sr_{2}MgWO_{6} (Achary et al., 2006) and (110)-cut perovskite LaTaO_{4} (Cordrey et al., 2015) which fall within this general paradigm. Substantial research opportunities exist in this area to more fully explore the NTE behavior of these classes of materials.

### 3.4. Dynamic Driving Force for NTE

So far we have not addressed the thermodynamic driving force for thermal expansion Φ(*T*), and it is rather more computationally expensive to calculate than κ as it requires the full phonon density of states to be computed for different strained structures within the QHA. Furthermore, this procedure is only possible if all phonons within the NTE structure simulated using DFT have real frequencies—suggesting that the phase must be stable at 0 K. For *n* > 1 structures in the Ca_{n+1}Ge_{n}O_{3n+1} series, this latter condition is not met, and therefore full computation of Φ(*T*) as we performed previously (Ablitt et al., 2017) would not be possible. However, we also previously identified that the most important modes driving NTE are octahedral tilts about an in-plane axis. Thus in Figure 7 we have plotted the frequency of the lowest frequency tilt mode to occur at a high-symmetry q-point in high-symmetry phases and phases with a frozen in-plane rotation. We note that these tilts are not the same as the vibrations of ions along the *z* direction predicted using a QHA-inspired method to cause NTE in a *A*2_{1}*am* Ca_{3}Ti_{2}O_{7} phase after 30 GPa hydrostatic pressure has been applied and in which these octahedral tilts are already frozen (Huang et al., 2016). As in previous studies (Senn et al., 2016, Ablitt et al., 2017) the Ca_{n+1}Ge_{n}O_{3n+1} series is being treated as an analog to Ca_{n+1}Mn_{n}O_{3n+1} for comparison against experimental data to avoid expensive magnetic calculations since Ge^{4+} and Mn^{4+} are known to have equal ionic radii (Shannon, 1976). The authors have previously shown that this substitution has little bearing on the properties (phonon frequencies, elastic constants) relevant to modeling thermal expansion within the QHA.

**Figure 7**. The lowest frequency harmonic phonon mode with octahedral tilt character found in DFT simulations at high-symmetry q-points in the Ca_{n+1}Ge_{n}O_{3n+1} series against interface fraction, 1/*n*. Tilt frequencies have been plotted in the high-symmetry and rotation phases with 0 GPa external pressure, and also for the rotation phase with a 4.3 GPa biaxial pressure applied to give better agreement of the Ca_{n+1}Ge_{n}O_{3n+1} structure to Ca_{n+1}Mn_{n}O_{3n+1}. DFT only simulates the system at 0 K whereas phonon frequencies may harden with increased temperature. This means that many octahedral tilts are predicted with imaginary frequencies (shown as negative) even if the mode has a real frequency and the structure is stable at some higher temperature. To give an idea of how the equilibrium phase changes with temperature, the inset graphs plot the experimentally observed phase diagram for each *n* in the Ca_{n+1}Mn_{n}O_{3n+1} system. The data to make these illustrations was taken from: Ca_{2}MnO_{4} (Takahashi and Kamegashira, 1993); Ca_{3}Mn_{2}O_{7} (Senn et al., 2015); Ca_{4}Mn_{3}O_{10} (Battle et al., 1998); CaMnO_{3} (Taguchi et al., 1989).

It can be seen from Figure 7 that in all phases, the frequency of the softest tilt mode stiffens with increased Ca_{2}GeO_{4} content, indicating that the inclusion of the CaGeO_{3}:CaO interface reduces the propensity of octahedra to tilt. The tilt mode is unstable (has an imaginary frequency) in all high-symmetry parent phases, which is unsurprising since these phases are not observed experimentally at low temperatures at any *n* for the Ca_{n+1}Mn_{n}O_{3n+1} series. The tilt stiffens between the parent and child phases for all *n*, showing that the dynamic tilt couples competitively with the frozen in-plane rotation. In both the high- and low-symmetry phases, the (imaginary) octahedral tilt frequency varies approximately linearly with Ca_{2}GeO_{4} mole fraction, although the stiffening effect going between the high- and low- symmetry phases becomes greater at higher *n* (lower Ca_{2}GeO_{4} fraction). The CaO rock salt layer also stiffens octahedral rotations and thus the angle of the frozen rotation increases with *n* in rotation phases from 10.55° in *I*4_{1}/*acd* Ca_{2}GeO_{4} (*n* = 1) up to the limit of 12.7° in *I*4/*mcm* CaGeO_{3} (*n* = ∞; see Figure S2). Hence, the greater amplitude of frozen rotation observed in CaGeO_{3} means that the competitive coupling between the frozen rotation and dynamic tilt is greatest for this largest *n* member, and leads to further hardening of the dynamic tilt frequency.

In the 0 GPa relaxed Ca_{n+1}Ge_{n}O_{3n+1} structures, the tilt is still unstable for all *n*. However, it was found previously (Ablitt et al., 2017) that a discrepancy arises between the Ca_{2}GeO_{4} and Ca_{2}MnO_{4} in-plane lattice parameters in the *I*4_{1}/*acd* phase—due to in-plane magneto-strictive coupling—but that applying a biaxial 4.3 GPa pressure corrects for this small difference, yielding very close agreement in the frequencies of the softest phonon modes. After applying this biaxial pressure, the *n* = 1 *I*4_{1}/*acd* Ca_{2}GeO_{4} rotation phase has all real mode frequencies, but the softest tilt in all structures with *n* > 1 is still unstable. This agrees with experimental observation of the low temperature stable phases, shown as inset figures for each composition on Figure 7.

Ca_{3}Mn_{2}O_{7} is found at low temperature in the improper ferroelectric *A*2_{1}*am* phase, but undergoes a wide temperature region of phase coexistence with the uniaxial NTE *Acaa* phase between 150 and 280 K on cooling and 300 and 360 K on heating (Senn et al., 2015). The strong first order nature of the phase transition arises because frozen octahedral rotations in the *A*2_{1}*am* and *Acaa* phases have opposite sense (in-phase vs. out-of-phase about *c* within each perovskite layer—see Figure A1 in Appendix 2) but the approximate transformation temperature gives an indication of the temperature at which octahedral tilts in a rotation phase develop real frequencies. Furthermore, the *n* = ∞ perovskite, CaMnO_{3}, transforms around 1166 K from a *Pnma* ground state with frozen octahedral tilts, to an *I*4/*mcm* phase with only an out-of-phase octahedral rotation—the *n* = ∞ analog of the *n* = 2 *Acaa*—remaining in this phase for only a brief temperature window before transforming again to the cubic parent structure at 1,184 K (Taguchi et al., 1989). For *n* = 3, the reported symmetry for Ca_{4}Mn_{3}O_{10} from 5 K up until room temperature is *Pbca* (Battle et al., 1998), in which a static rotation and tilt of the octahedra are present, and, although to the best of our knowledge a transformation to a higher-symmetry phase with only a frozen octahedral rotation has not yet been reported^{2}, interpolating the experimental observations in Figure 7 predicts a transformation temperature in the 250–1,100 K window. The magnitude of the imaginary tilt frequency computed in Figure 7 may therefore by interpreted as a crude indicator of the stability of the structure with condensed rotation and tilt and thus of the temperature required to transform to the higher-symmetry rotation phase.

As well as having optimal elastic anisotropy to facilitate uniaxial NTE, the *n* = 1 Ca_{2}MnO_{4} *I*4_{1}/*acd* phase has soft tilt modes at low temperatures providing a thermodynamic driving force for cooperative in-plane positive and out-of-plane NTE. This feature is not unique to the *n* = 1 structure however, since at some higher temperature all Ca_{n+1}Mn_{n}O_{3n+1} compounds should transform to a phase in which the tilt frequencies are real and soft, at least over some temperature range. Furthermore, we have demonstrated previously in the Ca_{3−x}Sr_{x}Mn_{2}O_{7} system, that for a given layer thickness *n*, chemical substitution (changing *x*) may be used to alter the Goldschmidt tolerance factor and thus the frequencies of these octahedral tilts, switching between positive and negative uniaxial thermal expansion (Senn et al., 2016). Therefore, Figure 7 shows that the RP structure, through the layer thickness *n*, determines a ballpark value for the 0 K tilt frequency in the uniaxial NTE phase, and thus a ballpark temperature window in which the NTE phase will be stable. Hence, given this structural constraint we may use chemical control to optimize the proximity of the structure to the phase transition and so enhance the NTE. In this respect, the *n* = 1 family member, since one may presume that it may be tuned chemically to arbitrary proximity to the competing phase transition, is the most promising candidate for exhibiting the largest NTE on account of the anisotropy with respect to its most and least compliant directions which we have shown is maximized for this system.

## 4. Conclusions

We have shown that the elastic anisotropy ratio, κ, found previously to be an essential ingredient for uniaxial NTE, increases linearly in the RP Ca_{n+1}Ge_{n}O_{3n+1} series (*n* = 1, 2, 3, 4 …∞) with the CaGeO_{3}:CaO content (expressed by the ratio 1/*n*), reaching a maximum in the structure with maximal interface (*n* = 1). By decomposing the components of the elastic compliance matrix for high-symmetry and NTE phases (with a frozen octahedral rotation about the layering axis) into different regimes that show similar trends with 1/*n*, we have been able to link these regimes with internal degrees of freedom in the structure that allow atomic mechanisms to operate that couple to cell strains. The most important of these is the “corkscrew” mechanism that operates locally at the CaGeO_{3}:CaO interface in phases with a frozen octahedral rotation about the layering axis and therefore explains the trend that anisotropic compliance correlates with the fraction of interface in these phases. This local atomic compliance mechanism is analogous in certain ways to the wine-rack mechanism that operates in many much softer framework materials. The compliance matrices can be rapidly calculated by DFT methods and diagonalized to assess them for cross coupling terms that promote pronounced uni or biaxial NTE. This makes them suitable descriptors for high throughput computational searching for novel NTE materials, especially when symmetry constraints may be employed to narrow the space of candidate phases.

We further investigated the trend in frequency of octahedral tilts with RP layer thickness and found that the 0 K tilt frequencies in NTE or analogous structures become softer with increasing *n*. This implies that a window of stability of the NTE phase with soft active tilt modes exists at increasingly higher temperatures with higher *n*. We had previously shown that the thermodynamic driver for NTE for a given *n* might be tuned with chemical substitution, and we now show that the anisotropic compliance necessary for NTE in these systems is maximized for a high fraction of CaGeO_{3}:CaO interface layers in the structure. On the basis of this analysis, we thus predict that the *n* = 1 systems, such as Ca_{2}MnO_{4}, will be the RP systems in which the maximum NTE can be achieved via chemical substitution.

## Data Availability Statement

Data underlying this article can be accessed on figshare at DOI:10.6084/m9.figshare.6729287, and used under the Creative Commons Attribution licence.

## Author Contributions

CA performed the calculations and data analysis. All authors contributed to the design of the study and the analysis and interpretation of the results. The paper was drafted by CA and MS, and all authors contributed to its development into final form.

## Funding

CA is supported by a studentship in the Centre for Doctoral Training on Theory and Simulation of Materials at Imperial College London funded by the EPSRC (EP/L015579/1). Calculations were performed on the Imperial College London high-performance computing facility. This work was supported by the Thomas Young Centre under grant TYC-101. We are grateful to the UK Materials and Molecular Modelling Hub for computational resources, which is partially funded by EPSRC (EP/P020194/1). MS acknowledges the Royal Society for a fellowship.

## Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2018.00455/full#supplementary-material

## Footnotes

1. ^However, we note that under these conditions one should use a Voigt average to interpolate *s*_{33}, whereas by drawing a straight line of compliance vs mole fraction we have actually performed a Reuss interpolation. In Figure S1 in the SI we show that a Voigt fit against CaO mole fraction is actually less close.

2. ^Although we do note that the high-symmetry *I*4/*mmm* parent structure has been stabilized at room temperature from solid state synthesis at high temperatures and pressures (Yu et al., 2001).

## References

Ablitt, C., Craddock, S., Senn, M. S., Mostofi, A. A., and Bristowe, N. C. (2017). The origin of uniaxial negative thermal expansion in layered perovskites. *NPJ Comput. Mater. * 3:44. doi: 10.1038/s41524-017-0040-0

Ablitt, C., Senn, M. S., Bristowe, N. C., and Mostofi, A. A. (2018). A corkscrew model for highly coupled anisotropic compliance in Ruddlesden-Popper oxides with frozen octahedral rotations. *arXiv:1810.02697*.

Achary, S. N., Chakraborty, K. R., Patwe, S. J., Shinde, A. B., Krishna, P. S. R., and Tyagi, A. K. (2006). Anisotropic thermal expansion behavior in tetragonal Sr_{2}MgWO_{6}. *Mater. Res. Bull. * 41, 674–682. doi: 10.1016/j.materresbull.2005.07.025

Azuma, M., Chen, W. T., Seki, H., Czapski, M., Olga, S., Olga, K., et al. (2011). Colossal negative thermal expansion in BiNiO_{3} induced by intermetallic charge transfer. *Nat. Commun.* 2:347. doi: 10.1038/ncomms1361

Barron, T. H. K., and Munn, R. W. (1967). Analysis of the thermal expansion of anisotropic solids: application to zinc. *Philos. Magaz.* 15:85

Battle, P. D., Green, M. A., Lago, J., Millburn, J. E., Rosseinsky, M. J., and Vente, J. F. (1998). Crystal and magnetic structures of Ca_{4}Mn_{3}O_{10}, an *n* = 3 Ruddlesden-Popper compound. *Chem. Mater. * 10, 658–664. doi: 10.1021/cm970647r

Benedek, N. A., and Fennie, C. J. (2011). Hybrid improper ferroelectricity: a mechanism for controllable polarization-magnetization coupling. *Phys. Rev. Lett.* 106:107204. doi: 10.1103/PhysRevLett.106.107204

Benedek, N. A., Rondinelli, J. M., Djani, H., Ghosez, P., and Lightfoot, P. (2015). Understanding ferroelectricity in layered perovskites: new ideas and insights from theory and experiments. *Dalton Trans.* 44, 10543–10558. doi: 10.1039/C5DT00010F

Borzi, R. A., Grigera, S. A., Farrell, J., Perry, R. S., Lister, S. J. S., Lee, S. L., et al. (2007). Formation of a nematic fluid at high fields in Sr_{3}Ru_{2}O_{7}. *Science* 315, 214–217. doi: 10.1126/science.1134796

Chen, J., Wang, F., Huang, Q., Hu, L., Song, X., Deng, J., et al. (2013). Effectively control negative thermal expansion of single-phase ferroelectrics of PbTiO_{3}-(Bi,La)FeO_{3} over a giant range. *Sci. Rep.* 3:2458. doi: 10.1038/srep02458

Clark, S. J., Segall, M. D., Pickard, C. J., Hasnip, P. J., Probert, M. J., Refson, K., et al. (2005). First principles methods using CASTEP. *Z. Kristall.* 220, 567–570. doi: 10.1524/zkri.220.5.567.65075

Cordrey, K. J., Stanczk, M., Dixon, C. A. L., Knight, K. S., Gardner, J., Morrison, F. D., et al. (2015). Structural and dielectric studies of the phase behaviour of the topological ferroelectric La_{1−x}Nd_{x}TaO_{4}. *Dalton Trans.* 44, 10673–10680. doi: 10.1039/C4DT03721A

Dove, M. T., Heine, V., and Hammond, K. D. (1995). Rigid unit modes in framework silicates. *Min. Mag.* 59, 629–639.

Dwivedi, A., and Cormack, A. N. (1991). Crystal chemistry of ruddlesden-popper type structures in high Tc ceramic superconductors. *Bull. Mater. Sci. * 14, 575–584. doi: 10.1007/BF02744688

Fortes, A. D., Suard, E., and Knight, K. S. (2011). Negative linear compressibility and massive anisotropic thermal expansion in methanol monohydrate. *Science* 331, 742–746. doi: 10.1126/science.1198640

Glasser, L. (2017). Systematic thermodynamics of layered perovskites: Ruddlesden–Popper phases. *Inorgan. Chem.* 56, 8920–8925. doi: 10.1021/acs.inorgchem.7b00884

Grima, J. N., Caruana-Gauci, R., Attard, D., and Gatt, R. (2012). Three-dimensional cellular structures with negative Poisson's ratio and negative compressibility properties. *Proc. R. Soc. Lond. A* 468, 3121–3138. doi: 10.1098/rspa.2011.0667

Grima, J. N., and Evans, K. E. (2000). Auxetic behavior from rotating squares. *J. Mater. Sci. Lett. * 19, 1563–1565. doi: 10.1023/A:1006781224002

Haeni, J. H., Theis, C. D., Schlom, D. G., Tian, W., Pan, X. Q., Chang, H., et al. (2001). Epitaxial growth of the first five members of the Srn+1TinO_{3}n+1 Ruddlesden-Popper homologous series. *Appl. Phys. Lett. * 78, 3292–3294. doi: 10.1063/1.1371788

Hatt, A. J., and Spaldin, N. A. (2010). Structural phases of strained LaAlO_{3} driven by octahedral tilt instabilities. *Phys. Rev. B.* 82:195402. doi: 10.1103/PhysRevB.82.195402

Heine, V., Welche, P. R. L., and Dove, M. T. (1999). Geometrical origin and theory of negative thermal expansion in framework structures. *J. Am. Ceramic Soc. * 82, 1793–1802. doi: 10.1111/j.1151-2916.1999.tb02001.x

Huang, L.-F., Lu, X.-Z., and Rondinelli, J. M. (2016). Tunable negative thermal expansion in layered perovskites from quasi-two-dimensional vibrations. *Phys. Rev. Lett.* 117:115901. doi: 10.1103/PhysRevLett.117.115901

Lee, C.-H., Podraza, N. J., Zhu, Y., Berger, R. F., Shen, S., Sestak, M., et al. (2013). Effect of reduced dimensionality on the optical band gap of SrTiO_{3}. *Appl. Phys. Lett. * 102:122901. doi: 10.1063/1.4798241

Mackenzie, A. P., and Maeno, Y. (2003). The superconductivity of Sr_{2}RuO_{4} and the physics of spin-triplet pairing. *Rev. Mod. Phys.* 75, 657–712. doi: 10.1103/RevModPhys.75.657

McCoy, M. A., Grimes, R. W., and Lee, W. E. (1997). Phase stability and interfacial structures in the SrO–SrTiO_{3} system. *Philos. Magaz. A* 75, 833–846. doi: 10.1080/01418619708207205

Palgrave, R. G., Borisov, P., Dyer, M. S., McMitchell, S. R. C., Darling, G. R., Claridge, J. B., et al. (2012). Artificial construction of the layered Ruddlesden-Popper manganite La_{2}Sr_{2}Mn_{3}O_{1}0 by reflection high energy electron diffraction monitored pulsed laser deposition. *J. Am. Chem. Soc.* 134, 7700–7714. doi: 10.1021/ja211138x

Perdew, J. P., Ruzsinszky, A., Csonka, G. I., Vydrov, O. A., Scuseria, G. E., Constantin, L. A., et al. (2008). Restoring the density-gradient expansion for exchange in solids and surfaces. *Phys. Rev. Lett.* 100:136406. doi: 10.1103/PhysRevLett.100.136406

Qi, T. F., Korneta, O. B., Parkin, S., De Long, L. E., Schlottmann, P., and Cao, G. (2010). Negative volume thermal expansion via orbital and magnetic orders in Ca_{2}Ru_{1−x}Cr_{x}O_{4} (0 < *x* < 0.13). *Phys. Rev. Lett.* 105:177203. doi: 10.1103/PhysRevLett.105.177203

Refson, K., Tulip, P. R., and Clark, S. J. (2006). Variational density-functional perturbation theory for dielectrics and lattice dynamics. *Phys. Rev. B* 73:155114. doi: 10.1103/PhysRevB.73.155114

Sarlosi, J., and Bocko, J. (2016). Effective material moduli for composites. *Am. J. Mech. Eng.* 4, 289–292. doi: 10.12691/ajme-4-7-11

Senn, M. S., Bombardi, A., Murray, C. A., Vecchini, C., Scherillo, A., Lui, X., et al. (2015). Negative thermal expansion in hybrid improper ferroelectric Ruddlesden-Popper perovskites by symmetry trapping. *Phys. Rev. Lett.* 114:035701. doi: 10.1103/PhysRevLett.114.035701

Senn, M. S., Murray, C. A., Luo, X., Wang, L., Huang, F.-T., Cheong, S.-W., et al. (2016). Symmetry switching of negative thermal expansion by chemical control. *J. Am. Chem. Soc.* 138, 5479–5482. doi: 10.1021/jacs.5b13192

Shannon, R. D. (1976). Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides. *Acta Crystallogr.* A32, 751–767. doi: 10.1107/S0567739476001551

Smith, C., Grima, J., and Evans, K. (2000). A novel mechanism for generating auxetic behaviour in reticulated foams: missing rib foam model. *Acta Mater.* 48, 4349–4356. doi: 10.1016/S1359-6454(00)00269-X

Taguchi, H., Nagao, M., Sato, T., and Shimada, M. (1989). High-temperature phase transition of CaMnO_{3}-δ. *J. Solid State Chem.* 78, 312–315. doi: 10.1016/0022-4596(89)90113-8

Takahashi, J., and Kamegashira, N. (1993). X-ray structural study of calcium manganese oxide by rietveld analysis at high temperatures. *MRS Bull.* 28, 565–573.

Wang, L., Luo, H., Deng, S., Sun, Y., and Wang, C. (2017). Uniaxial negative thermal expansion, negative linear compressibility, and negative Poisson's ratio induced by specific topology in Zn[Au(Cn)_{2}]_{2}. *Inorgan. Chem.* 56, 15101–15109. doi: 10.1021/acs.inorgchem.7b02416

Weber, M. C., Guennou, M., Dix, N., Pesquera, D., Sánchez, F., Herranz, G., et al. (2016). Multiple strain-induced phase transitions in LaNiO_{3} thin films. *Phys. Rev. B.* 94:014118. doi: 10.1103/PhysRevB.94.014118

Yan, L., Niu, H. J., Duong, G. V., Suchomel, M. R., Bacsa, J., Chalker, P. R., et al. (2011). Cation ordering within the perovskite block of a six-layer ruddlesden-popper oxide from layer-by-layer growth - artificial interfaces in complex unit cells. *Chem. Sci.* 2, 261–272. doi: 10.1039/C0SC00482K

Yoshida, Y., Ikeda, S.-I., Matsuhata, H., Shirakawa, N., Lee, C. H., and Katano, S. (2005). Crystal and magnetic structure of Ca_{3}Ru_{2}O_{7}. *Phys. Rev. B* 72:054412. doi: 10.1103/PhysRevB.72.054412

Yoshizawa, H., Kakeshita, T., Kajimoto, R., Tanabe, T., Katsufuji, T., and Tokura, Y. (2000). Stripe order at low temperatures in La_{2−x}Sr_{x}NiO_{4} with 0.289≲*x*≲0.5. *Phys. Rev. B* 61, R854–R857. doi: 10.1103/PhysRevB.61.R854

Keywords: NTE, perovskite, Ruddlesden-Popper, anisotropy, compliance, corkscrew, thermal expansion

Citation: Ablitt C, Mostofi AA, Bristowe NC and Senn MS (2018) Control of Uniaxial Negative Thermal Expansion in Layered Perovskites by Tuning Layer Thickness. *Front. Chem*. 6:455. doi: 10.3389/fchem.2018.00455

Received: 03 July 2018; Accepted: 13 September 2018;

Published: 18 October 2018.

Edited by:

Andrea Sanson, Università degli Studi di Padova, ItalyReviewed by:

Qiang Sun, Zhengzhou University, ChinaLei Hu, University of Science and Technology Beijing, China

Copyright © 2018 Ablitt, Mostofi, Bristowe and Senn. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Mark S. Senn, m.senn@warwick.ac.uk