Showing 1–50 of 89 results for author: Arras, P

Bayesian polarization calibration and imaging in very long baseline interferometry

Authors: Jong-Seo Kim, Jakob Roth, Jongho Park, Jack D. Livingston, Philipp Arras, Torsten A. Enßlin, Michael Janssen, J. Anton Zensus, Andrei P. Lobanov

Abstract: Extracting polarimetric information from very long baseline interferometry (VLBI) data is demanding but vital for understanding the synchrotron radiation process and the magnetic fields of celestial objects, such as active galactic nuclei (AGNs). However, conventional CLEAN-based calibration and imaging methods provide suboptimal resolution without uncertainty estimation of calibration solutions,… ▽ More Extracting polarimetric information from very long baseline interferometry (VLBI) data is demanding but vital for understanding the synchrotron radiation process and the magnetic fields of celestial objects, such as active galactic nuclei (AGNs). However, conventional CLEAN-based calibration and imaging methods provide suboptimal resolution without uncertainty estimation of calibration solutions, while requiring manual steering from an experienced user. We present a Bayesian polarization calibration and imaging method using Bayesian imaging software resolve for VLBI data sets, that explores the posterior distribution of antenna-based gains, polarization leakages, and polarimetric images jointly from pre-calibrated data. We demonstrate our calibration and imaging method with observations of the quasar 3C273 with the VLBA at 15 GHz and the blazar OJ287 with the GMVA+ALMA at 86 GHz. Compared to the CLEAN method, our approach provides physically realistic images that satisfy positivity of flux and polarization constraints and can reconstruct complex source structures composed of various spatial scales. Our method systematically accounts for calibration uncertainties in the final images and provides uncertainties of Stokes images and calibration solutions. The automated Bayesian approach for calibration and imaging will be able to obtain high-fidelity polarimetric images using high-quality data from next-generation radio arrays. The pipeline developed for this work is publicly available. △ Less

Submitted 20 November, 2025; originally announced November 2025.

Comments: submitted to A&A

Dynamical tide modified Roche limit in eccentric, asynchronous binaries

Authors: Hang Yu, Shu Yan Lau, Ethan Mckeever, Phil Arras, Nevin N. Weinberg

Abstract: The Roche limit, or the threshold separation within which a celestial object (the donor) M cannot remain in a stable configuration due to a companion's tidal field, has been well established when M is in hydrostatic equilibrium and has synchronous rotation in a circular orbit. However, limited analyses exist considering corrections to the Roche limit due to hydrodynamical effects. We fill in the g… ▽ More The Roche limit, or the threshold separation within which a celestial object (the donor) M cannot remain in a stable configuration due to a companion's tidal field, has been well established when M is in hydrostatic equilibrium and has synchronous rotation in a circular orbit. However, limited analyses exist considering corrections to the Roche limit due to hydrodynamical effects. We fill in the gap by providing a general theoretical framework involving nonlinear hydrodynamics. We consider both exact nonlinear equations derived from an affine model describing incompressible ellipsoids and series-expanded ones that can be calculated for realistic stars and planets. Our formulation addresses the Roche problem in generic orbits and synchronization levels of M, and fully accounts for the history-dependent hydrodynamical effects. We show that as the orbital eccentricity increases, fluid instability is more likely to develop at the pericenter due to the increased dynamical tide that accumulates over multiple orbits. When M moves in a highly eccentric orbit (with eccentricity around 0.9) and the damping of the fluid is small, the threshold pericenter separation at which mass loss from M can occur can be at least 30% higher than the value predicted for a circular orbit with hydrostatic equilibrium. If only a single passage is considered, however, the threshold separation is 20% smaller than the static limit. The nonlinear interaction at each pericenter passage can also trigger a chaotic fluid evolution inside M even with moderate eccentricities, complementing previous studies of chaotic tides caused by random propagation phases. Our work has broad implications for interacting binaries in eccentric orbits, including migrating gaseous exoplanets, repeated partial tidal disruption events, and more. △ Less

Submitted 21 January, 2026; v1 submitted 27 August, 2025; originally announced August 2025.

Comments: 31 pages, 8 figures, accepted by ApJ

Bayesian Imaging of Interferometric Data from Polarized Electromagnetic Signals

Authors: Philipp Arras, Jakob Roth, Martin Reinecke, Richard A. Perley, Andrei Frolov, Rüdiger Westermann, Torsten A. Enßlin

Abstract: We present an imaging algorithm for polarimetric interferometric data from radio telescopes. It is based on Bayesian statistics and thereby able to provide uncertainties and to incorporate prior information such as positivity of the total emission (Stokes I) or consistency constraints (polarized fraction can only be between 0% and 100%). By comparing our results to the output of the de-facto stand… ▽ More We present an imaging algorithm for polarimetric interferometric data from radio telescopes. It is based on Bayesian statistics and thereby able to provide uncertainties and to incorporate prior information such as positivity of the total emission (Stokes I) or consistency constraints (polarized fraction can only be between 0% and 100%). By comparing our results to the output of the de-facto standard algorithm called CLEAN, we show that these constraints paired with a consistent treatment of measurement uncertainties throughout the algorithm significantly improve image quality. In particular, our method reveals that depolarization canals in CLEAN images do not necessarily indicate a true absence of polarized emission, e.g., after frequency averaging, but can also stem from uncertainty in the polarization direction. This demonstrates that our Bayesian approach can distinguish between true depolarization and mere uncertainty, providing a more informative representation of polarization structures. △ Less

Submitted 31 March, 2025; originally announced April 2025.

Bayesian self-calibration and imaging in Very Long Baseline Interferometry

Authors: Jong-Seo Kim, Aleksei S. Nikonov, Jakob Roth, Torsten A. Ensslin, Michael Janssen, Philipp Arras, Hendrik Mueller, Andrei P. Lobanov

Abstract: Self-calibration methods with the CLEAN algorithm have been widely employed in Very Long Baseline Interferometry (VLBI) data processing in order to correct antenna-based amplitude and phase corruptions present in the data. However, human interaction during the conventional CLEAN self-calibration process can impose a strong effective prior, which in turn may produce artifacts within the final image… ▽ More Self-calibration methods with the CLEAN algorithm have been widely employed in Very Long Baseline Interferometry (VLBI) data processing in order to correct antenna-based amplitude and phase corruptions present in the data. However, human interaction during the conventional CLEAN self-calibration process can impose a strong effective prior, which in turn may produce artifacts within the final image and hinder the reproducibility of final results. In this work, we aim to demonstrate a combined self-calibration and imaging method for VLBI data in a Bayesian inference framework. The method corrects for amplitude and phase gains for each antenna and polarization mode by inferring the temporal correlation of the gain solutions. We use Stokes I data of M87 taken with the Very Long Baseline Array (VLBA) at 43GHz, pre-calibrated using the rPICARD CASA-based pipeline. For antenna-based gain calibration and imaging, we use the Bayesian imaging software resolve. To estimate gain and image uncertainties, we use a Variational Inference method. We obtain a high-resolution M87 Stokes I image at 43GHz in conjunction with antenna-based gain solutions using our Bayesian self-calibration and imaging method. The core with counter-jet structure is better resolved, and extended jet emission is better described compared to the CLEAN reconstruction. Furthermore, uncertainty estimation of the image and antenna-based gains allows us to quantify the reliability of the result. Our Bayesian self-calibration and imaging method is able to reconstruct robust and reproducible Stokes I images and gain solutions with uncertainty estimation by taking into account the uncertainty information in the data. △ Less

Submitted 20 July, 2024; originally announced July 2024.

Comments: accepted for publication in A&A

Journal ref: A&A 690, A129 (2024)

Resolving Horizon-Scale Dynamics of Sagittarius A*

Authors: Jakob Knollmüller, Philipp Arras, Torsten Enßlin

Abstract: Sagittarius A* (Sgr A*), the supermassive black hole at the heart of our galaxy, provides unique opportunities to study black hole accretion, jet formation, and gravitational physics. The rapid structural changes in Sgr A*'s emission pose a significant challenge for traditional imaging techniques. We present dynamic reconstructions of Sgr A* using Event Horizon Telescope (EHT) data from April 6th… ▽ More Sagittarius A* (Sgr A*), the supermassive black hole at the heart of our galaxy, provides unique opportunities to study black hole accretion, jet formation, and gravitational physics. The rapid structural changes in Sgr A*'s emission pose a significant challenge for traditional imaging techniques. We present dynamic reconstructions of Sgr A* using Event Horizon Telescope (EHT) data from April 6th and 7th, 2017, analyzed with a one-minute temporal resolution with the Resolve framework. This Bayesian approach employs adaptive Gaussian Processes and Variational Inference for data-driven self-regularization. Our results not only fully confirm the initial findings by the EHT Collaboration for a time-averaged source but also reveal intricate details about the temporal dynamics within the black hole environment. We find an intriguing dynamic feature on April 6th that propagates in a clock-wise direction. Geometric modelling with ray-tracing, although not fully conclusive, indicates compatibility with high-inclination configurations of about $θ_o = 160^\circ$, as seen in other studies. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Astrophysical Parameter Inference on Accreting White Dwarf Binaries using Gravitational Waves

Authors: Sophia Yi, Shu Yan Lau, Kent Yagi, Phil Arras

Abstract: Accreting binary white dwarf systems are among the sources expected to emanate gravitational waves that the Laser Interferometer Space Antenna (LISA) will detect. We investigate how accurately the binary parameters may be measured from LISA observations. We complement previous studies by performing our parameter estimation on binaries containing a low-mass donor with a thick, hydrogen-rich envelop… ▽ More Accreting binary white dwarf systems are among the sources expected to emanate gravitational waves that the Laser Interferometer Space Antenna (LISA) will detect. We investigate how accurately the binary parameters may be measured from LISA observations. We complement previous studies by performing our parameter estimation on binaries containing a low-mass donor with a thick, hydrogen-rich envelope. The evolution is followed from the early, pre-period minimum stage, in which the donor is non-degenerate, to a later, post-period minimum stage with a largely degenerate donor. We present expressions for the gravitational wave amplitude, frequency, and frequency derivative in terms of white dwarf parameters (masses, donor radius, etc.), where binary evolution is driven by gravitational wave radiation and accretion torques, and the donor radius and logarithmic change in radius ($η_{\rm d}$) due to mass loss are treated as model parameters. We then perform a Fisher analysis to reveal the accuracy of parameter measurements, using models from Modules for Experiments in Stellar Astrophysics (MESA) to estimate realistic fiducial values at which we evaluate the measurement errors. We find that the donor radius can be measured relatively well with LISA observations alone, while we can further measure the individual masses if we have an independent measurement of the luminosity distance from electromagnetic observations. When applied to the parameters of the recently-discovered white dwarf binary ZTF J0127+5258, our Fisher analysis suggests that we will be able to constrain the system's individual masses and donor radius using LISA's observations, given ZTF's measurement of the luminosity distance. △ Less

Submitted 28 June, 2024; v1 submitted 24 October, 2023; originally announced October 2023.

The Surface Mass Density of the Milky Way: Does the Traditional $K_Z$ Approach Work in the Context of New Surveys?

Authors: Xinlun Cheng, Borja Anguiano, Steven R. Majewski, Phil Arras

Abstract: We revisit the classical $K_Z$ problem -- determination of the vertical force and implied total mass density distribution of the Milky Way disk -- for a wide range of Galactocentric radius and vertical height using chemically selected thin and thick disk samples based on APOGEE spectroscopy combined with the Gaia astrometry. We derived the velocity dispersion profiles in Galactic cylindrical coord… ▽ More We revisit the classical $K_Z$ problem -- determination of the vertical force and implied total mass density distribution of the Milky Way disk -- for a wide range of Galactocentric radius and vertical height using chemically selected thin and thick disk samples based on APOGEE spectroscopy combined with the Gaia astrometry. We derived the velocity dispersion profiles in Galactic cylindrical coordinates, and solved the Jeans Equation for the two samples separately. The result is surprising that the total surface mass density as a function of vertical height as derived for these two chemically distinguished populations are different. The discrepancies are larger in the inner compared to the outer Galaxy, with the density calculated from thick disk being larger, independent of the Galactic radius. Furthermore, while there is an overall good agreement between the total mass density derived for the thick disk population and the Standard Halo Model for vertical heights larger than 1 kpc, close to the midplane the mass density observed using the thick disk population is larger than the predicted from the Standard Halo Model. We explore various implications of these discrepancies, and speculate their sources, including problems associated with the assumed density laws, velocity dispersion profiles, and the Galactic rotation curve, potential non-equilibrium of the Galactic disk, or a failure of the NFW dark matter halo profile for the Milky Way. We conclude that the growing detail in hand on the chemodynamical distributions of Milky Way stars challenges traditional analytical treatments of the $K_Z$ problem. △ Less

