Refinement Type Inference via Horn Constraint Optimization

We propose a novel method for inferring refinement types of higher-order functional programs. The main advantage of the proposed method is that it can infer maximally preferred (i.e., Pareto optimal) refinement types with respect to a user-specified preference order. The flexible optimization of refinement types enabled by the proposed method paves the way for interesting applications, such as inferring most-general characterization of inputs for which a given program satisfies (or violates) a given safety (or termination) property. Our method reduces such a type optimization problem to a Horn constraint optimization problem by using a new refinement type system that can flexibly reason about non-determinism in programs. Our method then solves the constraint optimization problem by repeatedly improving a current solution until convergence via template-based invariant generation. We have implemented a prototype inference system based on our method, and obtained promising results in preliminary experiments.Comment: 19 page

Theoretical models for classical Cepheids. VIII. Effects of helium and heavy elements abundance on the Cepheid distance scale

Previous nonlinear fundamental pulsation models for classical Cepheids with metal content Z <= 0.02 are implemented with new computations at super-solar metallicity (Z=0.03, 0.04) and selected choices of the helium-to-metal enrichment ratio DeltaY/Delta Z. On this basis, we show that the location into the HR diagram of the Cepheid instability strip is dependent on both metal and helium abundance, moving towards higher effective temperatures with decreasing the metal content (at fixed Y) or with increasing the helium content (at fixed Z). The contributions of helium and metals to the predicted Period-Luminosity and Period-Luminosity-Color relations are discussed, as well as the implications on the Cepheid distance scale. Based on these new results, we finally show that the empirical metallicity correction suggested by Cepheid observations in two fields of the galaxy M101 may be accounted for, provided that the adopted helium-to-metal enrichment ratio is reasonably high (Delta Y/Delta Z ~ 3.5).Comment: 23 pages, including 6 postscript figures, accepted for publication on Ap

Porto Oscillation Code (POSC)

The Porto Oscillation Code (POSC) has been developed in 1995 and improved over the years, with the main goal of calculating linear adiabatic oscillations for models of solar-type stars. It has also been used to estimate the frequencies and eigenfunctions of stars from the pre-main sequence up to the sub-giant phase, having a mass between 0.8 and 4 solar masses. The code solves the linearised perturbation equations of adiabatic pulsations for an equilibrium model using a second order numerical integration method. The possibility of using Richardson extrapolation is implemented. Several options for the surface boundary condition can be used. In this work we briefly review the key ingredients of the calculations, namely the equations, the numerical scheme and the output.Comment: Accepted for publication in Astrophysics and Space Science

Inclusion of turbulence in solar modeling

The general consensus is that in order to reproduce the observed solar p-mode oscillation frequencies, turbulence should be included in solar models. However, until now there has not been any well-tested efficient method to incorporate turbulence into solar modeling. We present here two methods to include turbulence in solar modeling within the framework of the mixing length theory, using the turbulent velocity obtained from numerical simulations of the highly superadiabatic layer of the sun at three stages of its evolution. The first approach is to include the turbulent pressure alone, and the second is to include both the turbulent pressure and the turbulent kinetic energy. The latter is achieved by introducing two variables: the turbulent kinetic energy per unit mass, and the effective ratio of specific heats due to the turbulent perturbation. These are treated as additions to the standard thermodynamic coordinates (e.g. pressure and temperature). We investigate the effects of both treatments of turbulence on the structure variables, the adiabatic sound speed, the structure of the highly superadiabatic layer, and the p-mode frequencies. We find that the second method reproduces the SAL structure obtained in 3D simulations, and produces a p-mode frequency correction an order of magnitude better than the first method.Comment: 10 pages, 12 figure

The underlying physical meaning of the νmax−νc\nu_{\rm max}-\nu_{\rm c} relation

Asteroseismology of stars that exhibit solar-like oscillations are enjoying a growing interest with the wealth of observational results obtained with the CoRoT and Kepler missions. In this framework, scaling laws between asteroseismic quantities and stellar parameters are becoming essential tools to study a rich variety of stars. However, the physical underlying mechanisms of those scaling laws are still poorly known. Our objective is to provide a theoretical basis for the scaling between the frequency of the maximum in the power spectrum ($\nu_{\rm max}$) of solar-like oscillations and the cut-off frequency ($\nu_{\rm c}$). Using the SoHO GOLF observations together with theoretical considerations, we first confirm that the maximum of the height in oscillation power spectrum is determined by the so-called \emph{plateau} of the damping rates. The physical origin of the plateau can be traced to the destabilizing effect of the Lagrangian perturbation of entropy in the upper-most layers which becomes important when the modal period and the local thermal relaxation time-scale are comparable. Based on this analysis, we then find a linear relation between $\nu_{\rm max}$ and $\nu_{\rm c}$, with a coefficient that depends on the ratio of the Mach number of the exciting turbulence to the third power to the mixing-length parameter.Comment: 8 pages, 11 figures. Accepted in A&

Theory of Stellar Oscillations

In recent years, astronomers have witnessed major progresses in the field of stellar physics. This was made possible thanks to the combination of a solid theoretical understanding of the phenomena of stellar pulsations and the availability of a tremendous amount of exquisite space-based asteroseismic data. In this context, this chapter reviews the basic theory of stellar pulsations, considering small, adiabatic perturbations to a static, spherically symmetric equilibrium. It starts with a brief discussion of the solar oscillation spectrum, followed by the setting of the theoretical problem, including the presentation of the equations of hydrodynamics, their perturbation, and a discussion of the functional form of the solutions. Emphasis is put on the physical properties of the different types of modes, in particular acoustic (p-) and gravity (g-) modes and their propagation cavities. The surface (f-) mode solutions are also discussed. While not attempting to be comprehensive, it is hoped that the summary presented in this chapter addresses the most important theoretical aspects that are required for a solid start in stellar pulsations research.Comment: Lecture presented at the IVth Azores International Advanced School in Space Sciences on "Asteroseismology and Exoplanets: Listening to the Stars and Searching for New Worlds" (arXiv:1709.00645), which took place in Horta, Azores Islands, Portugal in July 201

