Home | About | Journals | Submit | Contact Us | Français |

**|**Scientific Reports**|**PMC5556081

Formats

Article sections

Authors

Related links

Sci Rep. 2017; 7: 8023.

Published online 2017 August 14. doi: 10.1038/s41598-017-08605-7

PMCID: PMC5556081

Igor I. Smolyaninov, Email: ude.dmu@yloms.

Received 2017 January 23; Accepted 2017 July 12.

Copyright © The Author(s) 2017

We consider the recently suggested model of a multiverse based on a ferrofluid. When the ferrofluid is subjected to a modest external magnetic field, the nanoparticles inside the ferrofluid form small hyperbolic metamaterial domains, which from the electromagnetic standpoint behave as individual “Minkowski universes” exhibiting different “laws of physics”, such as different strength of effective gravity, different versions of modified Newtonian dynamics (MOND) and different radiation lifetimes. When the ferrofluid “multiverse” is populated with atomic or molecular species, and these species are excited using an external laser source, the radiation lifetimes of atoms and molecules in these “universes” depend strongly on the individual physical properties of each “universe” via the Purcell effect. Some “universes” are better fine-tuned than others to sustain the excited states of these species. Thus, the ferrofluid-based metamaterial “multiverse” may be used to study models of MOND and to illustrate the fine-tuning mechanism in cosmology.

Many physical properties of our universe, such as the relative strength of the fundamental interactions, the value of the cosmological constant, etc., appear to be fine-tuned for existence of human life. One possible explanation of this fine tuning assumes existence of a multiverse, which consists of a very large number of individual universes having different physical properties. Intelligent observers populate only a small subset of these universes, which are fine-tuned for life.

While this point of view may not be falsifiable based on astrophysical observations, one possible way to ascertain its viability may rely on condensed matter physics. In particular, it was suggested very recently^{1, 2} that magnetic nanoparticle-based ferrofluids share some common features with the class of cosmological models based on loop quantum gravity. This analogy relies on the fact that a modest external magnetic field aligns most of the individual magnetic nanoparticles in the ferrofluid into long parallel chains (see Fig. 1), so that the ferrofluid becomes a self-assembled hyperbolic metamaterial^{3} (an extremely anisotropic uniaxial material, which behaves like a metal in one direction and like a dielectric in the orthogonal direction). Photons in a hyperbolic metamaterial exhibit gravity-like nonlinear optical interactions^{4}, which may be traced to perturbations of an emergent effective Minkowski spacetime^{5} that describes light propagation in such metamaterials. It appears that both loop quantum gravity models and hyperbolic metamaterials may exhibit metric signature phase transitions^{6}, during which the spacetime metric used to describe the system changes its signature. Moreover, the metric signature transition in a ferrofluid leads to separation of the effective spacetime into a multitude of intermingled Minkowski and Euclidean domains, giving rise to a picture of “metamaterial multiverse”^{7}. Inflation-like behaviour appears to be generic within individual Minkowski domains^{8}. Thus, ferrofluid-based self-assembled metamaterial geometry captures many features of several cosmological models of the multiverse, such as metric signature transition scenario in loop quantum cosmology^{1, 2}, natural emergence of a large number of Minkowski universes, and inflation. Moreover, in this paper we will demonstrate that individual Minkowski domains in the ferrofluid exhibit different “laws of physics”, such as different versions of modified Newtonian dynamics (MOND), different strength of effective gravity, and different radiation lifetimes due to variations in the local Purcell factor, so that ferrofluids may be used to illustrate the fine-tuning mechanism in cosmology. On the other hand, it is also interesting to study the Minkowski domains which show the same (or very similar) effective laws of physics and properties. Such a study may illustrate potential existence of “parallel universes” and shed some light on the “measure problem” in a multiverse, which has to do with making probabilistic predictions of some particular measurement outcomes in a multiverse setting. All these effects may be studied in the described ferrofluid model via direct microscopic observations.

The metric signature transition in ferrofluids (shown in Fig. 1) may be understood based on the theory of superparamagnets^{9} and the Maxwell-Garnett approximation, which describes electromagnetic properties of wire array metamaterials^{10}. Our experiments use the cobalt magnetic fluid 27-0001 from Strem Chemicals, which is composed of 10nm cobalt nanoparticles in kerosene coated with sodium dioctylsulfosuccinate and a monolayer of LP4 fatty acid condensation polymer. The volume fraction of cobalt nanoparticles in this ferrofluid is *α*=8.2%. Magnetic interaction of these nanoparticles is rather weak, so in the absence of external magnetic field the nanoparticles are randomly distributed within the fluid, as illustrated in Fig. 1(a). On the other hand, application of a modest external magnetic field (of the order of 300 Gauss) leads to formation of nanocolumns aligned along the external field^{11}, as illustrated in Fig. 1(b). In addition, depending on the magnetic field magnitude, solvent used, and the nanoparticle concentration, the ferrofluid undergoes phase separation into nanoparticle rich and nanoparticle poor phases^{3}, which is easy to observe in microscopic images of the ferrofluid, as illustrated in Fig. 2(b): the typical spatial scale of phase separation is 0.1–1μm.

