The Latest Research Papers in Astrophysics

Astrophysics Of Galaxies

---

Disruption of satellite galaxies in simulated groups and clusters: the roles of accretion time, baryons, and pre-processing (1901.03336v2)

Yannick M. Bahé, Joop Schaye, David J. Barnes, Claudio Dalla Vecchia, Scott T. Kay, Richard G. Bower, Henk Hoekstra, Sean L. McGee, Tom Theuns

2019-01-10

We investigate the disruption of group and cluster satellite galaxies with total mass (dark matter plus baryons) above 10^10 M_sun in the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical zoom-in simulations based on the EAGLE model. The simulations predict that ~50 per cent of satellites survive to redshift z = 0, with higher survival fractions in massive clusters than in groups and only small differences between baryonic and pure N-body simulations. For clusters, up to 90 per cent of galaxy disruption occurs in lower-mass sub-groups (i.e., during pre-processing); 96 per cent of satellites in massive clusters that were accreted at z < 2 and have not been pre-processed survive. Of those satellites that are disrupted, only a few per cent merge with other satellites, even in low-mass groups. The survival fraction changes rapidly from less than 10 per cent of those accreted at high z to more than 90 per cent at low z. This shift, which reflects faster disruption of satellites accreted at higher z, happens at lower z for more massive galaxies and those accreted onto less massive haloes. The disruption of satellite galaxies is found to correlate only weakly with their pre-accretion baryon content, star formation rate, and size, so that surviving galaxies are nearly unbiased in these properties. These results suggest that satellite disruption in massive haloes is uncommon, and that it is predominantly the result of gravitational rather than baryonic processes.

Measuring the local matter density using Gaia DR2 (1811.07911v2)

Axel Widmark

2018-11-19

We determine the total dynamical matter density in the solar neighbourhood using the second Gaia data release (DR2). The dynamical matter density distribution is inferred in a framework of a Bayesian hierarchical model, which accounts for position and velocity of all individual stars, as well as the full error covariance matrix of astrometric observables, in a joint fit of the vertical velocity distribution and stellar number density distribution. This was done for eight separate data samples, with different cuts in observed absolute magnitude, each containing about 25,000 stars. The model for the total matter density does not rely on any underlying baryonic model, although we assumed that it is symmetrical, smooth, and monotonically decreasing with distance from the mid-plane. We infer a density distribution which is strongly peaked in the region close to the Galactic plane ( pc), for all eight stellar samples. Assuming a baryonic model and a dark matter halo of constant density, this corresponds to a surplus surface density of approximately 5--9 M/pc. For the Sun's position and vertical velocity with respect to the Galactic plane, we infer pc and km/s. These results suggest a surplus of matter close to the Galactic plane, possibly explained by an underestimated density of cold gas. We discuss possible systematic effects that could bias our result, for example unmodelled non-equilibrium effects, and how to account for such effects in future extensions of this work.

AREPO-RT: Radiation hydrodynamics on a moving mesh (1804.01987v2)

Rahul Kannan, Mark Vogelsberger, Federico Marinacci, Ryan McKinnon, Rüdiger Pakmor, Volker Springel

2018-04-05

We introduce AREPO-RT, a novel radiation hydrodynamic (RHD) solver for the unstructured moving-mesh code AREPO. Our method solves the moment-based radiative transfer equations using the M1 closure relation. We achieve second order convergence by using a slope limited linear spatial extrapolation and a first order time prediction step to obtain the values of the primitive variables on both sides of the cell interface. A Harten-Lax-Van Leer flux function, suitably modified for moving meshes, is then used to solve the Riemann problem at the interface. The implementation is fully conservative and compatible with the individual timestepping scheme of AREPO. It incorporates atomic Hydrogen (H) and Helium (He) thermochemistry, which is used to couple the ultra-violet (UV) radiation field to the gas. Additionally, infrared radiation is coupled to the gas under the assumption of local thermodynamic equilibrium between the gas and the dust. We successfully apply our code to a large number of test problems, including applications such as the expansion of regions, radiation pressure driven outflows and the levitation of optically thick layer of gas by trapped IR radiation. The new implementation is suitable for studying various important astrophysical phenomena, such as the effect of radiative feedback in driving galactic scale outflows, radiation driven dusty winds in high redshift quasars, or simulating the reionisation history of the Universe in a self consistent manner.