Submitted 29 September, 2023; originally announced September 2023.

Comments: Accepted by MNRAS

First spatio-spectral Bayesian imaging of SN1006 in X-ray

Authors: Margret Westerkamp, Vincent Eberle, Matteo Guardiani, Philipp Frank, Lukas Platz, Philipp Arras, Jakob Knollmüller, Julia Stadler, Torsten Enßlin

Abstract: Supernovae are an important source of energy in the interstellar medium. Young remnants of supernovae have a peak emission in the X-ray region, making them interesting objects for X-ray observations. In particular, the supernova remnant SN1006 is of great interest due to its historical record, proximity and brightness. It has therefore been studied by several X-ray telescopes. Improving the X-ray… ▽ More Supernovae are an important source of energy in the interstellar medium. Young remnants of supernovae have a peak emission in the X-ray region, making them interesting objects for X-ray observations. In particular, the supernova remnant SN1006 is of great interest due to its historical record, proximity and brightness. It has therefore been studied by several X-ray telescopes. Improving the X-ray imaging of this and other remnants is important but challenging as it requires to address a spatially varying instrument response in order to achieve a high signal-to-noise ratio. Here, we use Chandra observations to demonstrate the capabilities of Bayesian image reconstruction using information field theory. Our objective is to reconstruct denoised, deconvolved and spatio-spectral resolved images from X-ray observations and to decompose the emission into different morphologies, namely diffuse and point-like. Further, we aim to fuse data from different detectors and pointings into a mosaic and quantify the uncertainty of our result. Utilizing prior knowledge on the spatial and spectral correlation structure of the two components, diffuse emission and point sources, the presented method allows the effective decomposition of the signal into these. In order to accelerate the imaging process, we introduce a multi-step approach, in which the spatial reconstruction obtained for a single energy range is used to derive an informed starting point for the full spatio-spectral reconstruction. The method is applied to 11 Chandra observations of SN1006 from 2008 and 2012, providing a detailed, denoised and decomposed view of the remnant. In particular, the separated view of the diffuse emission should provide new insights into its complex small-scale structures in the center of the remnant and at the shock front profiles. △ Less

Submitted 18 December, 2023; v1 submitted 17 August, 2023; originally announced August 2023.

Large Dynamical Tide Amplitudes from Small Kicks at Pericenter

Authors: Phil Arras, Hang Yu, Nevin N. Weinberg

Abstract: The effect of dynamical tide ``kicks" on eccentric binary orbits is considered using the orbital mapping method. It is demonstrated that when mode damping is negligible the mode amplitude will generically grow in time for all values of orbital eccentricity and semi-major axis, even for small kicks outside the regime exhibiting diffusive growth. The origin of the small-kick growth is the change in… ▽ More The effect of dynamical tide ``kicks" on eccentric binary orbits is considered using the orbital mapping method. It is demonstrated that when mode damping is negligible the mode amplitude will generically grow in time for all values of orbital eccentricity and semi-major axis, even for small kicks outside the regime exhibiting diffusive growth. The origin of the small-kick growth is the change in kick size from orbit to orbit, an effect quadratic in the mode amplitude. When damping of the mode is included, the growth is shut off when the damping time is shorter than the growth time. Hence, in practice, kicks of sufficient size and long mode damping times are required for interesting levels of growth to occur. Application to the circularization of hot Jupiters is discussed. Previous investigations found that diffusive growth of the planetary f-mode in the large-kick regime would lead to rapid orbital shrinkage, but upon exiting the diffusive regime at $e \sim 0.9$ the theory would predict a large population of highly eccentric orbits. Simulations presented here show that subsequent orbital evolution relying on the small-kick regime may further decrease the eccentricity to $e \sim 0.2$ on timescales much less than the Gyrs ages of these systems. △ Less

Submitted 4 June, 2023; originally announced June 2023.

Comments: Submitted to ApJ. 18 pages, 10 figures

Orbital Decay of Hot Jupiters due to Weakly Nonlinear Tidal Dissipation

Authors: Nevin N. Weinberg, Niyousha Davachi, Reed Essick, Hang Yu, Phil Arras, Brent Belland

Abstract: We study tidal dissipation in hot Jupiter host stars due to the nonlinear damping of tidally driven $g$-modes, extending the calculations of Essick & Weinberg (2016) to a wide variety of non-solar type hosts. This process causes the planet's orbit to decay and has potentially important consequences for the evolution and fate of hot Jupiters. Previous studies either only accounted for linear dissip… ▽ More We study tidal dissipation in hot Jupiter host stars due to the nonlinear damping of tidally driven $g$-modes, extending the calculations of Essick & Weinberg (2016) to a wide variety of non-solar type hosts. This process causes the planet's orbit to decay and has potentially important consequences for the evolution and fate of hot Jupiters. Previous studies either only accounted for linear dissipation processes or assumed that the resonantly excited primary mode becomes strongly nonlinear and breaks as it approaches the stellar center. However, the great majority of hot Jupiter systems are in the weakly nonlinear regime in which the primary mode does not break but instead excites a sea of secondary modes via three-mode interactions. We simulate these nonlinear interactions and calculate the net mode dissipation for stars that range in mass from $0.5 M_\odot \le M_\star \le 2.0 M_\odot$ and in age from the early main sequence to the subgiant phase. For stars with $M_\star \lesssim 1.0 M_\odot$ of nearly any age, we find that the orbital decay time is $\lesssim 100 \textrm{ Myr}$ for orbital periods $P_{\rm orb} \lesssim 1 \textrm{ day}$. For $M_\star \gtrsim 1.2 M_\odot$, the orbital decay time only becomes short on the subgiant branch, where it can be $\lesssim 10 \textrm{ Myr}$ for $P_{\rm orb} \lesssim 2 \textrm{ days}$ and result in significant transit time shifts. We discuss these results in the context of known hot Jupiter systems and examine the prospects for detecting their orbital decay with transit timing measurements. △ Less

Submitted 17 May, 2023; originally announced May 2023.

Comments: 22 pages, 7 figures, submitted to ApJ

Bayesian radio interferometric imaging with direction-dependent calibration

Authors: Jakob Roth, Philipp Arras, Martin Reinecke, Richard A. Perley, Rüdiger Westermann, Torsten A. Enßlin

Abstract: Context: Radio interferometers measure frequency components of the sky brightness, modulated by the gains of the individual radio antennas. Due to atmospheric turbulence and variations in the operational conditions of the antennas these gains fluctuate. Thereby the gains do not only depend on time but also on the spatial direction on the sky. To recover high quality radio maps an accurate reconstr… ▽ More Context: Radio interferometers measure frequency components of the sky brightness, modulated by the gains of the individual radio antennas. Due to atmospheric turbulence and variations in the operational conditions of the antennas these gains fluctuate. Thereby the gains do not only depend on time but also on the spatial direction on the sky. To recover high quality radio maps an accurate reconstruction of the direction and time-dependent individual antenna gains is required. Aims: This paper aims to improve the reconstruction of radio images, by introducing a novel joint imaging and calibration algorithm including direction-dependent antenna gains. Methods: Building on the \texttt{resolve} framework, we designed a Bayesian imaging and calibration algorithm utilizing the image domain gridding method for numerically efficient application of direction-dependent antenna gains. Furthermore by approximating the posterior probability distribution with variational inference, our algorithm can provide reliable uncertainty maps. Results: We demonstrate the ability of the algorithm to recover high resolution high dynamic range radio maps from VLA data of the radio galaxy Cygnus A. We compare the quality of the recovered images with previous work relying on classically calibrated data. Furthermore we compare with a compressed sensing algorithm also incorporating direction-dependent gains. Conclusions: Including direction-dependent effects in the calibration model significantly improves the dynamic range of the reconstructed images compared to reconstructions from classically calibrated data. Compared to the compressed sensing reconstruction, the resulting sky images have a higher resolution and show fewer artifacts. For utilizing the full potential of radio interferometric data, it is essential to consider the direction dependence of the antenna gains. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 13 pages, 9 figures

Journal ref: A&A 678, A177 (2023)

The ngEHT Analysis Challenges

Authors: Freek Roelofs, Lindy Blackburn, Greg Lindahl, Sheperd S. Doeleman, Michael D. Johnson, Philipp Arras, Koushik Chatterjee, Razieh Emami, Christian Fromm, Antonio Fuentes, Jakob Knollmueller, Nikita Kosogorov, Hendrik Mueller, Nimesh Patel, Alexander Raymond, Paul Tiede, Thalia Traianou, Justin Vega

Abstract: The next-generation Event Horizon Telescope (ngEHT) will be a significant enhancement of the Event Horizon Telescope (EHT) array, with $\sim 10$ new antennas and instrumental upgrades of existing antennas. The increased $uv$-coverage, sensitivity, and frequency coverage allow a wide range of new science opportunities to be explored. The ngEHT Analysis Challenges have been launched to inform develo… ▽ More The next-generation Event Horizon Telescope (ngEHT) will be a significant enhancement of the Event Horizon Telescope (EHT) array, with $\sim 10$ new antennas and instrumental upgrades of existing antennas. The increased $uv$-coverage, sensitivity, and frequency coverage allow a wide range of new science opportunities to be explored. The ngEHT Analysis Challenges have been launched to inform development of the ngEHT array design, science objectives, and analysis pathways. For each challenge, synthetic EHT and ngEHT datasets are generated from theoretical source models and released to the challenge participants, who analyze the datasets using image reconstruction and other methods. The submitted analysis results are evaluated with quantitative metrics. In this work, we report on the first two ngEHT Analysis Challenges. These have focused on static and dynamical models of M87* and Sgr A*, and shown that high-quality movies of the extended jet structure of M87* and near-horizon hourly timescale variability of Sgr A* can be reconstructed by the reference ngEHT array in realistic observing conditions, using current analysis algorithms. We identify areas where there is still room for improvement of these algorithms and analysis strategies. Other science cases and arrays will be explored in future challenges. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: 32 pages, 14 figures, accepted for publication in Galaxies

Fast-Cadence High-Contrast Imaging with Information Field Theory

Authors: Jakob Roth, Gianluca Li Causi, Vincenzo Testa, Philipp Arras, Torsten A. Enßlin

Abstract: Although many exoplanets have been indirectly detected over the last years, direct imaging of them with ground-based telescopes remains challenging. In the presence of atmospheric fluctuations, it is ambitious to resolve the high brightness contrasts at the small angular separation between the star and its potential partners. Post-processing of telescope images has become an essential tool to impr… ▽ More Although many exoplanets have been indirectly detected over the last years, direct imaging of them with ground-based telescopes remains challenging. In the presence of atmospheric fluctuations, it is ambitious to resolve the high brightness contrasts at the small angular separation between the star and its potential partners. Post-processing of telescope images has become an essential tool to improve the resolvable contrast ratios. This paper contributes a post-processing algorithm for fast-cadence imaging, which deconvolves sequences of telescope images. The algorithm infers a Bayesian estimate of the astronomical object as well as the atmospheric optical path length, including its spatial and temporal structures. For this, we utilize physics-inspired models for the object, the atmosphere, and the telescope. The algorithm is computationally expensive but allows to resolve high contrast ratios despite short observation times and no field rotation. We test the performance of the algorithm with point-like companions synthetically injected into a real data set acquired with the SHARK-VIS pathfinder instrument at the LBT telescope. Sources with brightness ratios down to $6\cdot10^{-4}$ to the star are detected at $185$ mas separation with a short observation time of $0.6\,\text{s}$. △ Less

Submitted 7 February, 2023; v1 submitted 15 December, 2022; originally announced December 2022.

Comments: 12 pages, 6 figures

Journal ref: AJ 165 86 (2023)

Accretion Flow Morphology in Numerical Simulations of Black Holes from the ngEHT Model Library: The Impact of Radiation Physics

