# Universal Behaviour of the Superfluid Fraction and Tc of 3He in Aerogel

## Abstract

We have investigated the superfluid transition of <sup>3</sup>He in different samples of silica aerogel. By comparing new measurements on a 99.5% sample with previous observations on the behaviour of <sup>3</sup>He in 98% porous aerogel we have found evidence for universal behaviour of <sup>3</sup>He in aerogel. We relate both the transition temperature and superfluid density to the correlation length of the aerogel.

The properties of bulk <sup>3</sup>He are well understood. The extreme purity of <sup>3</sup>He at low temperatures makes it an ideal system to study the agreement between theoretical and experimental results on non-conventional Cooper pairing in the absence of disorder. Disorder plays a crucial role in suppressing the pairing interaction in high T<sub>c</sub> superconductors, the other well established non s-wave paired system. The superfluid transition of <sup>3</sup>He confined to a sample of very porous silica aerogel was first reported four years ago. The aerogel provides a structural disorder background to the liquid. <sup>3</sup>He is compressible, and the density can be continuously tuned by $``$30% while maintaining a fixed disorder. The <sup>3</sup>He zero temperature coherence length $`\xi _0`$, defined as $`\xi _0=\mathrm{}v_f/k_BT_c`$, varies from 180Å to over 700Å as a function of density. Because the Cooper pairs in <sup>3</sup>He form in a p-wave state, quasiparticle scattering from the aerogel strands is pairbreaking. Thus the <sup>3</sup>He-in-aerogel system is well suited to the exploration of the effect of impurity scattering and disorder on the superfluid transition and phase diagram.

The superfluidity of <sup>3</sup>He in silica aerogel has been studied using torsional oscillators, NMR and sound propagation techniques. These measurements show that both the superfluid transition temperature (T<sub>c</sub>) and superfluid density ($`\rho _s`$) of the <sup>3</sup>He are suppressed by the disorder, but that the transition remains sharp. This suppression is sensitive to both the density and the microstructure of the aerogel sample.

The simplest model for the effect of impurity scattering on the <sup>3</sup>He superfluid transition is the homogeneous scattering model (HSM) which is based on the Abrikosov-Gorkov model for a superconductor with magnetic impurities that induce pair-breaking via spinflip scattering. This mechanism is similar to that of diffuse scattering of Cooper paired <sup>3</sup>He from a surface, and is unable to explain the observed behaviour. Specifically, the observed suppression of the superfluid density is much greater than predicted by this model. More sophisticated models, such as the isotropic inhomogeneous scattering model (IISM) proposed by Thuneberg and co-workers are able to quantitatively predict the superfluid transition temperature of <sup>3</sup>He in aerogel (for small suppressions) and have had success at qualitatively explaining the observed superfluid densities.

In this Letter we present data from several different experiments on <sup>3</sup>He in aerogel, including new results on <sup>3</sup>He confined to a 99.5% porosity sample. This sample is a factor of four more dilute than any previously investigated and is crucial for understanding the evolution from bulk <sup>3</sup>He to a regime where impurity scattering dominates. In comparing these different samples we find evidence that the relation between superfluid density and the superfluid transition temperature of <sup>3</sup>He in aerogel follows a universal behaviour, independent of the aerogel sample. This is significant because both of these quantities are individually sensitive to the microstructure of the aerogel, and vary greatly from sample to sample. We also present evidence that the suppression of T<sub>c</sub> can be related to the correlation length ($`\xi _a`$) of the aerogel sample.

The aerogels used in the experiments discussed in this Letter were grown under basic conditions. Under these conditions gelation is the result of diffusion limited aggregation of small ($``$30Å diameter) primary silica particles. The aerogels are characterized by a fractal dimension (D<sub>f</sub>) related to the real space correlations, and a long length scale cutoff to these correlations ($`\xi _a`$) above which the sample appears homogeneous. The fractal exponent depends only on the gelation process, while the cutoff length is also dependent on the average density. Simulations based on the diffusion-limited cluster-cluster aggregation (DLCA) algorithm predict that the fractal exponent should lie between 1.7 and 1.9, which is in good agreement with small-angle X-ray scattering (SAXS) measurements (Table I). We note that $`\xi _a`$ in the most dilute sample, D, could not be inferred from the data as the SAXS did not extend to sufficiently small q.

Samples A and C have the same density, but were made under different conditions and have slightly different fractal dimension and significantly different $`\xi _a`$. Both samples A and C have been studied with SAXS. Sample B was made under conditions very similar to sample C. We do not have direct information on its microstructure from SAXS, but assume here that samples B and C are essentially identical. The correlation length for sample D is within the range obtained from simulations based on the DLCA algorithm and is consistent with the SAXS data of Mulders on a 99.6% sample ($`\xi _a`$2000Å). For a more extensive discussion the reader is referred to reference .

