The Latest Research Papers in Astrophysics

Astrophysics Of Galaxies

---

Stochastic transport of high-energy particles through a turbulent plasma (1808.04430v2)

L. E. Chen, A. F. A. Bott, P. Tzeferacos, A. Rigby, A. Bell, R. Bingham, C. Graziani, J. Katz, M. Koenig, C. K. Li, R. Petrasso, H. -S. Park, J. S. Ross, D. Ryu, D. Ryutov, T. G. White, B. Reville, J. Matthews, J. Meinecke, F. Miniati, E. G. Zweibel, S. Sarkar, A. A. Schekochihin, D. Q. Lamb, D. H. Froula, G. Gregori

2018-08-13

The interplay between charged particles and turbulent magnetic fields is crucial to understanding how cosmic rays propagate through space. A key parameter which controls this interplay is the ratio of the particle gyroradius to the correlation length of the magnetic turbulence. For the vast majority of cosmic rays detected at the Earth, this parameter is small, and the particles are well confined by the Galactic magnetic field. But for cosmic rays more energetic than about 30 EeV, this parameter is large. These highest energy particles are not confined to the Milky Way and are presumed to be extragalactic in origin. Identifying their sources requires understanding how they are deflected by the intergalactic magnetic field, which appears to be weak, turbulent with an unknown correlation length, and possibly spatially intermittent. This is particularly relevant given the recent detection by the Pierre Auger Observatory of a significant dipole anisotropy in the arrival directions of cosmic rays of energy above 8 EeV. Here we report measurements of energetic-particle propagation through a random magnetic field in a laser-produced plasma. We characterize the diffusive transport of these particles and recover experimentally pitch-angle scattering measurements and extrapolate to find their mean free path and the associated diffusion coefficient, which show scaling-relations consistent with theoretical studies. This experiment validates these theoretical tools for analyzing the propagation of ultra-high energy cosmic rays through the intergalactic medium.

Tidal disruption event rates in galaxy merger remnants (1903.09124v1)

Hugo Pfister, Ben Bar-Or, Marta Volonteri, Yohan Dubois, Pedro R. Capelo

2019-03-21

The rate of tidal disruption events (TDEs) depends sensitively on the stellar properties of the central galactic regions. Simulations show that galaxy mergers cause gas inflows that trigger nuclear starbursts that increase the central stellar density. Motivated by these numerical results, and by the observed over-representation of post-starburst galaxies among TDE hosts, we study the evolution of the TDE rate in high-resolution hydrodynamical simulations of a galaxy merger, in which we capture the evolution of the stellar density around the massive black holes (BHs) in the merging galaxies and in the final remnant. We post-process the simulation to obtain the TDE rate. At the second pericentre, a nuclear starburst enhances the stellar density around the BH in the least massive galaxy, leading to an enhancement of the TDE rate around the secondary BH, although the magnitude and the duration of the increase depend on the stochasticity of star formation on very small scales. The central stellar density around the primary BH remains instead fairly constant, and so is its TDE rate, until the two BHs bind and are surrounded by the same stellar environment. After the formation of the binary, the stellar density decreases, and so does the TDE rate.

Derivation and precision of mean field electrodynamics with mesoscale fluctuations (1710.04064v5)

Hongzhe Zhou, Eric G. Blackman, Luke Chamandy

2017-10-07

Mean field electrodynamics (MFE) facilitates practical modeling of secular, large scale properties of astrophysical or laboratory systems with fluctuations.Practitioners commonly assume wide scale separation between mean and fluctuating quantities, to justify equality of ensemble and spatial or temporal averages.Often however, real systems do not exhibit such scale separation. This raises two questions: (I) what are the appropriate generalized equations of MFE in the presence of mesoscale fluctuations? (II) how precise are theoretical predictions from MFE? We address both by first deriving the equations of MFE for different types of averaging, along with mesoscale correction terms that depend on the ratio of averaging scale to variation scale of the mean. We then show that even if these terms are small, predictions of MFE can still have a significant precision error. This error has an intrinsic contribution from the dynamo input parameters and a filtering contribution from differences in the way observations and theory are projected through the measurement kernel.Minimizing the sum of these contributions can produce an optimal scale of averaging that makes the theory maximally precise.The precision error is important to quantify when comparing to observations because it quantifies the resolution of predictive power. We exemplify these principles for galactic dynamos, comment on broader implications, and identify possibilities for further work.

