Researcher profile

F. S. Kitaura

F. S. Kitaura contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2016arXiv

Search for Extremely Metal-poor Galaxies in the Sloan Digital Sky Survey (II): high electron temperature objects

Extremely metal-poor (XMP) galaxies are defined to have gas-phase metallicity smaller than a tenth of the solar value (12 + log[O/H] < 7.69). They are uncommon, chemically and possibly dynamically primitive, with physical conditions characteristic of earlier phases of the Universe. We search for new XMPs in the Sloan Digital Sky Survey (SDSS) in a work that complements Paper I. This time high electron temperature objects are selected; since metals are a main coolant of the gas, metal- poor objects contain high-temperature gas. Using the algorithm k-means, we classify 788677 spectra to select 1281 galaxies having particularly intense [OIII]4363 with respect to [OIII]5007, which is a proxy for high electron temperature. The metallicity of these candidates was computed using a hybrid technique consistent with the direct method, rendering 196 XMPs. A less restrictive noise constraint provides a larger set with 332 candidates. Both lists are provided in electronic format. The selected XMP sample have mean stellar mass around 10^8Msun, with dust-mass \sim 10^3Msun for typical star-forming regions. In agreement with previous findings, XMPs show a tendency to be tadpole-like or cometary. Their underlying stellar continuum corresponds to a fairly young stellar population (< 1Gyr), although young and aged stellar populations co-exists at the low-metallicity starbursts. About 10% of the XMPs shows large N/O. Based on their location in constrained cosmological numerical simulations, XMPs have a strong tendency to appear in voids and to avoid galaxy clusters. The puzzling 2%-solar low-metallicity threshold exhibited by XMPs remains.

preprint2016arXiv

The Low Redshift survey at Calar Alto (LoRCA)

The Baryon Acoustic Oscillation (BAO) feature in the power spectrum of galaxies provides a standard ruler to measure the accelerated expansion of the Universe. To extract all available information about dark energy, it is necessary to measure a standard ruler in the local, z<0.2, universe where dark energy dominates most the energy density of the Universe. Though the volume available in the local universe is limited, it is just big enough to measure accurately the long 100 Mpc/h wave-mode of the BAO. Using cosmological N-body simulations and approximate methods based on Lagrangian perturbation theory, we construct a suite of a thousand light-cones to evaluate the precision at which one can measure the BAO standard ruler in the local universe. We find that using the most massive galaxies on the full sky (34,000 sq. deg.), i.e. a K(2MASS)<14 magnitude-limited sample, one can measure the BAO scale up to a precision of 4\% and 1.2\% using reconstruction). We also find that such a survey would help to detect the dynamics of dark energy.Therefore, we propose a 3-year long observational project, named the Low Redshift survey at Calar Alto (LoRCA), to observe spectroscopically about 200,000 galaxies in the northern sky to contribute to the construction of aforementioned galaxy sample. The suite of light-cones is made available to the public.

preprint2013arXiv

Detection of the cosmic γ-ray horizon from multiwavelength observations of blazars

The first statistically significant detection of the cosmic γ-ray horizon (CGRH) that is independent of any extragalactic background light (EBL) model is presented. The CGRH is a fundamental quantity in cosmology. It gives an estimate of the opacity of the Universe to very high energy (VHE) γ-ray photons due to photon-photon pair production with the EBL. The only estimations of the CGRH to date are predictions from EBL models and lower limits from γ-ray observations of cosmological blazars and γ-ray bursts. Here, we present homogeneous synchrotron/synchrotron self-Compton (SSC) models of the spectral energy distributions of 15 blazars based on (almost) simultaneous observations from radio up to the highest energy γ-rays taken with the Fermi satellite. These synchrotron/SSC models predict the unattenuated VHE fluxes, which are compared with the observations by imaging atmospheric Cherenkov telescopes. This comparison provides an estimate of the optical depth of the EBL, which allows a derivation of the CGRH through a maximum likelihood analysis that is EBL-model independent. We find that the observed CGRH is compatible with the current knowledge of the EBL.

preprint2012arXiv

Linearisation with Cosmological Perturbation Theory

