Relatively Complete Verification of Probabilistic Programs: An Expressive Language for Expectation-Based Reasoning

We study a syntax for specifying quantitative “assertions” - functions mapping program states to numbers - for probabilistic program verification. We prove that our syntax is expressive in the following sense: Given any probabilistic program C, if a function f is expressible in our syntax, then the function mapping each initial state σ to the expected value of f evaluated in the final states reached after termination C on σ (also called the weakest preexpectation wp[C](f)) is also expressible in our syntax. As a consequence, we obtain a relatively complete verification system for verifying expected values and probabilities in the sense of Cook: Apart from a single reasoning step about the inequality of two functions given as syntactic expressions in our language, given f, g, and C, we can check whether g ≤ wp[C](f)

Weighted programming: A programming paradigm for specifying mathematical models

We study weighted programming, a programming paradigm for specifying mathematical models. More specifically, the weighted programs we investigate are like usual imperative programs with two additional features: (1) nondeterministic branching and (2) weighting execution traces. Weights can be numbers but also other objects like words from an alphabet, polynomials, formal power series, or cardinal numbers. We argue that weighted programming as a paradigm can be used to specify mathematical models beyond probability distributions (as is done in probabilistic programming). We develop weakest-precondition- and weakest-liberal-precondition-style calculi à la Dijkstra for reasoning about mathematical models specified by weighted programs. We present several case studies. For instance, we use weighted programming to model the ski rental problem - an optimization problem. We model not only the optimization problem itself, but also the best deterministic online algorithm for solving this problem as weighted programs. By means of weakest-precondition-style reasoning, we can determine the competitive ratio of the online algorithm on source code level

Stellar sources in the ISOGAL intermediate bulge fields

We present a study of ISOGAL sources in the "intermediate" galactic bulge ($|$$l$$|$ $<$ 2$^\circ$, $|$$b$$|$ $\sim$ 1$^\circ$--4$^\circ$), observed by ISOCAM at 7 and 15 $\mu m$. In combination with near-infrared (I, J, K$_{\rm s}$) data of DENIS survey, complemented by 2MASS data, we discuss the nature of the ISOGAL sources, their luminosities, the interstellar extinction and the mass-loss rates. A large fraction of the 1464 detected sources at 15 $\mu m$ are AGB stars above the RGB tip, a number of them show an excess in ([7]-[15])$_{\rm 0}$ and (K$_{\rm s}$-[15])$_{\rm 0}$ colours, characteristic of mass-loss. The latter, especially (K$_{\rm s}$-[15])$_{\rm 0}$, provide estimates of the mass-loss rates and show their distribution in the range 10$^{-8}$ to 10$^{-5}$ M$_{\rm \odot}$/yr.Comment: 16 pages, accepted for publication in Astronomy and Astrophysic

Multisite spectroscopic seismic study of the beta Cep star V2052 Oph: inhibition of mixing by its magnetic field

We used extensive ground-based multisite and archival spectroscopy to derive observational constraints for a seismic modelling of the magnetic beta Cep star V2052 Ophiuchi. The line-profile variability is dominated by a radial mode (f_1=7.14846 d^{-1}) and by rotational modulation (P_rot=3.638833 d). Two non-radial low-amplitude modes (f_2=7.75603 d^{-1} and f_3=6.82308 d^{-1}) are also detected. The four periodicities that we found are the same as the ones discovered from a companion multisite photometric campaign (Handler et al. 2012) and known in the literature. Using the photometric constraints on the degrees l of the pulsation modes, we show that both f_2 and f_3 are prograde modes with (l,m)=(4,2) or (4,3). These results allowed us to deduce ranges for the mass (M \in [8.2,9.6] M_o) and central hydrogen abundance (X_c \in [0.25,0.32]) of V2052 Oph, to identify the radial orders n_1=1, n_2=-3 and n_3=-2, and to derive an equatorial rotation velocity v_eq \in [71,75] km s^{-1}. The model parameters are in full agreement with the effective temperature and surface gravity deduced from spectroscopy. Only models with no or mild core overshooting (alpha_ov \in [0,0.15] local pressure scale heights) can account for the observed properties. Such a low overshooting is opposite to our previous modelling results for the non-magnetic beta Cep star theta Oph having very similar parameters, except for a slower surface rotation rate. We discuss whether this result can be explained by the presence of a magnetic field in V2052 Oph that inhibits mixing in its interior.Comment: 12 pages, 6 figures and 5 tables; accepted for publication in MNRAS on 2012 August 1

ISOGAL: A deep survey of the obscured inner Milky Way with ISO at 7 and 15 micron and with DENIS in the near-infrared

The ISOGAL project is an infrared survey of specific regions sampling the Galactic Plane selected to provide information on Galactic structure,stellar populations,stellar mass-loss and the recent star formation history of the inner disk and Bulge of the Galaxy. ISOGAL combines 7 and 15 micron ISOCAM observations - with a resolution of 6'' at worst - with DENIS IJKs data to determine the nature of the sources and theinterstellar extinction. We have observed about 16 square degrees with a sensitivity approaching 10-20mJy, detecting ~10^5 sources,mostly AGB stars,red giants and young stars. The main features of the ISOGAL survey and the observations are summarized in this paper,together with a brief discussion of data processing and quality. The primary ISOGAL products are described briefly (a full description is given in Schuller et al. 2003, astro-ph/0304309): viz. the images and theISOGAL-DENIS five-wavelength point source catalogue. The main scientific results already derived or in progress are summarized. These include astrometrically calibrated 7 and 15um images,determining structures of resolved sources; identification and properties of interstellar dark clouds; quantification of the infrared extinction law and source dereddening; analysis of red giant and (especially) AGB stellar populations in the central Bulge,determining luminosity,presence of circumstellar dust and mass--loss rate,and source classification,supplemented in some cases by ISO/CVF spectroscopy; detection of young stellar objects of diverse types,especially in the inner Bulge with information about the present and recent star formation rate; identification of foreground sources with mid-IR excess. These results are the subject of about 25 refereed papers published or in preparation.Comment: A&A in press. 19 pages,10 Ps figures; problems with figures fixe