Testing Stellar Population Fitting Ingredients with Globular Clusters I: Stellar Libraries (1902.00420v1)

Lucimara P. Martins, Círia Lima-Dias, Paula R. T. Coelho, Tatiana F. Laganá

2019-02-01

The integrated spectra of stellar systems contain a wealth of information, and its analysis can reveal fundamental parameters such as metallicity, age and star formation history. Widely used methods to analyze these spectra are based on comparing the galaxy spectra to stellar population (SP) models. Despite being a powerful tool, SP models contain many ingredients, each with their assumptions and uncertainties. Among the several possible sources of uncertainties, it is not straightforward to identify which ingredient dominates the errors in the models. In this work we propose a study of one of the SP model ingredients -- the spectral stellar libraries -- independently of the other ingredients. To that aim, we will use the integrated spectra of globular clusters which have color-magnitude diagrams (CMDs) available. From these CMDs it is possible to model the integrated spectra of these objects without having to adopt -- or make assumptions -- on the other two main ingredients of SP models, evolutionary tracks or an IMF. Here we tested four widely used stellar libraries. We found that the libraries are able to reproduce the integrated spectra of 18 of the 30 cluster spectra inside a mean flux uncertainty of 5%. For the larger wavelength range tested, a theoretical library outperforms the empirical ones in the comparison. Without the blue part of the spectra, empirical libraries fare better than the theoretical, in particular when individual features are concerned. However, the results are promising for theoretical libraries, which are equally efficient to reproduce the whole spectrum.

The Role of Internal Photons on the Chemistry of the Circumstellar Envelopes of AGB Stars (1902.00416v1)

M Van de Sande, T J Millar

2019-02-01

Recent high spatial resolution observations of gas and dust in the circumstellar envelopes (CSEs) of AGB stars indicate morphologies much more complex than the smooth density distributions generated by spherically symmetric, constant mass loss rates. In particular, the observation of spiral arcs and disks indicate the likely presence of a binary companion which in some cases give rise to the UV photons detected by GALEX. In this Article, we extend our recent model of the chemistry in a clumpy, porous CSE around an AGB star to include the influence of stellar blackbody photons on the CSE chemistry. Our results indicate that internal photons, in a clumpy, porous CSE, can alter chemistry within a few stellar radii and, for some molecules, alter abundances out to several hundred stellar radii. They further suggest that harder radiation from companion stars or accretion disks will have a substantial impact on chemistry in the dust formation zones and inner CSEs of AGB stars.

Earth And Planetary Astrophysics

---

Stellar Spectroscopy in the Near-infrared with a Laser Frequency Comb (1902.00500v1)

Andrew J. Metcalf, Tyler Anderson, Chad F. Bender, Scott Blakeslee, Wesley Brand, David R. Carlson, William D. Cochran, Scott A. Diddams, Michael Endl, Connor Fredrick, Sam Halverson, Dan D. Hickstein, Fred Hearty, Jeff Jennings, Shubham Kanodia, Kyle F. Kaplan, Eric Levi, Emily Lubar, Suvrath Mahadevan, Andrew Monson, Joe P. Ninan, Colin Nitroy, Steve Osterman, Scott B. Papp, Franklyn Quinlan, Larry Ramsey, Paul Robertson, Arpita Roy, Christian Schwab, Steinn Sigurdsson, Kartik Srinivasan, Gudmundur Stefansson, David A. Sterner, Ryan Terrien, Alex Wolszczan, Jason T. Wright, Gabriel Ycas

2019-02-01

The discovery and characterization of exoplanets around nearby stars is driven by profound scientific questions about the uniqueness of Earth and our Solar System, and the conditions under which life could exist elsewhere in our Galaxy. Doppler spectroscopy, or the radial velocity (RV) technique, has been used extensively to identify hundreds of exoplanets, but with notable challenges in detecting terrestrial mass planets orbiting within habitable zones. We describe infrared RV spectroscopy at the 10 m Hobby-Eberly telescope that leverages a 30 GHz electro-optic laser frequency comb with nanophotonic supercontinuum to calibrate the Habitable Zone Planet Finder spectrograph. Demonstrated instrument precision <10 cm/s and stellar RVs approaching 1 m/s open the path to discovery and confirmation of habitable zone planets around M-dwarfs, the most ubiquitous type of stars in our Galaxy.