We propose a new method to linearise cosmological mass density fields using higher order Lagrangian perturbation theory (LPT). We demonstrate that a given density field can be expressed as the sum of a linear and a nonlinear component which are tightly coupled to each other by the tidal field tensor within the LPT framework. The linear component corresponds to the initial density field in Eulerian coordinates, and its mean relation with the total field can be approximated by a logarithm (giving theoretical support to recent attempts to find such component). We also propose to use a combination of the linearisation method and the continuity equation to find the mapping between Eulerian and Lagrangian coordinates. In addition, we note that this method opens the possibility of use directly higher order LPT on nonlinear fields. We test our linearization scheme by applying it to the z~0.5 density field from an N-body simulation. We find that the linearised version of the full density field can be successfully recovered on >~5 h^{-1}Mpc, reducing the skewness and kurtosis of the distribution by about one and two orders of magnitude, respectively. This component can also be successfully traced back in time, converging towards the initial unevolved density field at z~100. We anticipate a number of applications of our results, from predicting velocity fields to estimates of the initial conditions of the universe, passing by improved constraints on cosmological parameters derived from galaxy clustering via reconstruction methods.

preprint2012arXiv

Non-Gaussian gravitational clustering field statistics

In this work we investigate the multivariate statistical description of the matter distribution in the nonlinear regime. We introduce the multivariate Edgeworth expansion of the lognormal distribution to model the cosmological matter field. Such a technique could be useful to generate and reconstruct three-dimensional nonlinear cosmological density fields with the information of higher order correlation functions. We explicitly calculate the expansion up to third order in perturbation theory making use of the multivariate Hermite polynomials up to sixth order. The probability distribution function for the matter field includes at this level the two-point, the three-point and the four-point correlation functions. We use the hierarchical model to formulate the higher order correlation functions based on combinations of the two-point correlation function. This permits us to find compact expressions for the skewness and kurtosis terms of the expanded lognormal field which can be efficiently computed. The method is, however, flexible to incorporate arbitrary higher order correlation functions which have analytical expressions. The applications of such a technique can be especially useful to perform weak-lensing or neutral hydrogen 21 cm line tomography, as well as to directly use the galaxy distribution or the Lyman-alpha forest to study structure formation.

preprint2009arXiv

Bayesian non-linear large scale structure inference of the Sloan Digital Sky Survey data release 7

In this work we present the first non-linear, non-Gaussian full Bayesian large scale structure analysis of the cosmic density field conducted so far. The density inference is based on the Sloan Digital Sky Survey data release 7, which covers the northern galactic cap. We employ a novel Bayesian sampling algorithm, which enables us to explore the extremely high dimensional non-Gaussian, non-linear log-normal Poissonian posterior of the three dimensional density field conditional on the data. These techniques are efficiently implemented in the HADES computer algorithm and permit the precise recovery of poorly sampled objects and non-linear density fields. The non-linear density inference is performed on a 750 Mpc cube with roughly 3 Mpc grid-resolution, while accounting for systematic effects, introduced by survey geometry and selection function of the SDSS, and the correct treatment of a Poissonian shot noise contribution. Our high resolution results represent remarkably well the cosmic web structure of the cosmic density field. Filaments, voids and clusters are clearly visible. Further, we also conduct a dynamical web classification, and estimated the web type posterior distribution conditional on the SDSS data.

preprint2009arXiv

Bayesian power-spectrum inference for Large Scale Structure data

We describe an exact, flexible, and computationally efficient algorithm for a joint estimation of the large-scale structure and its power-spectrum, building on a Gibbs sampling framework and present its implementation ARES (Algorithm for REconstruction and Sampling). ARES is designed to reconstruct the 3D power-spectrum together with the underlying dark matter density field in a Bayesian framework, under the reasonable assumption that the long wavelength Fourier components are Gaussian distributed. As a result ARES does not only provide a single estimate but samples from the joint posterior of the power-spectrum and density field conditional on a set of observations. This enables us to calculate any desired statistical summary, in particular we are able to provide joint uncertainty estimates. We apply our method to mock catalogs, with highly structured observational masks and selection functions, in order to demonstrate its ability to reconstruct the power-spectrum from real data sets, while fully accounting for any mask induced mode coupling.

preprint2009arXiv

