Primordial Galaxy Formation and IGM Reionization

In this talk I will present a model for primordial galaxy formation. In particular, I will review the feedback effects that regulate the process: (i) radiative (i.e. ionizing and H_2-photodissociating photons) and (ii) stellar (i.e. SN explosions) feedback produced by massive stars. I will show the results of a model for galaxy formation and IGM reionization, which includes a self-consistent treatment of the above feedback effects. Finally, I will describe a Monte Carlo method for the radiative transfer of ionizing photons through the IGM and discuss its application to the IGM reionization problem.Comment: 4 pages, 3 figures. To appear in "A New Era in Cosmology", (ASP Conference Proceedings), eds. T. Shanks and N. Metcalf

Reionization History from Coupled CMB/21cm Line Data

We study CMB secondary anisotropies produced by inhomogeneous reionization by means of cosmological simulations coupled with the radiative transfer code CRASH. The reionization history is consistent with the WMAP Thomson optical depth determination. We find that the signal arising from this process dominates over the primary CMB component for l > 4000 and reaches a maximum amplitude of l(l+1)C_l/2\pi ~ 1.6 x 10^{-13} on arcmin scale, i.e. l as large as several thousands. We then cross-correlate secondary CMB anisotropy maps with neutral hydrogen 21cm line emission fluctuations obtained from the same simulations. The two signals are highly anti-correlated on angular scales corresponding to the typical size of HII regions (including overlapping) at the 21cm map redshift. We show how the CMB/21cm cross-correlation can be used to: (a) study the nature of the reionization sources, (b) reconstruct the cosmic reionization history, (c) infer the mean cosmic ionization level at any redshift. We discuss the feasibility of the proposed experiment with forthcoming facilities.Comment: 7 pages, 5 figures, MNRAS in pres

The Kepler-19 System: A Transiting 2.2 R_⊕ Planet and a Second Planet Detected via Transit Timing Variations

We present the discovery of the Kepler-19 planetary system, which we first identified from a 9.3 day periodic transit signal in the Kepler photometry. From high-resolution spectroscopy of the star, we find a stellar effective temperature T_(eff) = 5541 ± 60 K, a metallicity [Fe/H] = –0.13 ± 0.06, and a surface gravity log(g) = 4.59 ± 0.10. We combine the estimate of T_(eff) and [Fe/H] with an estimate of the stellar density derived from the photometric light curve to deduce a stellar mass of M_*= 0.936 ± 0.040 M_☉ and a stellar radius of R_* = 0.850 ± 0.018 R_☉ (these errors do not include uncertainties in the stellar models). We rule out the possibility that the transits result from an astrophysical false positive by first identifying the subset of stellar blends that reproduce the precise shape of the light curve. Using the additional constraints from the measured color of the system, the absence of a secondary source in the high-resolution spectrum, and the absence of a secondary source in the adaptive optics imaging, we conclude that the planetary scenario is more than three orders of magnitude more likely than a blend. The blend scenario is independently disfavored by the achromaticity of the transit: we measure a transit depth with Spitzer at 4.5 μm of 547^(+113)_(–110) ppm, consistent with the depth measured in the Kepler optical bandpass of 567 ± 6 ppm (corrected for stellar limb darkening). We determine a physical radius of the planet Kepler-19b of R_p = 2.209 ± 0.048 R_⊕; the uncertainty is dominated by uncertainty in the stellar parameters. From radial velocity observations of the star, we find an upper limit on the planet mass of 20.3 M_⊕, corresponding to a maximum density of 10.4 g cm^(–3). We report a significant sinusoidal deviation of the transit times from a predicted linear ephemeris, which we conclude is due to an additional perturbing body in the system. We cannot uniquely determine the orbital parameters of the perturber, as various dynamical mechanisms match the amplitude, period, and shape of the transit timing signal and satisfy the host star's radial velocity limits. However, the perturber in these mechanisms has a period ≾ 160 days and mass ≾ 6 M_(Jup), confirming its planetary nature as Kepler-19c. We place limits on the presence of transits of Kepler-19c in the available Kepler data

Modeling Kepler Transit Light Curves as False Positives: Rejection of Blend Scenarios for Kepler-9, and Validation of Kepler-9 d, A Super-earth-size Planet in a Multiple System