Bulge plus disc and Sérsic decomposition catalogues for 16,908 galaxies in the SDSS Stripe 82 co-adds: A detailed study of the structural measurements (1903.09068v1)

Connor Bottrell, Luc Simard, J. Trevor Mendel, Sara L. Ellison

2019-03-21

Quantitative characterization of galaxy morphology is vital in enabling comparison of observations to predictions from galaxy formation theory. However, without significant overlap between the observational footprints of deep and shallow galaxy surveys, the extent to which structural measurements for large galaxy samples are robust to image quality (e.g., depth, spatial resolution) cannot be established. Deep images from the Sloan Digital Sky Survey (SDSS) Stripe 82 co-adds provide a unique solution to this problem - offering magnitudes improvement in depth with respect to SDSS Legacy images. Having similar spatial resolution to Legacy, the co-adds make it possible to examine the sensitivity of parametric morphologies to depth alone. Using the Gim2D surface-brightness decomposition software, we provide public morphology catalogs for 16,908 galaxies in the Stripe 82 co-adds. Our methods and selection are completely consistent with the Simard et al. (2011) and Mendel et al. (2014) photometric decompositions. We rigorously compare measurements in the deep and shallow images. We find no systematics in total magnitudes and sizes except for faint galaxies in the -band and the brightest galaxies in each band. However, characterization of bulge-to-total fractions is significantly improved in the deep images. Furthermore, statistics used to determine whether single-S'ersic or two-component (e.g., bulge+disc) models are required become more bimodal in the deep images. Lastly, we show that asymmetries are enhanced in the deep images and that the enhancement is positively correlated with the asymmetries measured in Legacy images.

The shape alignment of satellite galaxies in galaxy pairs in SDSS (1805.06096v2)

Peng Wang, Quan Guo, Noam I. Libeskind, Elmo Tempel, Chengliang Wei, Xi Kang

2018-05-16

It has been shown, both in simulations and observationally, that the tidal field of a large galaxy can torque its satellites such that the major axis of satellite galaxies points towards their hosts. This so-called shape alignment' has been observed in isolated Milky Way-like galaxies but not inLocal Group'-like pairs. In this study, we investigate the shape alignment of satellite galaxies in galaxy pairs similar to the Local Group identified in the Sloan Digital Sky Survey Data Release 13 (SDSS DR13). By stacking tens of thousands of satellite galaxies around primary galaxy pairs, we find two statistically strong alignment signals. (1) The major axes of satellite galaxies located in the (projected) area between two primaries (the {\it facing} region) tend to be perpendicular to the line connecting the satellite to its host (tangential alignment), while (2) the major axes of satellite galaxies located in regions away from the other host (the {\it away} region) tend to be aligned with the line connecting the satellite to its host (radial alignment). These alignments are confirmed at levels. The alignment signal increases with increasing primary brightness, decreasing pair separation, and decreasing satellite distance. The alignment signal is also found to be stronger in filamentary environments. These findings will shed light on understanding the mechanisms of how satellite galaxies are affected by the tidal field in galaxy pairs and will be useful for investigating galaxy intrinsic alignment in the analyses of weak gravitational lensing.

Earth And Planetary Astrophysics

---

Mapping out the time-evolution of exoplanet processes (1903.09110v1)

Jessie L. Christiansen, Charles Beichman, David R. Ciardi, Daniel Huber, Mark S. Marley

2019-03-21

There are many competing theories and models describing the formation, migration and evolution of exoplanet systems. As both the precision with which we can characterize exoplanets and their host stars, and the number of systems for which we can make such a characterization increase, we begin to see pathways forward for validating these theories. In this white paper we identify predicted, observable correlations that are accessible in the near future, particularly trends in exoplanet populations, radii, orbits and atmospheres with host star age. By compiling a statistically significant sample of well-characterized exoplanets with precisely measured ages, we should be able to begin identifying the dominant processes governing the time-evolution of exoplanet systems.

A high performance and portable all-Mach regime flow solver code with well-balanced gravity. Application to compressible convection (1903.09043v1)

T. Padioleau, P. Tremblin, E. Audit, P. Kestener, S. Kokh

2019-03-21