The 99.5% porosity aerogel used for our experiment was grown inside the (on average) 100 $`\mu `$m large pores of a coarse silver sinter. Previous torsional oscillator experiments have been affected by the presence of spurious resonances resulting from composite modes of <sup>3</sup>He and aerogel whose frequency crosses the resonant frequency of the cell. The strength of these resonances grow as the porosity of the aerogel sample increases, affecting the quality of data. In our cell the aerogel is clamped to the silver sinter, and the effect of spurious resonances is strongly reduced.

We operated our torsional oscillator in self-resonant mode near 483 Hz. The temperature in the cell was measured using a lanthanum diluted cerous magnesium nitrate ac susceptibility thermometer, thermally connected to the sample through a shared reservoir of <sup>3</sup>He. Data was collected while the temperature increased at a rate of 20 $`\mu `$K per hour. The period shift of the oscillator as the superfluid <sup>3</sup>He decoupled from the torsion head provided both the transition temperature of the <sup>3</sup>He in aerogel and the superfluid density.

Figure 1 shows the superfluid transition temperature for several different aerogel samples. Three of these samples (A, B and C) have a nominal porosity of 98%, the fourth one (D) is our 99.5% aerogel. All of these measurements were done by monitoring the period shift in torsional oscillators filled with the aerogel and <sup>3</sup>He. The difference in transition temperature between samples A, B and C arises from differences in the microstructure of the aerogels (see table I and reference ). The relative suppression of the transition temperature (1-T<sub>c</sub>/T<sub>c0</sub>) (with T<sub>c0</sub> the transition temperature in bulk <sup>3</sup>He) is larger at lower pressures (larger $`\xi _0`$) than at higer pressures (smaller $`\xi _0`$). In view of the fact that all the aerogels used in these experiments have a similar fractal dimension and primary particle size, one would expect them to be mutually self similar on length scales shorter than the fractal cutoff length. Thus for temperature dependent <sup>3</sup>He coherence lengths $`\xi (T)`$ shorter than $`\xi _a`$, the ratio of density of silica sampled at two different $`\xi (T)`$ is independent of aerogel density. Experimentally one observes a strong dependence of the suppression of T<sub>c</sub> on the aerogel density.

In figure 2 we show the dependence of the relative suppression of T<sub>c</sub> as a function of $`\xi _0`$. At short coherence lengths, the relative suppression is small for all samples. However, when the pressure is reduced, increasing $`\xi _0`$, this suppression shows a marked dependence on the microstructure of the aerogel. The transition temperature shows an evolution from a strong impurity scattering regime when the <sup>3</sup>He is confined to a 98% porosity aerogel sample towards the behaviour of bulk <sup>3</sup>He in the 99.5% sample.

Motivated by figure 2 we plot the relative suppression of T<sub>c</sub> against $`\xi _0`$ scaled by the aerogel correlation length $`\xi _a`$. The aerogel-limited mean free path ($`l_g`$) would be another natural choice to compare against $`\xi _0`$. However, this length scale has not been independently measured. If we use the HSM to determine $`l_g`$ from T<sub>c</sub>/T<sub>c0</sub>, the values of $`l_g`$ show a strong pressure dependence. Figure 3 shows that the scaled transition temperature depends solely on the ratio of $`\xi _0`$ to $`\xi _a`$ *independent* of the aerogel density. The error bars for sample D result from the high and low estimates (2000Å and 3000Å respectively) for the correlation length for a 99.5% porosity sample from DLCA simulations.

We determined the temperature-dependent superfluid density of <sup>3</sup>He in 99.5% aerogel using the shift in resonant frequency of our torsional oscillator upon warming. A small, temperature dependent contribution due to bulk <sup>3</sup>He in the cell was subtracted. The remaining shift was scaled by the period change due to filling the cell with <sup>3</sup>He at 50 mK and by the tortuosity (measured with <sup>4</sup>He) to obtain the superfluid density. Figure 4 shows the bare superfluid density ($`\rho _s^b`$/$`\rho `$) plotted versus the reduced temperature T/T<sub>c</sub> for different aerogel samples and pressures. The bare superfluid density is obtained from $`\rho _s`$ by stripping away the Fermi liquid factor according to:

$$\frac{\rho _s^b}{\rho }=\frac{(1+\frac{1}{3}F_1)\frac{\rho _s}{\rho }}{1+\frac{1}{3}F_1\frac{\rho _s}{\rho }}$$
(1)   
and is equivalent to 1-Y(T) where Y(T) is the temperature dependent Yosida function for bulk <sup>3</sup>He.

As with the transition temperature T<sub>c</sub>, the $`\rho _s^b`$ for these two different samples are similar at high pressures and are both close to the bulk value. At the lowest pressure, there is a factor of five difference in $`\rho _s^b`$/$`\rho `$ at the same reduced temperature between the 99.5% sample and the 98% sample; $`\rho _s^b`$/$`\rho `$ is more strongly suppressed by the aerogel than T<sub>c</sub>. There is also a large difference between the suppression factors of T<sub>c</sub> and $`\rho _s^b`$ in the same aerogel sample. This large suppression in $`\rho _s`$ with respect to T<sub>c</sub> is consistent with measurements made on other samples, and cannot be explained with homogenous scattering models for <sup>3</sup>He in aerogel.

