## Abstract

The Mott metal-insulator transition, a paradigm of strong electron-electron correlations, has been considered as a source of intriguing phenomena. Despite its importance for a wide range of materials, fundamental aspects of the transition, such as its universal properties, are still under debate. We report detailed measurements of relative length changes Δ*L*/*L* as a function of continuously controlled helium-gas pressure *P* for the organic conductor κ-(BEDT-TTF)_{2}Cu[N(CN)_{2}]Cl across the pressure-induced Mott transition. We observe strongly nonlinear variations of Δ*L*/*L* with pressure around the Mott critical endpoint, highlighting a breakdown of Hooke’s law of elasticity. We assign these nonlinear strain-stress relations to an intimate, nonperturbative coupling of the critical electronic system to the lattice degrees of freedom. Our results are fully consistent with mean-field criticality, predicted for electrons in a compressible lattice with finite shear moduli. We argue that the Mott transition for all systems that are amenable to pressure tuning shows the universal properties of an isostructural solid-solid transition.

- Strongly correlated electron systems
- Mott metal-insulator transiton
- critical phenomena
- thermodynamic studies
- thermal expansion
- organic charge-transfer salts
- effects of hydrostatic pressure
- coupling of electrons to lattice degrees of freedom

## INTRODUCTION

The Mott metal-insulator transition, driven by strong Coulomb repulsion between electrons (*1*), is a key phenomenon in strongly correlated electron systems. Prime examples for materials where Mott physics is relevant include transition metal oxides (*2*), notably cuprates and manganites, as well as organic charge-transfer salts (*3*, *4*). The strong interest shown to these materials relates to the wealth of intriguing states that was observed to be close to the transition on either side, with spin-liquid phases (*5*, *6*), multiferroicity (*7*, *8*), and unconventional superconductivity (*9*–*11*) being among the most prominent ones. Likewise, the character of the transition itself is of special importance because it is the origin of emergent behavior. The basic notion is that for systems with approximately one electron per lattice site, the electrons become localized once the ratio of Coulomb repulsion *U*, for two electrons on the same site, over the kinetic energy *W* exceeds a critical value (*U*/*W*)_{c}. Experiments have demonstrated that the Mott transition can be controlled either by doping [as, for example, in the cuprates (*10*)] or by applying pressure *P* [organic charge-transfer salts (*3*)], giving rise to a first-order transition line. This line terminates in a second-order critical endpoint (Fig. 1A), which is expected to be characterized by universal properties.

The universal properties of the Mott transition continue to be a subject of great interest. For the idealized case of a purely electronic system, as, for example, simulated in cold-atom experiments (*12*), there is strong evidence from theory that the Mott transition shows the universal properties of a liquid-gas transition (*13*, *14*), represented by the Ising universality class, with the double occupancy playing the role of an order parameter. The crucial question is whether this simplified picture can adequately describe the properties of a real material where the critical electronic system is coupled to a compressible lattice. To date, the involvement of the elastic degrees of freedom in the Mott criticality has not been elucidated experimentally. This is remarkable, because the importance of the lattice for the Mott transition can be directly inferred from the fact that the transition can be tuned by pressure, that is, by lattice contraction, which increases the overlap of electronic wave functions on adjacent sites and thereby promotes metallicity. Conversely, the critical electronic system exerts an internal pressure on the elastic system to which the crystal lattice responds. Sufficiently away from the Mott endpoint, the lattice response is weak and acts as a small perturbation on the critical electronic system. However, upon approaching the critical endpoint, the lattice necessarily reacts in a nonperturbative way to the internal strain, implying a vanishing elastic modulus (*15*–*17*), as also shown for the case of piezoelectric ferroelectric criticality (*18*). The fingerprint of this effect would be a strongly nonlinear strain-stress relation, that is, a breakdown of Hooke’s law of elasticity. Because of the long-ranged shear forces of the crystal lattice, this breakdown is expected to be accompanied by a crossover to Landau criticality characterized by mean-field critical exponents (Fig. 1A). The remarkable consequence would be that, although the Mott transition is driven by strong electronic correlations, the critical behavior around the endpoint is eventually governed by the universal properties of a solid-solid transition (*17*, *19*) with the strain field becoming a primary order parameter. Measurements of the strain are thus an ideal tool to probe the criticality of the finite-temperature Mott endpoint.

## RESULTS

Here, we use measurements of relative length changes for directly probing the lattice effects around the Mott transition. By the combination of capacitive dilatometry with continuously controllable helium-gas pressure (*20*), we can perform high-resolution measurements of relative length changes, Δ*L*/*L*, both as a function of temperature at constant pressure and as a function of pressure at constant temperature. The material selected for our study is the quasi–two-dimensional (quasi-2D) organic charge-transfer salt κ-(BEDT-TTF)_{2}Cu[N(CN)_{2}]Cl (abbreviated as κ-Cl hereafter). The reason for choosing this material is twofold. First, this material is situated right on the verge of the Mott transition (*21*, *22*) so that moderate pressure *P* ≲ 30 MPa is sufficient to cover a wide-enough range around the critical endpoint. Second, for this material, a considerably large elasticity-governed critical regime around the critical endpoint has been predicted (*17*). The present study, which focuses on the critical elasticity, goes beyond previous investigations on the critical properties of this and other materials via conductivity (*23*–*25*) and nuclear magnetic resonance (*26*) measurements, where only electronic degrees of freedom have been considered.

The κ-Cl system crystallizes in the orthorhombic *Pnma* structure (*27*) in which the organic BEDT-TTF layers, where BEDT-TTF stands for bis(ethylenedithio)tetrathiafulvalene, alternate with the polymeric anion layers along the *b* axis (Fig. 1B); see section A in the Supplementary Materials for the effect of pressure on the low-temperature structure. This gives rise to a quasi-2D electronic structure with one hole carrier per *ac* plane (Fig. 1C), establishing a half-filled conduction band. Measurements have been performed on three different single crystals from different sources, yielding essentially the same results.

### Relative length changes across the Mott transition

In Fig. 2, we show results of relative length changes Δ*L*_{i}/*L*_{i} measured along the in-plane *i* = *a* axis (A) and the out-of-plane *b* axis (B) on the same single crystal as a function of pressure *P* at varying temperatures around the Mott transition. These two uniaxial effects represent the dominant contributions of the Mott transition to the volume expansivity, (Δ*V*/*V*)_{Mott} ≈ (Δ*L*_{a}/*L*_{a})_{Mott} + (Δ*L*_{b}/*L*_{b})_{Mott}, because there is practically no effect visible along the *c* axis, (Δ*L*_{c}/*L*_{c})_{Mott} ≈ 0 (see fig. S1). The data in Fig. 2 disclose two essential aspects. The first aspect concerns the character of the transition and how this character changes with increasing temperature.

For the lowest temperature studied at 30 K, which is well below the critical temperature *T*_{c} ≈ 36.5 K (see below for the determination of *T*_{c}), we observe an abrupt, slightly broadened jump of the sample length as a function of pressure along both axes, reflecting the first-order character of the Mott transition. The data reveal a distinct decrease of the *a* and *b* axes lattice parameters while tuning the system from the low-pressure insulating side of the transition to the high-pressure metallic side. This observation is consistent with the notion that the itinerant electrons on the metallic side provide an additional contribution to the cohesion of the solid and supports the notion that (Δ*V*/*V*)_{Mott} is proportional to the order parameter of the transition. For pressures distinctly below and above the jump, the sample length is changing in an approximately linear manner with pressure, reflecting normal elastic behavior following Hooke’s law of elasticity. This elastic background exhibits the usual negative slope for the *b* axis lattice parameter, that is, a shrinking of the lattice with increasing pressure, whereas it shows an unusual positive slope for the *a* axis lattice parameter. We assign this behavior to the Poisson effect as a result of a pressure-induced progressive inclination of the BEDT-TTF molecules with respect to the *b* axis (Fig. 1B). Upon increasing the temperature, the discontinuity gradually decreases and evolves into a continuous crossover behavior at higher temperatures.

