-
Astrophysical foreground cleanup using non-local means
Authors:
Guillermo F. Quispe Peña,
Andrei V. Frolov
Abstract:
To create high-fidelity cosmic microwave background maps, current component separation methods rely on availability of information on different foreground components, usually through multi-band frequency coverage of the instrument. Internal linear combination (ILC) methods provide an unbiased estimators for CMB which are easy to implement, but component separation quality crucially depends on the…
▽ More
To create high-fidelity cosmic microwave background maps, current component separation methods rely on availability of information on different foreground components, usually through multi-band frequency coverage of the instrument. Internal linear combination (ILC) methods provide an unbiased estimators for CMB which are easy to implement, but component separation quality crucially depends on the signal to noise ratio of the input maps. In the present paper, we develop an efficient non-linear filter along the lines of non-local means used in digital imaging research which significantly improves signal to noise ratio for astrophysical foreground maps, while having minimal signal attenuation, and evaluate it performance in map and spectral domains. Noise reduction is achieved by averaging ``similar'' pixels in the map. We construct the rotationally-invariant feature vector space and compute the similarity metric on it for the case of non-Gaussian signal contaminated by an additive Gaussian noise. The proposed filter has two tuneable parameters, and with minimal tweaking achieves a factor of two improvement in signal to noise spectral density in Planck dust maps. A particularly desirable feature is that signal loss is extremely small at all scales.
△ Less
Submitted 31 May, 2023;
originally announced June 2023.
-
A GPU-accelerated viewer for HEALPix maps
Authors:
Andrei V. Frolov
Abstract:
HEALPix by Górski et. al. (2005) is de-facto standard for Cosmic Microwave Background (CMB) data storage and analysis, and is widely used in current and upcoming CMB experiments. Almost all the datasets in Legacy Archive for Microwave Background Data Analysis (LAMBDA) use HEALPix as a format of choice. Visualizing the data plays important role in research, and several toolsets were developed to do…
▽ More
HEALPix by Górski et. al. (2005) is de-facto standard for Cosmic Microwave Background (CMB) data storage and analysis, and is widely used in current and upcoming CMB experiments. Almost all the datasets in Legacy Archive for Microwave Background Data Analysis (LAMBDA) use HEALPix as a format of choice. Visualizing the data plays important role in research, and several toolsets were developed to do that with HEALPix maps, most notably original Fortran facilities and Python integration with healpy.
With the current state of GPU performance, it is now possible to visualize extremely large maps in real time on a laptop or a tablet. HEALPix Viewer described here is developed for macOS, and takes full advantage of GPU acceleration to handle extremely large datasets in real time. It compiles natively on Intel and Arm64 architectures, and uses Metal framework for high-performance GPU computations. The aim of this project is to reduce the effort required for interactive data exploration, as well as time overhead for producing publication-quality maps. Drag and drop integration with Keynote and Powerpoint makes creating presentations easy.
The main codebase is written in Swift, a modern and efficient compiled language, with high-performance computing parts delegated entirely to GPU, and a few inserts in C interfacing to cfitsio library for I/O. Graphical user interface is written in SwiftUI, a new declarative UI framework based on Swift. Most common spherical projections and colormaps are supported out of the box, and the available source code makes it easy to customize the application and to add new features if desired. Overall performance improves on the current Python software stack by a factor of 3-180x depending on the task at hand.
△ Less
Submitted 23 September, 2023; v1 submitted 19 May, 2023;
originally announced May 2023.
-
The Simons Observatory: Galactic Science Goals and Forecasts
Authors:
Brandon S. Hensley,
Susan E. Clark,
Valentina Fanfani,
Nicoletta Krachmalnicoff,
Giulio Fabbian,
Davide Poletti,
Giuseppe Puglisi,
Gabriele Coppi,
Jacob Nibauer,
Roman Gerasimov,
Nicholas Galitzki,
Steve K. Choi,
Peter C. Ashton,
Carlo Baccigalupi,
Eric Baxter,
Blakesley Burkhart,
Erminia Calabrese,
Jens Chluba,
Josquin Errard,
Andrei V. Frolov,
Carlos Hervías-Caimapo,
Kevin M. Huffenberger,
Bradley R. Johnson,
Baptiste Jost,
Brian Keating
, et al. (9 additional authors not shown)
Abstract:
Observing in six frequency bands from 27 to 280 GHz over a large sky area, the Simons Observatory (SO) is poised to address many questions in Galactic astrophysics in addition to its principal cosmological goals. In this work, we provide quantitative forecasts on astrophysical parameters of interest for a range of Galactic science cases. We find that SO can: constrain the frequency spectrum of pol…
▽ More
Observing in six frequency bands from 27 to 280 GHz over a large sky area, the Simons Observatory (SO) is poised to address many questions in Galactic astrophysics in addition to its principal cosmological goals. In this work, we provide quantitative forecasts on astrophysical parameters of interest for a range of Galactic science cases. We find that SO can: constrain the frequency spectrum of polarized dust emission at a level of $Δβ_d \lesssim 0.01$ and thus test models of dust composition that predict that $β_d$ in polarization differs from that measured in total intensity; measure the correlation coefficient between polarized dust and synchrotron emission with a factor of two greater precision than current constraints; exclude the non-existence of exo-Oort clouds at roughly 2.9$σ$ if the true fraction is similar to the detection rate of giant planets; map more than 850 molecular clouds with at least 50 independent polarization measurements at 1 pc resolution; detect or place upper limits on the polarization fractions of CO(2-1) emission and anomalous microwave emission at the 0.1% level in select regions; and measure the correlation coefficient between optical starlight polarization and microwave polarized dust emission in $1^\circ$ patches for all lines of sight with $N_{\rm H} \gtrsim 2\times10^{20}$ cm$^{-2}$. The goals and forecasts outlined here provide a roadmap for other microwave polarization experiments to expand their scientific scope via Milky Way astrophysics.
△ Less
Submitted 3 November, 2021;
originally announced November 2021.
-
Observational Constraints on Constant Roll Inflation
Authors:
Jose T. Galvez Ghersi,
Alex Zucca,
Andrei V. Frolov
Abstract:
Constant-roll inflation was recently introduced by Motohashi, Starobinsky and Yokoyama as a phenomenological way to parametrize deviations from the slow-roll scenarios. In this paper, we investigate the dynamics of both the background and the perturbations in this model, without making any slow-roll assumptions. The perturbation spectra are computed with an efficient and accurate novel method that…
▽ More
Constant-roll inflation was recently introduced by Motohashi, Starobinsky and Yokoyama as a phenomenological way to parametrize deviations from the slow-roll scenarios. In this paper, we investigate the dynamics of both the background and the perturbations in this model, without making any slow-roll assumptions. The perturbation spectra are computed with an efficient and accurate novel method that allowed us to quickly scan the parameter space of constant-roll inflation. We derive the constraints on the model parameters from the cosmic microwave background anisotropy measurements provided by the joint analysis of the Planck Collaboration and the BICEP2/Keck Array data.
△ Less
Submitted 8 May, 2019; v1 submitted 3 August, 2018;
originally announced August 2018.
-
Unbraiding the Bounce: Superluminality around the Corner
Authors:
David A. Dobre,
Andrei V. Frolov,
José T. Gálvez Ghersi,
Sabir Ramazanov,
Alexander Vikman
Abstract:
We study a particular realization of the cosmological bounce scenario proposed recently by Ijjas and Steinhardt. First, we find that their bouncing solution starts from a divergent sound speed and ends with its vanishing. Thus, the solution connects two strongly coupled configurations. These pathologies are separated from the bouncing regime by only a few Planck times. We then reveal the exact str…
▽ More
We study a particular realization of the cosmological bounce scenario proposed recently by Ijjas and Steinhardt. First, we find that their bouncing solution starts from a divergent sound speed and ends with its vanishing. Thus, the solution connects two strongly coupled configurations. These pathologies are separated from the bouncing regime by only a few Planck times. We then reveal the exact structure of the Lagrangian, which reproduces this bouncing solution. This reconstruction allowed us to consider other cosmological solutions of the theory and analyze the phase space. In particular, we find other bouncing solutions and solutions with superluminal sound speed. These stable superluminal states can be continuously transformed into the solution constructed by Ijjas and Steinhardt. We discuss the consequences of this feature for a possible UV-completion.
△ Less
Submitted 13 March, 2018; v1 submitted 29 December, 2017;
originally announced December 2017.
-
Two-point correlators revisited: Fast and slow scales in multifield models of inflation
Authors:
José T. Gálvez Ghersi,
Andrei V. Frolov
Abstract:
We study the structure of two-point correlators of the inflationary field fluctuations in order to improve the accuracy and efficiency of the existing methods to calculate primordial spectra. We present a description motivated by the separation of the fast and slow evolving components of the spectrum which is based on Cholesky decomposition of the field correlator matrix. Our purpose is to rewrite…
▽ More
We study the structure of two-point correlators of the inflationary field fluctuations in order to improve the accuracy and efficiency of the existing methods to calculate primordial spectra. We present a description motivated by the separation of the fast and slow evolving components of the spectrum which is based on Cholesky decomposition of the field correlator matrix. Our purpose is to rewrite all the relevant equations of motion in terms of slowly varying quantities. This is important in order to consider the contribution from high-frequency modes to the spectrum without affecting computational performance. The slow-roll approximation is not required to reproduce the main distinctive features in the power spectrum for each specific model of inflation.
△ Less
Submitted 7 April, 2017; v1 submitted 15 September, 2016;
originally announced September 2016.
-
Cosmological dynamics in f(R) gravity
Authors:
Jun-Qi Guo,
Andrei V. Frolov
Abstract:
In this paper, we study the cosmological viability conditions, the phase-space dynamics, and the cosmological evolution of f(R) gravity. In contrast to most previous works in the literature, which analyzed the background dynamics of f(R) gravity by means of a dynamical system, we proceed by focusing on the equivalent scalar field description of the theory, which we believe is a more intuitive way…
▽ More
In this paper, we study the cosmological viability conditions, the phase-space dynamics, and the cosmological evolution of f(R) gravity. In contrast to most previous works in the literature, which analyzed the background dynamics of f(R) gravity by means of a dynamical system, we proceed by focusing on the equivalent scalar field description of the theory, which we believe is a more intuitive way of treating the problem. In order to study how the physical solutions evolve in f(R) cosmology, we explore the cosmological dynamics of a range of f(R) models, including models that yield a large hierarchy of scales and are singularity free. We present generic features of the phase-space dynamics in f(R) cosmology. We study the global structure of the phase space in f(R) gravity by compactifying the infinite phase space into a finite space via the Poincare transformation. On the expansion branch of the phase space, the constraint surface has a repeller and a de Sitter attractor; while on the contraction branch, the constraint surface has an attractor and a de Sitter repeller. Generally, the phase currents originate from the repeller and terminate at the corresponding attractor in each space. The trajectories between the repeller and the attractor in the presence of matter density are different from those in the vacuum case. The phase analysis techniques developed in this paper are very general, and can be applied to other similar dynamical systems.
△ Less
Submitted 9 February, 2014; v1 submitted 30 May, 2013;
originally announced May 2013.
-
Small Cosmological Constant from Running Gravitational Coupling
Authors:
Andrei V. Frolov,
Jun-Qi Guo
Abstract:
In this paper, we explore an idea of having Newton's constant change its value depending on the curvature scale involved. Such modification leads to a particular scalar-tensor gravity theory, with the Lagrangian derived from renormalization group (RG) flow arguments. Several of the well-known f(R) modified gravity models have remarkably simple description in terms of the infrared renormalization g…
▽ More
In this paper, we explore an idea of having Newton's constant change its value depending on the curvature scale involved. Such modification leads to a particular scalar-tensor gravity theory, with the Lagrangian derived from renormalization group (RG) flow arguments. Several of the well-known f(R) modified gravity models have remarkably simple description in terms of the infrared renormalization group, but not the "designer" types in general. We find that de Sitter-like accelerated expansion can be generated even in the absence of cosmological constant term, entirely due to running of the Newton's constant. In hopes of tackling the problem of cosmological constant's smallness, we explore the flows which are capable of generating exponential hierarchy between infrared and ultraviolet scales, and investigate cosmological evolution in the models thus derived.
△ Less
Submitted 4 September, 2012; v1 submitted 26 January, 2011;
originally announced January 2011.
-
Non-linear Dynamics and Primordial Curvature Perturbations from Preheating
Authors:
Andrei V. Frolov
Abstract:
In this paper I review the theory and numerical simulations of non-linear dynamics of preheating, a stage of dynamical instability at the end of inflation during which homogeneous inflaton explosively decays and deposits its energy into excitation of other matter fields. I focus on preheating in chaotic inflation models, which proceeds via broad parametric resonance. I describe a simple method to…
▽ More
In this paper I review the theory and numerical simulations of non-linear dynamics of preheating, a stage of dynamical instability at the end of inflation during which homogeneous inflaton explosively decays and deposits its energy into excitation of other matter fields. I focus on preheating in chaotic inflation models, which proceeds via broad parametric resonance. I describe a simple method to evaluate Floquet exponents, calculating stability diagrams of Mathieu and Lame equations describing development of instability in $m^2φ^2$ and $λφ^4$ preheating models. I discuss basic numerical methods and issues, and present simulation results highlighting non-equilibrium transitions, topological defect formation, late-time universality, turbulent scaling and approach to thermalization. I explain how preheating can generate large-scale primordial (non-Gaussian) curvature fluctuations manifest in cosmic microwave background anisotropy and large scale structure, and discuss potentially observable signatures of preheating.
△ Less
Submitted 13 May, 2010; v1 submitted 20 April, 2010;
originally announced April 2010.
-
Non-Gaussian Spikes from Chaotic Billiards in Inflation Preheating
Authors:
J. Richard Bond,
Andrei V. Frolov,
Zhiqi Huang,
Lev Kofman
Abstract:
A new class of non-Gaussian curvature fluctuations ζ_{pr} (\bx) \equiv δN(χ_i) arises from the post-inflation preheating behaviour of a non-inflaton field χ_i. Its billiard-like chaotic dynamics imprints regular log-spaced narrow spikes in the number of preheating e-folds N(χ_i). We perform highly accurate lattice simulations of SUSY-inspired quartic inflaton and coupling potentials in a separat…
▽ More
A new class of non-Gaussian curvature fluctuations ζ_{pr} (\bx) \equiv δN(χ_i) arises from the post-inflation preheating behaviour of a non-inflaton field χ_i. Its billiard-like chaotic dynamics imprints regular log-spaced narrow spikes in the number of preheating e-folds N(χ_i). We perform highly accurate lattice simulations of SUSY-inspired quartic inflaton and coupling potentials in a separate-universe approximation to compute N(χ_i) as a function of the (nearly homogeneous) initial condition χ_i. The super-horizon modes of χ_i(\bx) result in positive spiky excursions in ζ_{pr} and hence negative gravitational potential fluctuations added to the usual sign-independent inflaton-induced perturbations, observably manifested in large cosmic structures and as (polarized) temperature CMB cold spots.
△ Less
Submitted 20 March, 2009;
originally announced March 2009.
-
Non-Gaussianity as a Probe of the Physics of the Primordial Universe and the Astrophysics of the Low Redshift Universe
Authors:
E. Komatsu,
N. Afshordi,
N. Bartolo,
D. Baumann,
J. R. Bond,
E. I. Buchbinder,
C. T. Byrnes,
X. Chen,
D. J. H. Chung,
A. Cooray,
P. Creminelli,
N. Dalal,
O. Dore,
R. Easther,
A. V. Frolov,
K. M. Gorski,
M. G. Jackson,
J. Khoury,
W. H. Kinney,
L. Kofman,
K. Koyama,
L. Leblond,
J. -L. Lehners,
J. E. Lidsey,
M. Liguori
, et al. (36 additional authors not shown)
Abstract:
A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in f…
▽ More
A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.
△ Less
Submitted 19 April, 2009; v1 submitted 27 February, 2009;
originally announced February 2009.
-
DEFROST: A New Code for Simulating Preheating after Inflation
Authors:
Andrei V. Frolov
Abstract:
At the end of inflation, dynamical instability can rapidly deposit the energy of homogeneous cold inflaton into excitations of other fields. This process, known as preheating, is rather violent, inhomogeneous and non-linear, and has to be studied numerically. This paper presents a new code for simulating scalar field dynamics in expanding universe written for that purpose. Compared to available…
▽ More
At the end of inflation, dynamical instability can rapidly deposit the energy of homogeneous cold inflaton into excitations of other fields. This process, known as preheating, is rather violent, inhomogeneous and non-linear, and has to be studied numerically. This paper presents a new code for simulating scalar field dynamics in expanding universe written for that purpose. Compared to available alternatives, it significantly improves both the speed and the accuracy of calculations, and is fully instrumented for 3D visualization. We reproduce previously published results on preheating in simple chaotic inflation models, and further investigate non-linear dynamics of the inflaton decay. Surprisingly, we find that the fields do not want to thermalize quite the way one would think. Instead of directly reaching equilibrium, the evolution appears to be stuck in a rather simple but quite inhomogeneous state. In particular, one-point distribution function of total energy density appears to be universal among various two-field preheating models, and is exceedingly well described by a lognormal distribution. It is tempting to attribute this state to scalar field turbulence.
△ Less
Submitted 28 October, 2008; v1 submitted 29 September, 2008;
originally announced September 2008.
-
A Singularity Problem with f(R) Dark Energy
Authors:
Andrei V. Frolov
Abstract:
In this paper, I point out that there is a curvature singularity problem appearing on non-linear level that generally plagues f(R) models that modify Einstein gravity in the infrared. It is caused by the fact that for the effective scalar degree of freedom, the curvature singularity is at a finite field value and energy level, and can be easily accessed by the field dynamics in the presence of m…
▽ More
In this paper, I point out that there is a curvature singularity problem appearing on non-linear level that generally plagues f(R) models that modify Einstein gravity in the infrared. It is caused by the fact that for the effective scalar degree of freedom, the curvature singularity is at a finite field value and energy level, and can be easily accessed by the field dynamics in the presence of matter. This problem is invisible in linearized analysis, except for the tell-tale growing oscillatory modes it causes. In view of this, viability of many f(R) models in current literature will have to be re-evaluated.
△ Less
Submitted 21 June, 2008; v1 submitted 17 March, 2008;
originally announced March 2008.
-
Persistent Patterns in Accretion Disks
Authors:
Mustafa A. Amin,
Andrei V. Frolov
Abstract:
We present a set of new characteristic frequencies associated with accretion disks around compact objects. These frequencies arise from persistent rotating patterns in the disk that are finite in radial extent and driven purely by the gravity of the central body. Their existence depends on general relativistic corrections to orbital motion and, if observed, could be used to probe the strong grav…
▽ More
We present a set of new characteristic frequencies associated with accretion disks around compact objects. These frequencies arise from persistent rotating patterns in the disk that are finite in radial extent and driven purely by the gravity of the central body. Their existence depends on general relativistic corrections to orbital motion and, if observed, could be used to probe the strong gravity region around a black hole. We also discuss a possible connection to the puzzle of quasi-periodic oscillations.
△ Less
Submitted 24 March, 2006;
originally announced March 2006.
-
Accretion of Ghost Condensate by Black Holes
Authors:
Andrei V. Frolov
Abstract:
The intent of this letter is to point out that the accretion of a ghost condensate by black holes could be extremely efficient. We analyze steady-state spherically symmetric flows of the ghost fluid in the gravitational field of a Schwarzschild black hole and calculate the accretion rate. Unlike minimally coupled scalar field or quintessence, the accretion rate is set not by the cosmological ene…
▽ More
The intent of this letter is to point out that the accretion of a ghost condensate by black holes could be extremely efficient. We analyze steady-state spherically symmetric flows of the ghost fluid in the gravitational field of a Schwarzschild black hole and calculate the accretion rate. Unlike minimally coupled scalar field or quintessence, the accretion rate is set not by the cosmological energy density of the field, but by the energy scale of the ghost condensate theory. If hydrodynamical flow is established, it could be as high as tenth of a solar mass per second for 10MeV-scale ghost condensate accreting onto a stellar-sized black hole, which puts serious constraints on the parameters of the ghost condensate model.
△ Less
Submitted 29 July, 2004; v1 submitted 28 April, 2004;
originally announced April 2004.
-
BRANECODE: A Program for Simulations of Braneworld Dynamics
Authors:
Johannes Martin,
Gary N. Felder,
Andrei V. Frolov,
Lev Kofman,
Marco Peloso
Abstract:
We describe an algorithm and a C++ implementation that we have written and made available for calculating the fully nonlinear evolution of 5D braneworld models with scalar fields. Bulk fields allow for the stabilization of the extra space. However, they complicate the dynamics of the system, so that analytic calculations (performed within an effective 4D theory) are typically only reliable close…
▽ More
We describe an algorithm and a C++ implementation that we have written and made available for calculating the fully nonlinear evolution of 5D braneworld models with scalar fields. Bulk fields allow for the stabilization of the extra space. However, they complicate the dynamics of the system, so that analytic calculations (performed within an effective 4D theory) are typically only reliable close to stabilized configurations or when the evolution of the extra space is negligible. In the general case, a numerical study of the 5D equations is necessary, and the algorithm and code we describe are the first ones designed for this task. The program and its full documentation are available on the Web at http://www.cita.utoronto.ca/~jmartin/BRANECODE/. In this paper we provide a brief overview of what the program does and how to use it.
△ Less
Submitted 19 April, 2004;
originally announced April 2004.
-
Can Inflating Braneworlds be Stabilized?
Authors:
Andrei V. Frolov,
Lev Kofman
Abstract:
We investigate scalar perturbations from inflation in braneworld cosmologies with extra dimensions. For this we calculate scalar metric fluctuations around five dimensional warped geometry with four dimensional de Sitter slices. The background metric is determined self-consistently by the (arbitrary) bulk scalar field potential, supplemented by the boundary conditions at both orbifold branes. As…
▽ More
We investigate scalar perturbations from inflation in braneworld cosmologies with extra dimensions. For this we calculate scalar metric fluctuations around five dimensional warped geometry with four dimensional de Sitter slices. The background metric is determined self-consistently by the (arbitrary) bulk scalar field potential, supplemented by the boundary conditions at both orbifold branes. Assuming that the inflating branes are stabilized (by the brane scalar field potentials), we estimate the lowest eigenvalue of the scalar fluctuations - the radion mass. In the limit of flat branes, we reproduce well known estimates of the positive radion mass for stabilized branes. Surprisingly, however, we found that for de Sitter (inflating) branes the square of the radion mass is typically negative, which leads to a strong tachyonic instability. Thus, parameters of stabilized inflating braneworlds must be constrained to avoid this tachyonic instability. Instability of "stabilized" de Sitter branes is confirmed by the BraneCode numerical calculations in the accompanying paper hep-th/0309001. If the model's parameters are such that the radion mass is smaller than the Hubble parameter, we encounter a new mechanism of generation of primordial scalar fluctuations, which have a scale free spectrum and acceptable amplitude.
△ Less
Submitted 1 September, 2003;
originally announced September 2003.
-
Braneworld dynamics with the BraneCode
Authors:
Johannes Martin,
Gary N. Felder,
Andrei V. Frolov,
Marco Peloso,
Lev Kofman
Abstract:
We give a full nonlinear numerical treatment of time-dependent 5d braneworld geometry, which is determined self-consistently by potentials for the scalar field in the bulk and at two orbifold branes, supplemented by boundary conditions at the branes. We describe the BraneCode, an algorithm which we designed to solve the dynamical equations numerically. We applied the BraneCode to braneworld mode…
▽ More
We give a full nonlinear numerical treatment of time-dependent 5d braneworld geometry, which is determined self-consistently by potentials for the scalar field in the bulk and at two orbifold branes, supplemented by boundary conditions at the branes. We describe the BraneCode, an algorithm which we designed to solve the dynamical equations numerically. We applied the BraneCode to braneworld models and found several novel phenomena of the brane dynamics. Starting with static warped geometry with de Sitter branes, we found numerically that this configuration is often unstable due to a tachyonic mass of the radion during inflation. If the model admits other static configurations with lower values of de Sitter curvature, this effect causes a violent re-structuring towards them, flattening the branes, which appears as a lowering of the 4d effective cosmological constant. Braneworld dynamics can often lead to brane collisions. We found that in the presence of the bulk scalar field, the 5d geometry between colliding branes approaches a universal, homogeneous, anisotropic strong gravity Kasner-like asymptotic, irrespective of the bulk/brane potentials. The Kasner indices of the brane directions are equal to each other but different from that of the extra dimension.
△ Less
Submitted 1 September, 2003;
originally announced September 2003.
-
The naked singularity in the global structure of critical collapse spacetimes
Authors:
Andrei V. Frolov,
Ue-Li Pen
Abstract:
We examine the global structure of scalar field critical collapse spacetimes using a characteristic double-null code. It can integrate past the horizon without any coordinate problems, due to the careful choice of constraint equations used in the evolution. The limiting sequence of sub- and supercritical spacetimes presents an apparent paradox in the expected Penrose diagrams, which we address i…
▽ More
We examine the global structure of scalar field critical collapse spacetimes using a characteristic double-null code. It can integrate past the horizon without any coordinate problems, due to the careful choice of constraint equations used in the evolution. The limiting sequence of sub- and supercritical spacetimes presents an apparent paradox in the expected Penrose diagrams, which we address in this paper. We argue that the limiting spacetime converges pointwise to a unique limit for all r>0, but not uniformly. The r=0 line is different in the two limits. We interpret that the two different Penrose diagrams differ by a discontinuous gauge transformation. We conclude that the limiting spacetime possesses a singular event, with a future removable naked singularity.
△ Less
Submitted 17 July, 2003;
originally announced July 2003.