We investigate the nonlinear dynamics of cold atom systems that can in principle serve as quantum simulators of false vacuum decay. The analog false vacuum manifests as a metastable vacuum state for the relative phase in a two-species Bose-Einstein condensate (BEC), induced by a driven periodic coupling between the two species. In the appropriate low energy limit, the evolution of the relative phase is approximately governed by a relativistic wave equation exhibiting true and false vacuum configurations. In previous work, a linear stability analysis identified exponentially growing short-wavelength modes driven by the time-dependent coupling. These modes threaten to destabilize the analog false vacuum. Here, we employ numerical simulations of the coupled Gross-Pitaevski equations (GPEs) to determine the non-linear evolution of these linearly unstable modes. We find that unless a physical mechanism modifies the GPE on short length scales, the analog false vacuum is indeed destabilized. We briefly discuss various physically expected corrections to the GPEs that may act to remove the exponentially unstable modes. To investigate the resulting dynamics in cases where such a removal mechanism exists, we implement a hard UV cutoff that excludes the unstable modes as a simple model for these corrections. We use this to study the range of phenomena arising from such a system. In particular, we show that by modulating the strength of the time-dependent coupling, it is possible to observe the crossover between a second and first order phase transition out of the false vacuum.

We introduce a new picture of vacuum decay which, in contrast to existing semiclassical techniques, provides a real-time description and does not rely on classically-forbidden tunneling paths. Using lattice simulations, we observe vacuum decay via bubble formation by generating realizations of vacuum fluctuations and evolving with the classical equations of motion. The decay rate obtained from an ensemble of simulations is in excellent agreement with existing techniques. Future applications include bubble correlation functions, fast decay rates, and decay of non-vacuum states.

Analog condensed matter systems present an exciting opportunity to simulate early Universe models in table-top experiments. We consider a recent proposal for an analog condensed matter experiment to simulate the relativistic quantum decay of the false vacuum. In the proposed experiment, two ultra-cold condensates are coupled via a time-varying radio-frequency field. The relative phase of the two condensates in this system is approximately described by a relativistic scalar field with a potential possessing a series of false and true vacuum local minima. If the system is set up in a false vacuum, it would then decay to a true vacuum via quantum mechanical tunnelling. Should such an experiment be realized, it would be possible to answer a number of open questions regarding non-perturbative phenomena in quantum field theory and early Universe cosmology. In this paper, we illustrate a possible obstruction: the time-varying coupling that is invoked to create a false vacuum for the long-wavelength modes of the condensate leads to a destabilization of shorter wavelength modes within the system via parametric resonance. We focus on an idealized setup in which the two condensates have identical properties and identical background densities. Describing the system by the coupled Gross-Pitaevskii equations (GPE), we use the machinery of Floquet theory to perform a linear stability analysis, calculating the wavenumber associated with the first instability band for a variety of experimental parameters. However, we demonstrate that, by tuning the frequency of the time-varying coupling, it may be possible to push the first instability band outside the validity of the GPE, where dissipative effects are expected to damp any instabilities. This provides a viable range of experimental parameters to perform analog experiments of false vacuum decay.

Gravitational lensing in a weak but otherwise arbitrary gravitational field can be described in terms of a 3 x 3 tensor, the "effective refractive index". If the sources generating the gravitational field all have small internal fluxes, stresses, and pressures, then this tensor is automatically isotropic and the "effective refractive index" is simply a scalar that can be determined in terms of a classic result involving the Newtonian gravitational potential. In contrast if anisotropic stresses are ever important then the gravitational field acts similarly to an anisotropic crystal. We derive simple formulae for the refractive index tensor, and indicate some situations in which this will be important.