Authors: Koushik Chatterjee, Andrew Chael, Paul Tiede, Yosuke Mizuno, Razieh Emami, Christian Fromm, Angelo Ricarte, Lindy Blackburn, Freek Roelofs, Michael D. Johnson, Sheperd S. Doeleman, Philipp Arras, Antonio Fuentes, Jakob Knollmüller, Nikita Kosogorov, Greg Lindahl, Hendrik Müller, Nimesh Patel, Alexander Raymond, Efthalia Traianou, Justin Vega

Abstract: In the past few years, the Event Horizon Telescope (EHT) has provided the first-ever event horizon-scale images of the supermassive black holes (BHs) (M87*) and Sagittarius A$^*$ (Sgr A*). The next-generation EHT project is an extension of the EHT array that promises larger angular resolution and higher sensitivity to the dim, extended flux around the central ring-like structure, possibly connecti… ▽ More In the past few years, the Event Horizon Telescope (EHT) has provided the first-ever event horizon-scale images of the supermassive black holes (BHs) (M87*) and Sagittarius A$^*$ (Sgr A*). The next-generation EHT project is an extension of the EHT array that promises larger angular resolution and higher sensitivity to the dim, extended flux around the central ring-like structure, possibly connecting the accretion flow and the jet. The ngEHT Analysis Challenges aim to understand the science extractability from synthetic images and movies to inform the ngEHT array design and analysis algorithm development. In this work, we compare the accretion flow structure and dynamics in numerical fluid simulations that specifically target M87* and Sgr A*, and were used to construct the source models in the challenge set. We consider (1) a steady-state axisymmetric radiatively inefficient accretion flow model with a time-dependent shearing hotspot, (2) two time-dependent single fluid general relativistic magnetohydrodynamic (GRMHD) simulations from the H-AMR code, (3) a two-temperature GRMHD simulation from the BHAC code, and (4) a two-temperature radiative GRMHD simulation from the KORAL code. We find that the different models exhibit remarkably similar temporal and spatial properties, except for the electron temperature, since radiative losses substantially cool down electrons near the BH and the jet sheath, signaling the importance of radiative cooling even for slowly accreting BHs such as M87*. We restrict ourselves to standard torus accretion flows, and leave larger explorations of alternate accretion models to future work. △ Less

Submitted 7 March, 2023; v1 submitted 4 December, 2022; originally announced December 2022.

Comments: Accepted in Galaxies; 23 pages, 7 figures

Beyond the linear tide: impact of the non-linear tidal response of neutron stars on gravitational waveforms from binary inspirals

Authors: Hang Yu, Nevin N. Weinberg, Phil Arras, James Kwon, Tejaswi Venumadhav

Abstract: Tidal interactions in coalescing binary neutron stars modify the dynamics of the inspiral and hence imprint a signature on their gravitational wave (GW) signals in the form of an extra phase shift. We need accurate models for the tidal phase shift in order to constrain the supranuclear equation of state from observations. In previous studies, GW waveform models were typically constructed by treati… ▽ More Tidal interactions in coalescing binary neutron stars modify the dynamics of the inspiral and hence imprint a signature on their gravitational wave (GW) signals in the form of an extra phase shift. We need accurate models for the tidal phase shift in order to constrain the supranuclear equation of state from observations. In previous studies, GW waveform models were typically constructed by treating the tide as a linear response to a perturbing tidal field. In this work, we incorporate non-linear corrections due to hydrodynamic three- and four-mode interactions and show how they can improve the accuracy and explanatory power of waveform models. We set up and numerically solve the coupled differential equations for the orbit and the modes and analytically derive solutions of the system's equilibrium configuration. Our analytical solutions agree well with the numerical ones up to the merger and involve only algebraic relations, allowing for fast phase shift and waveform evaluations for different equations of state over a large parameter space. We find that, at Newtonian order, non-linear fluid effects can enhance the tidal phase shift by $\gtrsim 1\,{\rm radian}$ at a GW frequency of 1000 Hz, corresponding to a $10-20\%$ correction to the linear theory. The scale of the additional phase shift near the merger is consistent with the difference between numerical relativity and theoretical predictions that account only for the linear tide. Non-linear fluid effects are thus important when interpreting the results of numerical relativity and in the construction of waveform models for current and future GW detectors. △ Less

Submitted 6 January, 2023; v1 submitted 13 November, 2022; originally announced November 2022.

Comments: 20 pages, 9 figures. Accepted for publication in MNRAS

Bayesian statistics approach to imaging of aperture synthesis data: RESOLVE meets ALMA

Authors: Łukasz Tychoniec, Fabrizia Guglielmetti, Philipp Arras, Torsten Enßlin, Eric Villard

Abstract: The Atacama Large Millimeter/submillimeter Array (ALMA) is currently revolutionizing observational astrophysics. The aperture synthesis technique provides angular resolution otherwise unachievable with the conventional single-aperture telescope. However, recovering the image from the inherently undersampled data is a challenging task. The CLEAN algorithm has proven successful and reliable and is c… ▽ More The Atacama Large Millimeter/submillimeter Array (ALMA) is currently revolutionizing observational astrophysics. The aperture synthesis technique provides angular resolution otherwise unachievable with the conventional single-aperture telescope. However, recovering the image from the inherently undersampled data is a challenging task. The CLEAN algorithm has proven successful and reliable and is commonly used in imaging the interferometric observations. It is not, however, free of limitations. Point-source assumption, central to the CLEAN is not optimal for the extended structures of molecular gas recovered by ALMA. Additionally, negative fluxes recovered with CLEAN are not physical. This begs to search for alternatives that would be better suited for specific science cases. We present the recent developments in imaging ALMA data using Bayesian inference techniques, namely the RESOLVE algorithm This algorithm, based on information field theory \cite{Ensslin2013}, has been already successfully applied to image the Very Large Array data. We compare the capability of both CLEAN and RESOLVE to recover known sky signal, convoluted with the simulator of ALMA observation data and we investigate the problem with a set of actual ALMA observations. △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: Proceedings of International Workshop on Bayesian Inference and Maximum Entropy Methods in Science and Engineering, IHP, Paris, July 18-22, 2022

Bayesian and Machine Learning Methods in the Big Data era for astronomical imaging

Authors: Fabrizia Guglielmetti, Philipp Arras, Michele Delli Veneri, Torsten Enßlin, Giuseppe Longo, Łukasz Tychoniec, Eric Villard

Abstract: The Atacama Large Millimeter/submillimeter Array with the planned electronic upgrades will deliver an unprecedented amount of deep and high resolution observations. Wider fields of view are possible with the consequential cost of image reconstruction. Alternatives to commonly used applications in image processing have to be sought and tested. Advanced image reconstruction methods are critical to m… ▽ More The Atacama Large Millimeter/submillimeter Array with the planned electronic upgrades will deliver an unprecedented amount of deep and high resolution observations. Wider fields of view are possible with the consequential cost of image reconstruction. Alternatives to commonly used applications in image processing have to be sought and tested. Advanced image reconstruction methods are critical to meet the data requirements needed for operational purposes. Astrostatistics and astroinformatics techniques are employed. Evidence is given that these interdisciplinary fields of study applied to synthesis imaging meet the Big Data challenges and have the potentials to enable new scientific discoveries in radio astronomy and astrophysics. △ Less

Submitted 4 October, 2022; originally announced October 2022.

Comments: 8 pages, 5 figures, proceedings International Workshop on Bayesian Inference and Maximum Entropy Methods in Science and Engineering, IHP, Paris, July 18-22, 2022

Periastron precession effect of $f$-mode dynamical tides on gravitational waves from eccentric double white dwarfs

Authors: Shu Yan Lau, Kent Yagi, Phil Arras

Abstract: The dynamical tide can play an important role in the orbital motion of close eccentric double white dwarf binaries. As the launching of the space-based gravitational-wave detector, the Laser Interferometer Space Antenna (LISA), is just around the corner, detection of gravitational wave signals from such systems is anticipated. In this paper, we discuss the influence of the dynamical tide on eccent… ▽ More The dynamical tide can play an important role in the orbital motion of close eccentric double white dwarf binaries. As the launching of the space-based gravitational-wave detector, the Laser Interferometer Space Antenna (LISA), is just around the corner, detection of gravitational wave signals from such systems is anticipated. In this paper, we discuss the influence of the dynamical tide on eccentric orbits, focusing on the effect on orbital precession. We show that in orbits with a high eccentricity, resonance can cause a large precession when a harmonic of the orbital frequency matches the natural frequencies of the normal modes of the star. In contrast to the case with circular orbits, each mode can encounter multiple resonances with different harmonics and these resonant regions can cover about 10% of the frequency space for orbits with close separations. In this case, the tidal precession effect is distinct from the other contributions and can be identified with LISA if the signal-to-noise ratio is high enough. However, within the highly eccentric-small separation region, the dynamical tide causes chaotic motion and the gravitational wave signal becomes unpredictable. Even not at resonance, the dynamical tide can contribute up to 20% of the precession for orbits close to Roche-lobe filling separation with low eccentricities and LISA can resolve these off-resonant dynamical tide effects within the low eccentricity-small orbital separation region of the parameter space. For lower mass systems, the dynamical tide effect can degenerate with the uncertainties of the eccentricity, making it unmeasurable from the precession rate alone. For higher mass systems, the radiation reaction effect becomes significant enough to constrain the eccentricity, allowing the measurement of the dynamical tide. △ Less

Submitted 23 May, 2024; v1 submitted 14 September, 2022; originally announced September 2022.

Comments: 13 pages. 9 figures

Journal ref: Phys. Rev. D 106, 103038 (2022)

Dust Accretion onto Exoplanets

Authors: Phil Arras, Megan Wilson, Matthew Pryal, Jordan Baker

Abstract: Accretion of interplanetary dust onto gas giant exoplanets is considered. Poynting-Robertson drag causes dust particles from distant reservoirs to slowly inspiral toward the star. Orbital simulations for the three-body system of the star, planet, and dust particle show that a significant fraction of the dust may accrete onto massive planets in close orbits. The deceleration of the supersonic dust… ▽ More Accretion of interplanetary dust onto gas giant exoplanets is considered. Poynting-Robertson drag causes dust particles from distant reservoirs to slowly inspiral toward the star. Orbital simulations for the three-body system of the star, planet, and dust particle show that a significant fraction of the dust may accrete onto massive planets in close orbits. The deceleration of the supersonic dust in the planet's atmosphere is modeled, including ablation by thermal evaporation and sputtering. The fraction of the accreted dust mass deposited as gas-phase atoms is found to be large for close-in orbits and massive planets. If mass outflow and vertical mixing are sufficiently weak, the accreted dust produces a constant mixing ratio of atoms and remnant dust grains below the stopping layer. When vertical mixing is included along with settling, the solutions interpolate between the mixing ratio due to the meteoric source above the homopause, and that of the well-mixed deeper atmosphere below the homopause. The line opacity from atoms and continuum opacity from remnant dust may be observable in transmission spectra for sufficiently large dust accretion rates, a grain size distribution tilted toward the blowout size, and sufficiently weak vertical mixing. If mixing is strong, the meteoric source may still act to augment heavy elements mixed up from the deep atmosphere as well as provide nucleation sites for the formation of larger particles. The possible role of the Lorentz drag force in limiting the flow speeds and mixing coefficient for pressures $\la 1\, \rm mbar$ is discussed. △ Less

Submitted 17 June, 2022; originally announced June 2022.

Comments: 28 pages, 23 figures

Journal ref: The Astrophysical Journal, 932, 2022

A Novel Solution for Resonant Scattering Using Self-Consistent Boundary Conditions

Authors: B. Connor McClellan, Shane Davis, Phil Arras