On the 9:7 Mean Motion Resonance Capture in a System of Two Equal-mass Super-Earths (1902.00437v1)

Zijia Cui, John C. B. Papaloizou, Ewa Szuszkiewicz

2019-02-01

We study the formation of the 9:7 mean motion resonance in a system of two low-mass planets () embedded in a gaseous protoplanetary disk employing a full 2D hydrodynamic treatment of the disk-planet interactions. Our aim is to determine the disk properties that favor a capture of two equal-mass super-Earths into this second-order resonance. For this purpose, we have performed a series of numerical hydrodynamic simulations of the system of two super-Earths migrating in disks with a variety of different initial parameters and found conditions for the permanent or temporary locking in the 9:7 resonance. We observe that capture occurs during the convergent migration of planets if their resonance angle at the moment of arrival at the resonance assumes values in a certain range (inside a window of capture). The width of such a window depends on the relative migration and circularization rates that are determined by the disk parameters. The window is wide if the relative migration rate is slow, and it becomes narrower as the relative migration rate increases. The window will be closed if the migration rate is sufficiently high, and the capture will not take place. We illustrate also how the 9:7 resonance window of capture is affected by the initial eccentricities and the initial orbits of the planets.

Bolometric treatment of irradiation effects: general discussion and application to binary stars (1809.07963v2)

Martin Horvat, Kyle E. Conroy, David Jones, Andrej Prša

2018-09-21

A general framework for dealing with irradiation effects in the bolometric sense --- specifically, reflection with heat absorption and the consequent redistribution of the absorbed heat, for systems of astrophysical bodies where the boundaries are used as support for the description of the processes, is presented. Discussed are its mathematical and physical properties, as well as its implementation approximations, with a focus on three plausible redistribution processes (uniform, latitudinal, and local redistribution). These are tested by extending PHOEBE 2.1 ( http://phoebe-project.org ), the open-source package for modeling eclipsing binaries, and applied to a toy model of the known two-body eclipsing systems.

MASER: A Toolbox for Low Frequency Radio Astronomy (1902.00300v1)

Baptiste Cecconi, Alan Loh, Pierre Le Sidaner, Renaud Savalle, Xavier Bonnin, Quynh Nhu Nguyen, Sonny Lion, Albert Shih, Stéphane Aicardi, Philippe Zarka, Corentin Louis, Andrée Coffre, Laurent Lamy, Laurent Denis, Jean-Mathias Grießmeier, Jeremy Faden, Chris Piker, Nicolas André, Vincent Génot, Stéphane Erard, Joseph N Mafi, Todd A King, Mark Sharlow, Jim Sky, Markus Demleitner

2019-02-01

The MASER (Measurements, Analysis, and Simulation of Emission in the Radio range) project provides a comprehensive infrastructure dedicated to low frequency radio emissions (typically below 50 to 100 MHz). The main radio sources observed in this frequency are the Earth, the Sun, Jupiter, Saturn, and our Galaxy. They are observed either from ground (down to 10 MHz) or from space. Ground observatories are more sensitive than space observatories and capture high resolution data streams (up to a few TB per day for modern instruments). Conversely, space-borne instruments can observe below the ionospheric cut-off (10 MHz) and can be placed closer to the studied object. Several tools have been developed in the last decade for sharing space physics data. Data visualization tools developed by, e.g., The CDPP (http://cdpp.eu, Centre de Donn'ees de la Physique des Plasmas, in Toulouse, France) or the University of Iowa (Autoplot, http://autoplot.org) are available to display and analyze space physics time series and spectrograms. Other tools include ExPRES (Exoplanetary and Planetary Radio Emission Simulator) developed at LESIA. The VESPA (Virtual European Solar and Planetary Access) which provides a search interface that allows the discovery of data of interest for scientific users, and is based on IVOA standards (astronomical International Virtual Observatory Alliance). The University of Iowa has developed the Das2 server that allows the distribution of data with adjustable temporal resolution. MASER is making use of all these tools and standards to distribute datasets from space and ground radio instruments available from the Observatoire de Paris, the Station de Radioastronomie de Nancay and the CDPP deep archive. MASER also includes a Python software library for reading raw data from agency archives.

A generic frequency dependence for the atmospheric tidal torque of terrestrial planets (1902.00280v1)

Pierre Auclair-Desrotour, Jérémy Leconte, Cyril Mergny

2019-02-01

Thermal atmospheric tides have a strong impact on the rotation of terrestrial planets. They can lock these planets into an asynchronous rotation state of equilibrium. We aim at characterizing the dependence of the tidal torque resulting from the semidiurnal thermal tide on the tidal frequency, the planet orbital radius, and the atmospheric surface pressure. The tidal torque is computed from full 3D simulations of the atmospheric climate and mean flows using a generic version of the LMDZ general circulation model (GCM) in the case of a nitrogen-dominated atmosphere. Numerical results are discussed with the help of an updated linear analytical framework. Power scaling laws governing the evolution of the torque with the planet orbital radius and surface pressure are derived. The tidal torque exhibits i) a thermal peak in the vicinity of synchronization, ii) a resonant peak associated with the excitation of the Lamb mode in the high frequency range, and iii) well defined frequency slopes outside these resonances. These features are well explained by our linear theory. Whatever the star-planet distance and surface pressure, the torque frequency spectrum -- when rescaled with the relevant power laws -- always presents the same behaviour. This allows us to provide a single and easily usable empirical formula describing the atmospheric tidal torque over the whole parameter space. With such a formula, the effect of the atmospheric tidal torque can be implemented in evolutionary models of the rotational dynamics of a planet in a computationally efficient, and yet relatively accurate way.

