Proposal to Observe Transverse Sound in Normal Liquid $^3$He in Aerogel

In the Fermi liquid metallic state, a static local magnetic moment is induced on the application of a circularly polarized electromagnetic wave, via the inverse Fara-day effect (IFE). The direction of this moment is along the direction of propagation of light, and the magnitude of the moment depends on the frequency of light, the temperature and various material parameters characteristic of the metal. I propose an analogous effect in the Fermi liquid state of $^3$H. A static circulating current is induced when liquid $^3$H is driven by a circularly polarized transverse acoustic wave. For liquid $^3$H filled into aerogel, the coupled system supports a low-attenuation transverse sound mode. I estimate the magnitude of induced circulating currents for this system and find that these are within the range of experimental measurement in the low-attenuation regime. The axis of circulation is along the direction of propagation of the acoustic wave. I propose this analogue of the inverse Faraday effect as a scheme to experimentally demonstrate the propagation of transverse sound in $^3$H-aerogel.


Introduction
The propagation of transverse sound in pure liquid 3 He was first predicted by Landau.While the existence of shear waves demonstrates the Fermi liquid nature of liquid 3 He, shear waves have never been observed in the normal liquid state of 3 He, primarily due to the large attenuation of this sound mode [1].However, transverse sound and the acoustic Faraday effect have been observed in superfluid 3 He [2] confirming the existense of transverse sound, as predicted in this system. 3He in aerogel is a realisation of a disordered Fermi liquid state.This system has a rich superfluid phase diagram including new superfluid phases unseen in the bulk and has been studied extensively over the past two decades [3].The superfluid state realises a disordered p-wave order parameter and the effect of impurities, whose density and isotropy (or lack of it) can be controlled, has been a subject of continued interest in the context of unconventional superconductors and novel phenomena.In this paper, I focus on the Fermi liquid state of 3 He in aerogel.
Acoustic studies have conventionally been a powerful probe of liquid 3 He physics.For the 3 He-aerogel system, the first to zero sound transition, characteristic of the Fermi liquid state, is not observed [4,5].This is explained by the collision drag effect [6], viz., coupling of simultaneous oscillations of the liquid 3 He and the aerogel matrix.The effect of collision drag on transverse sound is discussed in [7] and the authors find that the attentuation in the 3 He-aerogel coupled system exhibits a different behaviour as a function of temperature, as opposed to the pure 3 He case.At low temperatures, the attenuation remains finite due to impurity scattering( as opposed to attenuation going to zero in the pure liquid case).At higher temperatures, the attenuation crosses over to a regime in which hydrodynamic transverse sound can propagate over a long distance with the aid of the elasticity of the aerogel.However, inspite of the existence of this low-attenuation mode, transverse sound has not been experimentally observed in 3 Heaerogel.Although shear impedance measurements have been attempted in this system, no conclusive evidence of a shear mode could be made from the response of an acoustic transducer to pressure modulation from a transverse sound mode in an acoustic cavity [8].In this paper, I present a proposal to demonstrate the existence of propagating transverse sound in the 3 He-aerogel system, by measuring static circulating currents induced when the 3 He-aerogel system is subject to a circularly polarized shear wave.These circulating currents are induced ONLY for transverse acoustic perturbations driven by a circularly polarized drive, thus providing an unambiguous demonstration of the propagation of transverse sound in normal-state 3 He-aerogel.The motivation for this proposal is derived from an analogous effect in the Fermi liquid state of metals, the Inverse Faraday Effect(IFE).
I introduce the IFE in metals in Sec 2. The induced circulation of orbital currents of electrons (which can be expressed in terms of a local magnetic moment) can be calculated within the Eilenberger formulation of quasiclassical theory.In Section 3, I make the analogy to the liquid 3 He case and show that the estimate for the induced circulating current is large enough to detect in experiments.

Inverse Faraday Effect in Metals
The Inverse Faraday Effect(IFE) was first proposed by Pitaevskii [9] as an effect of a time-varying electric field on the stress tensor in a dispersive medium.[10] coined the name and first reported an observed induced magnetization in various liquids and glasses.The effect has since been observed in several magnetic materials [11] and more recently, in metals [12].Theoretical studies of the IFE have been reported using semiclassical [13], ab initio [14] and Ginzburg Landau [15,16] schemes in a variety of materials.A microscopic theory for the IFE using quasiclassical Green's functions has been recently developed and used to predict the induced magnetic moments in metals [17].The primary response of electrons in a metal exposed to circularly polarized electromagnetic radiation of frequency ω, is to the oscillating electric field given by E = E 0 e iK.R−iωt .Ignoring the band structure of electrons and for an isotropic Fermi surface, the corrections to the Green's functions for electrons responding to the oscillating electric field, give a static induced current density via the IFE.The rotational part of this current density is the curl of a magnetic moment given by [17], in the K → 0 limit.Here, Γ is the rate at which quasiparticles scatter from impurities/disorder, ω p is the plasma frequency in the metal and ϵ 0 and µ B are the free-space permittivity and the Bohr magneton, respectively.k B is the Boltzmann constant and T is the ambient temperature.The function ζ 2 (x) ≤ 1 over the temperature range and is discussed in [17].The DC magnetization is induced via a nonlinear effect upon application of an AC external field and can be regarded as a "rectification" effect.

