Evolution of radial profiles in regular Lemaitre-Tolman-Bondi dust models

We undertake a comprehensive and rigorous analytic study of the evolution of radial profiles of covariant scalars in regular Lemaitre-Tolman-Bondi dust models. We consider specifically the phenomenon of "profile inversions" in which an initial clump profile of density, spatial curvature or the expansion scalar, might evolve into a void profile (and vice versa). Previous work in the literature on models with density void profiles and/or allowing for density profile inversions is given full generalization, with some erroneous results corrected. We prove rigorously that if an evolution without shell crossings is assumed, then only the 'clump to void' inversion can occur in density profiles, and only in hyperbolic models or regions with negative spatial curvature. The profiles of spatial curvature follow similar patterns as those of the density, with 'clump to void' inversions only possible for hyperbolic models or regions. However, profiles of the expansion scalar are less restrictive, with profile inversions necessarily taking place in elliptic models. We also examine radial profiles in special LTB configurations: closed elliptic models, models with a simultaneous big bang singularity, as well as a locally collapsing elliptic region surrounded by an expanding hyperbolic background. The general analytic statements that we obtain allow for setting up the right initial conditions to construct fully regular LTB models with any specific qualitative requirements for the profiles of all scalars and their time evolution. The results presented can be very useful in guiding future numerical work on these models and in revising previous analytic work on all their applications.Comment: Final version to appear in Classical and Quantum Gravity. Readers eager to know the results and implications without having to go through the technical detail are recommended to go directly to the summary and discussion in the final section (section 11). Typos have been corrected and an important reference has been adde

Towards a physical interpretation for the Stephani Universes

A physicaly reasonable interpretation is provided for the perfect fluid, sphericaly symmetric, conformally flat ``Stephani Universes''. The free parameters of this class of exact solutions are determined so that the ideal gas relation $p=n k T$ is identicaly fulfiled, while the full equation of state of a classical monatomic ideal gas and a matter-radiation mixture holds up to a good approximation in a near dust, matter dominated regime. Only the models having spacelike slices with positive curvature admit a regular evolution domain that avoids an unphysical singularity. In the matter dominated regime these models are dynamicaly and observationaly indistinguishable from ``standard'' FLRW cosmology with a dust source.Comment: 17 pages, 2 figures, LaTeX with revtex style, submitted to General Relativity and Gravitatio

Inhomogeneous models of interacting dark matter and dark energy

We derive and analyze a class of spherically symmetric cosmological models whose source is an interactive mixture of inhomogeneous cold dark matter (DM) and a generic homogeneous dark energy (DE) fluid. If the DE fluid corresponds to a quintessense scalar field, the interaction term can be associated with a well motivated non--minimal coupling to the DM component. By constructing a suitable volume average of the DM component we obtain a Friedman evolution equation relating this average density with an average Hubble scalar, with the DE component playing the role of a repulsive and time-dependent $\Lambda$ term. Once we select an ``equation of state'' linking the energy density ($\mu$) and pressure ($p$) of the DE fluid, as well as a free function governing the radial dependence, the models become fully determinate and can be applied to known specific DE sources, such as quintessense scalar fields or tachyonic fluids. Considering the simple equation of state $p= (\gamma-1) \mu$ with $0\leq\gamma <2/3$, we show that the free parameters and boundary conditions can be selected for an adequate description of a local DM overdensity evolving in a suitable cosmic background that accurately fits current observational data. While a DE dominated scenario emerges in the asymptotic future, with total $\Omega$ and $q$ tending respectively to 1 and -1/2 for all cosmic observers, the effects of inhomogeneity and anisotropy yield different local behavior and evolution rates for these parameters in the local overdense region. We suggest that the models presented can be directly applied to explore the effects of various DE formalisms on local DM cosmological inhomogeneities.Comment: 15 pages, revtex4, 10 eps figure

On the Thermodynamics of Simple Non-Isentropic Perfect Fluids in General Relativity

We examine the consistency of the thermodynamics of irrotational and non-isentropic perfect fluids complying with matter conservation by looking at the integrability conditions of the Gibbs-Duhem relation. We show that the latter is always integrable for fluids of the following types: (a) static, (b) isentropic (admits a barotropic equation of state), (c) the source of a spacetime for which $r\ge 2$, where $r$ is the dimension of the orbit of the isometry group. This consistency scheme is tested also in two large classes of known exact solutions for which $r< 2$, in general: perfect fluid Szekeres solutions (classes I and II). In none of these cases, the Gibbs-Duhem relation is integrable, in general, though specific particular cases of Szekeres class II (all complying with $r<2$) are identified for which the integrability of this relation can be achieved. We show that Szekeres class I solutions satisfy the integrability conditions only in two trivial cases, namely the spherically symmetric limiting case and the Friedman-Roberson-Walker (FRW) cosmology. Explicit forms of the state variables and equations of state linking them are given explicitly and discussed in relation to the FRW limits of the solutions. We show that fixing free parameters in these solutions by a formal identification with FRW parameters leads, in all cases examined, to unphysical temperature evolution laws, quite unrelated to those of their FRW limiting cosmologies.Comment: 29 pages, Plain.Te

Weighed scalar averaging in LTB dust models, part I: statistical fluctuations and gravitational entropy

