Numerical Simulations of Mass Transfer in Binaries with Bipolytropic Components

We present the first self-consistent, three dimensional study of hydrodynamic simulations of mass transfer in binary systems with bipolytropic (composite polytropic) components. In certain systems, such as contact binaries or during the common envelope phase, the core-envelope structure of the stars plays an important role in binary interactions. In this paper, we compare mass transfer simulations of bipolytropic binary systems in order to test the suitability of our numerical tools for investigating the dynamical behaviour of such systems. The initial, equilibrium binary models possess a core-envelope structure and are obtained using the bipolytropic self-consistent field technique. We conduct mass transfer simulations using two independent, fully three-dimensional, Eulerian codes - Flow-ER and Octo-tiger. These hydrodynamic codes are compared across binary systems undergoing unstable as well as stable mass transfer, and the former at two resolutions. The initial conditions for each simulation and for each code are chosen to match closely so that the simulations can be used as benchmarks. Although there are some key differences, the detailed comparison of the simulations suggests that there is remarkable agreement between the results obtained using the two codes. This study puts our numerical tools on a secure footing, and enables us to reliably simulate specific mass transfer scenarios of binary systems involving components with a core-envelope structure

VADER: A Flexible, Robust, Open-Source Code for Simulating Viscous Thin Accretion Disks

The evolution of thin axisymmetric viscous accretion disks is a classic problem in astrophysics. While models based on this simplified geometry provide only approximations to the true processes of instability-driven mass and angular momentum transport, their simplicity makes them invaluable tools for both semi-analytic modeling and simulations of long-term evolution where two- or three-dimensional calculations are too computationally costly. Despite the utility of these models, the only publicly-available frameworks for simulating them are rather specialized and non-general. Here we describe a highly flexible, general numerical method for simulating viscous thin disks with arbitrary rotation curves, viscosities, boundary conditions, grid spacings, equations of state, and rates of gain or loss of mass (e.g., through winds) and energy (e.g., through radiation). Our method is based on a conservative, finite-volume, second-order accurate discretization of the equations, which we solve using an unconditionally-stable implicit scheme. We implement Anderson acceleration to speed convergence of the scheme, and show that this leads to factor of $\sim 5$ speed gains over non-accelerated methods in realistic problems, though the amount of speedup is highly problem-dependent. We have implemented our method in the new code Viscous Accretion Disk Evolution Resource (VADER), which is freely available for download from https://bitbucket.org/krumholz/vader/ under the terms of the GNU General Public License.Comment: 58 pages, 13 figures, accepted to Astronomy & Computing; this version includes more discussion, but no other changes; code is available for download from https://bitbucket.org/krumholz/vader

The potential of the variable DA white dwarf G117-B15A as a tool for Fundamental Physics

White dwarfs are well studied objects. The relative simplicity of their physics allows to obtain very detailed models which can be ultimately compared with their observed properties. Among white dwarfs there is a specific class of stars, known as ZZ-Ceti objects, which have a hydrogen-rich envelope and show periodic variations in their light curves. G117-B15A belongs to this particular set of stars. The luminosity variations have been successfully explained as due to g-mode pulsations. G117-B15A has been recently claimed to be the most stable optical clock ever found, being the rate of change of its 215.2 s period very small: \dot{P}= (2.3 +- 1.4)x10^{-15} s s^-1, with a stability comparable to that of the most stable millisecond pulsars. The rate of change of the period is closely related to its cooling timescale, which can be accurately computed. In this paper we study the pulsational properties of G117-B15A and we use the observed rate of change of the period to impose constraints on the axion emissivity and, thus, to obtain a preliminary upper bound to the mass of the axion. This upper bound turns out to be 4cos^{2}{\beta} meV at the 95% confidence level. Although there are still several observational and theoretical uncertainties, we conclude that G117-B15A is a very promising stellar object to set up constraints on particle physics.Comment: 32 pages, 14 figures, accepted for publication in New Astronom

Evolution of galaxies due to self-excitation

These lectures will cover methods for studying the evolution of galaxies since their formation. Because the properties of a galaxy depend on its history, an understanding of galaxy evolution requires that we understand the dynamical interplay between all components. The first part will emphasize n-body simulation methods which minimize sampling noise. These techniques are based on harmonic expansions and scale linearly with the number of bodies, similar to Fourier transform solutions used in cosmological simulations. Although fast, until recently they were only efficiently used for small number of geometries and background profiles. These same techniques may be used to study the modes and response of a galaxy to an arbitrary perturbation. In particular, I will describe the modal spectra of stellar systems and role of damped modes which are generic to stellar systems in interactions and appear to play a significant role in determining the common structures that we see. The general development leads indirectly to guidelines for the number of particles necessary to adequately represent the gravitational field such that the modal spectrum is resolvable. I will then apply these same excitation to understanding the importance of noise to galaxy evolution.Comment: 24 pages, 7 figures, using Sussp.sty (included). Lectures presented at the NATO Advanced Study Institute, "The Restless Universe: Applications of Gravitational N-Body Dynamics to Planetary, Stellar and Galactic Systems," Blair Atholl, July 200