(**a**) Wavelength dependencies of the real parts of ε_{z}=ε_{2} and ε_{x,y}=ε_{1} for a cobalt nanoparticle-based ferrofluid at α=8.2% volume fraction of nanoparticles. While **...**

Ferrofluids subjected to external magnetic field are known to exhibit classical superparamagnetic behaviour, which is a form of magnetism that is exhibited by magnetic materials consisting of small ferromagnetic or ferrimagnetic nanoparticles. When an external magnetic field *H* is applied to an assembly of magnetic nanoparticles, their magnetic moments tend to align along the applied field, leading to a net magnetization^{9}:

$$M(H,T)=\alpha \mu \phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right),$$

1

where *α* is the nanoparticle concentration, *μ* is their magnetic moment, and *T* is the temperature. Ferrofluids are known from the experiment to obey the same relationship. However, as illustrated in Fig. 1(b) at small magnetic fields (100 Gauss^{3}) not all the nanoparticles in the ferrofluid are aligned into chains, and these free floating nanoparticles are much less aligned by the weak external field due to absence of the neighbors. Therefore, it is reasonable to assume that the volume fraction of the nanoparticles aligned into chains *α(Η, Τ)* is

$$\alpha (H,T)={\alpha}_{\infty}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)$$

2

where the limiting volume fraction *α*
_{∞}=0.082 is reached at a very large magnetic field when all the available nanoparticles in the ferrofluid participate in the chain formation.

Let us now consider the electromagnetic properties of the ferrofluid within the scope of this model. They may be understood based on the Maxwell-Garnett approximation via the dielectric permittivities *ε*
_{m} and *ε*
_{d} of cobalt and kerosene, respectively. At small magnetic fields the difference *α*
_{∞}−*α*(*H, T*) describes cobalt nanoparticles, which are not aligned and distributed homogeneously inside the ferrofluid. Dielectric polarizability of these nanoparticles may be included into *ε*
_{d}, leading to slight increase in its value. When the nanoparticles are aligned into chains, the diagonal components of the ferrofluid permittivity may be calculated using the Maxwell-Garnett approximation as derived in ref. ^{10}, where the volume fraction of the metal phase is assumed to be equal to *α(Η, Τ)*:

3

$${\epsilon}_{x}={\epsilon}_{y}={\epsilon}_{1}=\frac{2\alpha (H,T){\epsilon}_{m}{\epsilon}_{d}+(1-\alpha (H,T)){\epsilon}_{d}({\epsilon}_{d}+{\epsilon}_{m})}{(1-\alpha (H,T))({\epsilon}_{d}+{\epsilon}_{m})+2\alpha (H,T){\epsilon}_{d}}$$

4

where z direction is aligned with the direction of external magnetic field. As far as the magnetic permeability is concerned, at the visible and infrared frequencies the ferrofluid may be considered as a non-magnetic (μ=1) medium. Calculated wavelength dependencies of *ε*
_{2} and *ε*
_{1} at *α(H, T)*=*α*
_{∞}=8.2% are plotted in Fig. 2(a). These calculations are based on the optical properties of cobalt in the infrared range^{12}. While *ε*
_{1} stays positive and almost constant, *ε*
_{2} changes sign to negative around *λ*=1μm. So in a very strong magnetic field the ferrofluid becomes a hyperbolic metamaterial at *λ*>1μm.

On the other hand, Eq. (^{3}) demonstrates that if the volume fraction of cobalt nanoparticles varies depending on either magnetic field or temperature, the change of sign of *ε*
_{2} occurs at some critical value α_{H}:

$$\alpha (H,T)>{\alpha}_{H}=\frac{{\epsilon}_{d}}{{\epsilon}_{d}-{\epsilon}_{m}},$$

5

so that the ferrofluid becomes a hyperbolic metamaterial^{13, 14} above the critical volume fraction. Alignment of cobalt nanoparticle chains along the direction of external magnetic field is indeed clearly revealed by microscopic images of the ferrofluid, as shown in Fig. 2(b,c).

Propagation of monochromatic extraordinary photons inside the ferrofluid is described by the wave equation, which is formally equivalent to a 3D Klein-Gordon equation for a massive field
_{ω}=*E*
_{z} in a 3D Minkowski spacetime:

$$-\frac{{\partial}^{2}{\varphi}_{\omega}}{{\epsilon}_{1}\partial {z}^{2}}+\frac{1}{(\phantom{\rule{-.25em}{0ex}}-\phantom{\rule{-.25em}{0ex}}{\epsilon}_{2})}\left(\frac{{\partial}^{2}{\varphi}_{\omega}}{\partial {x}^{2}}+\frac{{\partial}^{2}{\varphi}_{\omega}}{\partial {y}^{2}}\right)=\frac{{{\omega}_{0}}^{2}}{{c}^{2}}{\varphi}_{\omega}$$

6

in which the spatial coordinate *z* behaves as a timelike variable, and *E*
_{z} is the z-component of electric field in the wave^{5}. Eq. (^{6}) exhibits effective Lorentz invariance under the coordinate transformation