The second important aspect contained in Fig. 2 (A and B) relates to the slopes of the Δ*L*_{i}(*P*)/*L*_{i} curves near the critical endpoint, reflecting the material’s uniaxial compressibilities κ_{i} = −*d*(Δ*L*_{i}/*L*_{i})/*dP* along the *i* = *a* and *b* axes. For temperatures around *T*_{c} ≈ 36.5 K up to at least 43 K, we observe highly nonlinear variations of Δ*L*_{i}/*L*_{i} with pressure for both axes over a wide pressure range. This is distinctly different from the usual elastic background behavior characterized by a smooth linear variation with pressure, in leading order. This approximate linear Δ*L*/*L* versus *P* behavior is revealed for the data set taken at 55 K (for the *b* axis, Fig. 2B), that is, way above *T*_{c}. As exemplarily indicated for the data at 43 K, which is about 6.5 K above *T*_{c}, strong variations in Δ*L*_{i}(*P*)/*L*_{i} with pressure, giving rise to two approximately linear regimes with distinctly different slopes, can be identified (broken lines in Fig. 2, A and B). The pressure-induced changes of the compressibilities revealed here are of the order of the compressibilities themselves, that is, *T*_{c}, *P*_{c}) to a critical elasticity as a result of the coupling of the critical electronic system to the lattice. The strongly nonlinear strain-stress relation observed here, representing the major result of our study, highlights a breakdown of Hooke’s law of elasticity over a considerably wide temperature-pressure range around the Mott transition. We note that indications for strong changes of the compressibility were also seen in previous ultrasonic investigations on the same compound (*28*).

### Phase diagram of κ-(BEDT-TTF)_{2}Cu[N(CN)_{2}]Cl

To classify the observed critical behavior and to compare it with theoretical predictions, precise knowledge about the system’s phase diagram (Fig. 3), including the location of the Mott transition line, the critical endpoint (*T*_{c}, *P*_{c}), and the various crossover lines is important. The inflection points of the Δ*L*_{i}(*P*)/*L*_{i} data in Fig. 2 (A and B) define a characteristic line, *P*_{c}(*T*) (Fig. 3), which is separated at the critical endpoint *P*_{c} = *P*_{c}(*T*_{c}) into a first-order transition line for *T* < *T*_{c} ≈ 36.5 K and the *Widom* line for *T* > *T*_{c}. We find a very good agreement of the Mott transition line obtained from our thermodynamic measurements with that of transport experiments (*23*) (see fig. S2). Our preliminary results on the effects of disorder, intentionally introduced into a crystal by x-ray irradiation (*29*), indicate a shift of the Mott transition line to lower pressures (see fig. S3), consistent with results from dynamic mean-field theory (*30*). We consider the good agreement of the phase transition lines, shown in fig. S2, as an indication that the crystals studied here are of similar high quality as those studied by Kagawa *et al.* (*23*). To determine *T*_{c} for the present crystal, we have taken two complementary approaches. First, we use the intersection point of the first-order transition line with a crossover line (Fig. 1A), which emanates from the critical point. As previously pointed out (*31*), the maximum response α_{max} in the coefficient of thermal expansion α(*T*) = *L*^{−1}*dL*/*dT* obtained from temperature-dependent measurements (see fig. S1) provides this crossover feature. Because of the finite slope of the first-order transition line at *T*_{c}, these measurements at varying *P* = const. also allow us to discriminate rounded crossover anomalies, expected for *P* > *P*_{c}, from jump-like features, characterizing the state at *P* < *P*_{c}. By including the positions of α_{max}(*T*, *P* = const.) into the phase diagram and applying the abovementioned criterion, we obtain a critical point at *T*_{c} = (36.5 ± 0.5) K and *P*_{c} = (23.4 ± 0.2) MPa.