Abstract: We present two novel additions to the semi-analytic solution of Lyman $α$ (Ly$α$) radiative transfer in spherical geometry: (1) implementation of the correct boundary condition for a steady source, and (2) solution of the time-dependent problem for an impulsive source. For the steady-state problem, the solution can be represented as a sum of two terms: a previously-known analytic solution of the e… ▽ More We present two novel additions to the semi-analytic solution of Lyman $α$ (Ly$α$) radiative transfer in spherical geometry: (1) implementation of the correct boundary condition for a steady source, and (2) solution of the time-dependent problem for an impulsive source. For the steady-state problem, the solution can be represented as a sum of two terms: a previously-known analytic solution of the equation with mean intensity $J=0$ at the surface, and a novel, semi-analytic solution which enforces the correct boundary condition of zero-ingoing intensity at the surface. This solution is compared to that of the Monte Carlo method, which is valid at arbitrary optical depth. It is shown that the size of the correction is of order unity when the spectral peaks approach the Doppler core and decreases slowly with line center optical depth, specifically as $(a τ_0)^{-1/3}$, which may explain discrepancies seen in previous studies. For the impulsive problem, the time, spatial, and frequency dependence of the solution are expressed using an eigenfunction expansion in order to characterize the escape time distribution and emergent spectra of photons. It is shown that the lowest-order eigenfrequency agrees well with the decay rate found in the Monte Carlo escape time distribution at sufficiently large line-center optical depths. The characterization of the escape-time distribution highlights the potential for a Monte Carlo acceleration method, which would sample photon escape properties from distributions rather than calculating every photon scattering, thereby reducing computational demand. △ Less

Submitted 10 May, 2022; originally announced May 2022.

Comments: 23 pages, 10 figures, 2 appendices, resubmitted to ApJ after review

Multi-Component Imaging of the Fermi Gamma-ray Sky in the Spatio-spectral Domain

Authors: Lukas I. Platz, Jakob Knollmüller, Philipp Arras, Philipp Frank, Martin Reinecke, Dominik Jüstel, Torsten A. Enßlin

Abstract: We perform two distinct spatio-spectral reconstructions of the gamma-ray sky in the range of 0.56-316 GeV based on Fermi Large Area Telescope (LAT) data. Both describe the sky brightness to be composed of a diffuse-emission and a point-source component. The first model requires minimal assumptions and provides a template-free reconstruction as a reference. It makes use of spatial and spectral corr… ▽ More We perform two distinct spatio-spectral reconstructions of the gamma-ray sky in the range of 0.56-316 GeV based on Fermi Large Area Telescope (LAT) data. Both describe the sky brightness to be composed of a diffuse-emission and a point-source component. The first model requires minimal assumptions and provides a template-free reconstruction as a reference. It makes use of spatial and spectral correlations to distinguish between the different components. The second model is physics-informed and further differentiates between diffuse emission of hadronic and leptonic origin. For this, we assume parametric, but spatially varying energy spectra to distinguish between the processes and use thermal Galactic dust observations to indicate the preferred sites of hadronic interactions. To account for instrumental effects we model the point-spread, the energy dispersion, and the exposure of the telescope throughout the observation. The reconstruction problem is formulated as a Bayesian inference task, that is solved by variational inference. We show decompositions of the Gamma-ray flux into diffuse and point-like emissions, and of the diffuse emissions into multiple physically motivated components. The diffuse decomposition provides an unprecedented view of the Galactic leptonic diffuse emission. It shows the Fermi bubbles and their spectral variations in high fidelity and other areas exhibiting strong cosmic ray electron contents, such as a thick disk in the inner Galaxy and outflow regions. Furthermore, we report a hard spectrum gamma ray arc in the northern outer bubble co-spatial with the reported X-ray arc by the eROSITA collaboration. All our spatio-spectral sky reconstructions and their uncertainty quantification are publicly available. △ Less

Submitted 6 December, 2023; v1 submitted 20 April, 2022; originally announced April 2022.

Comments: Please see the journal version for a language-edited release of this manuscript (open-access): https://doi.org/10.1051/0004-6361/202243819

Journal ref: A&A 680, A2 (2023)

Deep Images of the Galactic Center with GRAVITY

Authors: GRAVITY Collaboration, R. Abuter, N. Aimar, A. Amorim, P. Arras, M. Bauböck, J. P. Berger, H. Bonnet, W. Brandner, G. Bourdarot, V. Cardoso, Y. Clénet, R. Davies, P. T. de Zeeuw, J. Dexter, Y. Dallilar, A. Drescher, F. Eisenhauer, T. Enßlin, N. M. Förster Schreiber, P. Garcia, F. Gao, E. Gendron, R. Genzel, S. Gillessen , et al. (43 additional authors not shown)

Abstract: Stellar orbits at the Galactic Center provide a very clean probe of the gravitational potential of the supermassive black hole. They can be studied with unique precision, beyond the confusion limit of a single telescope, with the near-infrared interferometer GRAVITY. Imaging is essential to search the field for faint, unknown stars on short orbits which potentially could constrain the black hole s… ▽ More Stellar orbits at the Galactic Center provide a very clean probe of the gravitational potential of the supermassive black hole. They can be studied with unique precision, beyond the confusion limit of a single telescope, with the near-infrared interferometer GRAVITY. Imaging is essential to search the field for faint, unknown stars on short orbits which potentially could constrain the black hole spin. Furthermore, it provides the starting point for astrometric fitting to derive highly accurate stellar positions. Here, we present $\mathrm{G^R}$, a new imaging tool specifically designed for Galactic Center observations with GRAVITY. The algorithm is based on a Bayesian interpretation of the imaging problem, formulated in the framework of information field theory and building upon existing works in radio-interferometric imaging. Its application to GRAVITY observations from 2021 yields the deepest images to date of the Galactic Center on scales of a few milliarcseconds. The images reveal the complicated source structure within the central $100\,\mathrm{mas}$ around Sgr A*, where we detected the stars S29 and S55 and confirm S62 on its trajectory, slowly approaching Sgr A*. Furthermore, we were able to detect S38, S42, S60, and S63 in a series of exposures for which we offset the fiber from Sgr A*. We provide an update on the orbits of all aforementioned stars. In addition to these known sources, the images also reveal a faint star moving to the west at a high angular velocity. We cannot find any coincidence with any known source and, thus, we refer to the new star as S300. From the flux ratio with S29, we estimate its K-band magnitude as $m_\mathrm{K}\left(\mathrm{S300}\right)\simeq 19.0 - 19.3$. Images obtained with CLEAN confirm the detection. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: 24 pages, 14 figures

Journal ref: A&A 657, A82 (2022)

Tidal evolution and diffusive growth during high-eccentricity planet migration: revisiting the eccentricity distribution of hot Jupiters

Authors: Hang Yu, Nevin N. Weinberg, Phil Arras

Abstract: High-eccentricity tidal migration is a potential formation channel for hot Jupiters. During this process, the planetary f-mode may experience a phase of diffusive growth, allowing its energy to quickly build up to large values. In Yu et al. (2021, ApJ, 917, 31), we demonstrated that nonlinear mode interactions between a parent f-mode and daughter f- and p-modes expand the parameter space over whic… ▽ More High-eccentricity tidal migration is a potential formation channel for hot Jupiters. During this process, the planetary f-mode may experience a phase of diffusive growth, allowing its energy to quickly build up to large values. In Yu et al. (2021, ApJ, 917, 31), we demonstrated that nonlinear mode interactions between a parent f-mode and daughter f- and p-modes expand the parameter space over which the diffusive growth of the parent is triggered. We extend that study by incorporating (1) the angular momentum transfer between the orbit and the mode, and consequently the evolution of the pericenter distance, (2) a phenomenological correction to the nonlinear frequency shift at high parent mode energies, and (3) dissipation of the parent's energy due to both turbulent convective damping of the daughter modes and strongly nonlinear wave-breaking events. The new ingredients allow us to follow the coupled evolution of the mode and orbit over $\gtrsim 10^4$ years, covering the diffusive evolution from its onset to its termination. We find that the semi-major axis shrinks by a factor of nearly ten over $10^4$ years, corresponding to a tidal quality factor $\mathcal{Q}\sim10$. The f-mode's diffusive growth terminates while the eccentricity is still high, at around e=0.8-0.95. Using these results, we revisit the eccentricity distribution of proto-hot Jupiters. We estimate that less than 1 proto-HJ with eccentricity >0.9 should be expected in Kepler's data once the diffusive regime is accounted for, explaining the observed paucity of this population. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: 18 pages, 9 figures, submitted to ApJ

Damping of Oscillations in Red Giants by Resonant Mode Coupling

Authors: Nevin N. Weinberg, Phil Arras, Debaditya Pramanik

Abstract: Asteroseismic studies of red giants generally assume that the oscillation modes can be treated as linear perturbations to the background star. However, observations by the Kepler mission show that the oscillation amplitudes increase dramatically as stars ascend the red giant branch. The importance of nonlinear effects should therefore be assessed. In previous work, we found that mixed modes in red… ▽ More Asteroseismic studies of red giants generally assume that the oscillation modes can be treated as linear perturbations to the background star. However, observations by the Kepler mission show that the oscillation amplitudes increase dramatically as stars ascend the red giant branch. The importance of nonlinear effects should therefore be assessed. In previous work, we found that mixed modes in red giants are unstable to nonlinear three-wave interactions over a broad range of stellar mass and evolutionary state. Here we solve the amplitude equations that describe the mode dynamics for large networks of nonlinearly coupled modes. The networks consist of stochastically driven parent modes coupled to resonant secondary modes (daughters, granddaughters, etc.). We find that nonlinear interactions can lower the energy of gravity-dominated mixed modes by $\gtrsim 80\%$ compared to linear theory. However, they have only a mild influence on the energy of pressure-dominated mixed modes. Expressed in terms of the dipole mode visibility $V^2$, i.e., the summed amplitudes of dipole modes relative to radial modes, we find that $V^2$ can be suppressed by $50-80\%$ relative to the linear value for highly-evolved red giants whose frequency of maximum power $ν_{\rm max} \lesssim 100\,μ\textrm{Hz}$. However, for less evolved red giants with $150\lesssim ν_{\rm max} \lesssim 200\,μ\textrm{Hz}$, $V^2$ is suppressed by only $10-20\%$. We conclude that resonant mode coupling can have a potentially detectable effect on oscillations at $ν_{\rm max} \lesssim 100\,μ\textrm{Hz}$ but it cannot account for the population of red giants that exhibit dipole modes with unusually small amplitudes at high $ν_{\rm max}$. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Comments: Accepted for publication in ApJ; 19 pages, 9 figures

Tides in the high-eccentricity migration of hot Jupiters: Triggering diffusive growth by nonlinear mode interactions

Authors: Hang Yu, Nevin N. Weinberg, Phil Arras

Abstract: High eccentricity migration is a possible formation channel for hot Jupiters. However, in order for it to be consistent with the observed population of planets, tides must circularize the orbits in less than $\approx$ a Myr. A potential mechanism for such rapid circularization is the diffusive growth of the tidally driven planetary f-mode. Such growth occurs if the f-mode's phase at pericenter var… ▽ More High eccentricity migration is a possible formation channel for hot Jupiters. However, in order for it to be consistent with the observed population of planets, tides must circularize the orbits in less than $\approx$ a Myr. A potential mechanism for such rapid circularization is the diffusive growth of the tidally driven planetary f-mode. Such growth occurs if the f-mode's phase at pericenter varies chaotically from one pericenter passage to the next. Previous studies focused on the variation of the orbital period due to tidal back-reaction on the orbit as the source of chaos. Here we show that nonlinear mode interactions can also be an important source. Specifically, we show that nonlinear interactions between a parent f-mode and daughter f-/p-modes induce an energy-dependent shift in the oscillation frequency of the parent. This frequency shift varies randomly from orbit to orbit because the parent's energy varies. As a result, the parent's phase at pericenter varies randomly, which we find can trigger it to grow diffusively. We show that the phase shift induced by nonlinear mode interactions in fact dominates the shift induced by tidal back-reaction and significantly lowers the one-kick energy threshold for diffusive growth by about a factor of 5 compared to the linear theory's prediction. Nonlinear interactions could thus enhance the formation rate of hot Jupiters through the high-eccentricity migration channel and potentially mitigate the discrepancy between the observed and predicted occurrence rates for close-in gas giants as compared to those further from the star. △ Less

Submitted 11 April, 2021; originally announced April 2021.

Comments: 19 pages, 5 figures; submitted to ApJ

Exploring the Galactic Warp Through Asymmetries in the Kinematics of the Galactic Disk