Fast Hamiltonian sampling for large scale structure inference

In this work we present a new and efficient Bayesian method for nonlinear three dimensional large scale structure inference. We employ a Hamiltonian Monte Carlo (HMC) sampler to obtain samples from a multivariate highly non-Gaussian lognormal Poissonian density posterior given a set of observations. The HMC allows us to take into account the nonlinear relations between the observations and the underlying density field which we seek to recover. As the HMC provides a sampled representation of the density posterior any desired statistical summary, such as the mean, mode or variance, can be calculated from the set of samples. Further, it permits us to seamlessly propagate non-Gaussian uncertainty information to any final quantity inferred from the set of samples. The developed method is extensively tested in a variety of test scenarios, taking into account a highly structured survey geometry and selection effects. Tests with a mock galaxy catalog based on the millennium run show that the method is able to recover the filamentary structure of the nonlinear density field. The results further demonstrate the feasibility of non-Gaussian sampling in high dimensional spaces, as required for precision nonlinear large scale structure inference. The HMC is a flexible and efficient method, which permits for simple extension and incorporation of additional observational constraints. Thus, the method presented here provides an efficient and flexible basis for future high precision large scale structure inference.

preprint2009arXiv

Nucleosynthesis in Electron Capture Supernovae of AGB Stars

We examine nucleosynthesis in the electron capture supernovae of progenitor AGB stars with an O-Ne-Mg core (with the initial stellar mass of 8.8 M_\odot). Thermodynamic trajectories for the first 810 ms after core bounce are taken from a recent state-of-the-art hydrodynamic simulation. The presented nucleosynthesis results are characterized by a number of distinct features that are not shared with those of other supernovae from the collapse of stars with iron core (with initial stellar masses of more than 10 M_\odot). First is the small amount of 56Ni (= 0.002-0.004 M_\odot) in the ejecta, which can be an explanation for observed properties of faint supernovae such as SNe 2008S and 1997D. In addition, the large Ni/Fe ratio is in reasonable agreement with the spectroscopic result of the Crab nebula (the relic of SN 1054). Second is the large production of 64Zn, 70Ge, light p-nuclei (74Se, 78Kr, 84Sr, and 92Mo), and in particular, 90Zr, which originates from the low Y_e (= 0.46-0.49, the number of electrons per nucleon) ejecta. We find, however, that only a 1-2% increase of the minimum Y_e moderates the overproduction of 90Zr. In contrast, the production of 64Zn is fairly robust against a small variation of Y_e. This provides the upper limit of the occurrence of this type of events to be about 30% of all core-collapse supernovae.

preprint2006arXiv

Explosions of O-Ne-Mg Cores, the Crab Supernova, and Subluminous Type II-P Supernovae

We present results of simulations of stellar collapse and explosions in spherical symmetry for progenitor stars in the 8-10 solar mass range with an O-Ne-Mg core. The simulations were continued until nearly one second after core bounce and were performed with the Prometheus/Vertex code with a variable Eddington factor solver for the neutrino transport, including a state-of-the-art treatment of neutrino-matter interactions. Particular effort was made to implement nuclear burning and electron capture rates with sufficient accuracy to ensure a smooth continuation, without transients, from the progenitor evolution to core collapse. Using two different nuclear equations of state (EoSs), a soft version of the Lattimer & Swesty EoS and the significantly stiffer Wolff & Hillebrandt EoS, we found no prompt explosions, but instead delayed explosions, powered by neutrino heating and the neutrino-driven baryonic wind which sets in about 200 ms after bounce. The models eject little nickel (< 0.015 solar masses), explode with an energy of about or slightly more than 10**50 erg, and leave behind neutron stars (NSs) with a baryonic mass near 1.36 solar masses. Different from previous models of such explosions, the ejecta during the first second have a proton-to-baryon ratio of Ye > 0.46, which suggests a chemical composition that is not in conflict with galactic abundances. No low-entropy matter with Ye << 0.5 is ejected. This excludes such explosions as sites of a low-entropy r-process. The low explosion energy and nucleosynthetic implications are compatible with the observed properties of the Crab supernova, and the small nickel mass supports the possibility that our models explain some subluminous Type II-P supernovae.