A second approach for determining the critical point is provided by considering the width of the various features observed. This criterion also delivers insight into an additional scale that is introduced into the problem by the effects of disorder, relevant for any real material. Even far below the critical endpoint, we observe that the jump in Δ*L*_{i}/*L*_{i} is not infinitely sharp but broadened, giving rise to a finite width of the peak in κ_{i} at *P*_{c} that is practically temperature-independent. We assign this broadening to disorder, that is, spatial variations of internal stress induced by impurities and/or other crystal defects. However, upon increasing the temperature sufficiently above the endpoint, the width of the peak in κ_{i} becomes distinctly larger and strongly temperature-dependent and thus can be attributed to criticality. We estimate the width of the features in Δ*L*_{i}(*P*)/*L*_{i} by the full width at half maximum of the peak in the κ_{i} data (see fig. S4) and include this information [range delimited by red broken lines for *T* < *T*_{c} and blue broken lines for *T* > *T*_{c} in the phase diagram (Fig. 3)]. A linear extrapolation of these lines from high temperatures above 38 K, where the broadening is governed by criticality so that disorder effects are expected to play a minor role, to lower temperatures yields an intersection point of (36.7 ± 0.5) K, consistent with the abovementioned alternative estimate for *T*_{c}.

### Modeling of the relative length changes

After having pinpointed the critical endpoint (*T*_{c}, *P*_{c}) to a good degree of accuracy, we can proceed by further analyzing the critical behavior. Note that the disorder broadening of the critical signatures, implying an additional length scale, precludes a conventional one-parameter scaling collapse of the data (see section F, Supplementary Materials). Instead, we quantitatively compare our data to the theoretically expected behavior for a critical electronic system coupled to elastic degrees of freedom, yielding the universal properties of an isostructural solid-solid endpoint (*17*). For an orthorhombic crystal, the elastic constant tensor *C*_{ρν} in Voigt notation has six nondegenerate eigenvalues whose eigenvectors correspond to three shear and three compressive strain-singlet modes (*32*). One of the latter strain singlets, ε, can be identified with the order parameter of the solid-solid transition, which is measured by Δ*L*/*L*. This strain singlet exactly obeys a mean-field equation that, after appropriate rescaling, reads within Landau theory *r* ε + *u* ε^{3} = −σ (see section F and fig. S5, Supplementary Materials). Close to the endpoint, the dimensionless stress σ is given in terms of the hydrostatic pressure *P* by σ = −(*P* − *P*_{c}(*T*))/*P*_{c}(*T*_{c}) with the experimentally determined line *P*_{c}(*T*) and the rigidity *r* = (*T* − *T*_{c})/*T*_{c}, which is proportional to the associated eigenvalue of *C*_{ρν}. We neglect here a possible pressure dependence of the tuning parameter *r*, which would result in an additional but subleading nonlinearity in the strain-stress relation. Having experimentally determined the line *P*_{c}(*T*) as well as the critical temperature *T*_{c}, the mean-field equation is characterized by a single parameter, *u*, that quantifies the strength of the nonlinear cubic term ε^{3}. At *T*_{c}, the rigidity vanishes (*r* = 0), and the strain singlet responds in a nonlinear fashion to an applied stress, ε ~ (−σ)^{1/δ}, with the mean-field critical exponent δ = 3, signaling the breakdown of Hooke’s law.

To account for the disorder broadening in a phenomenological manner, we assume a multidomain state, caused by extrinsic defects of the crystal structure, where, for simplicity, each domain independently contributes to the mean strain. We average the mean-field solution ε(*r*, σ, *u*) with a Gaussian stress distribution *P*_{w}(*s*) with a constant variance, *w*, and zero mean, 〈ε(*r*, σ, *u*)〉_{w} = ∫*dsP*_{w}(*s*) ε(*r*, σ + *s*, *u*), that accounts for long-range stress fluctuations so that 〈ε(*r*, σ = 0, *u*)〉_{w} = 0. Close to the transition, slight variations of the stress can lead to the simultaneous compression and expansion of different parts of the crystal with the concomitant generation of dislocation defects and grain boundaries. This intrinsically generated disorder, in turn, may influence the stress distribution within the crystal (*33*, *34*), but a full self-consistent treatment of these effects is beyond the scope of the present paper. After normalizing the various data sets in Fig. 2 to their inflection point *P*_{c}(*T*) to eliminate offset contributions due to the temperature-dependent background expansion, the length changes are described by

Here, *A*_{i} is the proportionality constant between the singlet ε and the relative length change along the direction *i* = *a*, *b*. We determine the variance *w* from a Gaussian fit to the −*d(*Δ*L*_{i}/*L*_{i})/*dP* data at 30 K (see fig. S4 and section F, Supplementary Materials), that is, well below the critical point, which yields for the *a* axis *w*_{a} = 1.3 × 10^{−3} and the *b* axis *w*_{b} = 4.6 × 10^{−5}. Note that *w*^{1/2} measures the width of the broadening-affected pressure range. The fact that *a* axis than along the *b* axis. This implies strongly anisotropic domains, being elongated perpendicular to the planes, which can be rationalized by considering the spatial anisotropy inherent to the structure of κ-Cl (see fig. S6). As demonstrated in Fig. 4, fits based on Eq. 1 with only three parameters—*u*, *A*_{i} for *i* = *b*—provide an excellent description of the whole data set along the *b* axis covering a pressure range of 15 MPa ≲ *P* ≲ 30 MPa at 14 different temperatures of 30 K ≤ *T* ≤ 43 K (see fig. S7 for corresponding fits to the *a* axis data). For the dimensionless parameter *u* quantifying the nonlinearity of the Landau potential, we find *u* = 0.18. For the noncritical contribution, we obtain ^{−6} MPa^{−1}, which is of the same order as the background compressibility reported for a related material (*35*). The individual fits to the various *T* = const. data sets yield slightly different values for the proportionality constant *A*_{b}, which, however, follow to a good approximation a *T* linear variation *A*_{b} = (0.86 − 1.01 (*T* − *T*_{c})/*T*_{c}) × 10^{−4} (see fig. S8). We assign this effect to a smooth *T* dependence of the elastic tensor *C*_{ρν} affecting its eigenvalues and eigenvectors. We stress that a description of equally high quality is also obtained for the *a* axis data with identical *u* but different *A*_{a} *≠ A_{b}* (see fig. S7). Note that by changing the critical temperature within the indicated error margin, fits of equal quality are obtained with

*A*

_{i}values changing by 5% at maximum. Likewise, treating the critical temperature as an additional free parameter in the fitting procedure yields

*T*

_{c}= 36.2 K, consistent with the independently determined value of

*T*

_{c}= (36.5 ± 0.5) K. Taken together, the pronounced and strongly nonlinear lattice effects observed along the

*a*and

*b*axes of κ-Cl around its Mott transition can be well described on a quantitative level by mean-field criticality, as expected for a critical electronic system coupled to lattice degrees of freedom. We stress that mean-field criticality is also revealed in an independent approach based on the analysis of power-law dependencies in Δ

*L*

_{b}/

*L*

_{b}with respect to temperature and pressure upon approaching the critical endpoint (see section I, figs. S9 to S11, and table S1, Supplementary Materials).

## DISCUSSION

The raw data (Fig. 2) together with the fits based on Eq. 1 (Fig. 4) allow for a classification of the different scales involved. We find critically renormalized compressibilities κ_{i} around (*T*_{c}, *P*_{c}) = (36.5 K, 23.4 MPa) extending up to at least 43 K in temperature and up to about 25 MPa in pressure. This corresponds to a *T*-*P* range of critical elasticity of width Δ*T*/*T*_{c} ≈ 20% and Δ*P*/*P*_{c} ≈ 10% (yellow-shaded regime in Fig. 3). From this finding, we conclude that critical elasticity prevails in κ-Cl and that mean-field criticality governs its Mott endpoint. Note that critical elasticity can also be concluded (*36*) for the archetype Mott system V_{2}O_{3}, which is situated somewhat off the transition on the metallic side, where a substantial softening of one of the elastic constants has been observed (*37*).

At larger distances from the endpoint, where the lattice response becomes weak and acts as a small perturbation, a crossover to electronic criticality is expected with mean-field or nontrivial Ising exponents depending on the size of the underlying Ginzburg regime (blue-shaded regime in Fig. 1A). Whether this Ginzburg range exceeds the critical elastic regime (yellow-shaded regime in Fig. 1A) in the present system cannot be answered on the basis of the available data. However, the fact that thermal expansion measurements on a related system (*31*), the chemical pressure of which corresponds to (*P* − *P*_{c})/*P*_{c} ≈ 20%, were found to be consistent with 2D Ising criticality may indicate that the Ginzburg regime actually exceeds the mean-field regime for this material class. Moreover, we find that the width of the critical elasticity regime in pressure is a factor of 10 larger in comparison to the disorder-affected broadening regime, indicating that in this clean sample, disorder is only a minor effect and can be effectively treated within the framework of critical elasticity.

The critical elasticity revealed here for κ-Cl highlights the significance of strain measurements in probing the criticality of a finite-temperature Mott endpoint. At the same time, our result demonstrates that more theoretical work is needed to include the aspect of critical elasticity in the analysis of conductivity data (*23*), the relation of which to the order parameter is intricate (*38*) and whose classical critical signatures develop on a strongly varying background that has been assigned to quantum criticality (*39*, *40*).

Because of the long-range nature of the shear forces, mean-field criticality of an isostructural critical endpoint implies the absence of a diverging correlation length (*32*, *36*). As a consequence, one would expect that the critical slowing down observed by Hartmann *et al.* (*41*) saturates sufficiently close to the endpoint, provided that it is associated with a phonon mode. Alternatively, it could also arise from the slow dynamics associated with domain formation due to disorder. Further studies are required to identify the origin of the slowing down in κ-Cl close to its Mott endpoint.

In summary, our results presented for the organic conductor κ-(BEDT-TTF)_{2}Cu[N(CN)_{2}]Cl highlight the relevance of lattice degrees of freedom in the description of universal properties of the Mott transition in real materials, that is, in the presence of a compressible lattice. Whenever the Mott transition can be tuned by pressure, which is the case for many systems, such as various κ-(BEDT-TTF)_{2}*X* (*21*, *22*, *42*) or Et_{x}Me_{4–}* _{x}Z[Pd(dmit)_{2}]_{2} (43) organic charge-transfer salts, as well as transition metal oxides, such as (V_{1–x}M_{x}*)

_{2}O

_{3}(

*24*,

*44*–

*46*), NiO (

*47*), or PbCrO

_{3}(

*48*), the Mott endpoint will eventually be governed by critical elasticity. A different situation is encountered in VO

_{2}, where doping experiments suggest that the transition is insensitive to pressure so that the elastic coupling is expected to be weak. For this system, an extended regime was observed where the Mott criticality is governed by the random-field Ising model (

*49*). Generally, the steeper the first-order transition line is in the

*T*-

*P*phase diagram, the larger the size of the critical elastic regime should be. According to the Clausius-Clapeyron equation, the slope of this line

*d*/