Authors: Xinlun Cheng, Borja Anguiano, Steven R. Majewski, Christian Hayes, Phil Arras, Cristina Chiappini, Sten Hasselquist, Anna Bárbara de Andrade Queiroz, Christian Nitschelm, Domingo Anıbal Garcıa-Hernández, Richard R. Lane, Alexandre Roman-Lopes, Peter Frinchaboy

Abstract: Previous analyses of large databases of Milky Way stars have revealed the stellar disk of our Galaxy to be warped and that this imparts a strong signature on the kinematics of stars beyond the solar neighborhood. However, due to the limitation of accurate distance estimates, many attempts to explore the extent of these Galactic features have generally been restricted to a volume near the Sun. By c… ▽ More Previous analyses of large databases of Milky Way stars have revealed the stellar disk of our Galaxy to be warped and that this imparts a strong signature on the kinematics of stars beyond the solar neighborhood. However, due to the limitation of accurate distance estimates, many attempts to explore the extent of these Galactic features have generally been restricted to a volume near the Sun. By combining Gaia DR2 astrometric solution, StarHorse distance and stellar abundances from the APOGEE survey, we present the most detailed and radially expansive study yet of the vertical and radial motions of stars in the Galactic disk. We map stellar velocity with respect to their Galactocentric radius, angular momentum, and azimuthal angle and assess their relation to the warp. A decrease in vertical velocity is discovered at Galactocentric radius $R=13\ \text{kpc}$ and angular momentum $L_z=2800\ \text{kpc}\ \text{km}\ \text{s}^{-1}$. Smaller ripples in vertical and radial velocity are also discovered superposed on the main trend. We also discovered that trends in the vertical velocity with azimuthal angle are not symmetric about the peak, suggesting the warp to be lopsided. To explain the global trend in vertical velocity, we built a simple analytical model of the Galactic warp. Our best fit yields a starting radius of $8.87^{+0.08}_{-0.09}\ \text{kpc}$ and precession rate of $13.57^{+0.20}_{-0.18}\ \text{km}\ \text{s}^{-1}\ \text{kpc}^{-1}$. These parameters remain consistent across stellar age groups, a result that supports the notion that the warp is the result of an external, gravitationally induced phenomenon. △ Less

Submitted 20 October, 2020; originally announced October 2020.

Comments: Accepted for publication in ApJ

Towards Bayesian Data Compression

Authors: Johannes Harth-Kitzerow, Reimar Leike, Philipp Arras, Torsten A. Enßlin

Abstract: In order to handle large data sets omnipresent in modern science, efficient compression algorithms are necessary. Here, a Bayesian data compression (BDC) algorithm that adapts to the specific measurement situation is derived in the context of signal reconstruction. BDC compresses a data set under conservation of its posterior structure with minimal information loss given the prior knowledge on the… ▽ More In order to handle large data sets omnipresent in modern science, efficient compression algorithms are necessary. Here, a Bayesian data compression (BDC) algorithm that adapts to the specific measurement situation is derived in the context of signal reconstruction. BDC compresses a data set under conservation of its posterior structure with minimal information loss given the prior knowledge on the signal, the quantity of interest. Its basic form is valid for Gaussian priors and likelihoods. For constant noise standard deviation, basic BDC becomes equivalent to a Bayesian analog of principal component analysis. Using Metric Gaussian Variational Inference, BDC generalizes to non-linear settings. In its current form, BDC requires the storage of effective instrument response functions for the compressed data and corresponding noise encoding the posterior covariance structure. Their memory demand counteract the compression gain. In order to improve this, sparsity of the compressed responses can be obtained by separating the data into patches and compressing them separately. The applicability of BDC is demonstrated by applying it to synthetic data and radio astronomical data. Still the algorithm needs further improvement as the computation time of the compression and subsequent inference exceeds the time of the inference with the original data. △ Less

Submitted 29 December, 2020; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 39 pages, 15 figures, 1 table, for code, see https://gitlab.mpcdf.mpg.de/jharthki/bdc

Journal ref: ANNALEN DER PHYSIK 2021, 2000508

Efficient wide-field radio interferometry response

Authors: Philipp Arras, Martin Reinecke, Rüdiger Westermann, Torsten A. Enßlin

Abstract: Radio interferometers do not measure the sky brightness distribution directly but rather a modified Fourier transform of it. Imaging algorithms, thus, need a computational representation of the linear measurement operator and its adjoint, irrespective of the specific chosen imaging algorithm. In this paper, we present a C++ implementation of the radio interferometric measurement operator for wide-… ▽ More Radio interferometers do not measure the sky brightness distribution directly but rather a modified Fourier transform of it. Imaging algorithms, thus, need a computational representation of the linear measurement operator and its adjoint, irrespective of the specific chosen imaging algorithm. In this paper, we present a C++ implementation of the radio interferometric measurement operator for wide-field measurements which is based on "improved $w$-stacking". It can provide high accuracy (down to $\approx 10^{-12}$), is based on a new gridding kernel which allows smaller kernel support for given accuracy, dynamically chooses kernel, kernel support and oversampling factor for maximum performance, uses piece-wise polynomial approximation for cheap evaluations of the gridding kernel, treats the visibilities in cache-friendly order, uses explicit vectorisation if available and comes with a parallelisation scheme which scales well also in the adjoint direction (which is a problem for many previous implementations). The implementation has a small memory footprint in the sense that temporary internal data structures are much smaller than the respective input and output data, allowing in-memory processing of data sets which needed to be read from disk or distributed across several compute nodes before. △ Less

Submitted 28 January, 2021; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 13 pages, 8 figures

Journal ref: A&A 646, A58 (2021)

Comparison of classical and Bayesian imaging in radio interferometry

Authors: Philipp Arras, Hertzog L. Bester, Richard A. Perley, Reimar Leike, Oleg Smirnov, Rüdiger Westermann, Torsten A. Enßlin

Abstract: CLEAN, the commonly employed imaging algorithm in radio interferometry, suffers from a number of shortcomings: in its basic version it does not have the concept of diffuse flux, and the common practice of convolving the CLEAN components with the CLEAN beam erases the potential for super-resolution; it does not output uncertainty information; it produces images with unphysical negative flux regions… ▽ More CLEAN, the commonly employed imaging algorithm in radio interferometry, suffers from a number of shortcomings: in its basic version it does not have the concept of diffuse flux, and the common practice of convolving the CLEAN components with the CLEAN beam erases the potential for super-resolution; it does not output uncertainty information; it produces images with unphysical negative flux regions; and its results are highly dependent on the so-called weighting scheme as well as on any human choice of CLEAN masks to guiding the imaging. Here, we present the Bayesian imaging algorithm resolve which solves the above problems and naturally leads to super-resolution. We take a VLA observation of Cygnus~A at four different frequencies and image it with single-scale CLEAN, multi-scale CLEAN and resolve. Alongside the sky brightness distribution resolve estimates a baseline-dependent correction function for the noise budget, the Bayesian equivalent of weighting schemes. We report noise correction factors between 0.4 and 429. The enhancements achieved by resolve come at the cost of higher computational effort. △ Less

Submitted 25 January, 2021; v1 submitted 26 August, 2020; originally announced August 2020.

Comments: 23 pages, 16 figures, 4 tables, data published at https://doi.org/10.5281/zenodo.4267057

Journal ref: A&A 646, A84 (2021)

Geometry of the Draco C1 Symbiotic Binary

Authors: Hannah M. Lewis, Borja Anguiano, Keivan G. Stassun, Steven R. Majewski, Phil Arras, Craig L. Sarazin, Zhi-Yun Li, Nathan De Lee, Nicholas W. Troup, Carlos Allende Prieto, Carles Badenes, Katia Cunha, D. A. Garcia-Hernandez, David L. Nidever, Pedro A. Palicio, Joshua D. Simon, Verne V. Smith

Abstract: Draco C1 is a known symbiotic binary star system composed of a carbon red giant and a hot, compact companion -- likely a white dwarf -- belonging to the Draco dwarf spheroidal galaxy. From near-infrared spectroscopic observations taken by the Apache Point Observatory Galactic Evolution Experiment (APOGEE-2), part of Sloan Digital Sky Survey IV, we provide updated stellar parameters for the cool, g… ▽ More Draco C1 is a known symbiotic binary star system composed of a carbon red giant and a hot, compact companion -- likely a white dwarf -- belonging to the Draco dwarf spheroidal galaxy. From near-infrared spectroscopic observations taken by the Apache Point Observatory Galactic Evolution Experiment (APOGEE-2), part of Sloan Digital Sky Survey IV, we provide updated stellar parameters for the cool, giant component, and constrain the temperature and mass of the hot, compact companion. Prior measurements of the periodicity of the system, based on only a few epochs of radial velocity data or relatively short baseline photometric observations, were sufficient only to place lower limits on the orbital period ($P > 300$ days). For the first time, we report precise orbital parameters for the binary system: With 43 radial velocity measurements from APOGEE spanning an observational baseline of more than 3 years, we definitively derive the period of the system to be $1220.0^{+3.7}_{-3.5}$ days. Based on the newly derived orbital period and separation of the system, together with estimates of the radius of the red giant star, we find that the hot companion must be accreting matter from the dense wind of its evolved companion. △ Less

Submitted 18 September, 2020; v1 submitted 13 August, 2020; originally announced August 2020.

Comments: 8 pages, 4 figures, 1 table. Accepted for publication in ApJL

Atmospheric Escape Processes and Planetary Atmospheric Evolution

Authors: Guillaume Gronoff, Phil Arras, Suleiman M Baraka, Jared M Bell, Gaël Cessateur, Ofer Cohen, Shannon M. Curry, Jeremy J Drake, Meredith K Elrod, Justin T. Erwin, Katherine Garcia-Sage, Cecilia Garraffo, Alex Glocer, Nicholas Gray Heavens, Kylie Lovato, Romain Maggiolo, Christopher D. Parkinson, Cyril L. Simon Wedlund, Daniel R Weimer, William B. Moore

Abstract: The habitability of the surface of any planet is determined by a complex evolution of its interior, surface, and atmosphere. The electromagnetic and particle radiation of stars drive thermal, chemical and physical alteration of planetary atmospheres, including escape. Many known extrasolar planets experience vastly different stellar environments than those in our Solar system: it is crucial to und… ▽ More The habitability of the surface of any planet is determined by a complex evolution of its interior, surface, and atmosphere. The electromagnetic and particle radiation of stars drive thermal, chemical and physical alteration of planetary atmospheres, including escape. Many known extrasolar planets experience vastly different stellar environments than those in our Solar system: it is crucial to understand the broad range of processes that lead to atmospheric escape and evolution under a wide range of conditions if we are to assess the habitability of worlds around other stars. One problem encountered between the planetary and the astrophysics communities is a lack of common language for describing escape processes. Each community has customary approximations that may be questioned by the other, such as the hypothesis of H-dominated thermosphere for astrophysicists, or the Sun-like nature of the stars for planetary scientists. Since exoplanets are becoming one of the main targets for the detection of life, a common set of definitions and hypotheses are required. We review the different escape mechanisms proposed for the evolution of planetary and exoplanetary atmospheres. We propose a common definition for the different escape mechanisms, and we show the important parameters to take into account when evaluating the escape at a planet in time. We show that the paradigm of the magnetic field as an atmospheric shield should be changed and that recent work on the history of Xenon in Earth's atmosphere gives an elegant explanation to its enrichment in heavier isotopes: the so-called Xenon paradox. △ Less

Submitted 6 March, 2020; originally announced March 2020.

Variable structures in M87* from space, time and frequency resolved interferometry

Authors: Philipp Arras, Philipp Frank, Philipp Haim, Jakob Knollmüller, Reimar Leike, Martin Reinecke, Torsten Enßlin

