Accurate masses and radii of normal stars: modern results and applications

This paper presents and discusses a critical compilation of accurate, fundamental determinations of stellar masses and radii. We have identified 95 detached binary systems containing 190 stars (94 eclipsing systems, and alpha Centauri) that satisfy our criterion that the mass and radius of both stars be known to 3% or better. To these we add interstellar reddening, effective temperature, metal abundance, rotational velocity and apsidal motion determinations when available, and we compute a number of other physical parameters, notably luminosity and distance. We discuss the use of this information for testing models of stellar evolution. The amount and quality of the data also allow us to analyse the tidal evolution of the systems in considerable depth, testing prescriptions of rotational synchronisation and orbital circularisation in greater detail than possible before. The new data also enable us to derive empirical calibrations of M and R for single (post-) main-sequence stars above 0.6 M(Sun). Simple, polynomial functions of T(eff), log g and [Fe/H] yield M and R with errors of 6% and 3%, respectively. Excellent agreement is found with independent determinations for host stars of transiting extrasolar planets, and good agreement with determinations of M and R from stellar models as constrained by trigonometric parallaxes and spectroscopic values of T(eff) and [Fe/H]. Finally, we list a set of 23 interferometric binaries with masses known to better than 3%, but without fundamental radius determinations (except alpha Aur). We discuss the prospects for improving these and other stellar parameters in the near future.Comment: 56 pages including figures and tables. To appear in The Astronomy and Astrophysics Review. Ascii versions of the tables will appear in the online version of the articl

Order-of-magnitude speedup for steady states and traveling waves via Stokes preconditioning in Channelflow and Openpipeflow

Steady states and traveling waves play a fundamental role in understanding hydrodynamic problems. Even when unstable, these states provide the bifurcation-theoretic explanation for the origin of the observed states. In turbulent wall-bounded shear flows, these states have been hypothesized to be saddle points organizing the trajectories within a chaotic attractor. These states must be computed with Newton's method or one of its generalizations, since time-integration cannot converge to unstable equilibria. The bottleneck is the solution of linear systems involving the Jacobian of the Navier-Stokes or Boussinesq equations. Originally such computations were carried out by constructing and directly inverting the Jacobian, but this is unfeasible for the matrices arising from three-dimensional hydrodynamic configurations in large domains. A popular method is to seek states that are invariant under numerical time integration. Surprisingly, equilibria may also be found by seeking flows that are invariant under a single very large Backwards-Euler Forwards-Euler timestep. We show that this method, called Stokes preconditioning, is 10 to 50 times faster at computing steady states in plane Couette flow and traveling waves in pipe flow. Moreover, it can be carried out using Channelflow (by Gibson) and Openpipeflow (by Willis) without any changes to these popular spectral codes. We explain the convergence rate as a function of the integration period and Reynolds number by computing the full spectra of the operators corresponding to the Jacobians of both methods.Comment: in Computational Modelling of Bifurcations and Instabilities in Fluid Dynamics, ed. Alexander Gelfgat (Springer, 2018

Explosive higher-order Kuramoto dynamics on simplicial complexes

The higher-order interactions of complex systems, such as the brain, are captured by their simplicial complex structure and have a significant effect on dynamics. However the existing dynamical models defined on simplicial complexes make the strong assumption that the dynamics resides exclusively on the nodes. Here we formulate the higher-order Kuramoto model which describes the interactions between oscillators placed not only on nodes but also on links, triangles, and so on. We show that higher-order Kuramoto dynamics can lead to explosive synchronization transition by using an adaptive coupling dependent on the solenoidal and the irrotational component of the dynamics

Two-dimensional Vortices in Superconductors

Superconductors have two key characteristics. They expel magnetic field and they conduct electrical current with zero resistance. However, both properties are compromised in high magnetic fields which can penetrate the material and create a mixed state of quantized vortices. The vortices move in response to an electrical current dissipating energy which destroys the zero resistance state\cite{And64}. One of the central problems for applications of high temperature superconductivity is the stabilization of vortices to ensure zero electrical resistance. We find that vortices in the anisotropic superconductor Bi$_{2}$Sr$_{2}$CaCu$_{2}$O$_{8+\delta}$ (Bi-2212) have a phase transition from a liquid state, which is inherently unstable, to a two-dimensional vortex solid. We show that at high field the transition temperature is independent of magnetic field, as was predicted theoretically for the melting of an ideal two-dimensional vortex lattice\cite{Fis80,Gla91}. Our results indicate that the stable solid phase can be reached at any field as may be necessary for applications involving superconducting magnets\cite{Has04,Sca04,COHMAG}. The vortex solid is disordered, as suggested by previous studies at lower fields\cite{Lee93,Cub93}. But its evolution with increasing magnetic field displays unexpected threshold behavior that needs further investigation.Comment: 5 pages and 4 figures. submitted to Nature Physic

Planetary Dynamics and Habitable Planet Formation In Binary Star Systems