Analogue in He
To make the analogy to the charge neutral Fermi liquid 3 He, I follow closely along the lines of the Keldysh formulation [18] of quasiclassical theory [19,20] employed in [17] to calculate the IFE in metals.I consider the response of 3 He quasiparticles to an oscillating shear field of an impinging transverse wave, U = U 0 e iq.R−iωt .This is analogous to the oscillating electric field of light in the IFE.I estimate the second-order current response analogous to the induced current density via the IFE.
The quasiclassical Green's function, ĝ is the propagator for quasiparticles with effective mass m ⋆ , energy ε and Fermi momentum p = m ⋆ v F p, given by solutions to the Eilenberger equation [19].Observables such as the quasiparticle density are calculated from the Keldysh components of the Green's function, ĝK , viz., ρ = dε 2π d 3 p (2πℏ) 3 T (ĝ K ); here, T refers to a trace over Nambu and spin indices.I expand ĝK in the external field U, ĝK = ĝK 0 + ĝK 1 + ĝK 2 , with ĝK i ∝ O(U i ) being the i-th order correction to ĝK and evaluate the current response.The details of the definitions are in the Appendix.I am interested in the second-order O(U 2 ) corrections to the current density, j 2 , averaged over the period of the oscillating field.It has been shown [17] that where ⟨...⟩ p,(2π/ω) ≡ ( ω 2π ) dt d 3 p (2π) 3 refers to an average over the Fermi surface and an average over a time-period of the field; ρ1 and v1 are the first-order corrections to the quasiparticle density and velocity, respectively viz., where τ are Pauli matrices in Nambu space, ρ 0 is the equilibrium density of quasiparticles, N f is the density of states at the Fermi level and T refers to the trace over the Nambu and spin degrees of freedom.The first-order corrections ĝK 1 can be obtained by solving the Eilenberger equations.The Keldysh functions ĝK may also be expressed in terms of quasiparticle distribution functions, n and the Retarded and Advanced functions, The distribution function then satisfies the Landau-Boltzmann transport equation with driving terms given by the collision integrals for quasiparticles scattering via various scattering channels.The distribution functions may also be expanded in the external field as the quasiclassical propagators.Then the first order corrections may be related, Solving the Eilenberger equation for ĝR,A gives ĝR,A 1 = 0, as shown in the Appendix.Then equation( 5) reduces to [7] have calculated the first-order corrections to the quasiparticle distribution function by solving the Boltzmann-Landau transport equation for 3 He-aerogel, including elastic scattering of quasiparticles from the aerogel matrix, the motion of the aerogel driven by the acoustic field, and mutual collisions between quasiparticles that provide an inelastic scattering channel.They find that the corrections to the distribution function n 1 (p, ε) = f ′ (ε 0 p )ν p e iq.R−iωt , where f ′ (ε 0 p ) is the equilibrium distribution function and ν p = Σ 2 l=1 ν l Y 1 l (θ, ϕ); θ and ϕ being the polar and azimuthal angles of the quasiparticle momentum p with respect to the direction of q respectively.The coefficients, ν l depend on the Fermi liquid parameters in 3 He, the characteristic amplitude, frequency and wavelength of the driving sound wave and contain information about scattering parameters, including the relaxation times for impurity and mutual quasiparticle scattering.They represent the so-called collision drag effect due to the aerogel on the distribution function of quasiparticles, in response to an impinging transverse sound wave.
Using the result for n 1 from [7] in equation( 6) and the notation ĝX ≡ g where β = (2k B T ) −1 and the equilibrium distribution function is given by the Fermi function.I can now use equation (7) in equation(2-3) to find j 2 , which represents the static current response to an applied oscillating field analogous to the "rectification" represented by the IFE.
The objective is to calculate the time-averaged corrections to the current density j that are of second order in the field, j 2 ∝ O(U 2 ), If I ignore the collision drag effect, ρ 2 (p) is given by equation( 3) using the second-order corrections ĝK 2 .ĝK 2 is obtained from equation(A4) with the driving term ∝ (v F •q) ĝK 1 .