Abstract: Observing the dynamics of compact astrophysical objects provides insights into their inner workings, thereby probing physics under extreme conditions. The immediate vicinity of an active supermassive black hole with its event horizon, photon ring, accretion disk, and relativistic jets is a perfect place to study general relativity and magneto-hydrodynamics. The observations of M87* with Very Long… ▽ More Observing the dynamics of compact astrophysical objects provides insights into their inner workings, thereby probing physics under extreme conditions. The immediate vicinity of an active supermassive black hole with its event horizon, photon ring, accretion disk, and relativistic jets is a perfect place to study general relativity and magneto-hydrodynamics. The observations of M87* with Very Long Baseline Interferometry (VLBI) by the Event Horizon Telescope (EHT) allows to investigate its dynamical processes on time scales of days. Compared to regular radio interferometers, VLBI networks typically have fewer antennas and low signal to noise ratios (SNRs). Furthermore, the source is variable, prohibiting integration over time to improve SNR. Here, we present an imaging algorithm that copes with the data scarcity and temporal evolution, while providing uncertainty quantification. Our algorithm views the imaging task as a Bayesian inference problem of a time-varying brightness, exploits the correlation structure in time, and reconstructs a ${2+1+1}$ dimensional time-variable and spectrally resolved image at once. We apply this method to the EHT observation of M87* and validate our approach on synthetic data. The time- and frequency-resolved reconstruction of M87* confirms variable structures on the emission ring. The reconstruction indicates extended and time-variable emission structures outside the ring itself. △ Less

Submitted 5 June, 2022; v1 submitted 12 February, 2020; originally announced February 2020.

Comments: 32 pages, 17 figures, 6 tables

I-Love-Q Relations for Realistic White Dwarfs

Authors: Andrew Taylor, Kent Yagi, Phil Arras

Abstract: The space-borne gravitational wave interferometer, LISA, is expected to detect signals from numerous binary white dwarfs. At small orbital separation, rapid rotation and large tidal bulges may allow for the stellar internal structure to be probed through such observations. Finite-size effects are encoded in quantities like the moment of inertia ($I$), tidal Love number (Love), and quadrupole momen… ▽ More The space-borne gravitational wave interferometer, LISA, is expected to detect signals from numerous binary white dwarfs. At small orbital separation, rapid rotation and large tidal bulges may allow for the stellar internal structure to be probed through such observations. Finite-size effects are encoded in quantities like the moment of inertia ($I$), tidal Love number (Love), and quadrupole moment ($Q$). The universal relations among them (I-Love-Q relations) can be used to reduce the number of parameters in the gravitational-wave templates. We here study I-Love-Q relations for more realistic white dwarf models than used in previous studies. In particular, we extend previous works by including (i) differential rotation and (ii) internal temperature profiles taken from detailed stellar evolution calculations. We use the publicly available stellar evolution code MESA to generate cooling models of both low- and high-mass white dwarfs. We show that differential rotation causes the I-Q relation (and similarly the Love-Q relation) to deviate from that of constant rotation. We also find that the introduction of finite temperatures causes the white dwarf to move along the zero-temperature mass sequence of I-Q values, moving towards values that suggest a lower mass. We further find that after only a few Myrs, high-mass white dwarfs are well-described by the zero-temperature model, suggesting that the relations with zero-temperature may be good enough in most practical cases. Low-mass, He-core white dwarfs with thick hydrogen envelopes may undergo long periods of H burning which sustain the stellar temperature and allow deviations from the I-Love-Q relations for longer times. △ Less

Submitted 19 December, 2019; originally announced December 2019.

Comments: 15 pages (5 of which are appendices), 11 figures

The full $\mathit{Kepler}$ phase curve of the eclipsing hot white dwarf binary system KOI-964

Authors: Ian Wong, Avi Shporer, Juliette C. Becker, Benjamin J. Fulton, Travis A. Berger, Nevin N. Weinberg, Phil Arras, Andrew W. Howard, Björn Benneke

Abstract: We analyze the full $\mathit{Kepler}$ phase curve of KOI-964, a binary system consisting of a hot white dwarf on an eclipsing orbit around an A-type host star. Using all 18 quarters of long-cadence photometry, we carry out a joint light curve fit and obtain improved phase curve amplitudes, occultation depths, orbital parameters, and transit ephemeris over the previous results of Carter et al. 2011… ▽ More We analyze the full $\mathit{Kepler}$ phase curve of KOI-964, a binary system consisting of a hot white dwarf on an eclipsing orbit around an A-type host star. Using all 18 quarters of long-cadence photometry, we carry out a joint light curve fit and obtain improved phase curve amplitudes, occultation depths, orbital parameters, and transit ephemeris over the previous results of Carter et al. 2011. A periodogram of the residuals from the phase curve fit reveals an additional stellar variability signal from the host star with a characteristic period of $0.620276\pm0.000011$ days and a full amplitude of $24\pm2$ ppm. We also present new Keck/HIRES radial velocity observations which we use to measure the orbit and obtain a mass ratio of $q=0.106\pm0.012$. Combining this measurement with the results of a stellar isochrone analysis, we find that the masses of the host star and white dwarf companion are $2.23\pm0.12\,M_{\odot}$ and $0.236^{+0.028}_{-0.027}\,M_{\odot}$, respectively. The effective temperatures of the two components are $9940^{+260}_{-230}$ K and $15080\pm400$ K, respectively, and we determine the age of the system to be $0.21^{+0.11}_{-0.08}$ Gyr. We use the measured system properties to compute predicted phase curve amplitudes and find that while the measured Doppler boosting and mutual illumination components agree well with theory, the ellipsoidal distortion amplitude is significantly underestimated. We detail possible explanations for this discrepancy, including interactions between the dynamical tide of the host star and the tidal bulge and possible non-synchronous rotation of the host star. △ Less

Submitted 22 November, 2019; v1 submitted 18 October, 2019; originally announced October 2019.

Comments: Accepted for publication in AJ, 18 pages, 5 figures

Unified Radio Interferometric Calibration and Imaging with Joint Uncertainty Quantification

Authors: Philipp Arras, Philipp Frank, Reimar Leike, Rüdiger Westermann, Torsten Enßlin

Abstract: The data reduction procedure for radio interferometers can be viewed as a combined calibration and imaging problem. We present an algorithm that unifies cross-calibration, self-calibration, and imaging. Being a Bayesian method, that algorithm does not only calculate an estimate of the sky brightness distribution, but also provides an estimate of the joint uncertainty which entails both the uncerta… ▽ More The data reduction procedure for radio interferometers can be viewed as a combined calibration and imaging problem. We present an algorithm that unifies cross-calibration, self-calibration, and imaging. Being a Bayesian method, that algorithm does not only calculate an estimate of the sky brightness distribution, but also provides an estimate of the joint uncertainty which entails both the uncertainty of the calibration and the one of the actual observation. The algorithm is formulated in the language of information field theory and uses Metric Gaussian Variational Inference (MGVI) as the underlying statistical method. So far only direction-independent antenna-based calibration is considered. This restriction may be released in future work. An implementation of the algorithm is contributed as well. △ Less

Submitted 20 July, 2019; v1 submitted 26 March, 2019; originally announced March 2019.

Comments: 12 pages, 22 figures

Journal ref: A&A 627, A134 (2019)

Orbital Decay in Binaries Containing Post-Main Sequence Stars

Authors: Meng Sun, Phil Arras, Nevin N. Weinberg, Nicholas W. Troup, Steven R. Majewski

Abstract: The orbital decay of binaries containing a primary sub-giant or red giant star and a stellar or substellar companion is investigated. The tide raised in the primary by the companion leads to an exchange of angular momentum between the orbit and the stellar spin, causing the orbit to contract and the primary to spin up. The rate of orbital decay is computed including both the equilibrium tide, damp… ▽ More The orbital decay of binaries containing a primary sub-giant or red giant star and a stellar or substellar companion is investigated. The tide raised in the primary by the companion leads to an exchange of angular momentum between the orbit and the stellar spin, causing the orbit to contract and the primary to spin up. The rate of orbital decay is computed including both the equilibrium tide, damped by turbulent viscosity in the convective envelope, and the dynamical tide, assumed to be a traveling internal-gravity wave in the radiative core. For close binaries, the tidal forcing period is expected to be much shorter than the eddy turnover timescale in the convective envelope, and the prescription for "reduced" viscosity is an important consideration. The dynamical tide tends to dominate for the closest orbits, while the equilibrium tide dominates for more distant orbits, with the crossover point depending on the stellar mass. The spin up of the primary to synchronous rotation occurs for sufficiently massive secondaries, and this greatly slows the orbital decay until the Darwin instability occurs. A parameter survey is presented for orbital decay as a function of primary and secondary mass, as well as turbulent viscosity prescription. These results are summarized with analytic formulae and numerical results for the age-dependent critical separation, $a_{\rm crit}$, inside of which orbital decay is rapid, and few systems are expected to be observed. The calculations of $a_{\rm crit}$ are compared with APOGEE binaries, as well as solar mass exoplanet host stars. △ Less

Submitted 6 September, 2018; originally announced September 2018.

Comments: 18 pages, 19 figures, accepted to MNRAS

Nonlinear Mixed Modes in Red Giants

Authors: Nevin N. Weinberg, Phil Arras

Abstract: Turbulent motions in the convective envelope of red giants excite a rich spectrum of solar-like oscillation modes. Observations by CoRoT and Kepler have shown that the mode amplitudes increase dramatically as the stars ascend the red giant branch, i.e., as the frequency of maximum power, $ν_\mathrm{max}$, decreases. Most studies nonetheless assume that the modes are well described by the linearize… ▽ More Turbulent motions in the convective envelope of red giants excite a rich spectrum of solar-like oscillation modes. Observations by CoRoT and Kepler have shown that the mode amplitudes increase dramatically as the stars ascend the red giant branch, i.e., as the frequency of maximum power, $ν_\mathrm{max}$, decreases. Most studies nonetheless assume that the modes are well described by the linearized fluid equations. We investigate to what extent the linear approximation is justified as a function of stellar mass $M$ and $ν_\mathrm{max}$, focusing on dipole mixed modes with frequency near $ν_\mathrm{max}$. A useful measure of a mode's nonlinearity is the product of its radial wavenumber and its radial displacement, $k_r ξ_r$ (i.e., its shear). We show that $k_r ξ_r \propto ν_\mathrm{max}^{-9/2}$, implying that the nonlinearity of mixed modes increases significantly as a star evolves. The modes are weakly nonlinear ($k_r ξ_r > 10^{-3}$) for $ν_\mathrm{max} \lesssim 150 \, μ\mathrm{Hz}$ and strongly nonlinear ($k_r ξ_r > 1$) for $ν_\mathrm{max} \lesssim 30 \, μ\mathrm{Hz}$, with only a mild dependence on $M$ over the range we consider ($1.0 - 2.0 M_\odot$). A weakly nonlinear mixed mode can excite secondary waves in the stellar core through the parametric instability, resulting in enhanced, but partial, damping of the mode. By contrast, a strongly nonlinear mode breaks as it propagates through the core and is fully damped there. Evaluating the impact of nonlinear effects on observables such as mode amplitudes and linewidths requires large mode network simulations. We plan to carry out such calculations in the future and investigate whether nonlinear damping can explain why some red giants exhibit dipole modes with unusually small amplitudes, known as depressed modes. △ Less

Submitted 20 March, 2019; v1 submitted 5 September, 2018; originally announced September 2018.

Comments: 10 pages, 4 figures, matches version published in ApJ

Journal ref: ApJ, 873:67 (2019)

Urca reactions during neutron star inspiral

Authors: Phil Arras, Nevin N. Weinberg

Abstract: We study the impact of nonlinear bulk viscosity due to Urca reactions driven by tidally-induced fluid motion during binary neutron star inspiral. Fluid compression is computed for low radial order oscillation modes through an adiabatic, time-dependent solution of the mode amplitudes. Optically thin neutrino emission and heating rates are then computed from this adiabatic fluid motion. Calculations… ▽ More We study the impact of nonlinear bulk viscosity due to Urca reactions driven by tidally-induced fluid motion during binary neutron star inspiral. Fluid compression is computed for low radial order oscillation modes through an adiabatic, time-dependent solution of the mode amplitudes. Optically thin neutrino emission and heating rates are then computed from this adiabatic fluid motion. Calculations use direct and modified Urca reactions operating in a $M=1.4\, M_\odot$ neutron star, which is constructed using the Skyrme Rs equation of state. We find that the energy pumped into low order oscillation modes is not efficiently thermalized even by direct Urca reactions, with core temperatures reaching only $T \simeq 10^8\, {\rm K}$ during the inspiral. Although this is an order of magnitude larger than the heating due to shear viscosity considered by previous studies, it reinforces the result that the stars are quite cold at merger. Upon excitation of the lowest order g-mode, the chemical potential imbalance reaches $β> 1\, {\rm MeV}$ at orbital frequencies $ν_{\rm orb} > 200\, {\rm Hz}$, implying significant charged-current optical depths and Fermi blocking. To asses the importance of neutrino degeneracy effects, the neutrino transfer equation is solved in the static approximation for the three-dimensional density distribution, and the reaction rates are then computed including Fermi-blocking. We find that the heating rate is suppressed by factors of a few for $ν_{\rm orb} > 200\, {\rm Hz}$. The spectrum of emitted $ν_e$ and $\barν_e$, including radiation transfer effects, is presented for a range of orbital separations. △ Less