*T*_{c}*dP*is given by Δ

*V*/Δ

*S*, that is, the ratio of the volume change to the change of entropy across the first-order transition. Because Δ

*S*has to vanish for

*T*

_{c}→ 0, implying an increase of

*d*/

*T*_{c}*dP*, the range of critical elasticity is expected to grow for materials with lower critical temperature. However, irrespective of the width of this regime, our results demonstrate that the Mott transition in a real material with a finite shear modulus sufficiently close to the endpoint shows the universal properties of an isostructural solid-solid endpoint with mean-field critical exponents rather than a liquid-gas endpoint.

## MATERIALS AND METHODS

Measurements of the relative length changes were performed by using a capacitive dilatometer (*50*, *51*), enabling length changes of Δ*L* ≥ 5 × 10^{−3} nm to be resolved. By combining the dilatometer with a continuously controllable helium-gas pressure system, in which the pressure cell is connected by a capillary to a room-temperature pressure reservoir, we could perform measurements both as a function of temperature at constant pressure and as a function of pressure at constant temperature (*20*). In this capacitive technique, length changes were calculated from the measured changes of the capacitance of a plate capacitor; see the work by Pott and Schefzyk (*51*) for technical details. To account for effects that are related to temperature-dependent and pressure-dependent changes of the dielectric function ϵ_{r} of the pressure-transmitting medium helium, a quantity that is relevant for the data processing, ϵ_{r}(*T*, *P*) was determined from independent measurements on a standard material (aluminum) at the same temperatures where *T* = const. measurements on κ-Cl were performed. During a pressure sweep, the temperature is stabilized within ± 0.1 mK on short term (30 min) and |Δ*T*| < 15 mK on long term (24 hours). The pressure changes, with a typical rate of 1.5 MPa/hour, could be determined with a resolution of |Δ*P*| ≤ 0.03 MPa. For technical reasons, pressure sweeps could only be performed under decreasing pressure, precluding the measurements of hysteresis loops. During temperature sweeps, which had been performed upon warming at a slow rate of 1.5 K/hour, the pressure was kept constant within |Δ*P*| ≤ 0.05 MPa. To hold the crystal in the dilatometer cell, a small uniaxial force (typically 0.2 to 0.8 N) has to be applied to the crystal. The corresponding uniaxial pressure depends on the cross-sectional area of the crystal studied. For the measurements along the *a* and *b* axes shown in Fig. 2 (A and B), the uniaxial pressure components were estimated to 0.4 and 0.8 MPa, respectively (see fig. S12 and section J, Supplementary Materials).

