HPC4Mobilty w/ UCB

The purpose of this project is to examine the energy impact of urban-scale traffic for the Los Angeles Basin by developing and implementing a scalable traffic assignment model. An energy optimization function will be posed and when integrated into the optimization code for travel assignment it can be mathematically proven to converge. The energy optimization function can then be compared to the typical travel time optimization that is traditionally used in traffic assignment models. The analysis will begin with static traffic assignment models with the routing for all origin and destinations computed in parallel on high performance computing facilities. Convergence of the numerical methods rely on the solution of convex programs (or extensions of these). This step will mostly consist of demonstrating the ability to parallelize the Frank Wolfe algorithm on various platforms. This work will contribute to LBNL’s efforts to develop new processes, analytical tools, program designs, and business models to advance the state of the art in next-generation sustainable transportation solutions

A Generalization of Abel Inversion to non axisymmetric density distribution

Abel Inversion is currently used in laser-plasma studies in order to estimate the electronic density $n_e$ from the phase-shift map $\delta \phi$ obtained via interferometry. The main limitation of the Abel method is due to the assumption of axial symmetry of the electronic density, which is often hardly fulfilled. In this paper we present an improvement to the Abel inversion technique in which the axial symmetry condition is relaxed by means of a truncated Legendre Polinomial expansion in the azimutal angle. With the help of simulated interferograms, we will show that the generalized Abel inversion generates accurate densities maps when applied to non axisymmetric density sources

The volumetric rate of calcium-rich transients in the local universe

We present a measurement of the volumetric rate of `calcium-rich' optical transients in the local universe, using a sample of three events from the Palomar Transient Factory (PTF). This measurement builds on a detailed study of the PTF transient detection efficiencies, and uses a Monte Carlo simulation of the PTF survey. We measure the volumetric rate of calcium-rich transients to be higher than previous estimates: $1.21^{+1.13}_{-0.39}\times10^{-5}$ events yr$^{-1}$ Mpc$^{-3}$. This is equivalent to 33-94% of the local volumetric type Ia supernova rate. This calcium-rich transient rate is sufficient to reproduce the observed calcium abundances in galaxy clusters, assuming an asymptotic calcium yield per calcium-rich event of ~0.05$\mathrm{M}_{\odot}$. We also study the PTF detection efficiency of these transients as a function of position within their candidate host galaxies. We confirm as a real physical effect previous results that suggest calcium-rich transients prefer large physical offsets from their host galaxies.Comment: Accepted for publication in ApJ. 9 pages, 5 figure

Synthetic Spectra for Type Ia Supernovae at Early Epochs

We present the current status of our construction of synthetic spectra for type Ia supernovae. These properly take into account the effects of NLTE and an adequate representation of line blocking and blanketing. The models are based on a sophisticated atomic database. We show that the synthetic spectrum reproduces the observed spectrum of 'normal' SN-Ia near maximum light from the UV to the near-IR. However, further improvements are necessary before truly quantitative analyses of observed SN-Ia spectra can be performed. In particular, the inner boundary condition has to be fundamentally modified. This is due to the dominance of electron scattering over true absorption processes coupled with the flat density structure in these objectsComment: To appear in "Proceedings of the IAU Colloquium 192 - Supernovae (10 Years of SN1993J)", eds. J.M. Marcaide and K.W. Weile

Evidence for a Spectroscopic Sequence Among SNe Ia

In this Letter we present evidence for a spectral sequence among Type Ia supernovae (SNe Ia). The sequence is based on the systematic variation of several features seen in the near-maximum light spectrum. This sequence is analogous to the recently noted photometric sequence among SNe Ia which shows a relationship between the peak brightness of a SN Ia and the shape of its light curve. In addition to the observational evidence we present a partial theoretical explanation for the sequence. This has been achieved by producing a series of non-LTE synthetic spectra in which only the effective temperature is varied. The synthetic sequence nicely reproduces most of the differences seen in the observed one and presumably corresponds to the amount of 56Ni produced in the explosion.Comment: To appear in Astrophysical Journal Letters, uuencoded, gzipped postscript file, also available from http://www.nhn.uoknor.edu/~baron

Forward Modeling of Double Neutron Stars: Insights from Highly-Offset Short Gamma-Ray Bursts

We present a detailed analysis of two well-localized, highly offset short gamma-ray bursts---GRB~070809 and GRB~090515---investigating the kinematic evolution of their progenitors from compact object formation until merger. Calibrating to observations of their most probable host galaxies, we construct semi-analytic galactic models that account for star formation history and galaxy growth over time. We pair detailed kinematic evolution with compact binary population modeling to infer viable post-supernova velocities and inspiral times. By populating binary tracers according to the star formation history of the host and kinematically evolving their post-supernova trajectories through the time-dependent galactic potential, we find that systems matching the observed offsets of the bursts require post-supernova systemic velocities of hundreds of kilometers per second. Marginalizing over uncertainties in the stellar mass--halo mass relation, we find that the second-born neutron star in the GRB~070809 and GRB~090515 progenitor systems received a natal kick of $\gtrsim 200~\mathrm{km\,s}^{-1}$ at the 78\% and 91\% credible levels, respectively. Applying our analysis to the full catalog of localized short gamma-ray bursts will provide unique constraints on their progenitors and help unravel the selection effects inherent to observing transients that are highly offset with respect to their hosts.Comment: 18 pages, 7 figures, 1 table. ApJ, in pres

MAESTRO, CASTRO, and SEDONA -- Petascale Codes for Astrophysical Applications

Performing high-resolution, high-fidelity, three-dimensional simulations of Type Ia supernovae (SNe Ia) requires not only algorithms that accurately represent the correct physics, but also codes that effectively harness the resources of the most powerful supercomputers. We are developing a suite of codes that provide the capability to perform end-to-end simulations of SNe Ia, from the early convective phase leading up to ignition to the explosion phase in which deflagration/detonation waves explode the star to the computation of the light curves resulting from the explosion. In this paper we discuss these codes with an emphasis on the techniques needed to scale them to petascale architectures. We also demonstrate our ability to map data from a low Mach number formulation to a compressible solver.Comment: submitted to the Proceedings of the SciDAC 2010 meetin

MC generator TAUOLA: implementation of Resonance Chiral Theory for two and three meson modes. Comparison with experiment

We present a partial upgrade of the Monte Carlo event generator TAUOLA with the two and three hadron decay modes using the theoretical models based on Resonance Chiral Theory. These modes account for 88% of total hadronic width of the tau meson. First results of the model parameters have been obtained using BaBar data for three pion mode.Comment: 5 pages, 1 figure, contribution to the Proceedings of the QCD@Work12 Conferenc

Theoretical inputs and errors in the new hadronic currents in TAUOLA

The new hadronic currents implemented in the TAUOLA library are obtained in the unified and consistent framework of Resonance Chiral Theory: a Lagrangian approach in which the resonances exchanged in the hadronic tau decays are active degrees of freedom included in a way that reproduces the low-energy results of Chiral Perturbation Theory. The short-distance QCD constraints on the imaginary part of the spin-one correlators yield relations among the couplings that render the theory predictive. In this communication, the obtaining of the two- and three-meson form factors is sketched. One of the criticisms to our framework is that the error may be as large as 1/3, since it is a realization of the large-N_C limit of QCD in a meson theory. A number of arguments are given which disfavor that claim pointing to smaller errors, which would explain the phenomenological success of our description in these decays. Finally, other minor sources of error and current improvements of the code are discussed.Comment: 5 pages, no figures, contribution to the Proceedings of the QCD@Work12 Conferenc