Submitted 11 June, 2018; originally announced June 2018.

Comments: Submitted to MNRAS. 14 pages and 9 figures

Radio Imaging With Information Field Theory

Authors: Philipp Arras, Jakob Knollmüller, Henrik Junklewitz, Torsten A. Enßlin

Abstract: Data from radio interferometers provide a substantial challenge for statisticians. It is incomplete, noise-dominated and originates from a non-trivial measurement process. The signal is not only corrupted by imperfect measurement devices but also from effects like fluctuations in the ionosphere that act as a distortion screen. In this paper we focus on the imaging part of data reduction in radio a… ▽ More Data from radio interferometers provide a substantial challenge for statisticians. It is incomplete, noise-dominated and originates from a non-trivial measurement process. The signal is not only corrupted by imperfect measurement devices but also from effects like fluctuations in the ionosphere that act as a distortion screen. In this paper we focus on the imaging part of data reduction in radio astronomy and present RESOLVE, a Bayesian imaging algorithm for radio interferometry in its new incarnation. It is formulated in the language of information field theory. Solely by algorithmic advances the inference could be sped up significantly and behaves noticeably more stable now. This is one more step towards a fully user-friendly version of RESOLVE which can be applied routinely by astronomers. △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: 5 pages, 3 figures

Journal ref: 2018 26th European Signal Processing Conference (EUSIPCO)

A Model of the H$α$ and Na Transmission Spectrum of HD 189733b

Authors: Chenliang Huang, Phil Arras, Duncan Christie, Zhi-Yun Li

Abstract: This paper presents a detailed hydrostatic model of the upper atmosphere of HD 189733b, with the goal of constraining its temperature, particle densities, and radiation field over the pressure range $10^{-4}-10\, μ\rm bar$, where the observed H$α$ transmission spectrum is produced. The atomic hydrogen level population is computed including both collisional and radiative transition rates. The Ly… ▽ More This paper presents a detailed hydrostatic model of the upper atmosphere of HD 189733b, with the goal of constraining its temperature, particle densities, and radiation field over the pressure range $10^{-4}-10\, μ\rm bar$, where the observed H$α$ transmission spectrum is produced. The atomic hydrogen level population is computed including both collisional and radiative transition rates. The Ly$α$ resonant scattering is computed using a Monte-Carlo simulation. The model transmission spectra are in broad agreement with the data. Excitation of the H(2$\ell$) population is mainly by Ly$α$ radiative excitation due to the large Ly$α$ intensity. The density of H(2$\ell$) is nearly flat over two decades in pressure, and is optically thick to H$α$. Additional models computed for a range of the stellar Lyman continuum (LyC) flux suggest that the variability in H$α$ transit depth may be due to the variability in the stellar LyC. Since metal lines provide the dominant cooling of this part of the atmosphere, the atmosphere structure is sensitive to the density of species such as Mg and Na which may themselves be constrained by observations. Since the H$α$ and Na D lines have comparable absorption depths, we argue that the center of the Na D lines are also formed in the atomic layer where the H$α$ line is formed. △ Less

Submitted 15 November, 2017; originally announced November 2017.

Comments: 23 pages, 17 figures, 5 tables. Accepted by ApJ

Tidal Dissipation in WASP-12

Authors: Nevin N. Weinberg, Meng Sun, Phil Arras, Reed Essick

Abstract: WASP-12 is a hot Jupiter system with an orbital period of $P= 1.1\textrm{ day}$, making it one of the shortest-period giant planets known. Recent transit timing observations by Maciejewski et al. (2016) and Patra et al. (2017) find a decreasing period with $P/|\dot{P}| = 3.2\textrm{ Myr}$. This has been interpreted as evidence of either orbital decay due to tidal dissipation or a long term oscilla… ▽ More WASP-12 is a hot Jupiter system with an orbital period of $P= 1.1\textrm{ day}$, making it one of the shortest-period giant planets known. Recent transit timing observations by Maciejewski et al. (2016) and Patra et al. (2017) find a decreasing period with $P/|\dot{P}| = 3.2\textrm{ Myr}$. This has been interpreted as evidence of either orbital decay due to tidal dissipation or a long term oscillation of the apparent period due to apsidal precession. Here we consider the possibility that it is orbital decay. We show that the parameters of the host star are consistent with either a $M_\ast \simeq 1.3 M_\odot$ main sequence star or a $M_\ast \simeq 1.2 M_\odot$ subgiant. We find that if the star is on the main sequence, the tidal dissipation is too inefficient to explain the observed $\dot{P}$. However, if it is a subgiant, the tidal dissipation is significantly enhanced due to nonlinear wave breaking of the dynamical tide near the star's center. The subgiant models have a tidal quality factor $Q_\ast'\simeq 2\times10^5$ and an orbital decay rate that agrees well with the observed $\dot{P}$. It would also explain why the planet survived for $\simeq 3\textrm{ Gyr}$ while the star was on the main sequence and yet is now inspiraling on a 3 Myr timescale. Although this suggests that we are witnessing the last $\sim 0.1\%$ of the planet's life, the probability of such a detection is a few percent given the observed sample of $\simeq 30$ hot Jupiters in $P<3\textrm{ day}$ orbits around $M_\ast>1.2 M_\odot$ hosts. △ Less

Submitted 2 October, 2017; originally announced October 2017.

Comments: 6 pages, 3 figures, accepted to ApJ Letters

Formation of Extremely Low-mass White Dwarfs Binaries

Authors: Meng Sun, Phil Arras

Abstract: Motivated by the discovery of a handful of pulsating, extremely low mass white dwarfs (ELM WDs, mass $M \lesssim 0.17\, M_\odot$) which likely have WD companions, this paper discusses binary formation models for these systems. Formation of an ELM WD by unstable mass transfer (MT) or a common envelope (CE) event is unlikely. Stable Roche-lobe overflow with conservative MT produces only… ▽ More Motivated by the discovery of a handful of pulsating, extremely low mass white dwarfs (ELM WDs, mass $M \lesssim 0.17\, M_\odot$) which likely have WD companions, this paper discusses binary formation models for these systems. Formation of an ELM WD by unstable mass transfer (MT) or a common envelope (CE) event is unlikely. Stable Roche-lobe overflow with conservative MT produces only $M \gtrsim 0.2\, M_\odot$. This paper discusses the formation of ELM WD using angular momentum losses by magnetic braking. The initially more massive star becomes the companion WD through a CE event. The initially less massive star becomes the ELM WD as an evolved donor in a Cataclysmic Variable binary. Evolutionary models are constructed using the Modules for Experiments in Stellar Astrophysics (MESA), with ELM WD progenitors in the range $M_d=1.0-1.5\, M_\odot$ and WD companions in the range $M_{\rm a} \simeq 0.4-0.9\, M_\odot$. A prescription to reduce magnetic braking for thin surface convection zones is included. Upon the thinning of the evolved donor's envelope, the donor star shrinks out of contact and MT ceases, revealing the ELM WD. Systems with small masses have previously been suggested as possible AM CVN's. Systems with large masses, up to the limit $M \simeq 0.18\, M_\odot$ at which shell flashes occur on the WD cooling track, will tend to expand out to orbital periods $P_{\rm orb} \gtrsim 15\, {\rm hr}$. In between this range, ELM WD may become pulsators both as pre-WD and on the WD cooling track. Brickhill's criterion for convective mode driving is used to estimate the location of the blue edge of the g-mode instability strip. △ Less

Submitted 23 March, 2018; v1 submitted 5 March, 2017; originally announced March 2017.

Comments: 14 pages, 23 figures, accepted by ApJ

Using long-term millisecond pulsar timing to obtain physical characteristics of the bulge globular cluster Terzan 5

Authors: Brian Prager, Scott Ransom, Paulo Freire, Jason Hessels, Ingrid Stairs, Phil Arras, Mario Cadelano

Abstract: Over the past decade the discovery of three unique stellar populations and a large number of confirmed pulsars within the globular cluster Terzan 5 has raised questions over its classification. Using the long-term radio pulsar timing of 36 millisecond pulsars in the cluster core, we provide new measurements of key physical properties of the system. As Terzan 5 is located within the galactic bulge,… ▽ More Over the past decade the discovery of three unique stellar populations and a large number of confirmed pulsars within the globular cluster Terzan 5 has raised questions over its classification. Using the long-term radio pulsar timing of 36 millisecond pulsars in the cluster core, we provide new measurements of key physical properties of the system. As Terzan 5 is located within the galactic bulge, stellar crowding and reddening make optical and near infrared observations difficult. Pulsar accelerations, however, allow us to study the intrinsic characteristics of the cluster independent of reddening and stellar crowding and probe the mass density profile without needing to quantify the mass to light ratio. Relating the spin and orbital periods of each pulsar to the acceleration predicted by a King model, we find a core density of $1.58\times$10$^6$ M$_\odot$ pc$^{-3}$, a core radius of 0.16 pc, a pulsar density profile $n\propto r^{-3.14}$, and a total mass of M$_{\rm T}$($R_\perp<$1.0 pc)$\simeq3.0\times$10$^5$ M$_\odot$ assuming a cluster distance of 5.9 kpc. Using this information we argue against Terzan 5 being a disrupted dwarf galaxy and discuss the possibility of Terzan 5 being a fragment of the Milky Way's proto-bulge. We also discuss whether low-mass pulsars were formed via electron capture supernovae or exist in a core full of heavy white dwarfs and hard binaries. Finally we provide an upper limit for the mass of a possible black hole at the core of the cluster of 3.0$\times$10$^4$ M$_\odot$. △ Less

Submitted 12 July, 2017; v1 submitted 13 December, 2016; originally announced December 2016.

Comments: 27 pages, 20 figures, 5 tables, thesis research, accepted

A New Model of Roche-lobe Overflow for Short-Period Gaseous Planets and Binary Stars

Authors: Brian Jackson, Phil Arras, Kaloyan Penev, Sarah Peacock, Pablo Marchant

Abstract: Some close-in gaseous exoplanets are nearly in Roche-lobe contact, and previous studies show tidal decay can drive hot Jupiters into contact during the main sequence of their host stars. Improving upon a previous model, we present a revised model for mass transfer in a semi-detached binary system that incorporates an extended atmosphere around the donor and allows for an arbitrary mass ratio. We a… ▽ More Some close-in gaseous exoplanets are nearly in Roche-lobe contact, and previous studies show tidal decay can drive hot Jupiters into contact during the main sequence of their host stars. Improving upon a previous model, we present a revised model for mass transfer in a semi-detached binary system that incorporates an extended atmosphere around the donor and allows for an arbitrary mass ratio. We apply this new formalism to hypothetical, confirmed, and candidate planetary systems to estimate mass loss rates and compare with models of evaporative mass loss. Overflow may be significant for hot Neptunes out to periods of $\sim$ 2 days, while for hot Jupiters, it may only be important inward of 0.5 days. We find that CoRoT-24 b may be losing mass at a rate of more than an Earth mass in a Gyr. The hot Jupiter WASP-12 b may lose an Earth mass in a Myr, while the putative planet orbiting a T-Tauri star PTFO8-8695 might shed its atmosphere in a few Myrs. We point out that the orbital expansion that can accompany mass transfer may be less effective than previously considered because the gas accreted by the host star removes some of the system's angular momentum from the orbit, but simple scaling arguments suggest that the Roche-lobe overflow might remain stable. Consequently, the recently discovered small planets in ultra-short-periods ($<$ 1 day) may not be the remnants of hot Jupiters/Neptunes. The new model presented here has been incorporated into Modules for Experiments in Stellar Astrophysics (MESA). △ Less

Submitted 13 December, 2016; originally announced December 2016.