Convection is an important physical process in astrophysics well-studied using numerical simulations under the Boussinesq and/or anelastic approximations. However these approaches reach their limits when compressible effects are important in the high Mach flow regime, e.g. in stellar atmospheres or in the presence of accretion shocks. In order to tackle these issues, we propose a new high performance and portable code, called "ARK" with a numerical solver well-suited for the stratified compressible Navier-Stokes equations. We take a finite volume approach with machine precision conservation of mass, transverse momentum and total energy. Based on previous works in applied mathematics we propose the use of a low Mach correction to achieve a good precision in both low and high Mach regimes. The gravity source term is discretized using a well-balanced scheme in order to reach machine precision hydrostatic balance. This new solver is implemented using the Kokkos library in order to achieve high performance computing and portability across different architectures (e.g. multi-core, many-core, and GP-GPU). We show that the low-Mach correction allows to reach the low-Mach regime with a much better accuracy than a standard Godunov-type approach. The combined well-balanced property and the low-Mach correction allowed us to trigger Rayleigh-B'enard convective modes close to the critical Rayleigh number. Furthermore we present 3D turbulent Rayleigh-B'enard convection with low diffusion using the low-Mach correction leading to a higher kinetic energy power spectrum. These results are very promising for future studies of high Mach and highly stratified convective problems in astrophysics.

Surface evolution of the Anhur region on comet 67P from high-resolution OSIRIS images (1903.09017v1)

S. Fornasier, C. Feller, P. H. Hasselmann, M. A. Barucci, J. Sunshine, J. -B. Vincent, X. Shi, H. Sierks, G. Naletto, P. L. Lamy, R. Rodrigo, D. Koschny, B. Davidsson, J. -L. Bertaux, I. Bertini, D. Bodewits, G. Cremonese, V. Da Deppo, S. Debei, M. De Cecco, J. Deller, S. Ferrari, M. Fulle, P. J. Gutierrez, C. Güttler, W. -H. Ip, L. Jorda, H. U. Keller, M. L. Lara, M. Lazzarin, J. J. Lopez Moreno, A. Lucchetti, F. Marzari, S. Mottola, M. Pajola, I. Toth, C. Tubiana

2019-03-21

The southern hemisphere of comet 67P/Churyumov-Gerasimenko (67P) became observable by the Rosetta mission in March 2015, a few months before cometary southern vernal equinox. The Anhur region in the southern part of the comet's larger lobe was found to be highly eroded, enriched in volatiles, and highly active. We analyze high-resolution images of the Anhur region pre- and post-perihelion acquired by the OSIRIS imaging system on board the Rosetta mission. The Narrow Angle Camera is particularly useful for studying the evolution in Anhur in terms of morphological changes and color variations.}{Radiance factor images processed by the OSIRIS pipeline were coregistered, reprojected onto the 3D shape model of the comet, and corrected for the illumination conditions. We find a number of morphological changes in the Anhur region that are related to formation of new scarps; removal of dust coatings; localized resurfacing in some areas, including boulders displacements; and vanishing structures, which implies localized mass loss that we estimate to be higher than 50 million kg. The strongest changes took place in and nearby the Anhur canyon-like structure, where significant dust cover was removed, an entire structure vanished, and many boulders were rearranged. All such changes are potentially associated with one of the most intense outbursts registered by Rosetta during its observations, which occurred one day before perihelion passage. Moreover, in the niche at the foot of a new observed scarp, we also see evidence of water ice exposure that persisted for at least six months. The abundance of water ice, evaluated from a linear mixing model, is relatively high (> 20%). Our results confirm that the Anhur region is volatile-rich and probably is the area on 67P with the most pristine exposures near perihelion.

Analytic solutions to the maximum and average exoplanet transit depth for common stellar limb darkening laws (1901.01730v2)

René Heller

2019-01-07

The depth of an exoplanetary transit in the light curve of a distant star is commonly approximated as the squared planet-to-star radius ratio, (R_p/R_s)^2. Stellar limb darkening, however, results in significantly deeper transits. Here we derive analytical solutions to the overshoot of the mid-transit depth caused by stellar limb darkening compared to the (R_p/R_s)^2 estimate for arbitrary transit impact parameters. In turn, this allows us to compute the true planet-to-star radius ratio from the transit depth for a given parameterization of a limb darkening law and for a known transit impact parameter. We calculate the maximum emerging specific stellar intensity covered by the planet in transit and derive analytic solutions for the transit depth overshoot. Solutions are presented for the linear, quadratic, square-root, logarithmic, and non-linear stellar limb darkening with arbitrary transit impact parameters. We also derive formulae to calculate the average intensity along the transit chord, which allows us to estimate the actual transit depth (and therefore R_p/R_s) from the mean in-transit flux. The transit depth overshoot of exoplanets compared to the (R_p/R_s)^2 estimate increases from about 15% for A main-sequence stars to roughly 20% for sun-like stars and some 30% for K and M stars. The error in our analytical solutions for R_p/R_s from the small planet approximation is orders of magnitude smaller than the uncertainties arising from typical noise in real light curves and from the uncertain limb darkening. Our equations can be used to predict with high accuracy the expected transit depth of extrasolar planets. The actual planet radius can be calculated from the measured transit depth or from the mean in-transit flux if the stellar limb darkening can be properly parameterized and if the transit impact parameter is known. Light curve fitting is not required.