Solar And Stellar Astrophysics

---

Long-Term Cyclic Variability of YZ CMi in the Context of Solar and Stellar Physics (1902.00493v1)

N. I. Bondar, M. M. Katsova

2019-02-01

Manifestations of activity on the YZ CMi (M4.5Ve) flare star which has a rotation period of 2.77 days and belongs to a group of activity-saturated stars, are considered from long-term photometric data together with photographic measurements from 1926 to 2009. Long-term changes in yearly mean brightness of the star are found on a time scale of 27.5 years with an amplitude of that indicate variations in development of surface inhomogeneities and a large degree of surface spottedness in the maximum activity. Spots are distributed unevenly over the surface, their concentration is higher at certain longitudes spaced by intervals equal to 0.6 phases of the rotation period. The recovery of active longitude (a possible "flip-flop" effect) occurs over a period of about 6 years. We note differences in activity of this star from typical for the Sun and solar-like stars, which are associated with its dynamic characteristics and and its inner structure.

Resolved imaging of the AR Puppis circumbinary disk (1901.04109v3)

S. Ertel, D. Kamath, M. Hillen, H. van Winckel, J. Okumura, R. Manick, H. M. J. Boffin, J. Milli, G. H. -M. Bertrang, L. Guzman-Ramirez, J. Horner, J. P. Marshall, P. Scicluna, A. Vaz, E. Villaver, R. Wesson, S. Xu

2019-01-14

Circumbinary disks are common around post-asymptotic giant branch (post-AGB) stars with a stellar companion on orbital time scales of a few 100 to few 1000 days. The presence of a disk is usually inferred from the system's spectral energy distribution (SED), and confirmed, for a sub-sample, by interferometric observations. We used the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument on the Very Large Telescope to obtain extreme adaptive optics assisted scattered light images of the post-AGB binary system AR Puppis. Data have been obtained in the V, I, and H bands. Our observations have produced the first resolved images of AR Puppis' circumbinary disk and confirm its edge-on orientation. In our high angular-resolution and high dynamic-range images we identify several structural components such as a dark mid-plane, the disk surface, and arc-like features. We discuss the nature of these components and use complementary photometric monitoring to relate them to the orbital phase of the binary system. Because the star is completely obscured by the disk at visible wavelengths, we conclude that the long-term photometric variability of the system must be caused by variable scattering, not extinction, of star light by the disk over the binary orbit. Finally, we discuss how the short disk life times and fast evolution of the host stars compared to the ages at which protoplanetary disks are typically observed make systems like AR Puppis valuable extreme laboratories to study circumstellar disk evolution and constrain the time scale of dust grain growth during the planet formation process.

Carbon and oxygen in metal-poor halo stars (1901.03592v2)

A. M. Amarsi, P. E. Nissen, M. Asplund, K. Lind, P. S. Barklem

2019-01-11