In total, three different crystals from different sources were studied and found to show similar behavior. For crystal #063-1, as discussed here in the main text, data were taken along all three axes. The results were cross-checked by measurements along the *a* axis for crystals #063-2 and #5-1. These crystals, with typical dimensions *a* × *b* × *c* ≈ 0.7 × 1 × 0.5 mm^{3} were grown following the standard electrochemical procedure (*52*). Before the measurements, the crystals were cooled through the glass transition *T*_{g} ≈ 67 K by using a slow cooling rate of −3 K/hour to reduce disorder in the ethylene end groups of the BEDT-TTF molecules (*53*).

## SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/2/12/e1601646/DC1

section A. Structure information for κ-(BEDT-TTF)_{2}Cu[N(CN)_{2}]Cl

section B. Measurements of relative length changes performed as a function of temperature at constant pressure

section C. Phase transition and *Widom* line from thermodynamic versus conductivity measurements

section D. Effects of x-ray irradiation on the first-order Mott transition

section E. Measurements of relative length changes performed as a function of pressure at constant temperature

section F. Theory of critical elasticity for the solid-solid endpoint in κ-Cl

section G. Disorder affected broadening of the signatures below *T*_{c}

section H. Quantitative modeling of relative length changes along the *a* axis as a function of pressure around the critical endpoint