$$z\text{'}=\frac{1}{\sqrt{1-\frac{{\epsilon}_{1}}{(-{\epsilon}_{2})}\beta}}(z-\beta x)$$

7

$$x\text{'}=\frac{1}{\sqrt{1-\frac{{\epsilon}_{1}}{(-{\epsilon}_{2})}\beta}}\left(x-\beta \frac{{\epsilon}_{1}}{(\phantom{\rule{-.25em}{0ex}}-\phantom{\rule{-.25em}{0ex}}{\epsilon}_{2})}z\right),$$

where *β* is the effective boost. Similar to our own Minkowski spacetime, the effective Lorentz transformations in the *xz* and *yz* planes form the Poincare group together with translations along *x, y*, and *z* axis, and rotations in the *xy* plane. Thus, the wave Eq. (^{6}) describes world lines of massive particles which propagate in an effective 2+1 dimensional Minkowski spacetime. Components of the metamaterial dielectric tensor define the effective metric coefficients *g*
_{ik} of this spacetime: *g*
_{00}=*−ε*
_{1} and *g*
_{11}=*g*
_{22}=*−ε*
_{2}. Nonlinear optical Kerr effect bends this spacetime resulting in the effective gravitational interaction between extraordinary photons^{4}. Thus, transition to a hyperbolic state in the ferrofluid, which is described by Eq. (^{5}), may be understood as transition of the effective “electromagnetic spacetime” from Euclidean to Minkowski state. This transition may be driven by either temperature or external magnetic field. Since both parameters may experience variations across the volume of the ferrofluid, it is natural to expect that near the transition the ferrofluid separates into large number of individual Minkowski and Euclidean domains. Indeed, such domains may be clearly revealed in Fig. 2(c), which is obtained from the experimental image of light intensity $I(\overrightarrow{r})$ in Fig. 2(b) by numeric differentiation $\left|I(\overrightarrow{r})-I(\overrightarrow{r}+\overrightarrow{d}/2)\right|$, where $\overrightarrow{d}$ is the periodicity of stripes visible in Fig. 2(b). This numeric procedure reveals the stripe contrast, and emphasizes ferrofluid regions having larger volume fractions of nanoparticles *α(H,T)* aligned into the nanochains. According to Eq. (^{5}), these regions behave as “Minkowski domains” separated by regions of Euclidean spacetime, leading to a picture of “metamaterial multiverse”.

Here it is interesting to note that the effective metric near the Minkowski domain walls emulates cosmological inflation. Indeed, based on Eq. (^{6}) it is clear that the factor (−*ε*
_{2}) plays the role of a scale factor of the effective Minkowski spacetime

8

(since according to Eq. (^{4}) *ε*
_{1} is positive and almost constant, as plotted in Fig. 2(a)). In the case of *z*-dependent *ε*
_{1}=*ε*
_{x}=*ε*
_{y} and *ε*
_{2}=*ε*
_{z} the electromagnetic field separation into the ordinary and the extraordinary components remains well defined^{8, 15}. Taking into account *z* derivatives of *ε*
_{1} and *ε*
_{2}, Eq. (^{6}) becomes

$$\begin{array}{l}-\phantom{\rule{.25em}{0ex}}\frac{{\partial}^{2}{\varphi}_{\omega}}{{\epsilon}_{1}\partial {z}^{2}}+\frac{1}{(\phantom{\rule{-.25em}{0ex}}-\phantom{\rule{-.25em}{0ex}}{\epsilon}_{2})}\left(\frac{{\partial}^{2}{\varphi}_{\omega}}{\partial {x}^{2}}+\frac{{\partial}^{2}{\varphi}_{\omega}}{\partial {y}^{2}}\right)+\left(\frac{1}{{\epsilon}_{1}^{2}}\left(\frac{\partial {\epsilon}_{1}}{\partial z}\right)-\frac{2}{{\epsilon}_{1}{\epsilon}_{2}}\left(\frac{\partial {\epsilon}_{2}}{\partial z}\right)\right)\left(\frac{\partial {\varphi}_{\omega}}{\partial z}\right)\\ +\phantom{\rule{.25em}{0ex}}\frac{{\varphi}_{\omega}}{{\epsilon}_{1}{\epsilon}_{2}}\left(\frac{1}{{\epsilon}_{1}}\left(\frac{\partial {\epsilon}_{1}}{\partial z}\right)\left(\frac{\partial {\epsilon}_{2}}{\partial z}\right)-\left(\frac{{\partial}^{2}{\epsilon}_{2}}{\partial {z}^{2}}\right)\right)=\frac{{{\omega}_{0}}^{2}}{{c}^{2}}{\varphi}_{\omega}\end{array}$$

9

Since *ε*
_{1} is approximately equal to *ε*
_{d} (as evident from Fig. 2(a)) its derivatives may be neglected. If we also neglect the second derivative of *ε*
_{2} (thus assuming slow “adiabatic” changes of *ε*
_{2} as a function of spatial variables, which may be justified by small variations of temperature on microscopic scales) the wave equation for the extraordinary field
_{ω}=*E*
_{z} may be re-written as

$$-\frac{{\partial}^{2}{\varphi}_{\omega}}{{\epsilon}_{1}\partial {z}^{2}}+\frac{1}{(-\phantom{\rule{-.25em}{0ex}}{\epsilon}_{2})}\left(\frac{{\partial}^{2}{\varphi}_{\omega}}{\partial {x}^{2}}+\frac{{\partial}^{2}{\varphi}_{\omega}}{\partial {y}^{2}}\right)-\frac{2}{{\epsilon}_{1}{\epsilon}_{2}}\left(\frac{\partial {\epsilon}_{2}}{\partial z}\right)\left(\frac{\partial {\varphi}_{\omega}}{\partial z}\right)=\frac{{{\omega}_{0}}^{2}}{{c}^{2}}{\varphi}_{\omega}$$

10

It is easy to verify that the latter equation coincides with the Klein-Gordon equation

$$\frac{1}{\sqrt{-g}}\frac{\partial}{\partial {x}^{i}}\left({g}^{ik}\sqrt{-g}\frac{\partial \psi}{\partial {x}^{k}}\right)=\frac{{m}^{2}{c}^{2}}{{\hslash}^{2}}\psi $$

11

for a massive particle in a gravitational field^{16} described by the metric coefficients

*g*
_{00}=*−ε*
_{1} and *g*
_{11}=*g*
_{22}=*−ε*
_{2}. Indeed, for a wave field *ψ* = (−*ε*_{2})^{1/2}*φ* Eq. (^{10}) may be derived from Eq. (^{11}) if we also neglect the second derivatives of *ε*
_{2} in Eq. (^{11}).

Let us assume that the volume of ferrofluid is very large, and it has not reached thermodynamic equilibrium yet, while the external magnetic field was being changed during the experiment (this is possible, for example, in the presence of convective currents inside the ferrofluid, which are caused by changing external magnetic field). Under these conditions the local temperature distribution inside the ferrofluid may be described as

12

where *T*
_{c} is the temperature of the metric signature transition (the point where *ε*
_{2}=0). Since z coordinate plays the role of time in the effective spacetime described by Eq. (^{8}), the temperature gradient *η* will result in formation of a Minkowski domain, such as the one shown in the experimental image in Fig. 3(a). Based on Eq. (^{3}), the experimentally measured image contrast $|I(\overrightarrow{r})-I(\overrightarrow{r}+\overrightarrow{d}/2)|~\alpha ({\rm H},{\rm T})$ in Fig. 2(c) is proportional to the scale factor (−*ε*
_{2}) of the effective spacetime within the domain (see Eq. (^{8})). The cross section of Fig. 3(a) along the arrow indicating the effective time direction indeed shows very fast electromagnetic spacetime expansion near the Minkowski domain wall (we should also point out here that local magnetic field gradients may also cause formation of the Minkowski domains).

The behaviour of the experimentally measured scale factor in Fig. 3(b) may be compared to the theoretical prediction based on Eqs (^{2} and ^{4}). We may assume that −*ε*
_{m}*ε*
_{d} (since cobalt is a good metal in the infrared range), and that *α(Η,Τ)* is small (since it cannot be larger than *α*
_{∞}=8.2%). In such a case the effective metric coefficients are:

$${g}_{11}={g}_{22}=-\phantom{\rule{-.25em}{0ex}}{\epsilon}_{2}\approx -\phantom{\rule{-.25em}{0ex}}{\alpha}_{\infty}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right){\epsilon}_{m}-{\epsilon}_{d}$$

13

$${g}_{00}=-\phantom{\rule{-.25em}{0ex}}{\epsilon}_{1}\approx -\phantom{\rule{-.25em}{0ex}}{\epsilon}_{d}\left(1+2{\alpha}_{\infty}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)\right)$$