Carbon and oxygen are key tracers of the Galactic chemical evolution; in particular, a reported upturn in [C/O] towards decreasing [O/H] in metal-poor halo stars could be a signature of nucleosynthesis by massive Population III stars. We reanalyse carbon, oxygen, and iron abundances in thirty-nine metal-poor turn-off stars. For the first time, we take into account three-dimensional (3D) hydrodynamic effects together with departures from local thermodynamic equilibrium (LTE) when determining both the stellar parameters and the elemental abundances, by deriving effective temperatures from 3D non-LTE H profiles, surface gravities from Gaia parallaxes, iron abundances from 3D LTE Feii equivalent widths, and carbon and oxygen abundances from 3D non-LTE Ci and Oi equivalent widths. We find that [C/Fe] stays flat with [Fe/H], whereas [O/Fe] increases linearly up to dex with decreasing [Fe/H] down to dex. As such [C/O] monotonically decreases towards decreasing [O/H], in contrast to previous findings, mainly by virtue of less severe non-LTE effects for Oi at low [Fe/H] with our improved calculations.

STiC -- A multi-atom non-LTE PRD inversion code for full-Stokes solar observations (1810.08441v3)

J. de la Cruz Rodríguez, J. Leenaarts, S. Danilovic, H. Uitenbroek

2018-10-19

The inference of the underlying state of the plasma in the solar chromosphere remains extremely challenging because of the nonlocal character of the observed radiation and plasma conditions in this layer. Inversion methods allow us to derive a model atmosphere that can reproduce the observed spectra by undertaking several physical assumptions. The most advanced approaches involve a depth-stratified model atmosphere described by temperature, line-of-sight velocity, turbulent velocity, the three components of the magnetic field vector, and gas and electron pressure. The parameters of the radiative transfer equation are computed from a solid ground of physical principles. To apply these techniques to spectral lines that sample the chromosphere, NLTE effects must be included in the calculations. We developed a new inversion code STiC to study spectral lines that sample the upper chromosphere. The code is based the RH synthetis code, which we modified to make the inversions faster and more stable. For the first time, STiC facilitates the processing of lines from multiple atoms in non-LTE, also including partial redistribution effects. Furthermore, we include a regularization strategy that allows for model atmospheres with a complex stratification, without introducing artifacts in the reconstructed physical parameters, which are usually manifested in the form of oscillatory behavior. This approach takes steps toward a node-less inversion, in which the value of the physical parameters at each grid point can be considered a free parameter. In this paper we discuss the implementation of the aforementioned techniques, the description of the model atmosphere, and the optimizations that we applied to the code. We carry out some numerical experiments to show the performance of the code and the regularization techniques that we implemented. We made STiC publicly available to the community.

Testing Stellar Population Fitting Ingredients with Globular Clusters I: Stellar Libraries (1902.00420v1)

Lucimara P. Martins, Círia Lima-Dias, Paula R. T. Coelho, Tatiana F. Laganá

2019-02-01

The integrated spectra of stellar systems contain a wealth of information, and its analysis can reveal fundamental parameters such as metallicity, age and star formation history. Widely used methods to analyze these spectra are based on comparing the galaxy spectra to stellar population (SP) models. Despite being a powerful tool, SP models contain many ingredients, each with their assumptions and uncertainties. Among the several possible sources of uncertainties, it is not straightforward to identify which ingredient dominates the errors in the models. In this work we propose a study of one of the SP model ingredients -- the spectral stellar libraries -- independently of the other ingredients. To that aim, we will use the integrated spectra of globular clusters which have color-magnitude diagrams (CMDs) available. From these CMDs it is possible to model the integrated spectra of these objects without having to adopt -- or make assumptions -- on the other two main ingredients of SP models, evolutionary tracks or an IMF. Here we tested four widely used stellar libraries. We found that the libraries are able to reproduce the integrated spectra of 18 of the 30 cluster spectra inside a mean flux uncertainty of 5%. For the larger wavelength range tested, a theoretical library outperforms the empirical ones in the comparison. Without the blue part of the spectra, empirical libraries fare better than the theoretical, in particular when individual features are concerned. However, the results are promising for theoretical libraries, which are equally efficient to reproduce the whole spectrum.

---

Thank you for reading!

Don't forget to Follow and Resteem. @astrophysics
Keeping everyone inform.