BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Eugene A. Lim appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
GRChombo is an open-source code for performing Numerical Relativity time evolutions, built on top of the publicly available Chombo software for the solution of PDEs. Whilst GRChombo uses standard techniques in NR, it focusses on applications in theoretical physics where adaptability, both in terms of grid structure, and in terms of code modification, are key drivers.
We demonstrate the flexibility and utility of the Berger-Rigoutsos Adaptive Mesh Refinement (AMR) algorithm used in the open-source numerical relativity code GRChombo for generating gravitational waveforms from binary black-hole inspirals, and for studying other problems involving non-trivial matter configurations. We show that GRChombo can produce high quality binary black-hole waveforms through a code comparison with the established numerical relativity code Lean. We also discuss some of the technical challenges involved in making use of full AMR (as opposed to, e.g. moving box mesh refinement), including the numerical effects caused by using various refinement criteria when regridding. We suggest several "rules of thumb" for when to use different tagging criteria for simulating a variety of physical phenomena. We demonstrate the use of these different criteria through example evolutions of a scalar field theory. Finally, we also review the current status and general capabilities of GRChombo.
We study the formation of black holes from subhorizon and superhorizon perturbations in a matter dominated universe with 3+1D numerical relativity simulations. We find that there are two primary mechanisms of formation depending on the initial perturbation's mass and geometry -- via $\textit{direct collapse}$ of the initial overdensity and via $\textit{post-collapse accretion}$ of the ambient dark matter. In particular, for the latter case, the initial perturbation does not have to satisfy the hoop conjecture for a black hole to form. In both cases, the duration of the formation the process is around a Hubble time, and the initial mass of the black hole is $M_\mathrm{BH} \sim 10^{-2} H^{-1} M_\mathrm{Pl}^2$. Post formation, we find that the PBH undergoes rapid mass growth beyond the self-similar limit $M_\mathrm{BH}\propto H^{-1}$, at least initially. We argue that this implies that most of the final mass of the PBH is accreted from its ambient surroundings post formation.
Through numerical simulations of boson-star head-on collisions, we explore the quality of binary initial data obtained from the superposition of single-star spacetimes. Our results demonstrate that evolutions starting from a plain superposition of individual boosted boson-star spacetimes are vulnerable to significant unphysical artefacts. These difficulties can be overcome with a simple modification of the initial data suggested in [PRD 99 (2018) 044046] for collisions of oscillatons. While we specifically consider massive complex scalar field boson star models up to a 6th-order-polynomial potential, we argue that this vulnerability is universal and present in other kinds of exotic compact systems and hence needs to be addressed.
We show that the so-called "post-Born" effects of weak lensing at 4th order are equivalent to lens-lens couplings in the Born Approximation. We demonstrate this by explicitly showing the equivalence of the canonical weak lensing approach at 4th order using the anisotropy remapping method, to that of the 4th order calculation of the lens-lens coupling effects using the Boltzmann equation approach that was first developed in [Phys. Rev. D89, 123006]. Furthermore, we argue that to incorporate true "post-Born" effects, i.e. taking into account non-straight photon paths, require the addition of a photon deflection term which has not been taken into account in the canonical formalism nor the Boltzmann method.
We study the robustness of single-field inflation against inhomogeneities. We derive a simple analytic criterion on the shape of the potential for successful inflation in the presence of inhomogeneities, and demonstrate its validity using full 3+1 dimensional numerical relativity simulations on several classes of popular models of single-field inflation. We find that models with convex potentials are more robust to inhomogeneities than those with concave potentials, and that concave potentials that vary on super-Planckian scales are significantly more robust than those that vary on sub-Planckian scales.
We investigate the physics of black hole formation from the head-on collisions of boosted equal mass Oscillatons (OS) in full numerical relativity, for both the cases where the OS have equal phases or are maximally off-phase (anti-phase). While unboosted OS collisions will form a BH as long as their initial compactness $\mathcal{C}\equiv GM/R$ is above a numerically determined critical value $\mathcal{C}>0.035$, we find that imparting a small initial boost counter-intuitively \emph{prevents} the formation of black holes even if $\mathcal{C}> 0.035$. If the boost is further increased, at very high boosts $γ>1/12\mathcal{C}$, BH formation occurs as predicted by the hoop conjecture. These two limits combine to form a "stability band" where collisions result in either the OS "passing through" (equal phase) or "bouncing back" (anti-phase), with a critical point occurring around ${\cal C}\approx 0.07$. We argue that the existence of this stability band can be explained by the competition between the free fall and the interaction timescales of the collision.
We investigate numerically the critical behaviour which occurs in the collapse of both spherically symmetric and asymmetric scalar field bubbles with full general relativity. We use a minimally coupled scalar field subject to a "double well" interaction potential, with the bubble wall spanning the barrier between two degenerate minima. We find that the spherically symmetric case exhibits Type 2 critical behaviour with the critical index consistent with a value of gamma = 0.37 as expected. In the asymmetric case, we find that again our results are consistent with a value of gamma = 0.37 for the dominant unstable mode. We do not see strong evidence of echoing in the solutions, which could be due to the fact that the coordinates are not well adapted to the echoing behaviour, or due to being still too far from the critical point to properly observe the critical solution.
In this work, we introduce GRChombo: a new numerical relativity code which incorporates full adaptive mesh refinement (AMR) using block structured Berger-Rigoutsos grid generation. The code supports non-trivial "many-boxes-in-many-boxes" mesh hierarchies and massive parallelism through the Message Passing Interface (MPI). GRChombo evolves the Einstein equation using the standard BSSN formalism, with an option to turn on CCZ4 constraint damping if required. The AMR capability permits the study of a range of new physics which has previously been computationally infeasible in a full 3+1 setting, whilst also significantly simplifying the process of setting up the mesh for these problems. We show that GRChombo can stably and accurately evolve standard spacetimes such as binary black hole mergers and scalar collapses into black holes, demonstrate the performance characteristics of our code, and discuss various physics problems which stand to benefit from the AMR technique.
Two of the most commonly used tools to constrain the primordial non-Gaussianity are the bispectrum and the Minkowski functionals of CMB temperature anisotropies. These two measures of non-Gaussianity in principle provide distinct (though correlated) information, but in the past constraints from them have only been loosely compared, and not statistically combined. In this work we evaluate, for the first time, the covariance matrix between the local non-Gaussianity coefficient fnl estimated through the bispectrum and Minkowski functionals. We find that the estimators are positively correlated, with corerlation coefficient r ~ 0.3. Using the WMAP7 data to combine the two measures and accounting for the point-source systematics, we find the combined constraint fnl=37+/-28, which has a ~20% smaller error than either of the individual constraints.
It has been shown that the primordial perturbations sourced by inflation are driven to classicality by unitary evolution alone. However, their coupling with the environment such as photons and subsequent decoherence renders the cosmological correlations quantum, losing primordial information in the process. We argue that the quantumness of the resulting cosmological correlations is given by quantum discord, which captures non-classical behavior beyond quantum entanglement. By considering the environment as a quantum channel in which primordial information contained in the perturbations is transmitted to us, we can then ask how much of this information is inaccessible. We show that this amount of information is given by the discord of the joint primordial perturbations-environment system. To illustrate these points, we model the joint system as a mixed bi-modal Gaussian state, and show that quantum discord is dependent on the basis which decoherence occurs.
The Planck mission has conclusively detected lensing of the Cosmic Microwave Background (CMB) radiation from foreground sources to an overall significance of greater than $25σ$. The high precision of this measurement motivates the development of a more complete formulation of the calculation of this effect. While most effects on the CMB anisotropies are widely studied through direct solutions of the Boltzmann equation, the non-linear effect of CMB lensing is formulated through the solutions of the geodesic equation. In this paper, we present a new formalism to the calculation of the lensing effect by \emph{directly solving the Boltzmann equation}, as we did in the calculation of the CMB anisotropies at recombination. In particular, we developed a diagrammatic approach to efficiently keep track of all the interaction terms and calculate all possible non-trivial correlations to arbitrary high orders. Using this formalism, we explicitly articulate the approximations required to recover the usual remapping approach used in current studies of the weak lensing. In addition, we point out additional unexplored corrections that are manifest in our formalism to which experiments may be sensitive. As an example, we calculate the correction to the CMB temperature power spectrum for the \emph{lens-lens} coupling effects which are neglected in standard calculations. We find that the correction is $\lesssim 0.1\%$ of the CMB temperature power spectrum for $\ell$ up to 3000 and thus is comparable to the cosmic variance.
We derive a closed-form expression for the phase shift experienced by 1+1 dimensional kinks colliding at ultra-relativistic velocities (gamma v >> 1), valid for arbitrary periodic potentials. Our closed-form expression is the leading order result of a more general scattering theory of solitary waves described in a companion paper [Amin, Lim and Yang, arXiv:1308.0606]. This theory relies on a small kinematic parameter 1/(gamma v) << 1 rather than a small parameter in the Lagrangian. Our analytic results can be directly extracted from the Lagrangian without solving the equation of motion. Based on our closed-form expression, we prove that kink-kink and kink-antikink collisions have identical phase shifts at leading order.
We construct a perturbative framework for understanding the collision of solitons (more precisely, solitary waves) in relativistic scalar field theories. Our perturbative framework is based on the suppression of the space-time interaction area proportional to $1/(γv)$, where $v$ is the relative velocity of an incoming solitary wave and $γ= 1/\sqrt{1-v^2} \gg 1$. We calculate the leading order results for collisions of (1+1) dimensional kinks in periodic potentials, and provide explicit, closed form expressions for the phase shift and the velocity change after the collisions. We find excellent agreement between our results and detailed numerical simulations. Crucially, our perturbation series is controlled by a kinematic parameter, and hence not restricted to small deviations around integrable cases such as the Sine-Gordon model.
We study cosmological perturbations by solving the governing Boltzmann and Einstein Field equations up to second order, and calculate the corresponding CMB bispectrum during recombination. We include all the second-order Liouville and collision terms, truncating the multipole hierarchy at $l=10$, consistently including all $m\neq 0$ terms when calculating the bispectrum in the flat-sky limit. At this stage, we focus on contributions at recombination and we neglect 2nd-order vector and tensor terms, lensing effects, and late time non-linear ISW. We find that the signal-to-noise for the bispectrum is 0.69 for $l_{max}=2000$, yielding an overall signal $F_{\rm NL} = 3.19$ (normalised relative to the local model). We find that the effective $f_{NL}$'s of the equilateral and local type are 5.11 and 0.88 respectively. This recombination signal will have to be taken into account in a quantitative analysis of the Planck data.
We prove that the linear instability in a non-degenerate higher derivative theory, the Ostrogradski instability, can only be removed by the addition of constraints if the original theory's phase space is reduced.
We show that the instabilities of higher derivative gravity models with quadratic curvature invariant $αR^2+βR_{μν}R^{μν}$ can be removed by judicious addition of constraints at the quadratic level of metric fluctuations around Minkowski/de Sitter background. With a suitable parameter choice, we find that the instabilities of helicity-0, 1, 2 modes can be removed while reducing the dimensionality of the original phase space. To retain the renormalization properties of higher derivative gravity, Lorentz symmetry in the constrained theory is explicitly broken.
We investigate the effect of modifed gravity with screening mechanisms, such as the chameleon or symmetron models, upon the structure of main sequence stars. We find that unscreened stars can be significantly more luminous and ephemeral than their screened doppelgangers. By embedding these stars into dwarf galaxies, which can be unscreened for values of the parameters not yet ruled out observationally, we show that the cumulative effect of their increased luminosity can enhance the total galactic luminosity. We estimate this enhancement and find that it can be considerable given model parameters that are still under experimental scrutiny. By looking for systematic offsets between screened dwarf galaxies in clusters and unscreened galaxies in voids, these effects could form the basis of an independent observational test that can potentially lower the current experimental bounds on the model independent parameters of these theories by and order of magnitude or more.
We provide analytic solutions for the power spectrum and bispectrum of curvature fluctuations produced by a step feature in the inflaton potential, valid in the limit that the step is short and sharp. In this limit, the bispectrum is strongly scale dependent and its effective non-linearity attains a large oscillatory amplitude. The perturbations to the curvature power spectrum, on the other hand, remain a small component on top of the usual spectrum of fluctuations generated by slow roll. We utilize our analytic solutions to assess the observability of the predicted non-Gaussian signatures and show that, if present, only very sharp steps on scales larger than ~ 2 Gpc are likely to be able to be detected by Planck. Such features are not only consistent with WMAP7 data, but can also improve its likelihood by 2 Delta ln L ~ 12 for two extra parameters, the step location and height. If this improvement were due to a slow roll violating step as considered here, a bispectrum or corresponding polarization power spectrum detection would provide definitive checks as to its primordial origin.
In this paper, we investigate the utility of Minkowski Functionals as a probe of cold/hot disk-like structures in the CMB. In order to construct an accurate estimator, we resolve a long-standing issue with the use of Minkowski Functionals as probes of the CMB sky -- namely that of systematic differences ("residuals") when numerical and analytical MF are compared. We show that such residuals are in fact by-products of binning, and not caused by pixelation or masking as originally thought. We then derive a map-independent estimator that encodes the effects of binning, applicable to beyond our present work. Using this residual-free estimator, we show that small disk-like effects (as claimed by Vielva et al.) can be detected only when a large sample of such maps are averaged over. In other words, our estimator is noise-dominated for small disk sizes at WMAP resolution. To confirm our suspicion, we apply our estimator to the WMAP7 data to obtain a null result.
We investigate flux vacua on a variety of one-parameter Calabi-Yau compactifications, and find many examples that are connected through continuous monodromy transformations. For these, we undertake a detailed analysis of the tunneling dynamics and find that tunneling trajectories typically graze the conifold point---particular 3-cycles are forced to contract during such vacuum transitions. Physically, these transitions arise from the competing effects of minimizing the energy for brane nucleation (facilitating a change in flux), versus the energy cost associated with dynamical changes in the periods of certain Calabi-Yau 3-cycles. We find that tunneling only occurs when warping due to back-reaction from the flux through the shrinking cycle is properly taken into account.
Recent studies have shown that the number counts of convergence peaks N(kappa) in weak lensing (WL) maps, expected from large forthcoming surveys, can be a useful probe of cosmology. We follow up on this finding, and use a suite of WL convergence maps, obtained from ray-tracing N-body simulations, to study (i) the physical origin of WL peaks with different heights, and (ii) whether the peaks contain information beyond the convergence power spectrum P_ell. In agreement with earlier work, we find that high peaks (with amplitudes >~ 3.5 sigma, where sigma is the r.m.s. of the convergence kappa) are typically dominated by a single massive halo. In contrast, medium-height peaks (~0.5-1.5 sigma) cannot be attributed to a single collapsed dark matter halo, and are instead created by the projection of multiple (typically, 4-8) halos along the line of sight, and by random galaxy shape noise. Nevertheless, these peaks dominate the sensitivity to the cosmological parameters w, sigma_8, and Omega_m. We find that the peak height distribution and its dependence on cosmology differ significantly from predictions in a Gaussian random field. We directly compute the marginalized errors on w, sigma_8, and Omega_m from the N(kappa) + P_ell combination, including redshift tomography with source galaxies at z_s=1 and z_s=2. We find that the N(kappa) + P_ell combination has approximately twice the cosmological sensitivity compared to P_ell alone. These results demonstrate that N(kappa) contains non-Gaussian information complementary to the power spectrum.
In this paper, we show that Minkowski Functionals (MFs) of weak gravitational lensing (WL) convergence maps contain significant non-Gaussian, cosmology-dependent information. To do this, we use a large suite of cosmological ray-tracing N-body simulations to create mock WL convergence maps, and study the cosmological information content of MFs derived from these maps. Our suite consists of 80 independent 512^3 N-body runs, covering seven different cosmologies, varying three cosmological parameters Omega_m, w, and sigma_8 one at a time, around a fiducial LambdaCDM model. In each cosmology, we use ray-tracing to create a thousand pseudo-independent 12 deg^2 convergence maps, and use these in a Monte Carlo procedure to estimate the joint confidence contours on the above three parameters. We include redshift tomography at three different source redshifts z_s=1, 1.5, 2, explore five different smoothing scales theta_G=1, 2, 3, 5, 10 arcmin, and explicitly compare and combine the MFs with the WL power spectrum. We find that the MFs capture a substantial amount of information from non-Gaussian features of convergence maps, i.e. beyond the power spectrum. The MFs are particularly well suited to break degeneracies and to constrain the dark energy equation of state parameter w (by a factor of ~ three better than from the power spectrum alone). The non-Gaussian information derives partly from the one-point function of the convergence (through V_0, the "area" MF), and partly through non-linear spatial information (through combining different smoothing scales for V_0, and through V_1 and V_2, the boundary length and genus MFs, respectively). In contrast to the power spectrum, the best constraints from the MFs are obtained only when multiple smoothing scales are combined.
We introduce a novel class of field theories where energy always flows along timelike geodesics, mimicking in that respect dust, yet which possess non-zero pressure. This theory comprises two scalar fields, one of which is a Lagrange multiplier enforcing a constraint between the other's field value and derivative. We show that this system possesses no wave-like modes but retains a single dynamical degree of freedom. Thus, the sound speed is always identically zero on all backgrounds. In particular, cosmological perturbations reproduce the standard behaviour for hydrodynamics with vanishing sound speed. Using all these properties we propose a model unifying Dark Matter and Dark Energy in a single degree of freedom. In a certain limit this model exactly reproduces the evolution history of Lambda-CDM, while deviations away from the standard expansion history produce a potentially measurable difference in the evolution of structure.
It has recently been shown in high resolution numerical simulations that relativistic collisions of bubbles in the context of a multi-vacua potential may lead to the creation of bubbles in a new vacuum. In this paper, we show that scalar fields with only potential interactions behave like free fields during high-speed collisions; the kick received by them in a collision can be deduced simply by a linear superposition of the bubble wall profiles. This process is equivalent to the scattering of solitons in 1+1 dimensions. We deduce an expression for the field excursion (shortly after a collision), which is related simply to the field difference between the parent and bubble vacua, i.e. contrary to expectations, the excursion cannot be made arbitrarily large by raising the collision energy. There is however a minimum energy threshold for this excursion to be realized. We verify these predictions using a number of 3+1 and 1+1 numerical simulations. A rich phenomenology follows from these collision induced excursions - they provide a new mechanism for scanning the landscape, they might end/begin inflation, and they might constitute our very own big bang, leaving behind a potentially observable anisotropy.
We consider the 3-pt function (i.e. the bispectrum or non-Gaussianity) for stochastic backgrounds of gravitational waves. We estimate the amplitude of this signal for the primordial inflationary background, gravitational waves generated during preheating, and for gravitational waves produced by self-ordering scalar fields following a global phase transition. To assess detectability, we describe how to extract the 3-pt signal from an idealized interferometric experiment and compute the signal to noise ratio as a function of integration time. The 3-pt signal for the stochastic gravitational wave background generated by inflation is unsurprisingly tiny. For gravitational radiation generated by purely causal, classical mechanisms we find that, no matter how non-linear the process is, the 3-pt correlations produced vanish in direct detection experiments. On the other hand, we show that in scenarios where the B-mode of the CMB is sourced by gravitational waves generated by a global phase transition, a strong 3-pt signal among the polarization modes could also be produced. This may provide another method of distinguishing inflationary B-modes. To carry out this computation, we have developed a diagrammatic approach to the calculation of stochastic gravitational waves sourced by scalar fluids, which has applications beyond the present scenario.
We compute the 3-point correlation function for a general model of inflation driven by a single, minimally coupled scalar field. Our approach is based on the numerical evaluation of both the perturbation equations and the integrals which contribute to the 3-point function. Consequently, we can analyze models where the potential has a "feature", in the vicinity of which the slow roll parameters may take on large, transient values. This introduces both scale and shape dependent non-Gaussianities into the primordial perturbations. As an example of our methodology, we examine the ``step'' potentials which have been invoked to improve the fit to the glitch in the $<TT>$ $C_l$ for $l \sim 30$, present in both the one and three year WMAP data sets. We show that for the typical parameter values, the non-Gaussianities associated with the step are far larger than those in standard slow roll inflation, and may even be within reach of a next generation CMB experiment such as Planck. More generally, we use this example to explain that while adding features to potential can improve the fit to the 2-point function, these are generically associated with a greatly enhanced signal at the 3-point level. Moreover, this 3-point signal will have a very nontrivial shape and scale dependence, which is correlated with the form of the 2-point function, and may thus lead to a consistency check on the models of inflation with non-smooth potentials.
We consider the long standing puzzle of how to obtain meaningful probabilities in eternal inflation. We demonstrate a new algorithm to compute the probability distribution of pocket universe types, given a multivacua inflationary potential. The computed probability distribution is finite and manifestly gauge-independent. We argue that in some scenarios this technique can be applied to disfavor some eternally inflating potentials.
We consider the gravitational effects of a single, fixed-norm, Lorentz-violating timelike vector field. In a cosmological background, such a vector field acts to rescale the effective value of Newton's constant. The energy density of this vector field precisely tracks the energy density of the rest of the universe, but with the opposite sign, so that the universe experiences a slower rate of expansion for a given matter content. This vector field similarly rescales Newton's constant in the Newtonian limit, although by a different factor. We put constraints on the parameters of the theory using the predictions of primordial nucleosynthesis, demonstrating that the norm of the vector field should be less than the Planck scale by an order of magnitude or more.
We investigate the perturbation theory of a fixed-norm, timelike Lorentz-violating vector field. After consistently quantizing the vector field to put constraints on its parameters, we compute the primordial spectra of perturbations generated by inflation in the presence of this vector field. We find that its perturbations are sourced by the perturbations of the inflaton; without the inflaton perturbation the vector field perturbations decay away leaving no primordial spectra of perturbations. Since the inflaton perturbation does not have a spin-1 component, the vector field generically does not generate any spin-1 ``vector-type'' perturbations. Nevertheless, it will modify the amplitude of both the spin-0 ``scalar-type'' and spin-2 ``tensor-type'' perturbation spectra, leading to violations of the inflationary consistency relationship.
We investigate the ability to constrain oscillatory features in the primordial power spectrum using current and future cosmic microwave background observations. In particular, we study the observability of an oscillation arising from imprints of physics at the cut-off energy scale. We perform a likelihood analysis on the WMAP data set, and find that the current data set constrains the amplitude of the oscillations to be less than 0.77 at 2-sigma, consistent with a power spectrum without oscillations. In addition, we investigate the fundamental limitations in the measurement of oscillation parameters by studying the constraints from a cosmic variance limited experiment. We find that such an experiment is capable of constraining the amplitude of such oscillations to be below 0.005, implying that reasonable models with cut-off energy scales Lambda>200 H_infl are unobservable through the microwave background.
In the presence of a short-distance cutoff, the choice of a vacuum state in an inflating, non-de Sitter universe is unavoidably ambiguous. The ambiguity is related to the time at which initial conditions for the mode functions are specified and to the way the expansion of the universe affects those initial conditions. In this paper we study the imprint of these uncertainties on the predictions of inflation. We parametrize the most general set of possible vacuum initial conditions by two phenomenological variables. We find that the generated power spectrum receives oscillatory corrections whose amplitude is proportional to the Hubble parameter over the cutoff scale. In order to further constrain the phenomenological parameters that characterize the vacuum definition, we study gravitational particle production during different cosmological epochs.