14

The effective spacetime appears to be a Minkowski one if the temperature is low enough, or the magnetic field is strong enough, so that

$${\alpha}_{\infty}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)>\frac{{\epsilon}_{d}}{(\phantom{\rule{-.25em}{0ex}}-\phantom{\rule{-.25em}{0ex}}{\epsilon}_{m})}$$

15

As evident from Fig. 3(b), the measured data for several domain walls show good agreement with a theoretical fit obtained based on Eq. (^{13}) using the optical properties of cobalt^{12} at *λ*=1.5μm. As the ferrofluid temperature falls away from the *T*
_{c} boundary at z=0, Eq. (^{13}) predicts that the scale factor of the effective Minkowski spacetime increases sharply as a function of z. It is remarkable that this “cosmological” spacetime expansion may be visualized directly using an optical microscope. It is also interesting to note that the physical vacuum appears to exhibit hyperbolic metamaterial properties when subjected to a very strong magnetic field^{17, 18}. Thus, the ferrofluid in an external magnetic field provides us with an easily accessible experimental system, which may be used to illustrate some properties of vacuum and the hypothesized physics of multiverse. Below we will demonstrate that it also lets us trace the emergence of analogue gravity and the modified Newtonian dynamics (MOND) from the well-understood microscopic degrees of freedom.