We introduce a weighed scalar average formalism ("q-average") for the study of the theoretical properties and the dynamics of spherically symmetric Lemaitre-Tolman-Bondi (LTB) dust models models. The "q-scalars" that emerge by applying the q-averages to the density, Hubble expansion and spatial curvature (which are common to FLRW models) are directly expressible in terms of curvature and kinematic invariants and identically satisfy FLRW evolution laws without the back-reaction terms that characterize Buchert's average. The local and non-local fluctuations and perturbations with respect to the q-average convey the effects of inhomogeneity through the ratio of curvature and kinematic invariants and the magnitude of radial gradients. All curvature and kinematic proper tensors that characterize the models are expressible as irreducible algebraic expansions on the metric and 4-velocity, whose coefficients are the q-scalars and their linear and quadratic local fluctuations. All invariant contractions of these tensors are quadratic fluctuations, whose q-averages are directly and exactly related to statistical correlation moments of the density and Hubble expansion scalar. We explore the application of this formalism to a definition of a gravitational entropy functional proposed by Hosoya et al (2004 Phys. Rev. Lett. 92 141302). We show that a positive entropy production follows from a negative correlation between fluctuations of the density and Hubble scalar, providing a brief outline on its fulfillment in various LTB models and regions. While the q-average formalism is specially suited for LTB and Szekeres models, it may provide a valuable theoretical insight on the properties of scalar averaging in inhomogeneous spacetimes in general.Comment: 27 pages in IOP format, 1 figure. Matches version accepted for publication in Classical and Quantum Gravit

Lattice Boltzmann Model for Axisymmetric Multiphase Flows

In this paper, a lattice Boltzmann (LB) model is presented for axisymmetric multiphase flows. Source terms are added to a two-dimensional standard lattice Boltzmann equation (LBE) for multiphase flows such that the emergent dynamics can be transformed into the axisymmetric cylindrical coordinate system. The source terms are temporally and spatially dependent and represent the axisymmetric contribution of the order parameter of fluid phases and inertial, viscous and surface tension forces. A model which is effectively explicit and second order is obtained. This is achieved by taking into account the discrete lattice effects in the Chapman-Enskog multiscale analysis, so that the macroscopic axisymmetric mass and momentum equations for multiphase flows are recovered self-consistently. The model is extended to incorporate reduced compressibility effects. Axisymmetric equilibrium drop formation and oscillations, breakup and formation of satellite droplets from viscous liquid cylindrical jets through Rayleigh capillary instability and drop collisions are presented. Comparisons of the computed results with available data show satisfactory agreement.Comment: 17 pages, 11 figures, to be published in Physical Review

Specific IgE Response to Purified and Recombinant Allergens in Latex Allergy

Background In recent years, allergy to natural rubber latex has emerged as a major allergy among certain occupational groups and patients with underlying diseases. The sensitization and development of latex allergy has been attributed to exposure to products containing residual latex proteins. Although improved manufacturing procedures resulted in a considerable reduction of new cases, the potential risk for some patient groups is still great. In addition the prevalent cross-reactivity of latex proteins with other food allergens poses a major concern. A number of purified allergens and a few commercial kits are currently available, but no concerted effort was undertaken to evaluate them. Methods We studied 11 purified latex allergens, Hev b 1 to Hev b 10, and Hev b 13 along with several crude allergen extracts and two commercial ImmunoCAP assays to evaluate specific IgE antibody in the sera from latex allergic patients and controls. Health care workers and spina bifida patients with clinical symptoms of latex allergy, spina bifida patients without latex allergy, and non-atopic health care workers have been studied. Results The results suggest that Hev b 2, 5, 6, and 13 together identified over 80 percent health care workers with latex allergy, while Hev b 6 along with Hev b 1 or 3 detected specific IgE antibody in all sera studied from patients with spina bifida and latex allergy. The ImmunoCAP results using both Hev b 5 amplified and non-amplified closely agreed with the clinical diagnosis of latex allergy in health care workers and in spina bifida. Conclusion Although the purified allergens and crude extracts reacted diversely with IgE from different patient groups, the results indicated that use of certain combinations of purified recombinant antigens will be useful in commercial kits or in in-house assays for detecting specific IgE antibody in the sera. The results suggest that a combination of Hev b 2, 3, 5, 6, and 13 together detected specific IgE in 80% of the sera from latex allergic patients. Both ImmunoCAPs correctly identified over 95% of latex allergic patients, however, showed reactivity with a few normal control subject

Back-reaction and effective acceleration in generic LTB dust models

We provide a thorough examination of the conditions for the existence of back-reaction and an "effective" acceleration (in the context of Buchert's averaging formalism) in regular generic spherically symmetric Lemaitre-Tolman-Bondi (LTB) dust models. By considering arbitrary spherical comoving domains, we verify rigorously the fulfillment of these conditions expressed in terms of suitable scalar variables that are evaluated at the boundary of every domain. Effective deceleration necessarily occurs in all domains in: (a) the asymptotic radial range of models converging to a FLRW background, (b) the asymptotic time range of non-vacuum hyperbolic models, (c) LTB self-similar solutions and (d) near a simultaneous big bang. Accelerating domains are proven to exist in the following scenarios: (i) central vacuum regions, (ii) central (non-vacuum) density voids, (iii) the intermediate radial range of models converging to a FLRW background, (iv) the asymptotic radial range of models converging to a Minkowski vacuum and (v) domains near and/or intersecting a non-simultaneous big bang. All these scenarios occur in hyperbolic models with negative averaged and local spatial curvature, though scenarios (iv) and (v) are also possible in low density regions of a class of elliptic models in which local spatial curvature is negative but its average is positive. Rough numerical estimates between -0.003 and -0.5 were found for the effective deceleration parameter. While the existence of accelerating domains cannot be ruled out in models converging to an Einstein de Sitter background and in domains undergoing gravitational collapse, the conditions for this are very restrictive. The results obtained may provide important theoretical clues on the effects of back-reaction and averaging in more general non-spherical models.Comment: Final version accepted for publication in Classical and Quantum Gravity. 47 pages in IOP LaTeX macros, 12 pdf figure