In order to better understand the superfluid transition of <sup>3</sup>He in aerogel, in figure 5 we plot $`\rho _s^b`$/$`\rho _s`$ at 0.8T<sub>c</sub> against (T<sub>c</sub>/T<sub>c0</sub>)<sup>2</sup> for several aerogel samples. The error bars shown for cell B arise from a spurious resonance that introduces uncertainty in the determination of $`\rho _s`$. As in figure 3, the data collapses on a universal curve. The dashed line is the prediction for a homogeneously scattering model based on the Abrikosov-Gorkov equation. This plot compares aerogels with *different* densities—there is a factor of four difference in the average impurity density between the 98% samples and the 99.5% aerogel sample. Furthermore, the coherence length of the Cooper pairs varies from 180 Å to 600 Å over this data, yet the very strong pressure dependence shown in figures 2 and 4 has been factored out in this plot.

The correlations in the aerogel will affect the suppression of T<sub>c</sub> and the evolution of $`\rho _s^b`$ relative to a homogeous disorder, but this suppression apparently depends only on the correlation length of the sample and possibly the fractal exponent. The steep slope of the data in figure 5 is evidence that the fractal nature of the aerogel plays an important role in the development of $`\rho _s`$. Since all of these samples were base-catalyzed, the behaviour of <sup>3</sup>He in each aerogel is determined mainly by $`\xi _a`$. As long as $`\xi _0`$ is much less than $`\xi _a`$, the disorder sampled by the ensemble of Cooper pairs should be insensitive to changes in the temperature-dependent coherence length $`\xi `$(T), until $`\xi `$(T)$`\xi _a`$. That is, the system has a conformal symmetry normally absent in disordered systems. The evidence of this one-parameter scaling is displayed in figure 5; $`\rho _s^b`$ is a function of T<sub>c</sub>/T<sub>c0</sub> (or equivalently $`\xi _0`$/$`\xi _a`$) only. This behaviour is reminiscent of the compilation of data from disordered high T<sub>c</sub> materials by Franz *et al* with the exception that the relative suppression of T<sub>c</sub>/T<sub>c0</sub> and $`\rho _s/\rho _{s0}`$ in the high T<sub>c</sub> materials do not follow a universal behaviour, presumably because the impurities are not fractally correlated.

In order to compare the behaviour of T<sub>c</sub> and $`\rho _s^b`$ of <sup>3</sup>He in aerogel with <sup>3</sup>He in bulk it will be necessary to understand precisely *how* the fractal disorder affects the superfluid pairing mechanism. One test for models of non-conventional Cooper pairing in the presence of disorder would be to predict the functional form of the universal curve for $`\rho _s^b`$ versus (T<sub>c</sub>)<sup>2</sup> for <sup>3</sup>He in (base-catalyzed) aerogel. The IISM of Thuneberg *et al.* predicts a relation between $`\rho _s^b/\rho _{s0}^b`$ and T<sub>c</sub>/T<sub>c0</sub> similar to the trend illustrated in figure 5, showing behaviour very different than Abrikosov-Gorkov model. This model does not explicitly consider the fractal nature of the aerogel, but shows how inhomogeneities in the disorder can lead to a large suppression of $`\rho _s^b`$ relative to T<sub>c</sub>.

In this letter we have presented data from our measurements on <sup>3</sup>He in a very dilute 99.5% porous silica aerogel. The values of T<sub>c</sub> and $`\rho _s^b`$ fall between those of <sup>3</sup>He in bulk and <sup>3</sup>He in denser aerogel samples. We also present T<sub>c</sub> and $`\rho _s^b`$ data for <sup>3</sup>He in aerogel experiments performed at Cornell that show universal behaviour that can be traced to the fractal structure of the aerogel. In order to more fully understand this exciting physical system, more attention must be devoted to understanding the microstructure of the aerogel. Specifically, the universal scaling discussed above depends strongly on the fact that the fractal exponent real-space correlations is similar for all the aerogel samples. It would be interesting to study neutrally catalyzed silica aerogels, which have a different fractal exponent than the base-catalyzed samples. <sup>3</sup>He would be expected to follow *different* universal behaviour when confined to base-catalyzed and neutrally catalyzed aerogels. As yet, there is no clear theoretical picture for why the transition temperature T<sub>c</sub>should depend on the ratio $`\xi _0/\xi _a`$. The explanation for this behaviour will provide insight into how correlated disorder affects non-conventional Cooper pairing.

We would like to acknowledge helpful conversations with T.L. Ho, S. Yip, E. Thuneberg, J. Beamish, A. Golov and M.H.W. Chan. This research was supported by the NSF under DMR-9705295.