Since the language of effective spacetime appears to be quite efficient at describing light propagation inside the hyperbolic ferrofluid, it is also natural to express nonlinear optical interaction of photons via the Kerr effect, as “bending” of this spacetime, which results in effective gravitational interaction between extraordinary photons^{4}. Indeed, when the nonlinear optical effects become important, the dielectric tensor of the ferrofluid may be written as

$${\epsilon}_{ij}={\chi}_{ij}^{(1)}+{\chi}_{ijlm}^{(3)}{E}_{l}{E}_{m}+\mathrm{...},$$

16

where the second term provides gravity-like coupling between the effective metric and the energy-momentum tensor (note that due to axial symmetry of the wire array geometry, the term proportional to *χ*
^{(2)} is missing in Eq. (^{16})). It was demonstrated in ref. ^{4}. that in the weak field (Newtonian) limit the effective gravitational constant may be expressed as

$$\gamma =\frac{{c}^{2}{\omega}^{2}{\epsilon}_{2}}{4}{\chi}^{(3)}$$

17

Eq. (^{17}) establishes connection between the effective gravitational constant *γ* and the third order nonlinear susceptibility *χ*
^{(3)} of the ferrofluid. Since *ε*
_{2}<0, the sign of *χ*
^{(3)} must be negative for the effective gravity to be attractive. This condition is satisfied naturally in most liquids, and in particular, in kerosene. Because of the large and negative thermo-optic coefficient inherent to most liquids, heating produced by partial absorption of the propagating beam translates into a significant decrease of the refractive index at higher light intensity. The thermal origin of this effective gravity looks interesting in light of the modern advances in gravitation theory^{19, 20}, which strongly indicate that the classic general relativity description of gravity results from thermodynamic effects. It is also interesting to note that variations of the dielectric permittivity component *ε*
_{2} within the individual Minkowski domains (or “universes”) described by Eqs (^{3} and ^{13}) lead to different strength of the effective gravity in different domains.

Let us also analyze variations of the effective modified Newtonian dynamics (MOND)^{21}, which arise naturally in different Minkowski domains inside the ferrofluid. The effective “second law of Newtonian mechanics” for photons inside a hyperbolic ferrofluid domain may be derived based on the extraordinary photon dispersion law

$$\frac{{\omega}^{2}}{{c}^{2}}=\frac{{k}_{z}^{2}}{{\epsilon}_{1}}-\frac{{k}_{x}^{2}+{k}_{y}^{2}}{(\phantom{\rule{-.25em}{0ex}}-\phantom{\rule{-.25em}{0ex}}{\epsilon}_{2})}$$

18

which follows from Eq. (^{6}), and on the explicit expressions for *ε*
_{1} and *ε*
_{2} given by Eqs (^{13} and ^{14}). These equations result in the following expression for *k*
_{z}
^{2}:

$${k}_{z}^{2}=\frac{{\omega}^{2}}{{c}^{2}}{\epsilon}_{d}+2\frac{{\omega}^{2}}{{c}^{2}}{\epsilon}_{d}{\alpha}_{\infty}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)+({k}_{x}^{2}+{k}_{y}^{2}){\epsilon}_{d}\left(\frac{1+2{\alpha}_{\infty}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)}{-{\alpha}_{\infty}{\epsilon}_{m}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)-{\epsilon}_{d}}\right)$$

19

Note that *E* = *ℏ**c**k*_{z} plays the role of energy, while ${\overrightarrow{P}}^{2}={\hslash}^{2}({k}_{x}^{2}+{k}_{y}^{2})$ plays the role of the momentum squared in an effective 2+1 dimensional Minkowski spacetime introduced in Section 2, while the effective mass of the extraordinary photon appears to be $m=\hslash \omega {\epsilon}_{d}^{1/2}/{c}^{2}$. Using these newly introduced variables, Eq. (^{19}) may be re-written as

$${E}^{2}={m}^{2}{c}^{4}\left(1+2{\alpha}_{\infty}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)\right)+{\overrightarrow{P}}^{2}{c}^{2}\left(\frac{1+2{\alpha}_{\infty}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)}{-{\alpha}_{\infty}{\epsilon}_{m}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)-{\epsilon}_{d}}\right)$$

20

Further introducing a “non-relativistic” effective energy *E**−*mc*
^{2}, Eq. (^{20}) may be re-written in the “non-relativistic form” as

$${E}^{\u204e}=m{c}^{2}{\alpha}_{\infty}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)+\frac{{\overrightarrow{P}}^{2}}{2m}\left(\frac{1+2{\alpha}_{\infty}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)}{-{\alpha}_{\infty}{\epsilon}_{m}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)-{\epsilon}_{d}}\right)$$

21

The first term in this equation plays the role of the “gravitational potential energy” *mϕ*, which makes it apparent that the spatial gradients of *T* and *H* behave as an effective gravitational field in this model. The second term may be identified as the “kinetic energy” of the photon. The “Newton’s second law” may be obtained using the Hamilton’s equation as