"Auxiliary" Science with the WFIRST Microlensing Survey (1903.08986v1)

B. Scott Gaudi, Rachel Akeson, Jay Anderson, Etienne Bachelet, David P. Bennett, Aparna Bhattacharya, Valerio Bozza, Sebastiano Calchi Novati, Calen B. Henderson, Samson A. Johnson, Jeffrey Kruk, Jessica R. Lu, Shude Mao, Benjamin T. Montet, David M. Nataf, Matthew T. Penny, Radoslaw Poleski, Clément Ranc, Kailash Sahu, Yossi Shvartzvald, David N. Spergel, Daisuke Suzuki, Keivan G. Stassun, Rachel A. Street

2019-03-21

The Wide Field Infrared Survey Telescope (WFIRST) will monitor deg toward the Galactic bulge in a wide (m) W149 filter at 15-minute cadence with exposure times of 50s for 6 seasons of 72 days each, for a total 41,000 exposures taken over 432 days, spread over the 5-year prime mission. This will be one of the deepest exposures of the sky ever taken, reaching a photon-noise photometric precision of 0.01 mag per exposure and collecting a total of photons over the course of the survey for a W149 star. Of order stars will be monitored with W149<21, and 10 stars with W145<23. The WFIRST microlensing survey will detect 54,000 microlensing events, of which roughly 1% (500) will be due to isolated black holes, and 3% (1600) will be due to isolated neutron stars. It will be sensitive to (effectively) isolated compact objects with masses as low as the mass of Pluto, thereby enabling a measurement of the compact object mass function over 10 orders of magnitude. Assuming photon-noise limited precision, it will detect transiting planets with sizes as small as , perform asteroseismology of giant stars, measure the proper motions to and parallaxes to for the disk and bulge stars in the survey area, and directly detect Trans-Neptunian objects (TNOs) with diameters down to km, as well as detect occulations of stars by TNOs during the survey. All of this science will completely serendipitous, i.e., it will not require modifications of the WFIRST optimal microlensing survey design. Allowing for some minor deviation from the optimal design, such as monitoring the Galactic center, would enable an even broader range of transformational science.

Solar And Stellar Astrophysics

---

Two is a Crowd? Characterising the effect of photometric contamination on the extraction of the global asteroseismic parameter in Red-Giant Binaries (1903.09146v1)

Sanjay Sekaran, Cole Johnston, Paul G Beck, Andrej Prsa, Kelly M Hambleton

2019-03-21

Theoretical scaling relations for solar-like oscillators and red giants are widely used to estimate fundamental stellar parameters. The accuracy and precision of these relations have often been questioned in the literature, with studies often utilising binarity for model-independent validation. However, it has not been tested if the photometric effects of binarity introduce a systematic effect on the extraction of the seismic properties of the pulsating component(s). In this paper, we present an estimation of the impact of a contaminating photometric signal with a distinct background profile on the global asteroseismic parameter through the analysis of synthetic red-giant binary light curves. We generate the pulsational and granulation parameters for single red giants with different masses, radii and effective temperatures from theoretical scaling relations and use them to simulate single red-giant light curves with the characteristics of Kepler long-cadence photometric data. These are subsequently blended together according to their light ratio to generate binary red-giant light curves of various configurations. We then perform a differential analysis to characterise the systematic effects of binarity on the extraction of . We find that the extraction for red-giant power spectra featuring overlapping power excesses is unreliable if unconstrained priors are used. Outside of this scenario, we obtain results that are nearly identical to single-star case. We conclude that i) the photometric effects of binarity on the extraction of are largely negligible as long as the power excesses of the individual components do not overlap, and that ii) there is minimal advantage to using more than two super-Lorentzian components to model the granulation signal of a binary red-giant.

Inferring Explosion Properties from Type II-Plateau Supernova Light Curves (1903.09114v1)

Jared A. Goldberg, Lars Bildsten, Bill Paxton