Tomography of the Solar Interior

Solar oscillations consist of a rich spectrum of internal acoustic waves and surface gravity waves, stochastically excited by turbulent convection. They have been monitored almost continuously over the last ten years with high-precision Doppler images of the solar surface. The purpose of helioseismology is to retrieve information about the structure and the dynamics of the solar interior from the frequencies, phases, and amplitudes of solar waves. Methods of analysis are being developed to make three-dimensional images of subsurface motions and temperature inhomogeneities in order to study convective structures and regions of magnetic activity, like sunspots.Comment: Brief Revie

Probing the properties of convective cores through g modes: high-order g modes in SPB and gamma Doradus stars

In main sequence stars the periods of high-order gravity modes are sensitive probes of stellar cores and, in particular, of the chemical composition gradient that develops near the outer edge of the convective core. We present an analytical approximation of high-order g modes that takes into account the effect of the mu gradient near the core. We show that in main-sequence models, similarly to the case of white dwarfs, the periods of high-order gravity modes are accurately described by a uniform period spacing superposed to an oscillatory component. The periodicity and amplitude of such component are related, respectively, to the location and sharpness of the mu gradient. We investigate the properties of high-order gravity modes for stellar models in a mass domain between 1 and 10 Msun, and the effects of the stellar mass, evolutionary state, and extra-mixing processes on period spacing features. In particular, we show that for models of a typical SPB star, a chemical mixing that could likely be induced by the slow rotation observed in these stars, is able to significantly change the g-mode spectra of the equilibrium model. Prospects and challenges for the asteroseismology of gamma Doradus and SPB stars are also discussed.Comment: 18 pages, 29 figures, accepted for publication in MNRA

Correlated Gravitational Wave and Neutrino Signals from General-Relativistic Rapidly Rotating Iron Core Collapse

We present results from a new set of 3D general-relativistic hydrodynamic simulations of rotating iron core collapse. We assume octant symmetry and focus on axisymmetric collapse, bounce, the early postbounce evolution, and the associated gravitational wave (GW) and neutrino signals. We employ a finite-temperature nuclear equation of state, parameterized electron capture in the collapse phase, and a multi-species neutrino leakage scheme after bounce. The latter captures the important effects of deleptonization, neutrino cooling and heating and enables approximate predictions for the neutrino luminosities in the early evolution after core bounce. We consider 12-solar-mass and 40-solar-mass presupernova models and systematically study the effects of (i) rotation, (ii) progenitor structure, and (iii) postbounce neutrino leakage on dynamics, GW, and, neutrino signals. We demonstrate, that the GW signal of rapidly rotating core collapse is practically independent of progenitor mass and precollapse structure. Moreover, we show that the effects of neutrino leakage on the GW signal are strong only in nonrotating or slowly rotating models in which GW emission is not dominated by inner core dynamics. In rapidly rotating cores, core bounce of the centrifugally-deformed inner core excites the fundamental quadrupole pulsation mode of the nascent protoneutron star. The ensuing global oscillations (f~700-800 Hz) lead to pronounced oscillations in the GW signal and correlated strong variations in the rising luminosities of antineutrino and heavy-lepton neutrinos. We find these features in cores that collapse to protoneutron stars with spin periods <~ 2.5 ms and rotational energies sufficient to drive hyper-energetic core-collapse supernova explosions. Hence, joint GW + neutrino observations of a core collapse event could deliver strong evidence for or against rapid core rotation. [abridged]Comment: 29 pages, 14 figures. Replaced with version matching published versio

Whole Earth Telescope observations of the hot helium atmosphere pulsating white dwarf EC 20058-5234

We present the analysis of a total of 177h of high-quality optical time-series photometry of the helium atmosphere pulsating white dwarf (DBV) EC 20058-5234. The bulk of the observations (135h) were obtained during a WET campaign (XCOV15) in July 1997 that featured coordinated observing from 4 southern observatory sites over an 8-day period. The remaining data (42h) were obtained in June 2004 at Mt John Observatory in NZ over a one-week observing period. This work significantly extends the discovery observations of this low-amplitude (few percent) pulsator by increasing the number of detected frequencies from 8 to 18, and employs a simulation procedure to confirm the reality of these frequencies to a high level of significance (1 in 1000). The nature of the observed pulsation spectrum precludes identification of unique pulsation mode properties using any clearly discernable trends. However, we have used a global modelling procedure employing genetic algorithm techniques to identify the n, l values of 8 pulsation modes, and thereby obtain asteroseismic measurements of several model parameters, including the stellar mass (0.55 M_sun) and T_eff (~28200 K). These values are consistent with those derived from published spectral fitting: T_eff ~ 28400 K and log g ~ 7.86. We also present persuasive evidence from apparent rotational mode splitting for two of the modes that indicates this compact object is a relatively rapid rotator with a period of 2h. In direct analogy with the corresponding properties of the hydrogen (DAV) atmosphere pulsators, the stable low-amplitude pulsation behaviour of EC 20058 is entirely consistent with its inferred effective temperature, which indicates it is close to the blue edge of the DBV instability strip. (abridged)Comment: 19 pages, 8 figures, 5 tables, MNRAS accepte