$$\overrightarrow{a}=\frac{d}{d{t}^{\u204e}}\left(\frac{d{E}^{\u204e}}{d\overrightarrow{P}}\right)=\frac{d}{md{t}^{\u204e}}\left(\overrightarrow{P}\left(\frac{1+2{\alpha}_{\infty}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)}{-{\alpha}_{\infty}{\epsilon}_{m}\phantom{\rule{.25em}{0ex}}\mathrm{tanh}\left(\frac{\mu H}{kT}\right)-{\epsilon}_{d}}\right)\right)$$

22

where the role of effective time *t** is played by the z coordinate. Since *T* and *H* may depend on z within the individual Minkowski domains, different versions of MOND in these Minkowski domains (or “universes”) will be different from each other, and will not coincide with the conventional *F*=*ma* form of the second law. Moreover if we re-write Eq. (^{22}) via the effective gravitational potential *ϕ* as

$$\frac{d}{d{t}^{\u204e}}\left(\overrightarrow{P}\left(\frac{1+\frac{2\phi}{{c}^{2}}}{-{\epsilon}_{m}\frac{\phi}{{c}^{2}}-{\epsilon}_{d}}\right)\right)=m\overrightarrow{a}$$

23

and recall that according to the virial theorem the average potential energy of a bound system is proportional to its average kinetic energy, it is clear that the gravitational dynamics of bound systems in the ferrofluid would differ considerably from the Newtonian limit if *d**ϕ*/*d**t*^{⁎}~*d*(tanh(*μ**H*/*k**T*))/*d**z* ≠ 0. The latter condition is always satisfied in the presence of “cosmological expansion” near the Minkowski domain boundaries.

In particular, let us consider the rotation curves of test objects (extraordinary light beams) in the effective gravitational field of a central massive body (an intense extraordinary light beam) in an experimental geometry shown in Fig. 4. As has been noted above, in our case the role of the “central mass” *M* and the “test bodies” *m* are played by the light beams, which trajectories behave as world lines in the 2+1 dimensional effective spacetime inside the ferrofluid. In Newtonian 2D gravity the rotation curves are “flat”:

$$\frac{m{v}^{2}}{r}=\gamma \frac{mM}{r}$$

24

so that

$$v=\sqrt{\gamma M}$$

25

and does not depend on *r*. Variations of *ϕ* in Eq. (^{23}) as a function of *t**=*z* and *r* near the Minkowski domain boundaries will obviously change these flat rotation curves. However, the 2+1 dimensional character of the effective Minkowski spacetime in the ferrofluid makes the MOND analogies limited, and does not allow realistic modeling and study of such important issues as galaxy rotation curves in the 3+1 dimensional spacetime.

Let us now consider the difference in radiation lifetimes of excited atoms and molecules inside different Minkowski domains, which is another example of different “laws of physics” which operate in different regions of the ferrofluid-based “metamaterial multiverse”. When the ferrofluid “multiverse” is populated with atomic or molecular species, and these species are excited using an external laser source, the radiation lifetimes of atoms and molecules in these “universes” depend strongly on the individual physical properties of each “universe” via the Purcell effect^{22, 23}. Some “universes” are better fine-tuned than others to sustain the excited states of these species. Thus, the ferrofluid-based metamaterial “multiverse” may be used to illustrate the fine-tuning mechanism in cosmology.

The Purcell effect is the enhancement of a fluorescent molecule’s spontaneous emission rate by its environment, which is determined mainly by the enhancement of the photonic density of states. According to the Fermi’s golden rule, the transition rate between the excited and the ground state of an atom or molecule is proportional to the density of final states:

$${\Gamma}_{i\to f}=\frac{2\pi}{\hslash}{\left|\langle f\left|{H}_{\mathrm{int}}\right|i\rangle \right|}^{2}\rho ,$$

26

where *H*
_{int} is the perturbing Hamiltonian, and *ρ* is the density of final states. As has been shown in^{6}, the density of photonic states diverges in the hyperbolic metamaterials. This follows from the extraordinary photon dispersion law (Eq. (^{18})), which describes a hyperboloid in the phase space. Since the absolute value of the photon k-vector is not limited, the phase space volume between two hyperboloids corresponding to different values of frequency is formally infinite in the effective medium approximation. This divergence leads to a very large density of photonic states

$${\Gamma}^{meta}\approx \frac{{\mu}^{2}{k}_{\mathrm{max}}^{3}}{2\hslash}\frac{2\sqrt{-{\epsilon}_{1}{\epsilon}_{2}}}{(1-{\epsilon}_{1}{\epsilon}_{2})}$$

27

for every frequency where different components of the dielectric permittivity have opposite signs^{22}. Since, the effective medium description eventually fails at the point when the wavelength of the propagating mode becomes comparable to the distance between the nanocolumns *d* (see Fig. 1(b)), this introduces a wave number cut-off at

28

of this broadband divergence. According to Fig. 2(b) this cut-off occurs at d~0.2μmμm.

Since its theoretical prediction in ref. ^{22}, the broadband Purcell effect has been indeed observed in multiple experiments (see for example ref. ^{23}). As demonstrated in Fig. 5, narrow absorption lines of the ferrofluid in the long wavelength infrared (LWIR) range also exhibit the broadband Purcell effect. The inset in Fig. 5 shows FTIR transmission spectrum of the diluted (*α*
_{∞}=0.8%) ferrofluid in zero magnetic field. It reveals strong kerosene absorption lines, which separate low absorption regions of the LWIR spectrum where the Purcell effect measurements have been performed. Comparison of the FTIR transmission spectra of the ferrofluid with and without external magnetic field shows clear influence of the Purcell effect on the narrow absorption lines of the ferrofluid, which are marked by green arrows. It is clear from Fig. 5 that application of the external magnetic field, which leads to formation of the hyperbolic metamaterial structure, makes the radiation lifetime shorter, leading to observation of broader and less pronounced absorption lines. For example, the broadening of the 8.5μm absorption line corresponds to ~2 times shorter radiation lifetime. The more pronounced Purcell effect at longer wavelengths may be explained by better plasmonic properties of the cobalt nanoparticles at these wavelengths.

Comparison of the FTIR transmission spectra of the ferrofluid with and without external magnetic field shows clear influence of the Purcell effect on the narrow absorption lines of the ferrofluid. Application of the external magnetic field, which leads **...**

Eqs (^{27}, ^{28}) demonstrate that the Purcell effect in different Minkowski domains inside the ferrofluid is affected by such parameters as local variations of cobalt nanoparticle concentration *α*, external magnetic field *H* and temperature *T* (via their influence on *ε*
_{1} and *ε*
_{2} – see Eqs (^{3}, ^{4})), as well as local nanocolumn periodicity via the wave vector cutoff. The Minkowski domain size also affects the Purcell effect via the effective cavity dimensions. Considerable variations of the radiation lifetime due to variations in Purcell effect in different Minkowski domains lead to variations of IR absorption, and hence transmission of the domains, which can be easily observed using IR microscopy.

For experimental realization of such microscopic observations we have chosen to operate our setup at 1550nm using a readily available laser source in transmission mode. A FLIR infrared camera has been used to capture the microscopic images as a function of external magnetic field. This spectral range appears to be quite convenient since it is located above the metric signature transition in cobalt ferrofluid at large magnetic fields (see Fig. 2a), it guarantees sufficient spatial resolution for Minkowski domain observation, and because of availability of common C-H, O-H, and N-H absorption bands near 1550nm. The common C-H absorption bands are located around 1600nm, N-H bands are located around 1550nm, and O-H bands are located around 1450nm, so that IR spectroscopy at these bands is widely used to detect fatty acids^{24}, which are used as the cobalt nanoparticle coating. As shown in Fig. 6, experimental microscopic images observed at 1550nm indeed exhibit considerable contrast between different Minkowski domains. For example, cross section of the microscopic image in Fig. 6b along the yellow line, which is shown in Fig. 6c shows more than two fold variation in Minkowski domain transmission. The observed difference in infrared absorption between different domains indicates at least two fold difference in radiation lifetime between the darker and brighter looking domains. The excited states of the fatty acids appear to have considerably longer lifetime in brighter domains, so that these brighter “universes” appear to be better fine-tuned for their existence. This effect appears to be somewhat analogous to the proposed fine tuning mechanism in multiverse cosmology.

In conclusion, we have demonstrated that ferrofluid-based self-assembled metamaterial geometry captures many features of several cosmological models of the multiverse, such as metric signature transition scenario in loop quantum cosmology^{1, 2}, natural emergence of a large number of Minkowski universes, and inflation. Moreover, we have demonstrated that individual Minkowski domains in the ferrofluid exhibit different “laws of physics”, such as different versions of modified Newtonian dynamics (MOND), different strength of effective gravity, and different radiation lifetimes due to variations in the local Purcell effect, so that ferrofluids may be used to illustrate the fine-tuning mechanism in cosmology. All these effects may be studied via direct microscopic observations. While these effects were observed in cobalt nanoparticle-based ferrofluids, other ferrofluids (for example magnetite-based) may also be used leading to qualitatively similar results. However, since conductivity of magnetite is smaller than the conductivity of cobalt, a magnetite-based ferrofluid would become hyperbolic at lower frequencies in the long wavelength infrared range, making experimental observations more difficult.

Unfortunately, the described analogy between the extraordinary light propagation inside the ferrofluid and the dynamics of massive particles in Minkowski spacetime is far from being perfect. The main difficulty comes from the cross-coupling between extraordinary and ordinary light inside the ferrofluid, which may be caused by domain interfaces and internal defects. Since ordinary light does not obey the same wave Eq. (^{6}), such a cross-coupling breaks the effective Lorentz symmetry (7) of the system. In addition, such a model is necessarily limited to 2+1 spacetime dimensions. Nevertheless, despite these limitations the developed condensed matter model of the cosmological multiverse appears to be quite interesting, since it is able to replicate many of its hypothesized features in the laboratory setting.

This work was supported in part by NSF grant DMR-1104676 and the School of Emerging Technologies at Towson University.

Author Contributions

I.S. and V.S. collected experimental data and wrote the manuscript text.

The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

1. Mielczarek, J. Signature change in loop quantum cosmology, arXiv:1207.4657 [gr-qc].

2. Bojowald M, Mielczarek J. Some implications of signature-change in cosmological models of loop quantum gravity. J. of Cosmology and Astroparticle Phys. 2015;08 doi: 10.1088/1475-7516/2015/08/052. [Cross Ref]

3. Smolyaninova VN, Yost B, Lahneman D, Narimanov E, Smolyaninov II. Self-assembled tunable photonic hyper-crystals. Scientific Reports. 2014;4 doi: 10.1038/srep05706. [PMC free article] [PubMed] [Cross Ref]

4. Smolyaninov II. Analogue gravity in hyperbolic metamaterials. Phys. Rev. A. 2013;88 doi: 10.1103/PhysRevA.88.033843. [Cross Ref]

5. Smolyaninov II, Hung YJ. Modeling of time with metamaterials. JOSA B. 2011;28:1591–1595. doi: 10.1364/JOSAB.28.001591. [Cross Ref]

6. Smolyaninov II, Narimanov EE. Metric signature transitions in optical metamaterials. Phys. Rev. Letters. 2010;105 doi: 10.1103/PhysRevLett.105.067402. [PubMed] [Cross Ref]

7. Smolyaninov II, Yost B, Bates E, Smolyaninova VN. Experimental demonstration of metamaterial “multiverse” in a ferrofluid. Optics Express. 2013;21:14918–14925. doi: 10.1364/OE.21.014918. [PubMed] [Cross Ref]

8. Smolyaninov II, Hung YJ, Hwang E. Experimental modeling of cosmological inflation with metamaterials. Physics Letters A. 2012;376:2575–2579. doi: 10.1016/j.physleta.2012.07.010. [Cross Ref]

9. Kittel, C. *Introduction to Solid State Physics* (Wiley, New York, 2004).

10. Wangberg R, Elser J, Narimanov EE, Podolskiy VA. Nonmagnetic nanocomposites for optical and infrared negative-refractive-index media. J. Opt. Soc. Am. B. 2006;23:498–505. doi: 10.1364/JOSAB.23.000498. [Cross Ref]

11. Gao Y, et al. Optical negative refraction in ferrofluids with magnetocontrollability. Phys. Rev. Letters. 2010;104 doi: 10.1103/PhysRevLett.104.034501. [PubMed] [Cross Ref]

12. Johnson PB, Christy RW. Optical constants of transition metals: Ti, V, Cr, Mn, Fe, Co, Ni, and Pd. Phys. Rev. B. 1974;9 doi: 10.1103/PhysRevB.9.5056. [Cross Ref]

13. Jacob Z, Alekseyev LV, Narimanov E. Optical hyperlens: Far-field imaging beyond the diffraction limit. Optics Express. 2006;14:8247–8256. doi: 10.1364/OE.14.008247. [PubMed] [Cross Ref]

14. Smolyaninov II, Hung YJ, Davis CC. Magnifying superlens in the visible frequency range. Science. 2007;315:1699–1701. doi: 10.1126/science.1138746. [PubMed] [Cross Ref]

15. Smolyaninov II, Hwang E, Narimanov EE. Hyperbolic metamaterial interfaces: Hawking radiation from Rindler horizons and spacetime signature transtions. Phys. Rev. B. 2012;85 doi: 10.1103/PhysRevB.85.235122. [Cross Ref]

16. Davydov, A. S. *Quantum Mechanics* (Pergamon, 1976).

17. Smolyaninov II. Vacuum in strong magnetic field as a hyperbolic metamaterial. Phys. Rev. Letters. 2011;107 doi: 10.1103/PhysRevLett.107.253903. [PubMed] [Cross Ref]

18. Chernodub MN, Van Doorsselaere J, Verschelde H. Spontaneous electromagnetic superconductivity and superfluidity of QCDxQED vacuum in strong magnetic field. Proceedings of Science. 2012;QNP2012

19. Jacobson T. Thermodynamics of spacetime: The Einstein equation of state. Phys. Rev. Letters. 1995;75 doi: 10.1103/PhysRevLett.75.1260. [PubMed] [Cross Ref]

20. Verlinde EP. On the origin of gravity and the laws of Newton. J of High Energy Phys. 2011;2011 doi: 10.1007/JHEP04(2011)029. [Cross Ref]

21. Milgrom M. A modification of the Newtonian dynamics as a possible alternative to the hidden mass hypothesis. Astrophysical Journal. 1983;270:365–370. doi: 10.1086/161130. [Cross Ref]

22. Jacob Z, Smolyaninov II, Narimanov EE. Broadband Purcell effect: Radiative decay engineering with metamaterials. Appl. Phys. Lett. 2012;100 doi: 10.1063/1.4710548. [Cross Ref]

23. Jacob Z, et al. Engineering photonic density of states using metamaterials. App. Phys. B. 2010;100 doi: 10.1007/s00340-010-4096-5. [Cross Ref]

24. Robinson, J. W., Skelly, E. M. & Frame, G. M. *Undergraduate Instrumental Analysis* (CRC Press, 2004).

Articles from Scientific Reports are provided here courtesy of **Nature Publishing Group**