As ĝK
1 ∝ (v F • q) (analogous to the IFE case, ignoring collision drag), this gives ⟨ρ 2 v 0 ⟩ p ∝ ⟨v F (v F • q)(v F • q)⟩ p = 0. Similarly, the second term in equation( 8) also gives a zero contribution.This leaves only the third term non-zero in j 2 above.ρ 1 and v 1 are given by ĝK 1 , as in equation (3).Ignoring the effects of collision drag beyond linear order O(U), j 2 is given by equation(2-3) with ĝK 1 given by equation (7)[21].Further, the dynamics of quasiparticles in 3 He-aerogel is in the hydrodynamic limit over a large frequency range [7].Assuming local equilibrium, I use the gradient representation derived in the Appendix and evaluate ρ1 and v1 , To estimate the magnitude of the induced current, I estimate the magnitude of Here U 0 is the amplitude of the shear wave which is circularly polarized, i.e., U 0 ̸ = U ⋆ 0 ; and ν 1a contains details about aerogel parameters including quasiparticle relaxation times and the geometric mean free path of the aerogel.j 2 has both rotational and irrotational terms.The irrotational part corresponds to currents in the plane of polarization of the transverse acoustic drive.Writing out the rotational part of the current, I estimate the time-averaged (j 2 ) rot thus, For a shear wave of frequency ω = 20M Hz and q ∼ 10 6 m −1 , used in [7], and for an aerogel sample of size ∼ 1cm on each side, I estimate the magnitude of this current density to be ∼ 10 −6 kg m −2 sec −1 .The curl in the expression above reveals the circulating nature of currents.Clearly, this vanishes for linearly polarized drive where U 0 = U ⋆ 0 .The currents circulate around an axis along the direction of U 0 × U ⋆ 0 .For a circularly polarized shear wave impinging along the axis of a cm-sized cylindrical sample of aerogel, the currents circulate about an axis that is parallel (or antiparallel) to the axis of the cylinder, the specific direction (parallel or antiparallel) being determined by the handedness of polarization of the shear wave.A schematic is shown in Fig. 1.I estimate the total current for a 1cm tall sample to be ∼ 10 −10 kg/sec, which is within the range of mass current measurement by experiment [22].U 0 LCP transverse acoustic wave j rot 3 He -Aerogel Fig. 1 An illustration circulating currents jrot ∝ U 2 0 induced by an impinging left circularly polarized (LCP) transverse acoustic wave of amplitude U 0 .The direction of the induced current will be reversed for right circularly polarized drive.
The collision drag effect couples simultaneous oscillations of liquid 3 He with that of the aerogel matrix leading to a decrease in the attenuation of transverse sound at higher temperatures [7].This decrease in attenuation depends on the velocity of transverse sound (not reported thus far) and can be up to two orders of magnitude for higher sound velocities, for temperatures ∼ 50 − 100mK [7].Thus, in this temperature range, the aerogel supports the propagation of transverse sound over longer distances enabling the detection of circulating currents produced as a result.While these currents would also exist in the bulk via the mechanisms described in this work, the attenuation of transverse sound in pure liquid 3 He is rather high (of the order of 500cm −1 ) rendering the induced currents weak and undetectable.
Linearly polarized transverse acoustic waves have been successfully realised experimentally and used in acoustic studies of liquid 3 He [1,4].Circularly polarized shear waves have, however, not been favourably generated experimentally for such studies yet.The generation of circularly polarized shear waves for applications ranging from solid-state to geophysical are based on the use of two transverse transducers that generate linearly polarized shear waves in orthogonal directions, locked in phase to generate a resultant circularly polarized wave.While a perfect circularly polarized shear wave might be difficult to generate due to intrinsic differences between transducers and (lack of) precision of phase locking [23], a "good" circularly polarized wave is not necessary to observe the induced current predicted in this work.An elliptically polarized shear wave or in particular, any shear drive, U with U × U ⋆ ̸ = 0 is sufficient to induce circulating currents described above (U × U ⋆ = 0 for linearly polarized shear drive).

Summary
I consider the effect in normal-state Fermi liquid 3 He analogous to the IFE in metals.The response of 3 He quasiparticles to an impinging circularly polarized shear wave is to transfer the angular momentum of the wave to quasiparticles resulting in an induced circulating current.As 3 He in aerogel is predicted to support a low-attenuation transverse sound mode, I estimate the magnitude of the induced circulating current for 3 He-aerogel.For high porosity aerogels, I estimate this current to be within the range of experimental measurement for sound waves in the low-attenuation M Hz frequency range for cm-sized samples.As this circulating current is induces only for shear drives, I propose the observation of this current as a demonstration of transverse sound propagation in liquid 3 He-aerogel.be calculated using equations(A4).Expand ĝ in the external field U, with ĝX i ∝ O(U i ) being the i-th order correction to ĝX and evaluate the current response.The current density operator at a space-time point (R, t) is given by where ρ(R, t) and v(R, t) are the local density and velocity of quasiparticles, respectively.(R and t are the centre of mass space and time coordinates).
where β ≡ (2k B T ) −1 and ε R,A = ε ± i0+.I use the notation for the Nambu elements, ĝX ≡ g X f X f X g X .g X and g X (and analogously f X and f X ) are related by symmetries, [ĝ R ] † = −τ 3 ĝA τ3 and [ĝ K ] † = τ3 ĝK τ3 .ĝX 1 are given by solutions to equations(A4).Using the notation in equation(A5) and the definition(A3), to O(U), ), using the notation defined above, with the symmetry g R (p, ε) = −g R (−p, −ε) ⋆ and g K (p, ε) = g K (−p, −ε) ⋆ .Using the ansatz, equation(A11) in equation(A10), I get equations for g R ± (g ± can be obtained using symmetries), Similarly, ĝA 1 = 0. Assuming local equilibrium, to lowest order in spatial gradients and to zeroth-order in applied fields, the Eilenberger equation (A4) for ĝK may be written as I use this gradient representation (to zeroth order in applied field) in equation( 3) to reveal the circulating nature of the induced current density.