2019-03-21

We present advances in modeling Type IIP supernovae using \texttt{MESA} for evolution to shock breakout coupled with \texttt{STELLA} for generating light and radial velocity curves. Explosion models and synthetic light curves can be used to translate observable properties of supernovae (such as the luminosity at day 50 and the duration of the plateau, as well as the observable quantity ) into families of explosions which produce the same light curve and velocities on the plateau. These predicted families of explosions provide a useful guide towards modeling observed SNe, and can constrain explosion properties when coupled with other observational or theoretical constraints. For an observed supernova with a measured mass, breaking the degeneracies within these families of explosions (ejecta mass, explosion energy, and progenitor radius) requires independent knowledge of one parameter. We expect the most common case to be a progenitor radius measurement for a nearby supernova. We show that ejecta velocities inferred from the Fe II 5169 \AA\ line measured during the majority of the plateau phase provide little additional information about explosion characteristics. Only during the initial shock cooling phase can photospheric velocity measurements potentially aid in unraveling light curve degeneracies.

Mapping out the time-evolution of exoplanet processes (1903.09110v1)

Jessie L. Christiansen, Charles Beichman, David R. Ciardi, Daniel Huber, Mark S. Marley

2019-03-21

There are many competing theories and models describing the formation, migration and evolution of exoplanet systems. As both the precision with which we can characterize exoplanets and their host stars, and the number of systems for which we can make such a characterization increase, we begin to see pathways forward for validating these theories. In this white paper we identify predicted, observable correlations that are accessible in the near future, particularly trends in exoplanet populations, radii, orbits and atmospheres with host star age. By compiling a statistically significant sample of well-characterized exoplanets with precisely measured ages, we should be able to begin identifying the dominant processes governing the time-evolution of exoplanet systems.

Observations of red giants with suspected massive companions (1902.03145v2)

Valeri V. Makarov, Andrei Tokovinin

2019-02-08

Motivated by the existence of binary systems where a stellar-mass black hole is bound to a normal star, we selected four red giants with large radial velocity (RV) variation from the survey of SIM grid stars and monitored their RVs for several months. None turned out to contain a massive companion above 2.5 solar masses. The red giant TYC 9299-1080-1 with a large RV and a large proper motion is a single-lined spectroscopic binary with a period of 81 days. It is an extreme halo object moving at 350 km s almost directly toward the Galactic center. HD 206092 is a double-lined binary with a short period of 4.37 days. It belongs to the rare class of active RS CVn-type binaries with evolved primary components, apparently undergoing mass transfer. The X-ray luminosity of HD 206092 is about twice as high as the most luminous coronal X-ray emitters observed by ROSAT, including II Peg and the prototype star RS CVn. HD~318347 has a variable double-peaked emission-line spectrum (not a giant), while HD~324668 has a constant RV. Despite the overall good quality of the SIM survey data confirmed by a comparison with Gaia DR2 mean radial velocities, the few large RV variations are explained, mostly, by erroneous data. We discuss the significance of the non-detection of massive companions in the SIM grid sample and the associated work.

The O star hinterland of the Galactic starburst, NGC 3603 (1903.09053v1)

J. E. Drew, M. Monguió, N. J. Wright

2019-03-21

The very bright and compact massive young cluster, NGC 3603, has been cited as an example of a starburst in the Milky Way and compared with the much-studied R136/30 Doradus region in the Large Magellanic Cloud. Here we build on the discovery by Mohr-Smith et al. (2017) of a large number of reddened O stars around this cluster. We construct a list of 288 candidate O stars with proper motions, in a region of sky spanning 1.5x1.5 square degrees centered on NGC 3603, by cross-matching the Mohr-Smith et al. (2017) catalogue with Gaia DR2 (Gaia Collaboration et al. 2018). This provides the basis for a first comprehensive examination of the proper motions of these massive stars in the halo of NGC 3603, relative to the much better studied central region. We identify up to 11 likely O star ejections -- 8 of which would have been ejected between 0.60 and 0.95 Myr ago (supporting the age of c.1 Myr that has been attributed to the bright cluster centre). Seven candidate ejections are arranged in a partial ring to the south of the cluster core spanning radii of 9-18 arcmin (18-36 pc if the cluster is 7 kpc away). We also show that the cluster has a halo of a further 100 O stars extending to a radius of at least 5 arcmin, adding to the picture of NGC 3603 as a scaled down version of the R136/30 Dor region.

---

Thank you for reading!

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