section I. Independent estimate of critical exponents

section J. Effect of a small uniaxial pressure exerted by the dilatometer cell

fig. S1. Uniaxial coefficients of thermal expansion.

fig. S2. Phase transition and *Widom* line from thermodynamic versus conductivity measurements.

fig. S3. Effects of x-ray irradiation on the first-order Mott transition.

fig. S4. Uniaxial compressibilities.

fig. S5. Effect of higher-order couplings.

fig. S6. Schematic representation of the domain formation in κ-Cl.

fig. S7. Modeling of the relative length changes along the *a* axis.

fig. S8. Variation of the fit parameter *A*_{i} as a function of temperature.

fig. S9. Scheme of the scans in the phase diagram unraveling the different critical exponents.

fig. S10. Estimate of critical exponent δ.

fig. S11. Estimate of critical exponents β and γ.

fig. S12. Effect of a small uniaxial pressure exerted by the dilatometer cell.

table S1. Set of critical exponents β, γ, and δ of various universality classes.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is **not** for commercial advantage and provided the original work is properly cited.

## REFERENCES AND NOTES

**Acknowledgments:**We acknowledge useful discussions with J. Schmalian, P. Kopietz, I. Paul, M. Zacharias, and S. Winter.

**Funding:**The research was supported by the German Science Foundation via the Transregional Collaborative Research Center SFB/TR49 “Condensed Matter Systems with Variable Many-Body Interactions” and the Collaborative Research Center SFB 1143 “Correlated Magnetism: From Frustration to Topology.” Work in Sendai was partly supported by the Japan Society for the Promotion of Science KAKENHI grant #JP25287080. Work at Argonne, in a U.S. Department of Energy Office of Science laboratory, is operated under contract #DE-AC02-06CH11357. J.A.S. acknowledges support from the Independent Research and Development program while serving at the NSF.

**Author contributions:**M.L. conceived and supervised the project. Measurements of the relative length changes were performed by E.G. with contributions from R.S.M., U.T., B.W., and M.L. Theoretical description of the data was developed by M.G. with contributions from L.B. M.G. and E.G. analyzed the data. H.S., T.S., and J.A.S. grew the high-quality single crystals. M.L. wrote the paper with contributions from E.G. and M.G. All authors discussed the results and commented on the manuscript.

**Competing interests:**The authors declare that they have no competing interests.

**Data and materials availability:**All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

- Copyright © 2016, The Authors