Whether binaries can harbor potentially habitable planets depends on several factors including the physical properties and the orbital characteristics of the binary system. While the former determines the location of the habitable zone (HZ), the latter affects the dynamics of the material from which terrestrial planets are formed (i.e., planetesimals and planetary embryos), and drives the final architecture of the planets assembly. In order for a habitable planet to form in a binary star system, these two factors have to work in harmony. That is, the orbital dynamics of the two stars and their interactions with the planet-forming material have to allow terrestrial planet formation in the habitable zone, and ensure that the orbit of a potentially habitable planet will be stable for long times. We have organized this chapter with the same order in mind. We begin by presenting a general discussion on the motion of planets in binary stars and their stability. We then discuss the stability of terrestrial planets, and the formation of potentially habitable planets in a binary-planetary system.Comment: 56 pages, 29 figures, chapter to appear in the book: Planets in Binary Star Systems (Ed. N. Haghighipour, Springer publishing company

Fermi Gamma-ray Imaging of a Radio Galaxy

The Fermi Gamma-ray Space Telescope has detected the gamma-ray glow emanating from the giant radio lobes of the radio galaxy Centaurus A. The resolved gamma-ray image shows the lobes clearly separated from the central active source. In contrast to all other active galaxies detected so far in high-energy gamma-rays, the lobe flux constitutes a considerable portion (>1/2) of the total source emission. The gamma-ray emission from the lobes is interpreted as inverse Compton scattered relic radiation from the cosmic microwave background (CMB), with additional contribution at higher energies from the infrared-to-optical extragalactic background light (EBL). These measurements provide gamma-ray constraints on the magnetic field and particle energy content in radio galaxy lobes, and a promising method to probe the cosmic relic photon fields.Comment: 27 pages, includes Supplementary Online Material; corresponding authors: C.C. Cheung, Y. Fukazawa, J. Knodlseder, L. Stawar

Detection of Gamma-Ray Emission from the Starburst Galaxies M82 and NGC 253 with the Large Area Telescope on Fermi

We report the detection of high-energy gamma-ray emission from two starburst galaxies using data obtained with the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. Steady point-like emission above 200 MeV has been detected at significance levels of 6.8 sigma and 4.8 sigma respectively, from sources positionally coincident with locations of the starburst galaxies M82 and NGC 253. The total fluxes of the sources are consistent with gamma-ray emission originating from the interaction of cosmic rays with local interstellar gas and radiation fields and constitute evidence for a link between massive star formation and gamma-ray emission in star-forming galaxies.Comment: Submitted to ApJ Letter

The stellar halo of the Galaxy

Stellar halos may hold some of the best preserved fossils of the formation history of galaxies. They are a natural product of the merging processes that probably take place during the assembly of a galaxy, and hence may well be the most ubiquitous component of galaxies, independently of their Hubble type. This review focuses on our current understanding of the spatial structure, the kinematics and chemistry of halo stars in the Milky Way. In recent years, we have experienced a change in paradigm thanks to the discovery of large amounts of substructure, especially in the outer halo. I discuss the implications of the currently available observational constraints and fold them into several possible formation scenarios. Unraveling the formation of the Galactic halo will be possible in the near future through a combination of large wide field photometric and spectroscopic surveys, and especially in the era of Gaia.Comment: 46 pages, 16 figures. References updated and some minor changes. Full-resolution version available at http://www.astro.rug.nl/~ahelmi/stellar-halo-review.pd

Triadic percolation induces dynamical topological patterns in higher-order networks

Triadic interactions are higher-order interactions which occur when a set of nodes affects the interaction between two other nodes. Examples of triadic interactions are present in the brain when glia modulate the synaptic signals among neuron pairs or when interneuron axo-axonic synapses enable presynaptic inhibition and facilitation, and in ecosystems when one or more species can affect the interaction among two other species. On random graphs, triadic percolation has been recently shown to turn percolation into a fully fledged dynamical process in which the size of the giant component undergoes a route to chaos. However, in many real cases, triadic interactions are local and occur on spatially embedded networks. Here, we show that triadic interactions in spatial networks induce a very complex spatio-temporal modulation of the giant component which gives rise to triadic percolation patterns with significantly different topology. We classify the observed patterns (stripes, octopus, and small clusters) with topological data analysis and we assess their information content (entropy and complexity). Moreover, we illustrate the multistability of the dynamics of the triadic percolation patterns, and we provide a comprehensive phase diagram of the model. These results open new perspectives in percolation as they demonstrate that in presence of spatial triadic interactions, the giant component can acquire a time-varying topology. Hence, this work provides a theoretical framework that can be applied to model realistic scenarios in which the giant component is time dependent as in neuroscience

The stellar and sub-stellar IMF of simple and composite populations

The current knowledge on the stellar IMF is documented. It appears to become top-heavy when the star-formation rate density surpasses about 0.1Msun/(yr pc^3) on a pc scale and it may become increasingly bottom-heavy with increasing metallicity and in increasingly massive early-type galaxies. It declines quite steeply below about 0.07Msun with brown dwarfs (BDs) and very low mass stars having their own IMF. The most massive star of mass mmax formed in an embedded cluster with stellar mass Mecl correlates strongly with Mecl being a result of gravitation-driven but resource-limited growth and fragmentation induced starvation. There is no convincing evidence whatsoever that massive stars do form in isolation. Various methods of discretising a stellar population are introduced: optimal sampling leads to a mass distribution that perfectly represents the exact form of the desired IMF and the mmax-to-Mecl relation, while random sampling results in statistical variations of the shape of the IMF. The observed mmax-to-Mecl correlation and the small spread of IMF power-law indices together suggest that optimally sampling the IMF may be the more realistic description of star formation than random sampling from a universal IMF with a constant upper mass limit. Composite populations on galaxy scales, which are formed from many pc scale star formation events, need to be described by the integrated galactic IMF. This IGIMF varies systematically from top-light to top-heavy in dependence of galaxy type and star formation rate, with dramatic implications for theories of galaxy formation and evolution.Comment: 167 pages, 37 figures, 3 tables, published in Stellar Systems and Galactic Structure, Vol.5, Springer. This revised version is consistent with the published version and includes additional references and minor additions to the text as well as a recomputed Table 1. ISBN 978-90-481-8817-