CONCERTO: High-fidelity simulation of millimeter line emissions of galaxies and [CII] intensity mapping

The intensity mapping of the [CII] 158-μm line redshifted to the submillimeter window is a promising probe of the za&gt;4 star formation and its spatial distribution into large-scale structures. To prepare the first-generation experiments (e.g., CONCERTO), we need realistic simulations of the submillimeter extragalactic sky in spectroscopy. We present a new version of the simulated infrared dusty extragalactic sky (SIDES) model including the main submillimeter lines around 1 mm (CO, [CII], [CI]). This approach successfully reproduces the observed line luminosity functions. We then use our simulation to generate CONCERTO-like cubes (125-305 GHz) and forecast the power spectra of the fluctuations caused by the various astrophysical components at those frequencies. Depending on our assumptions on the relation between the star formation rate and [CII] luminosity, and the star formation history, our predictions of the za∼6 [CII] power spectrum vary by two orders of magnitude. This highlights how uncertain the predictions are and how important future measurements will be to improve our understanding of this early epoch. SIDES can reproduce the CO shot noise recently measured at a4;100 GHz by the millimeter-wavelength intensity mapping experiment (mmIME). Finally, we compare the contribution of the different astrophysical components at various redshifts to the power spectra. The continuum is by far the brightest, by a factor of three to 100, depending on the frequency. At 300 GHz, the CO foreground power spectrum is higher than the [CII] one for our base scenario. At lower frequencies, the contrast between [CII] and extragalactic foregrounds is even worse. Masking the known galaxies from deep surveys should allow us to reduce the foregrounds to 20% of the [CII] power spectrum up to z∼ 6.5. However, this masking method will not be sufficient at higher redshifts. The code and the products of our simulation are released publicly, and can be used for both intensity mapping experiments and submillimeter continuum and line surveys

Measurement of the Ar(e,e′^\prime p) and Ti(e,e′^\prime p) cross sections in Jefferson Lab Hall A

The E12-14-012 experiment, performed in Jefferson Lab Hall A, has collected exclusive electron-scattering data (e,e$^\prime$p) in parallel kinematics using natural argon and natural titanium targets. Here, we report the first results of the analysis of the data set corresponding to beam energy of 2,222 MeV, electron scattering angle 21.5 deg, and proton emission angle -50 deg. The differential cross sections, measured with $\sim$4% uncertainty, have been studied as a function of missing energy and missing momentum, and compared to the results of Monte Carlo simulations, obtained from a model based on the Distorted Wave Impulse Approximation.Comment: 14 pages, 8 figures (submitted to PRC

Multiple strategies to prevent oxidative stress in Arabidopsis plants lacking the malate valve enzyme NADP-malate dehydrogenase

The nuclear-encoded chloroplast NADP-dependent malate dehydrogenase (NADP-MDH) is a key enzyme controlling the malate valve, to allow the indirect export of reducing equivalents. Arabidopsis thaliana (L.) Heynh. T-DNA insertion mutants of NADP-MDH were used to assess the role of the light-activated NADP-MDH in a typical C3 plant. Surprisingly, even when exposed to high-light conditions in short days, nadp-mdh knockout mutants were phenotypically indistinguishable from the wild type. The photosynthetic performance and typical antioxidative systems, such as the Beck–Halliwell–Asada pathway, were barely affected in the mutants in response to high-light treatment. The reactive oxygen species levels remained low, indicating the apparent absence of oxidative stress, in the mutants. Further analysis revealed a novel combination of compensatory mechanisms in order to maintain redox homeostasis in the nadp-mdh plants under high-light conditions, particularly an increase in the NTRC/2-Cys peroxiredoxin (Prx) system in chloroplasts. There were indications of adjustments in extra-chloroplastic components of photorespiration and proline levels, which all could dissipate excess reducing equivalents, sustain photosynthesis, and prevent photoinhibition in nadp-mdh knockout plants. Such metabolic flexibility suggests that the malate valve acts in concert with other NADPH-consuming reactions to maintain a balanced redox state during photosynthesis under high-light stress in wild-type plants

Determination of the Titanium Spectral Function From (e, e\u27p) Data

The E12-14-012 experiment, performed in Jefferson Lab Hall A, has measured the (e,e′p) cross section in parallel kinematics using a natural titanium target. In this paper, we report the analysis of the dataset obtained in different kinematics for our solid natural titanium target. Data were obtained in a range of missing momentum and missing energy between 15 ≲ pm ≲ 250  MeV/c and 12 ≲ Em ≲ 80  MeV, respectively, and using an electron beam energy of 2.2 GeV. We measured the reduced cross section with ∼7% accuracy as a function of both missing momentum and missing energy. Our Monte Carlo simulation, including both a model spectral function and the effects of final-state interactions, satisfactorily reproduces the data