Light curves from the Kepler Mission contain valuable information on the nature of the phenomena producing the transit-like signals. To assist in exploring the possibility that they are due to an astrophysical false positive, we describe a procedure (BLENDER) to model the photometry in terms of a "blend" rather than a planet orbiting a star. A blend may consist of a background or foreground eclipsing binary (or star-planet pair) whose eclipses are attenuated by the light of the candidate and possibly other stars within the photometric aperture. We apply BLENDER to the case of Kepler-9 (KIC 3323887), a target harboring two previously confirmed Saturn-size planets (Kepler-9 b and Kepler-9 c) showing transit timing variations, and an additional shallower signal with a 1.59 day period suggesting the presence of a super-Earth-size planet. Using BLENDER together with constraints from other follow-up observations we are able to rule out all blends for the two deeper signals and provide independent validation of their planetary nature. For the shallower signal, we rule out a large fraction of the false positives that might mimic the transits. The false alarm rate for remaining blends depends in part (and inversely) on the unknown frequency of small-size planets. Based on several realistic estimates of this frequency, we conclude with very high confidence that this small signal is due to a super-Earth-size planet (Kepler-9 d) in a multiple system, rather than a false positive. The radius is determined to be 1.64^(+0.19)_(–0.14) R_⊕, and current spectroscopic observations are as yet insufficient to establish its mass

Social Norms, Cognitive Dissonance and the Timing of Marriage

We present a model of courtship in which the timing of marriage is affected by the cognitive dissonance between perceived norms and personal aims. We argue that as long as the family has been the main provider of social protection, marriage has been favoured by strongly felt social norms, and thus people accepted less-than-ideal partners early on in their search in order to minimize the dissonance caused by the non-adherence to the custom. Once the Welfare state has replaced the family, these norms have lost their strength, so that agents can afford the luxury of searching their preferred partners at length without feeling at odds with their social duties. The model yields predictions in line with relevant stylised facts: the raising age of marriage, the prevalence of assortative mating and the common occurrence of divorce in the early years of marriage. We finally discuss the impact of late marriages on fertility, and argue that there need not be negative consequences if the declining role of the family becomes socially accepted, and alternative arrangements are made possible and indeed encouraged by means of an appropriate family policy.marriage, cognitive dissonance, fertility

Lyalpha heating and its impact on early structure formation

In this paper we have calculated the effect of Lyalpha photons emitted by the first stars on the evolution of the IGM temperature. We have considered both a standard Salpeter IMF and a delta-function IMF for very massive stars with mass 300 M_sun. We find that the Lyalpha photons produced by the stellar populations considered here are able to heat the IGM at z<25, although never above ~100 K. Stars with a Salpeter IMF are more effective as, due to the contribution from small-mass long-living stars, they produce a higher Lyalpha background. Lyalpha heating can affect the subsequent formation of small mass objects by producing an entropy floor that may limit the amount of gas able to collapse and reduce the gas clumping.We find that the gas fraction in halos of mass below ~ 5 x 10^6 M_sun is less than 50% (for the smallest masses this fraction drops to 1% or less) compared to a case without Lyalpha heating. Finally, Lyalpha photons heat the IGM temperature above the CMB temperature and render the 21cm line from neutral hydrogen visible in emission at z<15.Comment: 7 pages, 5 figures, to be printed in MNRA

Radiative transfer of ionizing radiation through gas and dust: stellar source case

We present a new dust extension to the Monte Carlo radiative transfer code crash, which enables it to simulate the propagation of ionizing radiation through mixtures of gas and dust. The new code is applied to study the impact of dust absorption on idealized galactic H II regions and on small scale reionization. We find that H II regions are reduced in size by the presence of dust, while their inner temperature and ionization structure remain largely unaffected. In the small scale reionization simulation, dust hardens ionization fronts and delays the overlap of ionized bubbles. This effect is found to depend only weakly on the assumed abundance of dust in underdense regions.Comment: 17 pages, 14 figures. Accepted for publication in MNRA

AMADA-Analysis of Multidimensional Astronomical Datasets

We present AMADA, an interactive web application to analyse multidimensional datasets. The user uploads a simple ASCII file and AMADA performs a number of exploratory analysis together with contemporary visualizations diagnostics. The package performs a hierarchical clustering in the parameter space, and the user can choose among linear, monotonic or non-linear correlation analysis. AMADA provides a number of clustering visualization diagnostics such as heatmaps, dendrograms, chord diagrams, and graphs. In addition, AMADA has the option to run a standard or robust principal components analysis, displaying the results as polar bar plots. The code is written in R and the web interface was created using the Shiny framework. AMADA source-code is freely available at https://goo.gl/KeSPue, and the shiny-app at http://goo.gl/UTnU7I.Comment: Accepted for publication in Astronomy & Computin