Comments: Submitted to ApJ, 2016 April Accepted by ApJ on 2016 Dec 13; 43 pages, 9 figures, 2 tables; All code, figures, data, and manuscript files available here - https://github.com/decaelus/Jacksonetal_2016_RLO

Investigation of the environment around close-in transiting exoplanets using CLOUDY

Authors: Jake D. Turner, Duncan Christie, Phil Arras, Robert E. Johnson, Carl Schmidt

Abstract: It has been suggested that hot stellar wind gas in a bow shock around an exoplanet is sufficiently opaque to absorb stellar photons and give rise to an observable transit depth at optical and UV wavelengths. In the first part of this paper, we use the CLOUDY plasma simulation code to model the absorption from X-ray to radio wavelengths by 1-D slabs of gas in coronal equilibrium with varying densit… ▽ More It has been suggested that hot stellar wind gas in a bow shock around an exoplanet is sufficiently opaque to absorb stellar photons and give rise to an observable transit depth at optical and UV wavelengths. In the first part of this paper, we use the CLOUDY plasma simulation code to model the absorption from X-ray to radio wavelengths by 1-D slabs of gas in coronal equilibrium with varying densities ($10^{4}-10^{8} \, {\rm cm^{-3}}$) and temperatures ($2000-10^{6} \ {\rm K}$) illuminated by a solar spectrum. For slabs at coronal temperatures ($10^{6} \ {\rm K}$) and densities even orders of magnitude larger than expected for the compressed stellar wind ($10^{4}-10^{5} \, {\rm cm^{-3}}$), we find optical depths orders of magnitude too small ($> 3\times10^{-7}$) to explain the $\sim3\%$ UV transit depths seen with Hubble. Using this result and our modeling of slabs with lower temperatures ($2000-10^4 {\rm K}$), the conclusion is that the UV transits of WASP-12b and HD 189733b are likely due to atoms originating in the planet, as the stellar wind is too highly ionized. A corollary of this result is that transport of neutral atoms from the denser planetary atmosphere outward must be a primary consideration when constructing physical models. In the second part of this paper, additional calculations using CLOUDY are carried out to model a slab of planetary gas in radiative and thermal equilibrium with the stellar radiation field. Promising sources of opacity from the X-ray to radio wavelengths are discussed, some of which are not yet observed. △ Less

Submitted 3 March, 2016; originally announced March 2016.

Comments: 13 pages, 7 figures, 4 Tables. Accepted for publication in MNRAS (March 3, 2016)

Tidal Decay and Stable Roche-Lobe Overflow of Short-Period Gaseous Exoplanets

Authors: Brian Jackson, Emily Jensen, Sarah Peacock, Phil Arras, Kaloyan Penev

Abstract: Many gaseous exoplanets in short-period orbits are on the verge or are in the process of Roche-lobe overflow (RLO). Moreover, orbital stability analysis shows tides can drive many hot Jupiters to spiral inevitably toward their host stars. Thus, the coupled processes of orbital evolution and RLO likely shape the observed distribution of close-in exoplanets and may even be responsible for producing… ▽ More Many gaseous exoplanets in short-period orbits are on the verge or are in the process of Roche-lobe overflow (RLO). Moreover, orbital stability analysis shows tides can drive many hot Jupiters to spiral inevitably toward their host stars. Thus, the coupled processes of orbital evolution and RLO likely shape the observed distribution of close-in exoplanets and may even be responsible for producing some of the short-period rocky planets. However, the exact outcome for an overflowing planet depends on its internal response to mass loss, and the accompanying orbital evolution can act to enhance or inhibit RLO. In this study, we apply the fully-featured and robust Modules for Experiments in Stellar Astrophysics (MESA) suite to model RLO of short-period gaseous planets. We show that, although the detailed evolution may depend on several properties of the planetary system, it is largely determined by the core mass of the overflowing gas giant. In particular, we find that the orbital expansion that accompanies RLO often stops and reverses at a specific maximum period that depends on the core mass. We suggest that RLO may often strand the remnant of a gas giant near this orbital period, which provides an observational prediction that can corroborate the hypothesis that short-period gas giants undergo RLO. We conduct a preliminary comparison of this prediction to the observed population of small, short-period planets and find some planets in orbits that may be consistent with this picture. To the extent that we can establish some short-period planets are indeed the remnants of gas giants, that population can elucidate the properties of gas giant cores, the properties of which remain largely unconstrained. △ Less

Submitted 23 May, 2016; v1 submitted 1 March, 2016; originally announced March 2016.

Comments: Accepted to "Celestial Mechanics and Dynamical Astronomy" special issue on tides. Several changes based on referee comments, including to the title of the paper. Some new analysis of non-conservative (but still stable) mass transfer as well. Article repository and data files linked to here -- http://www.astrojack.com/research/

Axisymmetric Simulations of Hot Jupiter-Stellar Wind Hydrodynamic Interaction

Authors: Duncan Christie, Phil Arras, Zhi-Yun Li

Abstract: Gas giant exoplanets orbiting at close distances to the parent star are subjected to large radiation and stellar wind fluxes. In this paper, hydrodynamic simulations of the planetary upper atmosphere and its interaction with the stellar wind are carried out to understand the possible flow regimes and how they affect the Lyman-alpha transmission spectrum. Following Tremblin and Chiang, charge excha… ▽ More Gas giant exoplanets orbiting at close distances to the parent star are subjected to large radiation and stellar wind fluxes. In this paper, hydrodynamic simulations of the planetary upper atmosphere and its interaction with the stellar wind are carried out to understand the possible flow regimes and how they affect the Lyman-alpha transmission spectrum. Following Tremblin and Chiang, charge exchange reactions are included to explore the role of energetic atoms as compared to thermal particles. In order to understand the role of the tail as compared to the leading edge of the planetary gas, the simulations were carried out under axisymmetry, and photoionization and stellar wind electron impact ionization reactions were included to limit the extent of the neutrals away from the planet. By varying the planetary gas temperature, two regimes are found. At high temperature, a supersonic planetary wind is found, which is turned around by the stellar wind and forms a tail behind the planet. At lower temperatures, the planetary wind is shut off when the stellar wind penetrates inside where the sonic point would have been. In this regime mass is lost by viscous interaction at the boundary between planetary and stellar wind gases. Absorption by cold hydrogen atoms is large near the planetary surface, and decreases away from the planet as expected. The hot hydrogen absorption is in an annulus and typically dominated by the tail, at large impact parameter, rather than by the thin leading edge of the mixing layer near the substellar point. △ Less

Submitted 20 January, 2016; originally announced January 2016.

Comments: Accepted to ApJ

Seismology of Rapidly Rotating Accreting White Dwarfs

Authors: Dean M. Townsley, Phil Arras, Lars Bildsten

Abstract: A number of White Dwarfs (WDs) in cataclysmic binaries have shown brightness variations consistent with non-radial oscillations as observed in isolated WDs. A few objects have been well-characterized with photometric campaigns in the hopes of gleaning information about the mass, spin, and possibly internal structural characteristics. The novel aspect of this work is the possiblity to measure or co… ▽ More A number of White Dwarfs (WDs) in cataclysmic binaries have shown brightness variations consistent with non-radial oscillations as observed in isolated WDs. A few objects have been well-characterized with photometric campaigns in the hopes of gleaning information about the mass, spin, and possibly internal structural characteristics. The novel aspect of this work is the possiblity to measure or constrain the interior structure and spin rate of WDs which have spent gigayears accreting material from their companion, undergoing thousands of nova outbursts in the process. In addition, variations in the surface temperature affect the site of mode driving, and provide unique and challenging tests for mode driving theories previously applied to isolated WD's. Having undergone long-term accretion, these WDs are expected to have been spun up. Spin periods in the range 60-100 seconds have been measured by other means for two objects, GW Lib and V455 And. Compared to typical mode frequencies, the spin frequency may be similar or higher, and the Coriolis force can no longer be treated as a small perturbation on the fluid motions. We present the results of a non-perturbative calculation of the normal modes of these WDs, using interior thermal structures appropriate to accreting systems. This includes a discussion of the surface brightness distributions, which are strongly modified from the non-rotating case. Using the measured spin period of approximately 100 seconds, we show that the observed pulsations from GW Lib are consistent with the three lowest azimuthal order rotationally modified modes that have the highest frequency in the stellar frame. The high frequencies are needed for the convective driving, but are then apparently shifted to lower frequencies by a combination of their pattern motion and the WD rotation. △ Less

Submitted 8 January, 2016; originally announced January 2016.

Comments: 6 pages, 4 figures, proceedings from 2010 conference "The Physics of Accreting Compact Binaries" Kyoto, Japan. Submitted to proceedings in 2011

Magnetohydrodynamic Simulations of Hot Jupiter Upper Atmospheres

Authors: George B. Trammell, Zhi-Yun Li, Phil Arras

Abstract: Two-dimensional simulations of hot Jupiter upper atmospheres including the planet's magnetic field are presented. The goal is to explore magnetic effects on the layer of the atmosphere that is ionized and heated by stellar EUV radiation, and the imprint of these effects on the Ly$α$ transmission spectrum. The simulations are axisymmetric, isothermal, and include both rotation and azimuth-averaged… ▽ More Two-dimensional simulations of hot Jupiter upper atmospheres including the planet's magnetic field are presented. The goal is to explore magnetic effects on the layer of the atmosphere that is ionized and heated by stellar EUV radiation, and the imprint of these effects on the Ly$α$ transmission spectrum. The simulations are axisymmetric, isothermal, and include both rotation and azimuth-averaged stellar tides. Mass density is converted to atomic hydrogen density through the assumption of ionization equilibrium. The three-zone structure -- polar dead zone, mid-latitude wind zone, and equatorial dead zone -- found in previous analytic calculations is confirmed. For a magnetic field comparable to that of Jupiter, the equatorial dead zone, which is confined by the magnetic field and corotates with the planet, contributes at least half of the transit signal. For even stronger fields, the gas escaping in the mid-latitude wind zone is found to have a smaller contribution to the transit depth than the equatorial dead zone. Transmission spectra computed from the simulations are compared to HST STIS and ACS data for HD 209458b and HD 189733b, and the range of model parameters consistent with the data is found. The central result of this paper is that the transit depth increases strongly with magnetic field strength when the hydrogen ionization layer is magnetically dominated, for dipole magnetic field $B_0 > 10\ {\rm G}$. Hence transit depth is sensitive to magnetic field strength, in addition to standard quantities such as the ratio of thermal to gravitational binding energies. △ Less

Submitted 23 April, 2014; originally announced April 2014.

Comments: accepted at ApJ

Dynamical resonance locking in tidally interacting binary systems

Authors: Joshua Burkart, Eliot Quataert, Phil Arras

Abstract: We examine the dynamics of resonance locking in detached, tidally interacting binary systems. In a resonance lock, a given stellar or planetary mode is trapped in a highly resonant state for an extended period of time, during which the spin and orbital frequencies vary in concert to maintain the resonance. This phenomenon is qualitatively similar to resonance capture in planetary dynamics. We show… ▽ More We examine the dynamics of resonance locking in detached, tidally interacting binary systems. In a resonance lock, a given stellar or planetary mode is trapped in a highly resonant state for an extended period of time, during which the spin and orbital frequencies vary in concert to maintain the resonance. This phenomenon is qualitatively similar to resonance capture in planetary dynamics. We show that resonance locks can accelerate the course of tidal evolution in eccentric systems and also efficiently couple spin and orbital evolution in circular binaries. Previous analyses of resonance locking have not treated the mode amplitude as a fully dynamical variable, but rather assumed the adiabatic (i.e. Lorentzian) approximation valid only in the limit of relatively strong mode damping. We relax this approximation, analytically derive conditions under which the fixed point associated with resonance locking is stable, and further check these analytic results using numerical integrations of the coupled mode, spin, and orbital evolution equations. These show that resonance locking can sometimes take the form of complex limit cycles or even chaotic trajectories. We provide simple analytic formulae that define the binary and mode parameter regimes in which resonance locks of some kind occur (stable, limit cycle, or chaotic). We briefly discuss the astrophysical implications of our results for white dwarf and neutron star binaries as well as eccentric stellar binaries. △ Less

Submitted 25 October, 2014; v1 submitted 17 December, 2013; originally announced December 2013.

Comments: 16 pages, 11 figures

Journal ref: MNRAS, 443, 2957