Abstract visibility
Abstract text size

Papers for Tuesday, Mar 16 2021

Papers with local authors

Allyson A. Sheffield, Aidan Z. Subrahimovic, Mohammad Refat, Rachael L. Beaton, Sten Hasselquist, Christian R. Hayes, Adrian. M. Price-Whelan, Danny Horta, Steven R. Majewski, Katia Cunha, Verne V. Smith, Jose G. Fernandez-Trincado, Jennifer S. Sobeck, Ricardo R. Munoz, D. A. Garcia-Hernandez, Richard R. Lane, Christian Nitschelm, Alexandre Roman-Lopes

Accepted to ApJ

0 votes
Paper 6 — arXiv:2103.07488
0 votes
Paper 6 — arXiv:2103.07488

We present the kinematic and chemical profiles of red giant stars observed by the APOGEE-2 survey in the direction of the Jhelum stellar stream, a Milky Way substructure located in the inner halo of the Milky Way at a distance from the Sun of $\approx$ 13 kpc. From the six APOGEE-2 Jhelum pointings, we isolate stars with log($g$) $<$ 3.5, leaving a sample of 289 red giant stars. From this sample of APOGEE giants, we identified seven stars that are consistent with the astrometric signal from $Gaia$ DR2 for this stream. Of these seven, one falls onto the RGB along the same sequence as the Jhelum stars presented by \cite{ji20}. This new Jhelum member has [Fe/H]=-2.2 and is at the tip of the red giant branch. By selecting high orbital eccentricity, metal-rich stars, we identify red giants in our APOGEE sample that are likely associated with the $Gaia$-Enceladus-Sausage (GES) merger. We compare the abundance profiles of the Jhelum stars and GES stars and find similar trends in $\alpha$-elements, as expected for low-metallicity populations. However, we find that the orbits for GES and Jhelum stars are not generally consistent with a shared origin. The chemical abundances for the APOGEE Jhelum star and other confirmed members of the stream are similar to stars in known stellar streams and thus are consistent with an accreted dwarf galaxy origin for the progenitor of the stream, although we cannot rule out a globular cluster origin.

Ryan C Terrien, Joe P Ninan, Scott A Diddams, Suvrath Mahadevan, Samuel Halverson, Chad Bender, Connor Fredrick, Fred Hearty, Jeff Jennings, Andrew J. Metcalf, Andrew Monson, Arpita Roy, Christian Schwab, Gudmundur Stefansson

Accepted for publication in The Astronomical Journal

0 votes
Paper 61 — arXiv:2103.08456
0 votes
Paper 61 — arXiv:2103.08456

The comb-like spectrum of a white light-illuminated Fabry-P\'{e}rot etalon can serve as a cost-effective and stable reference for precise Doppler measurements. Understanding the stability of these devices across their broad (100's of nm) spectral bandwidths is essential to realize their full potential as Doppler calibrators. However, published descriptions remain limited to small bandwidths or short timespans. We present a $\sim6$ month broadband stability monitoring campaign of the Fabry-P\'{e}rot etalon system deployed with the near-infrared Habitable Zone Planet Finder spectrograph (HPF). We monitor the wavelengths of each of $\sim3500$ resonant modes measured in HPF spectra of this Fabry-P\'{e}rot etalon (free spectral range = 30 GHz, bandwidth = 820 - 1280 nanometers), leveraging the accuracy and precision of an electro-optic frequency comb reference. These results reveal chromatic structure in the Fabry-P\'{e}rot mode locations and in their evolution with time. We measure an average drift on the order of 2 cm s $^{-1}$ d$^{-1}$, with local departures up to $\pm5$ cm s $^{-1}$ d$^{-1}$. We discuss these behaviors in the context of the Fabry-P\'{e}rot etalon mirror dispersion and other optical properties of the system, and the implications for the use of similar systems for precise Doppler measurements. Our results show that this system supports the wavelength calibration of HPF at the $\lesssim10$ cm s $^{-1}$ level over a night and at the $\lesssim30$ cm s $^{-1}$ level over $\sim10$ d. Our results also highlight the need for long-term and spectrally-resolved study of similar systems that will be deployed to support Doppler measurement precision approaching $\sim10$ cm s $^{-1}$.

All other papers

A. Drlica-Wagner, J. L. Carlin, D. L. Nidever, P. S. Ferguson, N. Kuropatkin, M. Adamów, W. Cerny, Y. Choi, J. H. Esteves, C. E. Martínez-Vázquez, S. Mau, A. E. Miller, B. Mutlu-Pakdil, E. H. Neilsen, K. A. G. Olsen, A. B. Pace, A. H. Riley, J. D. Sakowska, D. J. Sand, L. Santana-Silva, E. J. Tollerud, D. L. Tucker, A. K. Vivas, E. Zaborowski, A. Zenteno, T. M. C. Abbott, S. Allam, K. Bechtol, C. P. M. Bell, E. F. Bell, P. Bilaji, C. R. Bom, J. A. Carballo-Bello, M.-R. L. Cioni, A. Diaz-Ocampo, T. J. L. de Boer, D. Erkal, R. A. Gruendl, D. Hernandez-Lang, A. K. Hughes, D. J. James, L. C. Johnson, T. S. Li, Y.-Y. Mao, D. Martínez-Delgado, P. Massana, M. McNanna, R. Morgan, E. O. Nadler, N. E. D. Noël, A. Palmese, A. H. G. Peter, E. S. Rykoff, J. Sánchez, N. Shipp, et al. (15 additional authors not shown)

30 pages, 14 figures, 4 tables; submitted to AAS Journals; public data release at this https URL

The DECam Local Volume Exploration survey (DELVE) is a 126-night survey program on the 4-m Blanco Telescope at the Cerro Tololo Inter-American Observatory in Chile. DELVE seeks to understand the characteristics of faint satellite galaxies and other resolved stellar substructures over a range of environments in the Local Volume. DELVE will combine new DECam observations with archival DECam data to cover ~15000 deg$^2$ of high-Galactic-latitude (|b| > 10 deg) southern sky to a 5$\sigma$ depth of g,r,i,z ~ 23.5 mag. In addition, DELVE will cover a region of ~2200 deg$^2$ around the Magellanic Clouds to a depth of g,r,i ~ 24.5 mag and an area of ~135 deg$^2$ around four Magellanic analogs to a depth of g,i ~ 25.5 mag. Here, we present an overview of the DELVE program and progress to date. We also summarize the first DELVE public data release (DELVE DR1), which provides point-source and automatic aperture photometry for ~520 million astronomical sources covering ~5000 deg$^2$ of the southern sky to a 5$\sigma$ point-source depth of g=24.3, r=23.9, i=23.3, and z=22.8 mag. DELVE DR1 is publicly available via the NOIRLab Astro Data Lab science platform.

Marilena Spavone, Davor Krajnović, Eric Emsellem, Enrichetta Iodice, Mark den Brok

Accepted for publication in A&A. 12 pages, 7 figures

In this paper we present the new deep images from the VEGAS survey of three massive ($M_{*} \simeq 10^{12}$~M$_\odot$) galaxies from the MUSE Most Massive Galaxies (M3G) project, with distances in the range $151\leq D \leq 183$ Mpc: PGC007748, PGC015524 and PGC049940. The long integration time and the wide field of view of OmegaCam@VST allowed us to map the light and color distributions down to $\mu_g\simeq30$~mag/arcsec$^2$ and out to $\sim 2R_e$. The deep data are crucial to estimate the contribution of the different galaxy's components, in particular the accreted fraction in the stellar halo. The available integral-field observations with MUSE cover a limited portion of each galaxy (out to $\sim 1R_e$), but, from the imaging analysis we find that they map the kinematics and stellar population beyond the first transition radius, where the contribution of the accreted component starts to dominate. The main goal of this work is to correlate the scales of the different components derived from the image analysis with the kinematics and stellar population profiles from the MUSE data. Results were used to address the assembly history of the three galaxies with the help of the theoretical predictions. Our results suggest that PGC049940 has the lowest accreted mass fraction of 77%. The higher accreted mass fraction estimated for PGC007748 and PGC015524 (86% and 89%, respectively), combined with the flat $\lambda_R$ profiles suggest that a great majority of the mass has been acquired through major mergers, which have also shaped the shallower metallicity profiles observed at larger radii.

James H. Leftley, Konrad R. W. Tristram, Sebastian F. Hönig, Daniel Asmus, Makoto Kishimoto, Poshak Gandhi

12 pages, 11 figures, and 5 tables in main text. Accepted for publication in ApJ

Infrared interferometry has fuelled a paradigm shift in our understanding of the dusty structure in the central parsecs of Active Galactic Nuclei (AGN). The dust is now thought to comprise of a hot ($\sim1000\,$K) equatorial disk, some of which is blown into a cooler ($\sim300\,$K) polar dusty wind by radiation pressure. In this paper, we utilise the new near-IR interferometer GRAVITY on the Very Large Telescope Interferometer (VLTI) to study a Type 1.2 AGN hosted in the nearby Seyfert galaxy ESO323-G77. By modelling the squared visibility and closure phase, we find that the hot dust is equatorially extended, consistent with the idea of a disk, and shows signs of asymmetry in the same direction. Furthermore, the data is fully consistent with the hot dust size determined by K band reverberation mapping as well as the predicted size from a CAT3D-WIND model created in previous work using the SED of ESO323-G77 and observations in the mid-IR from VLTI/MIDI.

Raghunath Ghara, Sambit K. Giri, Benedetta Ciardi, Garrelt Mellema, Saleem Zaroubi

13 Pages, 6 Figures, 5 Tables, Accepted for publication in MNRAS

The Murchison Widefield Array (MWA) team has derived new upper limits on the spherically averaged power spectrum of the 21-cm signal at six redshifts in the range $z \approx 6.5-8.7$. We use these upper limits and a Bayesian inference framework to derive constraints on the ionization and thermal state of the intergalactic medium (IGM) as well as on the strength of a possible additional radio background. We do not find any constraints on the state of the IGM for $z\gtrsim 7.8$ if no additional radio background is present. In the presence of such a radio background, the 95 per cent credible intervals of the disfavoured models at redshift $\gtrsim 6.5 $ correspond to an IGM with a volume averaged fraction of ionized regions below 0.6 and an average gas temperature $\lesssim 10^3$ K. In these models, the heated regions are characterised by a temperature larger than that of the radio background, and by a distribution with characteristic size $\lesssim 10$ $h^{-1}$ Mpc and a full width at half maximum (FWHM) of $\lesssim 30$ $h^{-1}$ Mpc. Within the same credible interval limits, we exclude an additional radio background of at least $0.008\%$ of the CMB at 1.42 GHz.

Elisa Bortolas, Matteo Bonetti, Massimo Dotti, Alessandro Lupi, Pedro R. Capelo, Lucio Mayer, Alberto Sesana

16 pages, 12 figures, 2 tables. Submitted to MNRAS. Comments welcome!

In this paper, we explore the impact of a galactic bar on the inspiral time-scale of a massive perturber (MP) within a Milky Way-like galaxy. We integrate the orbit of MPs in a multi-component galaxy model via a semi-analytical approach including an accurate treatment for dynamical friction generalized to rotationally supported backgrounds. We compare the MP evolution in a galaxy featuring a Milky Way-like rotating bar to the evolution within an analogous axisymmetric galaxy without the bar. We find that the bar presence may significantly affect the inspiral, sometimes making it shorter by a factor of a few, sometimes hindering it for a Hubble time, implying that dynamical friction alone is greatly insufficient to fully characterize the orbital decay. The effect of the bar is more prominent for initially in-plane, prograde MPs, especially those crossing the bar co-rotation radius or outer Lindblad resonance. In the barred galaxy, we find the sinking of the most massive MPs (>~10^7.5 Msun) approaching the galaxy from large separations (>~8 kpc) to be most efficiently hampered. Neglecting the effect of global torques associated to the non-symmetric mass distribution is thus not advisable even within our idealized, smooth Milky Way model, and it should be avoided when dealing with more complex and realistic galaxy systems. This has important implications for the orbital decay of massive black holes in late-type spirals, the natural candidate sources to be detected with the Laser Interferometer Space Antenna (LISA).

Brian L. Flores, D. John Hillier

Accepted for publication in MNRAS

Hot massive stars exhibit strong stellar winds that enrich the surrounding interstellar medium and affect the stars' evolution. However, the winds are inhomogeneous (clumped) and are difficult to model in radiative transfer codes. To produce more realistic spectra many codes use a volume-filling factor approach to incorporate the effects of clumping. While this approach is convenient it is simplistic. We introduce an alternative approach to incorporate clumping by assuming the wind is composed of dense spherical shells. Using this approach in the radiative transfer code CMFGEN we produce synthetic spectra for AzV83, an O7Iaf+ supergiant located in the Small Magellanic Cloud. The spectrum of AzV83 is rich in both photospheric and wind features, making it an ideal candidate with which to investigate the physical characteristics of stellar winds. Synthetic spectra are compared to the star's observed spectrum to better characterize the influence of clumped winds on spectral features, and to better understand the limitations of the volume-filling factor approach. The approach using spherical shells yields similar wind parameters to those obtained using the volume-filling factor approach although a slightly higher mass-loss rate is required to fit H$\alpha$. As expected, the interclump medium in the model with shells allows the high ionisation resonance transitions of N V and O VI to be fitted using $L_{\rm X-ray}/L_{\rm Bol} \approx 10^{-7}$ which is typically observed for O stars, and which is a factor of ten lower than needed with the volume-filling factor approach.

Alex Drlica-Wagner, Edgar Marrufo Villalpando, Judah O'Neil, Juan Estrada, Stephen Holland, Noah Kurinsky, Ting S. Li, Guillermo Fernandez Moroni, Javier Tiffenberg, Sho Uemura

14 pages, 6 figures, 2 tables

We characterize the response of a novel 250 $\mu$m thick, fully-depleted Skipper Charged-Coupled Device (CCD) to visible/near-infrared light with a focus on potential applications for astronomical observations. We achieve stable, single-electron resolution with readout noise $\sigma \sim 0.18$ e$^{-}$ rms/pix from 400 non-destructive measurements of the charge in each pixel. We verify that the gain derived from photon transfer curve measurements agrees with the gain calculated from the quantized charge of individual electrons to within < 1%. We also perform relative quantum efficiency measurements and demonstrate high relative quantum efficiency at optical/near-infrared wavelengths, as is expected for a thick, fully depleted detector. Finally, we demonstrate the ability to perform multiple non-destructive measurements and achieve sub-electron readout noise over configurable subregions of the detector. This work is the first step toward demonstrating the utility of Skipper CCDs for future astronomical and cosmological applications.

Azim Ahmadzadeh, Berkay Aydin, Manolis K. Georgoulis, Dustin J. Kempton, Sushant S. Mahajan, Rafal A. Angryk

The Astrophysical Journal Supplement series (2021)

We present a case study of solar flare forecasting by means of metadata feature time series, by treating it as a prominent class-imbalance and temporally coherent problem. Taking full advantage of pre-flare time series in solar active regions is made possible via the Space Weather Analytics for Solar Flares (SWAN-SF) benchmark dataset; a partitioned collection of multivariate time series of active region properties comprising 4075 regions and spanning over 9 years of the Solar Dynamics Observatory (SDO) period of operations. We showcase the general concept of temporal coherence triggered by the demand of continuity in time series forecasting and show that lack of proper understanding of this effect may spuriously enhance models' performance. We further address another well-known challenge in rare event prediction, namely, the class-imbalance issue. The SWAN-SF is an appropriate dataset for this, with a 60:1 imbalance ratio for GOES M- and X-class flares and a 800:1 for X-class flares against flare-quiet instances. We revisit the main remedies for these challenges and present several experiments to illustrate the exact impact that each of these remedies may have on performance. Moreover, we acknowledge that some basic data manipulation tasks such as data normalization and cross validation may also impact the performance -- we discuss these problems as well. In this framework we also review the primary advantages and disadvantages of using true skill statistic and Heidke skill score, as two widely used performance verification metrics for the flare forecasting task. In conclusion, we show and advocate for the benefits of time series vs. point-in-time forecasting, provided that the above challenges are measurably and quantitatively addressed.

Felipe Andrade-Santos, Gabriel W. Pratt, Jean-Baptiste Melin, Monique Arnaud, Christine Jones, William R. Forman, Etienne Pointecouteau, Iacopo Bartalucci, Alexey Vikhlinin, Stephen S. Murray, Pasquale Mazzotta, Stefano Borgani, Lorenzo Lovisari, Reinout J. van Weeren, Ralph P. Kraft, Laurence P. David, Simona Giacintucci

28 pages, 15 figures, and 6 tables. Submitted to ApJ

Using Chandra observations, we derive the $Y_{\rm X}$ proxy and associated total mass measurement, $M_{500}^{\rm Y_X}$, for 147 clusters with $z \leq 0.35$ from the Planck Early Sunyaev-Zel'dovich catalog, and for 80 clusters with $z \leq 0.30$ from an X-ray flux-limited sample. We re-extract the Planck $Y_{\rm SZ}$ measurements and obtain the corresponding mass proxy, $M_{500}^{\rm SZ}$, from the full Planck mission maps, minimizing the Malmquist bias due to observational scatter. The masses re-extracted using the more precise X-ray position and characteristic size agree with the published PSZ2 values, but yield a significant reduction in the scatter (by a factor of two) in the $M_{500}^{\rm SZ}$-$M_{500}^{\rm X}$ relation. The slope is $0.93\pm0.03$, and the median ratio, $M_{500}^{\rm SZ}/M_{500}^{\rm X}= 0.91\pm0.01$, is within the expectations from known X-ray calibration systematics. The $Y_{\rm SZ}/Y_{\rm X}$ ratio is $0.88\pm0.02$, in good agreement with predictions from cluster structure, and implying a low level of clumpiness. In agreement with the findings of the Planck Collaboration, the slope of the $Y_{\rm SZ}$-$D_{\rm A}^{-2} Y_{X}$ flux relation is significantly less than unity ($0.89\pm0.01$). Using extensive simulations, we show that this result is not due to selection effects, intrinsic scatter, or covariance between quantities. We demonstrate analytically that changing the $Y_{\rm SZ}$-$Y_{X}$ relation from apparent flux to intrinsic properties results in a best-fit slope that is closer to unity and increases the dispersion about the relation. The redistribution resulting from this transformation implies that the best fit parameters of the $M_{500}^{\rm SZ}$-$M_{500}^{\rm X}$ relation will be sample-dependent.

Thai Pham, Opher Ganel, Azita Valinia, Nicholas Siegler, Brendan Crill, Mario R. Perez

17 pages, 7 tables, 5 figures

NASA Astrophysics Division funds development of cutting-edge technology to enable its missions to achieve ambitious and groundbreaking science goals. These technology development efforts are managed by the Physics of the Cosmos, Cosmic Origins, and Exoplanet Exploration Programs. The NASA Strategic Astrophysics Technology Program (SAT) was established in 2009 as a new technology maturation program to fill the gap in the Technology Readiness Level range from 3 to 6. Since program inception, 100 SAT grants have been openly competed and awarded, along with dozens of direct-funded projects, leading to a host of technologies advancing their Technology Readiness Levels and/or being infused into space and suborbital missions and ground-based projects. We present the portfolio distribution in terms of specific technology areas addressed, including optics, detectors, coatings, corona graphs, star shades, lasers, electronics, and cooling subsystems. We show an analysis of the rate of Technology Readiness Level advances, infusion success stories, and other benefits such as training the future astrophysics workforce, including students and postdoctoral fellows hired by projects. Finally, we present the Astrophysics Division current strategic technology maturation priorities for investment, enabling a range of future strategic astrophysics missions.

Jonathan Mackey, Samuel Green, Maria Moutzouri, Thomas J. Haworth, Robert D. Kavanagh, Davit Zargaryan, Maggie Celeste

27 pages, accepted for publication in MNRAS

Expanding nebulae are produced by mass loss from stars, especially during late stages of evolution. Multi-dimensional simulation of these nebulae requires high resolution near the star and permits resolution that decreases with distance from the star, ideally with adaptive timesteps. We report the implementation and testing of static mesh-refinement in the radiation-magnetohydrodynamics code PION, and document its performance for 2D and 3D calculations. The bow shock produced by a hot, magnetized, slowly rotating star as it moves through the magnetized ISM is simulated in 3D, highlighting differences compared with 2D calculations. Latitude-dependent, time-varying magnetized winds are modelled and compared with simulations of ring nebulae around blue supergiants from the literature. A 3D simulation of the expansion of a fast wind from a Wolf-Rayet star into the slow wind from a previous red supergiant phase of evolution is presented, with results compared with results in the literature and analytic theory. Finally the wind-wind collision from a binary star system is modelled with 3D MHD, and the results compared with previous 2D hydrodynamic calculations. A python library is provided for reading and plotting simulation snapshots, and the generation of synthetic infrared emission maps using torus is also demonstrated. It is shown that state-of-the-art 3D MHD simulations of wind-driven nebulae can be performed using PION with reasonable computational resources. The source code and user documentation is made available for the community under a BSD3 licence.

Sebastiaan Y. Haffert, Jared R. Males, Laird M. Close, Kyle Van Gorkom, Joseph D. Long, Alexander D. Hedglen, Olivier Guyon, Lauren Schatz, Maggie Kautz, Jennifer Lumbres, Alex Rodack, Justin M. Knight, He Sun, Kevin Fogarty

Accepted for publication in JATIS

The search for exoplanets is pushing adaptive optics systems on ground-based telescopes to their limits. One of the major limitations at small angular separations, exactly where exoplanets are predicted to be, is the servo-lag of the adaptive optics systems. The servo-lag error can be reduced with predictive control where the control is based on the future state of the atmospheric disturbance. We propose to use a linear data-driven integral predictive controller based on subspace methods that is updated in real time. The new controller only uses the measured wavefront errors and the changes in the deformable mirror commands, which allows for closed-loop operation without requiring pseudo-open loop reconstruction. This enables operation with non-linear wavefront sensors such as the pyramid wavefront sensor. We show that the proposed controller performs near-optimal control in simulations for both stationary and non-stationary disturbances and that we are able to gain several orders of magnitude in raw contrast. The algorithm has been demonstrated in the lab with MagAO-X, where we gain more than two orders of magnitude in contrast.

Joss Bland-Hawthorn (U Sydney), Matt Sellars (ANU), John Bartholomew (U Sydney)

14 pages, 5 figures. Paper written for special double issue of Astrophotonics in JOSA-B and Applied Optics

Thomas Young's slit experiment lies at the heart of classical interference and quantum mechanics. Over the last fifty years, it has been shown that particles (e.g. photons, electrons, large molecules), even individual particles, generate an interference pattern at a distant screen after passage through a double slit, thereby demonstrating wave-particle duality. We revisit this famous experiment by replacing both slits with single-mode fibre inputs to two independent quantum memories that are capable of storing the incident electromagnetic field's amplitude and phase as a function of time. At a later time, the action is reversed: the quantum memories are read out in synchrony and the single-mode fibre outputs are allowed to interact consistent with the original observation. In contrast to any classical memory device, the write and read processes of a quantum memory are non-destructive and hence, preserve the photonic quantum states. In principle, with sufficiently long storage times and sufficiently high photonic storage capacity, quantum memories operating at widely separated telescopes can be brought together to achieve optical interferometry over arbitrarily long baselines.

Ashutosh Tripathi, Biao Zhou, Askar B. Abdikamalov, Dimitry Ayzenberg, Cosimo Bambi

8 pages, 5 figures

Einstein-Maxwell dilaton-axion gravity is a string-inspired model arising from the low energy effective action of heterotic string theory and an important candidate as alternative to General Relativity. Recently, some authors have explored its astrophysical implications in the spectra of accreting black holes and inferred the constraint $r_2 < 0.1$, where $r_2 \ge 0$ is the black hole dilaton charge and General Relativity is recovered for $r_2 = 0$. In the present paper, we study the impact of a non-vanishing black hole dilaton charge on the reflection spectrum of the disk. From the analysis of a NuSTAR spectrum of the black hole binary EXO 1846-031, we find the constraint $r_2 < 0.011$ (90% CL), which is an order of magnitude more stringent.

Hiroya Yamaguchi, Fabio Acero, Chuan-Jui Li, You-Hua Chu

12 pages, 12 figures, accepted for publication in ApJ Letters (originally submitted to ApJ, but offered a transfer to ApJL)

The geometric structure of supernova remnants (SNR) provides a clue to unveiling the pre-explosion evolution of their progenitors. Here we present an X-ray study of N103B (0509-68.7), a Type Ia SNR in the Large Magellanic Cloud, that is known to be interacting with dense circumstellar matter (CSM). Applying our novel method for feature extraction to deep Chandra observations, we have successfully resolved the CSM, Fe-rich ejecta, and intermediate-mass element (IME) ejecta components, and revealed each of their spatial distribution. Remarkably, the IME ejecta component exhibits a double-ring structure, implying that the SNR expands into an hourglass-shape cavity and thus forms bipolar bubbles of the ejecta. This interpretation is supported by more quantitative spectroscopy that reveals a clear bimodality in the distribution of the ionization state of the IME ejecta. These observational results can be naturally explained if the progenitor binary system had formed a dense CSM torus on the orbital plane prior to the explosion, providing further evidence that the SNR N103B originates from a single-degenerate progenitor.

Satoru Katsuda, Hitoshi Fujiwara, Yoshitaka Ishisaki, Yoshitomo Maeda, Koji Mori, Yuko Motizuki, Kosuke Sato, Makoto S. Tashiro, Yukikatsu Terada

29 pages, 11 figures, 4 tables, accepted for publication in Journal of Geophysical Research - Space Physics

We present new measurements of the vertical density profile of the Earth's atmosphere at altitudes between 70 and 200 km, based on Earth occultations of the Crab Nebula observed with the X-ray Imaging Spectrometer onboard Suzaku and the Hard X-ray Imager onboard Hitomi. X-ray spectral variation due to the atmospheric absorption is used to derive tangential column densities of the absorbing species, i.e., N and O including atoms and molecules, along the line of sight. The tangential column densities are then inverted to obtain the atmospheric number density. The data from 219 occultation scans at low latitudes in both hemispheres from September 15, 2005 to March 26, 2016 are analyzed to generate a single, highly-averaged (in both space and time) vertical density profile. The density profile is in good agreement with the NRLMSISE-00 model, except for the altitude range of 70-110 km, where the measured density is about 50% smaller than the model. Such a deviation is consistent with the recent measurement with the SABER aboard the TIMED satellite (Cheng et al. 2020). Given that the NRLMSISE-00 model was constructed some time ago, the density decline could be due to the radiative cooling/contracting of the upper atmosphere as a result of greenhouse warming in the troposphere. However, we cannot rule out a possibility that the NRL model is simply imperfect in this region. We also present future prospects for the upcoming Japan-US X-ray astronomy satellite, XRISM, which will allow us to measure atmospheric composition with unprecedented spectral resolution of dE ~ 5 eV in 0.3-12 keV.

Guilherme Limberg, Rafael M. Santucci, Silvia Rossi, Derek Shank, Vinicius M. Placco, Timothy C. Beers, Kevin C. Schlaufman, Andrew R. Casey, Hélio D. Perottoni, Young Sun Lee

Accepted for ApJ. Long tables at the end of the manuscript

We present the results of spectroscopic follow-up for 1897 low-metallicity star candidates, selected from the Best & Brightest (B&B) Survey, carried out with the GMOS-N/S (Gemini North/South telescopes) and Goodman (SOAR Telescope) spectrographs. From these low-resolution ($R \sim 2000$) spectra, we estimate stellar atmospheric parameters, as well as carbon and magnesium (representative of $\alpha$ elements) abundance ratios. We confirm that $56\%$ of our program stars are metal-poor ([Fe/H] $< -1.0$), $30\%$ are very metal-poor (VMP; [Fe/H] $< -2.0$) and $2\%$ are extremely metal-poor (EMP; [Fe/H] $< -3.0$). There are 191 carbon-enhanced metal-poor (CEMP) stars, resulting in CEMP fractions of $19\%$ and $43\%$ for the VMP and EMP regimes, respectively. A total of 94 confirmed CEMP stars belong to Group I ($A({\rm C}) \gtrsim 7.25$) and 97 to Group II ($A({\rm C}) \lesssim 7.25$) in the Yoon-Beers $A$(C)$-$[Fe/H] diagram. Moreover, we combine these data with Gaia EDR3 astrometric information to delineate new target-selection criteria, which have been applied to the Goodman/SOAR candidates, to more than double the efficiency for identification of bona-fide VMP and EMP stars in comparison to random draws from the B&B catalog. We demonstrate that this target-selection approach can achieve success rates of $96\%$, $76\%$, $28\%$ and $4\%$ for [Fe/H] $\leq -1.5$, $\leq -2.0$, $\leq -2.5$ and $\leq -3.0$, respectively. Finally, we investigate the presence of dynamically interesting stars in our sample. We find that several VMP/EMP ([Fe/H] $\leq -2.5$) stars can be associated with either the disk system or halo substructures like Gaia-Sausage/Enceladus and Sequoia.

Masanobu Kunitomo, Shigeru Ida, Taku Takeuchi, Olja Panić, James M. Miley, Takeru K. Suzuki

21 pages, 15 figures, published in ApJ. Evolutionary models of young 0.5-5 Msun stars are available at this https URL

We aim to understand the effect of stellar evolution on the evolution of protoplanetary disks. We focus in particular on the disk evolution around intermediate-mass (IM) stars, which evolve more rapidly than low-mass ones. We numerically solve the long-term evolution of disks around 0.5-5 solar-mass stars considering viscous accretion and photoevaporation (PE) driven by stellar far-ultraviolet (FUV), extreme-ultraviolet (EUV), and X-ray emission. We also take stellar evolution into account and consider the time evolution of the PE rate. We find that the FUV, EUV, and X-ray luminosities of IM stars evolve by orders of magnitude within a few Myr along with the time evolution of stellar structure, stellar effective temperature, or accretion rate. Therefore, the PE rate also evolves with time by orders of magnitude, and we conclude that stellar evolution is crucial for the disk evolution around IM stars.

Filip Morawski, Michał Bejger, Elena Cuoco, Luigia Petre

Submitted to the Machine Learning: Science and Technology; 29 pages, 18 figures

As of this moment, fifty gravitational waves (GW) detections have been announced, thanks to the observational efforts of the LIGO-Virgo Collaboration, working with the Advanced LIGO and the Advanced Virgo interferometers. The detection of signals is complicated by the noise-dominated nature of the data. Conventional approaches in GW detection procedures require either precise knowledge of the GW waveform in the context of matched filtering searches or coincident analysis of data from multiple detectors. Furthermore, the analysis is prone to contamination by instrumental or environmental artifacts called glitches which either mimic astrophysical signals or reduce the overall quality of data. In this paper, we propose an alternative generic method of studying GW data based on detecting anomalies. The anomalies we study are transient signals, different from the slow non-stationary noise of the detector. Presented in the manuscript anomalies are mostly based on the GW emitted by the mergers of binary black hole systems. However, the presented study of anomalies is not limited only to GW alone, but also includes glitches occurring in the real LIGO/Virgo dataset available at the Gravitational Waves Open Science Center.

Jian-Min Wang, Jun-Rong Liu, Luis C. Ho, Pu Du

10 pages, 3 figures, accepted

Compact objects are expected to exist in the accretion disks of supermassive black holes (SMBHs) in active galactic nuclei (AGNs), and in the presence of such a dense environment ($\sim 10^{14}\,{\rm cm^{-3}}$), they will form Thorne-\.Zytkow objects (TZOs). This hypothesis is supported by recent LIGO/Virgo detection of the mergers of very high-mass stellar binary black holes (BHs). We show that the TZOs will be trapped by the SMBH-disk within a typical AGN lifetime. In the context of SMBH-disks, the rates of Bondi accretion onto BHs are $\sim 10^{9}L_{\rm Edd}/c^{2}$, where $L_{\rm Edd}$ is the Eddington luminosity and $c$ is the speed of light. Outflows developed from the hyper-Eddington accretion strongly impact the Bondi sphere and induce episodic accretion. We show that the hyper-Eddington accretion will be halted after an accretion interval of $t_{\rm a}\sim 10^{5}m_{1}\,$s, where $m_{1}=\bhm/10\sunm$ is the BH mass. The kinetic energy of the outflows accumulated during $t_{\rm a}$ is equivalent to 10 supernovae driving an explosion of the Bondi sphere and developing blast waves. We demonstrate that a synchrotron flare from relativistic electrons accelerated by the blast waves peaks in the soft X-ray band ($\sim 0.1\,$keV), significantly contributing to the radio, optical, UV, and soft X-ray emission of typical radio-quiet quasars. External inverse Compton scattering of the electrons peaks around $40\,$GeV and is detectable through {\it Fermi}-LAT. The flare, decaying with $t^{-6/5}$ with a few months, will appear as a slowly varying transient. The flares, occurring at a rate of a few per year in radio-quiet quasars, provide a new mechanism for explaining AGN variability.

Edwin S. Kite, Laura Schaefer

Astrophysical Journal Letters, in press

Data suggest that most rocky exoplanets with orbital period $p$ $<$ 100 d ("hot" rocky exoplanets) formed as gas-rich sub-Neptunes that subsequently lost most of their envelopes, but whether these rocky exoplanets still have atmospheres is unknown. We identify a pathway by which 1-1.7 $R_{Earth}$ (1-10 $M_{Earth}$) rocky exoplanets with orbital periods of 10-100 days can acquire long-lived 10-2000 bar atmospheres that are H$_2$O-dominated, with mean molecular weight $>$10. These atmospheres form during the planets' evolution from sub-Neptunes into rocky exoplanets. H$_2$O that is made by reduction of iron oxides in the silicate magma is highly soluble in the magma, forming a dissolved reservoir that is protected from loss so long as the H$_2$-dominated atmosphere persists. The large size of the dissolved reservoir buffers the H$_2$O atmosphere against loss after the H$_2$ has dispersed. Within our model, a long-lived, water-dominated atmosphere is a common outcome for efficient interaction between a nebula-derived atmosphere (peak atmosphere mass fraction 0.1-0.6 wt%) and oxidized magma ($>$5 wt% FeO), followed by atmospheric loss. This idea predicts that most rocky planets that have orbital periods of 10-100 days and that have radii within 0.1-0.2 $R_{Earth}$ of the lower edge of the radius valley still retain H$_2$O atmospheres. This prediction is imminently testable with JWST and has implications for the interpretation of data for transiting super-Earths.

Stefan J. Kautsch, Dmitry Bizyaev, Dimitry I. Makarov, Vladimir P. Reshetnikov, Alexander V. Mosenkov, Alexandra V. Antipova

Four pages, one figure. Published in 2021, Research Notes of the Astronomical Society, 5, 43

We present spectroscopic observations of superthin galaxies. Superthin galaxies have the thinnest stellar disks among disk galaxies. A sample of 138 superthins was observed in visible light with the 3.5 m telescope at Apache Point Observatory in New Mexico to obtain the rotation curves of the ionized gas in the galaxies. The sample represents the largest survey of superthin galaxies so far and provides a database to investigate the kinematic and dynamic properties of this special type of extragalactic objects. Here we present the rotation curves of our sample objects.

Luciano Nicastro, Cristiano Guidorzi, Eliana Palazzi, Luca Zampieri, Massimo Turatto, Angela Gardini

47 pages, 8 figures, Accepted for publication on Universe

The origin and phenomenology of the Fast Radio Burst (FRB) remains unknown in spite of more than a decade of efforts (see arXiv:astro-ph.HE/1804.09092, arXiv:astro-ph.HE/1904.07947, arXiv:astro-ph.HE/1906.05878, arXiv:astro-ph.HE/2011.03500, arXiv:astro-ph.HE/2101.04907 for a review). Though several models have been proposed to explain the observed data, none is able to explain alone the variety of events so far recorded. The leading models consider magnetars as potential FRB sources. The recent detection of FRBs from the galactic magnetar SGR J1935+2154 seems to support them. Still emission duration and energetic budget challenge all these models. Like for other classes of objects initially detected in a single band, it appeared clear that any solution to the FRB enigma could only come from a coordinated observational and theoretical effort in a as wider as possible energy band. In particular the detection and localisation of optical/NIR or/and high-energy counterparts seemed an unavoidable starting point that could shed light on the FRB physics. Multiwavelength (MWL) search campaigns were conducted for several FRBs, in particular for repeaters. Here we summarize the observational and theoretical results and the perspectives in view of the several new sources accurately localised that will likely be identified by various radio facilities worldwide. We conclude that more dedicated MWL campaigns sensitive to the millisecond--minute timescale transients are needed to address the various aspects involved in the identification of FRB counterparts. Dedicated instrumentation could be one of the key points in this respect. In the optical/NIR band, fast photometry looks to be the only viable strategy. Additionally, small/medium size radiotelescopes co-pointing higher energies telescopes look a very interesting and cheap complementary observational strategy.

I. Heywood, E. J. Murphy, E. F. Jiménez-Andrade, L. Armus, W. D. Cotton, C. DeCoursey, M. Dickinson, T. J. W. Lazio, E. Momjian, K. Penner, I. Smail, O. M. Smirnov

29 pages, 13 figures, 6 tables, accepted for publication in The Astrophysical Journal

The Frontier Fields project is an observational campaign targeting six galaxy clusters, with the intention of using the magnification provided by gravitational lensing to study galaxies that are extremely faint or distant. We used the Karl G. Jansky Very Large Array (VLA) at 3 and 6 GHz to observe three Frontier Fields: MACSJ0416.1$-$2403 ($z$ = 0.396), MACSJ0717.5+3745 ($z$ = 0.545), and MACSJ1149.5+2223 ($z$ = 0.543). The images reach noise levels of $\sim$1 $\mu$Jy beam$^{-1}$ with sub-arcsecond resolution ($\sim$2.5 kpc at $z$ = 3), providing a high-resolution view of high-$z$ star-forming galaxies that is unbiased by dust obscuration. We generate dual-frequency continuum images at two different resolutions per band, per cluster, and derive catalogs totalling 1966 compact radio sources. Components within the areas of Hubble Space Telescope and Subaru observations are cross-matched, providing host galaxy identifications for 1296 of them. We detect 13 moderately-lensed (2.1 $<$ $\mu$ $<$ 6.5) sources, one of which has a demagnified peak brightness of 0.9 $\mu$Jy beam$^{-1}$, making it a candidate for the faintest radio source ever detected. There are 66 radio sources exhibiting complex morphologies, and 58 of these have host galaxy identifications. We reveal that MACSJ1149.5+2223 is not a cluster with a double relic, as the western candidate relic is resolved as a double-lobed radio galaxy associated with a foreground elliptical at $z$ = 0.24. The VLA Frontier Fields project is a public legacy survey. The image and catalog products from this work are freely available.

E. F. Jiménez-Andrade, E. J. Murphy, I. Heywood, I. Smail, K. Penner, E. Momjian, M. Dickinson, L. Armus, T. J. W. Lazio

30 pages, 15 figures, 2 tables, accepted for publication in The Astrophysical Journal

To investigate the growth history of galaxies, we measure the rest-frame radio, ultraviolet (UV), and optical sizes of 98 radio-selected, star-forming galaxies (SFGs) distributed over $0.3 \lesssim z \lesssim 3$ and median stellar mass of $\log(M_\star/ \rm M_\odot)\approx10.4$. We compare the size of galaxy stellar disks, traced by rest-frame optical emission, relative to the overall extent of star formation activity that is traced by radio continuum emission. Galaxies in our sample are identified in three Hubble Frontier Fields: MACSJ0416.1$-$2403, MACSJ0717.5+3745, and MACSJ1149.5+2223. Radio continuum sizes are derived from 3 GHz and 6 GHz radio images ($\lesssim 0$''$.6$ resolution, $\approx0.9\, \rm \mu Jy\, beam^{-1}$ noise level) from the Karl G. Jansky Very Large Array. Rest-frame UV and optical sizes are derived using observations from the Hubble Space Telescope and the ACS and WFC3 instruments. We find no clear dependence between the 3 GHz radio size and stellar mass of SFGs, which contrasts with the positive correlation between the UV/optical size and stellar mass of galaxies. Focusing on SFGs with $\log(M_\star/\rm M_\odot)>10$, we find that the radio/UV/optical emission tends to be more compact in galaxies with high star-formation rates ($\rm SFR\gtrsim 100\,M_\odot\,yr^{-1}$), suggesting that a central, compact starburst (and/or an Active Galactic Nucleus) resides in the most luminous galaxies of our sample. We also find that the physical radio/UV/optical size of radio-selected SFGs with $\log(M_\star/\rm M_\odot)>10$ increases by a factor of $1.5-2$ from $z\approx 3$ to $z\approx0.3$, yet the radio emission remains two-to-three times more compact than that from the UV/optical. These findings indicate that these massive, {radio-selected} SFGs at $0.3 \lesssim z \lesssim 3$ tend to harbor centrally enhanced star formation activity relative to their outer-disks.

Abigail R. Harden, Navdeep K. Panesar, Ronald L. Moore, Alphonse C. Sterling, Mitzi L. Adams

Using EUV images and line-of-sight magnetograms from Solar Dynamics Observatory, we examine eight emerging bipolar magnetic regions (BMRs) in central-disk coronal holes for whether the emerging magnetic arch made any noticeable coronal jets directly, via reconnection with ambient open field as modeled by Yokoyama and Shibata (1995). During emergence, each BMR produced no obvious EUV coronal jet of normal brightness, but each produced one or more faint EUV coronal jets that are discernible in AIA 193 {\AA} images. The spires of these jets are much fainter and usually narrower than for typical EUV jets that have been observed to be produced by minifilament eruptions in quiet regions and coronal holes. For each of 26 faint jets from the eight emerging BMRs, we examine whether the faint spire was evidently made a la Yokoyama and Shibata (1995). We find: (1) 16 of these faint spires evidently originate from sites of converging opposite-polarity magnetic flux and show base brightenings like those in minifilament-eruption-driven coronal jets, (2) the 10 other faint spires maybe were made by a burst of the external-magnetic-arcade-building reconnection of the emerging magnetic arch with the ambient open field, reconnection directly driven by the arch's emergence, but (3) none were unambiguously made by such emergence-driven reconnection. Thus, for these eight emerging BMRs the observations indicate that emergence-driven external reconnection of the emerging magnetic arch with ambient open field at most produces a jet spire that is much fainter than in previously-reported, much more obvious coronal jets driven by minifilament eruptions.

We first report GeV $\gamma$-ray emission from SNR G15.9+0.2 in this work. We find that its power-law spectral index is 2.13$\pm$0.05 with 13.29$\sigma$ significance level, and the $\gamma$-ray emission can be characterized by a 2D Gauss spatial distribution, which has a better improvement than the case of a point source. Moreover, we find that its likely counterparts from radio, X-ray, and TeV energy bands are well coincident with its spatial location. Analyzing the variability from 12.4 years of the light curve (LC), we find that this LC exists weak variability with a 3.30$\sigma$ variability significance level. Investigating the 2$\sigma$ error region of its best-fit position, we do not find certified active galactic nuclei (AGNs) and AGNs' candidates from the region of this SNR, so we suggest that the new $\gamma$-ray emission is likely to originate from SNR G15.9+0.2. On this basis, we discussed the likely origins of its $\gamma$-ray radiation combined with the distribution of surrounding molecular clouds.

Yilin Guo, Yijun Hou, Ting Li, Jun Zhang

10 pages, 5 figures, 2 movies, accepted for publication in ApJL

Bubbles, the semi-circular voids below quiescent prominences (filaments), have been extensively investigated in the past decade. However, hitherto the magnetic nature of bubbles cannot be verified due to the lack of on-disk photospheric magnetic field observations. Here for the first time, we find and investigate an on-disk prominence bubble around a filament barb on 2019 March 18 based on stereoscopic observations from NVST, SDO, and STEREO-A. In high-resolution NVST H$\alpha$ images, this bubble has a sharp arch-like boundary and a projected width of $\thicksim$26 Mm. Combining SDO and STEREO-A images, we further reconstruct 3D structure of the bubble boundary, whose maximum height is $\thicksim$15.6 Mm. The squashing factor Q map deduced from extrapolated 3D magnetic fields around the bubble depicts a distinct arch-shaped interface with a height of $\thicksim$11 Mm, which agrees well with the reconstructed 3D structure of the observed bubble boundary. Under the interface lies a set of magnetic loops, which is rooted on a surrounding photospheric magnetic patch. To be more persuasive, another on-disk bubble on 2019 June 10 is presented as a supplement. According to these results obtained from on-disk bubble observations, we suggest that the bubble boundary corresponds to the interface between the prominence dips (barb) and the underlying magnetic loops rooted nearby. It is thus reasonable to speculate that the bubble can form around a filament barb below which there is a photospheric magnetic patch.

Ealeal Bear, Noam Soker (Technion, Israel)

Will be summitted in two days to allow comments

We study a white dwarf (WD) - neutron star (NS) reverse evolution that might lead to the rare explosion of the WD as a type Ia peculiar supernova (peculiar SN Ia) few months to several years before the core, which is a remnant of a red supergiant (RSG) star, explodes as a core collapse supernova (CCSN). Using the evolutionary code \textsc{mesa-binary} we simulate evolution of binary systems with stars of initial masses of 6-7.5 Mo. The more massive star, the primary, transfers mass to the secondary star and leaves a CO WD remnant of mass ~ 1 Mo. The secondary becomes massive enough to end in a CCSN. As the secondary evolves to the RSG phase it engulfs the WD and the system experience a common envelope evolution that ends with a WD-core binary system at an orbital separation of a_f ~ 1-5 Ro. Our simulations show that the core explodes as a CCSN at t_{CEE-CCSN} ~ 3000 - 10^5 yr after the CEE. The rest of the scenario is speculative. We assume that if the WD accretes helium-rich gas from the core it might explode as a SN Ia in the frame of the double detonation scenario for SNe Ia and peculiar SNe Ia. The peculiar SN Ia explosion might occur at the end of the CEE evolution or just few years before the CCSN. We predict the very rare occurrence of a peculiar SN Ia followed within months to years by a CCSN.

Llibert Aresté Saló, David Benisty, Eduardo I. Guendelman, Jaume de Haro

11 pages; 4 figures. arXiv admin note: substantial text overlap with arXiv:2102.09514

An exponential kind of quintessential inflation potential motivated by supergravity is studied. This type belongs to the class of $\alpha$-attractor models. The model was studied for the first time in [1] where the authors introduced a Cosmological Constant in order to ensure quintessence at late times. However, in this note we disregard this Cosmological Constant, showing that the obtained results are very close to the ones obtained recently in the context of Lorenzian Quintessential Inflation, and thus, depicting with great accuracy the early and late time acceleration of our universe. The model is compatible with the recent observations. Finally, we review the treatment of $\alpha$-attractor and we show that other forms of the potential does not depict the late time accelerated universe as the simple model.

Michael Shull, Jeremy Darling, Charles Danforth (Univ. of Colorado)

19 pages, 2 figures, 3 tables (submitted to Astrophysical Journal)

Using offset-corrected Gaia-EDR3 parallax measurements and spectrophotometric methods, we have determined distances for 69 massive stars in the Carina OB1 association and associated clusters: Trumpler 16 (21 stars), Trumpler 14 (20 stars), Trumpler 15 (3 stars), Bochum 11 (5 stars), and South Pillars region (20 stars). Past distance estimates to the Carina Nebula range from 2.2 to 3.6 kpc, with uncertainties arising from photometry and anomalous dust extinction. The EDR3 parallax solutions show considerable improvement over DR2, with typical errors $\sigma_{\varpi}/\varpi \approx$~3-5%. The O-type stars in the Great Carina Nebula lie at essentially the same distance ($2.35\pm0.08$ kpc), quoting mean and rms variance. The clusters have distances of $2.32\pm0.12$ kpc (Tr 16), $2.37\pm0.15$ kpc (Tr 14), $2.36\pm0.09$ kpc (Tr 15), and $2.33\pm0.12$ kpc (Bochum 11) in good agreement with the $\eta$ Car distance of around 2.3 kpc. O-star proper motions suggest internal (2D) velocity dispersions $\sim4$ km/s for Tr 14 and Tr 16. Reliable distances allow estimates of cluster sizes, stellar dynamics, luminosities, and fluxes of photoionizing radiation incident on photodissociation regions in the region. We estimate that Tr 14 and Tr 16 have half-mass radii $r_h = 1.5-1.8$ pc, stellar crossing times $t_{\rm cr} = r_h/v_m \approx 0.7-0.8$ Myr, and two-body relaxation times $t_{rh} \approx 40-80$ Myr. The underlying velocity dispersion for Tr 14, if a bound cluster, would be $v_m \approx 2.1^{+0.7}_{-0.4}$ km/s for $N = 7600^{+5800}_{-2600}$ stars. With the higher dispersions of the O-stars, mass segregation might occur slowly, on times scales of 3-6~Myr.

S. E. Woosley (1), Alexander Heger (2 and 3 and 4 and 5) ((1) UCSC, (2) Monash, (3) OzGrav, (4) ASTRO-3D, (5) JINA-CEE)

16 pages, 4 figure panels, 6 tables, submitted to ApJ

Stellar evolution theory predicts a "gap" in the black hole birth function caused by the pair instability. Presupernova stars that have a core mass below some limiting value, Mlo, after all pulsational activity is finished, collapse to black holes, whereas more massive ones, up to some limiting value, Mhi, explode, promptly and completely, as pair-instability supernovae. Previous work has suggested Mlo is approximately 50 solar masses and Mhi is approximately 130 solar masses. These calculations have been challenged by recent LIGO observations that show many black holes merging with individual masses, Mlo is least some 65 solar masses. Here we explore four factors affecting the theoretical estimates for the boundaries of this mass gap: nuclear reaction rates, evolution in detached binaries, rotation, and hyper-Eddington accretion after black hole birth. Current uncertainties in reaction rates by themselves allow Mlo to rise to 64 solar masses and Mhi as large as 161 solar masses. Rapid rotation could further increase Mlo to about 70 solar masses, depending on the treatment of magnetic torques. Evolution in detached binaries and super-Eddington accretion can, with great uncertainty, increase Mlo still further. Dimensionless Kerr parameters close to unity are allowed for the more massive black holes produced in close binaries, though they are generally smaller.

Rahul Datta, David T. Chuss, Joseph Eimer, Thomas Essinger-Hileman, Natalie N. Gandilo, Kyle Helson, Alan J. Kogut, Luke Lowe, Paul Mirel, Karwan Rostem, Marco Sagliocca, Danielle Sponseller, Eric R. Switzer, Peter A. Taraschi, Edward J. Wollack

Measuring the faint polarization signal of the cosmic microwave background (CMB) not only requires high optical throughput and instrument sensitivity but also control over systematic effects. Polarimetric cameras or receivers used in this setting often employ dielectric vacuum windows, filters, or lenses to appropriately prepare light for detection by cooled sensor arrays. These elements in the optical chain are typically designed to minimize reflective losses and hence improve sensitivity while minimizing potential imaging artifacts such as glint and ghosting. The Primordial Inflation Polarization ExploreR (PIPER) is a balloon-borne instrument designed to measure the polarization of the CMB radiation at the largest angular scales and characterize astrophysical dust foregrounds. PIPER's twin telescopes and detector systems are submerged in an open-aperture liquid helium bucket dewar. A fused-silica window anti-reflection (AR) coated with polytetrafluoroethylene (PTFE) is installed on the vacuum cryostat that houses the cryogenic detector arrays. Light passes from the skyward portions of the telescope to the detector arrays though this window, which utilizes an indium seal to prevent superfluid helium leaks into the vacuum cryostat volume. The AR coating implemented reduces reflections from each interface to <1% compared to ~10% from an uncoated window surface. The AR coating procedure and room temperature optical measurements of the window are presented. The indium vacuum sealing process is also described in detail and test results characterizing its integrity to superfluid helium leaks are provided.

Tien-Hao Hsieh, Michihiro Takami, Michael S. Connelley, Sheng-Yuan Liu, Yu-Nung Su, Naomi Hirano, Motohide Tamura, Masaaki Otsuka, Jennifer L. Karr, Tae-Soo Pyo

18 page, 7 figures, accepted to ApJ

A classical paradox in high-mass star formation is that powerful radiation pressure can halt accretion, preventing further growth of a central star. Disk accretion has been proposed to solve this problem, but the disks and the accretion process in high-mass star formation are poorly understood. We executed high-resolution ($R$=35,000-70,000) iSHELL spectroscopy in $K$-band for eleven high-mass protostars. Br-$\gamma$ emission was observed toward eight sources, and the line profiles for most of these sources are similar to those of low-mass PMS stars. Using an empirical relationship between the Br-$\gamma$ and accretion luminosities, we tentatively estimate disk accretion rates ranging from $\lesssim$10$^{-8}$ and $\sim$10$^{-4}$ $M_\odot$ yr$^{-1}$. These low-mass-accretion rates suggest that high-mass protostars gain more mass via episodic accretion as proposed for low-mass protostars. Given the detection limits, CO overtone emission ($v$=2-0 and 3-1), likely associated with the inner disk region ($r \ll 100$ au), was found towards two sources. This low-detection rate compared with Br-$\gamma$ emission is consistent with previous observations. Ten out of the eleven sources show absorption at the $v$=0-2 ${\rm R(7)-R(14)}$ CO R-branch. Most of them are either blueshifted or redshifted, indicating that the absorption is associated with an outflow or an inflow with a velocity of up to $\sim50$ km s$^{-1}$. Our analysis indicates that the absorption layer is well thermalized (and therefore $n_{\mathrm H_2} \gtrsim 10^6$ cm$^{-3}$) at a single temperature of typically 100-200 K, and located within 200-600 au of the star.

Taylor Jacovich, Daniel Patnaude, Pat Slane, Carles Badenes, Shiu-Hang Lee, Shigehiro Nagataki, Dan Milisavljevic

19 pages, 25 figures, 1 table, submitted to ApJ

Massive stars can shed material via steady, line-driven winds, eruptive outflows, or mass-transfer onto a binary companion. In the case of single stars, the mass is deposited by the stellar wind into the nearby environment. After the massive star explodes, the stellar ejecta interact with this circumstellar material (CSM), often-times resulting in bright X-ray line emission from both the shock-heated CSM and ejecta. The amount of material lost by the progenitor, the mass of ejecta, and its energetics all impact the bulk spectral characteristics of this X-ray emission. Here we present a grid of core-collapse supernova remnant models derived from models for massive stars with zero age main sequence masses of $\sim$ 10 - 30 M$_\odot$ evolved from the pre-main sequence stage with wind-driven mass-loss. Evolution is handled by a multi-stage pipeline of software packages. First, we use mesa (Modules for Experiments in Stellar Astrophysics) to evolve the progenitors from pre-main sequence to iron core collapse. We then use the Supernova Explosion Code (snec) to explode the mesa models, and follow them for the first 100 days following core-collapse. Finally, we couple the snec output, along with the CSM generated from mesa mass-loss rates, into the Cosmic-Ray Hydrodynamics code (ChN) to model the remnant phase to 7000 years post core-collapse. At the end of each stage, we compare our outputs with those found in the literature, and we examine any qualitative and quantitative differences in the bulk properties of the remnants and their spectra based on the initial progenitor mass, as well as mass-loss history.

Andreas Efstathiou, Katarzyna Malek, Denis Burgarella, Peter Hurley, Seb Oliver, Veronique Buat, Raphael Shirley, Steven Duivenvoorden, Vicky Papadopoulou Lesta, Duncan Farrah, Kenneth J. Duncan, Maria del Carmen Campos Varillas

Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters. Models used in this work are available at this http URL

In this work we report the discovery of the hyperluminous galaxy HELP_J100156.75+022344.7 at the photometric redshift of z ~ 4.3. The galaxy was discovered in the Cosmological Evolution Survey (COSMOS) field, one of the fields studied by the Herschel Extragalactic Legacy Project (HELP). We present the spectral energy distribution (SED) of the galaxy and fit it with the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) multi-component radiative transfer models. We find that its emission is dominated by an obscured quasar with a predicted total 1-1000um luminosity of $3.91^{+1.69}_{-0.55} \times 10^{13} L_\odot$ and an active galactic nucleus (AGN) fraction of ~89%. We also fit HELP_J100156.75+022344.7 with the Code Investigating GALaxy Emission (CIGALE) code and find a similar result. This is only the second z > 4 hyperluminous obscured quasar discovered to date. The discovery of HELP_J100156.75+022344.7 in the ~ 2deg^2 COSMOS field implies that a large number of obscured hyperluminous quasars may lie in the HELP fields which cover ~ 1300deg^2. If this is confirmed, tension between supermassive black hole evolution models and observations will be alleviated. We estimate the space density of objects like HELP_J100156.75+022344.7 at z ~ 4.5 to be $\sim 1.8 \times 10^{-8}$Mpc$^{-3}$. This is slightly higher than the space density of coeval hyperluminous optically selected quasars suggesting that the obscuring torus in z > 4 quasars may have a covering factor $\gtrsim 50\%$.

Elena Petrova, Tatiana Podladchikova, Astrid M. Veronig, Stijn Lemmens, Benjamin Bastida Virgili, Tim Flohrer

28 pages, 15 figures, accepted for publication in the Astrophysical Journal Supplement Series

The solar radio flux at F10.7 cm and F30 cm is required by most models characterizing the state of the Earth's upper atmosphere, such as the thermosphere and ionosphere to specify satellite orbits, re-entry services, collision avoidance maneuvers and modeling of space debris evolution. We develop a method called RESONANCE ("Radio Emissions from the Sun: ONline ANalytical Computer-aided Estimator") for the prediction of the 13-month smoothed monthly mean F10.7 and F30 indices 1-24 months ahead. The prediction algorithm includes three steps. First, we apply a 13-month optimized running mean technique to effectively reduce the noise in the radio flux data. Second, we provide initial predictions of the F10.7 and F30 indices using the McNish-Lincoln method. Finally, we improve these initial predictions by developing an adaptive Kalman filter with the error statistics identification. The root-mean-square-error of predictions with lead times from 1 to 24 months is 5-27 sfu for the F10.7 and 3-16 sfu for F30 index, which statistically outperforms current algorithms in use. The proposed approach based on Kalman filter is universal and can be applied to improve the initial predictions of a process under study provided by any other forecasting method. Furthermore, we present a systematic evaluation of re-entry forecast as an application to test the performance of F10.7 predictions on past ESA re-entry campaigns for payloads, rocket bodies, and space debris that re-entered from June 2006 to June 2019. The test results demonstrate that the predictions obtained by RESONANCE in general also lead to improvements in the forecasts of re-entry epochs.

Pascale Garaud

Proceedings of the summer school "Multi-Dimensional Processes In Stellar Physics". Edited by Michel Rieutord, Isabelle Baraffe and Yveline Lebreton

The past 20 years have witnessed a renewal of interest in the subject of double-diffusive processes in astrophysics, and their impact on stellar evolution. This lecture aims to summarize the state of the field as of early 2019, although the reader should bear in mind that it is rapidly evolving. An Annual Review of Fluid Mechanics article entitled "Double-diffusive convection at low Prandtl number" (Garaud, 2018) contains a reasonably comprehensive review of the topic, up to the summer of 2017. I focus here on presenting what I hope are clear derivations of some of the most important results with an astrophysical audience in mind, and discuss their implications for stellar evolution, both in an observational context, and in relation to previous work on the subject.

Yanxia Zhang, Yongheng Zhao, Xue-Bing Wu

11 pages, 4 figures, 10 tables, accepted by Monthly Notices of the Royal Astronomical Society

The ESA's X-ray Multi-Mirror Mission (XMM-Newton) created a new, high quality version of the XMM-Newton serendipitous source catalogue, 4XMM-DR9, which provides a wealth of information for observed sources. The 4XMM-DR9 catalogue is correlated with the Sloan Digital Sky Survey (SDSS) DR12 photometric database and the ALLWISE database, then we get the X-ray sources with information from X-ray, optical and/or infrared bands, and obtain the XMM-WISE sample, the XMM-SDSS sample and the XMM-WISE-SDSS sample. Based on the large spectroscopic surveys of SDSS and the Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST), we cross-match the XMM-WISE-SDSS sample with those sources of known spectral classes, and obtain the known samples of stars, galaxies and quasars. The distribution of stars, galaxies and quasars as well as all spectral classes of stars in 2-d parameter spaces is presented. Various machine learning methods are applied on different samples from different bands. The better classified results are retained. For the sample from X-ray band, rotation forest classifier performs the best. For the sample from X-ray and infrared bands, a random forest algorithm outperforms all other methods. For the samples from X-ray, optical and/or infrared bands, LogitBoost classifier shows its superiority. Thus, all X-ray sources in the 4XMM-DR9 catalogue with different input patterns are classified by their respective models which are created by these best methods. Their membership and membership probabilities to individual X-ray sources are assigned. The classified result will be of great value for the further research of X-ray sources in greater detail.

Maresuke Shiraishi, Kazuyuki Akitsu, Teppei Okumura

6 pages, 1 figure

The Alcock-Paczynski (AP) effect is a geometrical distortion in three-dimensional observed galaxy statistics. In anticipation of precision cosmology based on ongoing and upcoming all-sky galaxy surveys, we build an efficient method to compute the AP-distorted correlations of galaxy number density and peculiar velocity fields for any larger angular scale not relying on the conventionally used plane-parallel (PP) approximation. Here, instead of the usual Legendre polynomial basis, the correlation functions are decomposed using tripolar spherical harmonic basis; hence, characteristic angular dependence due to the wide-angle AP effect can rigorously captured. By means of this, we demonstrate the computation of the AP-distorted correlations over the various scales. Comparing our results with the PP-limit ones, we confirm that the errors due to the PP approximation become more remarkable as the visual angle of separation between target galaxies, $\Theta$, enlarges, and especially for the density auto correlation, the error exceeds $10\%$ when $\Theta \gtrsim 30^\circ$. This highlights the importance of the analysis beyond the PP approximation.

Spatial correlations between spiral arms and other galactic components such as giant molecular clouds and massive OB stars suggest that spiral arms can play vital roles in various aspects of disk galaxy evolution. Segmentation of spiral arms in disk galaxies is therefore a key task to investigate these correlations. We here try to decompose disk galaxies into spiral and non-spiral regions by using U-net, which is based on deep learning algorithms and has been invented for segmentation tasks in biology.

Kai-Xing Lu, Zhi-Xiang Zhang, Ying-Ke Huang, An-Bing Ren, Liang Xu, Hai-Cheng Feng, Yu-Xin Xin, Xu Ding, Xiao-Guang Yu, Jin-Ming Bai

15 pages, 10 figure, accepted for publication in RAA (Research in Astronomy and Astrophysics)

Observing a telluric standard star for correcting the telluric absorption lines of spectrum will take a significant amount of precious telescope time, especially in the long-term spectral monitoring project. Beyond that, it's difficult to select a suitable telluric standard star near in both time and airmass to the scientific object. In this paper, we present a method of correcting the telluric absorption lines by combining the advantages of long slit spectroscopy. By rotating the slit, we observed the scientific object and a nearby comparison star in one exposure, so that the spectra of both objects should have the same telluric transmission spectrum. The telluric transmission spectrum was constructed by dividing the observed spectrum of comparison star by its stellar template, and was used to correct the telluric absorption lines of the scientific object. Using the long slit spectrograph of Lijiang 2.4-meter telescope, we designed a long-term spectroscopic observation strategy, and finished a four-year spectroscopic monitoring for a pair of objects (an active galactic nuclei and an non-varying comparison star). We applied this method to correct the telluric absorption lines of the long-term monitored spectra by Lijiang 2.4-meter telescope, and investigated the variation of the telluric absorptions at Lijiang Observatory. We found that the telluric absorption transparency is mainly modulated by the seasonal variability of the relative humidity, airmass and seeing. Using the scatter of the [O~III] $\lambda$5007 fluxes emitted from the narrow-line region of active galactic nuclei as an indicator, we found that the correction accuracy of the telluric absorption lines is 1%.

Sunquakes are helioseismic power enhancements initiated by solar flares, but not all flares generate sunquakes. It is curious why some flares cause sunquakes while others do not. Here we propose a hypothesis to explain the disproportionate occurrence of sunquakes: during a flare's impulsive phase when the flare's impulse acts upon the photosphere, delivered by shock waves, energetic particles from higher atmosphere, or by downward Lorentz Force, a sunquake tends to occur if the background oscillation at the flare footpoint happens to oscillate downward in the same direction with the impulse from above. To verify this hypothesis, we select 60 strong flares in Solar Cycle 24, and examine the background oscillatory velocity at the sunquake sources during the flares' impulsive phases. Since the Doppler velocity observations at sunquake sources are usually corrupted during the flares, we reconstruct the oscillatory velocity in the flare sites using helioseismic holography method with an observation-based Green's function. A total of 24 flares are found to be sunquake active, giving a total of 41 sunquakes. It is also found that in 3-5 mHz frequency band, 25 out of 31 sunquakes show net downward oscillatory velocities during the flares' impulsive phases, and in 5-7 mHz frequency band, 33 out of 38 sunquakes show net downward velocities. These results support the hypothesis that a sunquake more likely occurs when a flare impacts a photospheric area with a downward background oscillation.

D. Minniti, V. Ripepi, J. G. Fernández-Trincado, J. Alonso-García, L. C. Smith, P. W. Lucas, M. Gómez, J. B. Pullen, E. R. Garro, F. Vivanco Cádiz, M. Hempel, M. Rejkuba, R. K. Saito, T. Palma, J. J. Clariá, M. Gregg, D. Majaess

12 pages, 6 Figures

Context. Globular clusters (GCs) are witnesses of the past accretion events onto the Milky Way (MW). In particular, the GCs of the Sagittarius (Sgr) dwarf galaxy are important probes of an on-going merger. Aims. Our main goal is to search for new GC members of this dwarf galaxy using the VISTA Variables in the Via Lactea Extended Survey (VVVX) near-infrared database combined with the Gaia Early Data Release 3 (EDR3) optical database. Methods. We investigated all VVVX-enabled discoveries of GC candidates in a region covering about 180 sq. deg. toward the bulge and the Sgr dwarf galaxy. We used multiband point-spread function photometry to obtain deep color-magnitude diagrams (CMDs) and luminosity functions (LFs) for all GC candidates, complemented by accurate Gaia-EDR3 proper motions (PMs) to select Sgr members and variability information to select RR Lyrae which are potential GC members. Results. After applying a strict PM cut to discard foreground bulge and disk stars, the CMDs and LFs for some of the GC candidates exhibit well defined red giant branches and red clump giant star peaks. We selected the best Sgr GCs, estimating their distances, reddenings, and associated RR Lyrae. Conclusions. We discover 12 new Sgr GC members, more than doubling the number of GCs known in this dwarf galaxy. In addition, there are 11 other GC candidates identified that are uncertain, awaiting better data for confirmation.

K.F. Cheng, X.H. Zhao, B.B. Zhang, J.M. Bai

16 pages, 5 figures, 1 table. Accepted for publication in RAA

The afterglow of GRB 170817A has been detected for more than three years, but the origin of the multi-band afterglow light curves remains under debate. A classical top-hat jet model is faced with difficulties in producing a shallow rise of the afterglow light curves as observed $(F_{\nu} \propto T^{0.8})$. Here we reconsider the model of stratified ejecta with energy profile of $E(>\Gamma \beta)=E_0(\Gamma \beta)^{-k}$ as the origin of the afterglow light curves of the burst, where $\Gamma$ and $\beta$ are the Lorentz factor and speed of the ejecta, respectively. $k$ is the power-law slope of the energy profile. We consider the ejecta are collimated into jets. Two kinds of jet evolutions are investigated, including a lateral-spreading jet and a non-lateral-spreading jet. We fit the multi-band afterglow light curves, including the X-ray data at one thousand days post-burst, and find that both the models of the spreading and non-spreading jets can fit the light curves well, but the observed angular size of the source and the apparent velocity of the flux centroid for the spreading jet model are beyond the observation limits, while the non-spreading jet model meets the observation limits. Some of the best-fit parameters for the non-spreading jet model, such as the number density of the circumburst medium $\sim10^{-2}$ cm$^{-3}$ and the total jet kinetic energy $E \sim 4.8\times 10^{51}$ erg, also appear plausible. The best-fit slope of the jet energy profile is $k \sim 7.1$. Our results suggest that the afterglow of GRB 170817A may arise from the stratified jet and that the lateral spreading of the jet is not significant.

G. Valle, M. Dell'Omodarme, E. Tognelli

Accepted for publication in A&A. Abstract abridged

[...] This paper presents a rigorous derivation of a goodness-of-fit statistics for colour-magnitude diagrams (CMD). We discuss the reliability of the underlying assumptions and their validity. We derived the distribution of the sum of squared Mahalanobis distances of stellar data and theoretical isochrone for a generic set of data and models. We applied this to the case of synthetic CMDs constructed to mimic real data of open clusters in the GAIA sample. We analysed the capability of distinguishing among different sets of input physics and parameters that were used to compute the stellar models. We generated synthetic clusters from isochrones computed with these perturbed quantities, and we evaluated the goodness-of-fit with respect to the unperturbed isochrone. We show that when $r$ magnitudes are available for each of the $N$ observational objects and $p$ hyperparameters are estimated in the fit, the error distribution follows a $\chi^2$ distribution with $(r-1)N - p$ degrees of freedom. We show that the linearisation of the isochrone causes negligible deviation from this result. We investigated the possibility of detecting the effects on stellar models that are induced when varying convective core overshooting efficiency, $^{14}$N$(p,\gamma)^{15}$O reaction rate, microscopic diffusion velocities, outer BCs, and colour transformation. The results suggest that it is possible to detect the effect induced by only some of the perturbed quantities. [...] A variation in the convective core overshooting efficiency was detectable only for photometric errors of 0.003 mag and only for the 1 Gyr case. The effects induced by the outer boundary conditions and the bolometric corrections are the largest. [...] As a last exercise, we addressed the validity of the goodness-of-fit statistics for real-world open cluster CMDs, contaminated by field stars or unresolved binaries. [...]

The motion of our solar system relative to the CMB rest frame leads to subtle distortions in the observed CMB sky map due to the aberration effect. Usually the corresponding peculiar velocity is determined from the CMB dipole but neglecting intrinsic dipole contributions. Here it is investigated whether certain invariant scalar measures, which are derived from first and second order covariant derivatives on the sphere, can detect the distortions caused by the aberration effect at high multipoles. This would in principle allow to disentangle the Doppler from intrinsic dipole contributions providing an independent method for the determination of our peculiar velocity. It is found that the eigenvalues of the Hessian matrix of the temperature field are well suited for that task.

I. Gallardo Cava, M. Gómez-Garrido, V. Bujarrabal, A. Castro-Carrizo, J. Alcolea, H. Van Winckel

Accepted 4 March 2021 in Astronomy & Astrophysics. 26 pages

There is a class of binary post-AGB stars with a remarkable near-infrared excess that are surrounded by Keplerian or quasi-Keplerian disks and extended outflows composed of gas escaping from the disk. The Keplerian dynamics had been well identified in four cases, namely the Red Rectangle, AC Her, IW Car, and IRAS 08544-4431. In these objects, the mass of the outflow represents ~ 10 % of the nebular mass, the disk being the dominant component of the nebula. We present interferometric NOEMA maps of 12CO and 13CO J=2-1 in 89 Her and 12CO J=2-1 in AC Her, IRAS 19125+0343, and R Sct. Several properties of the nebula are obtained from the data and model fitting, including the structure, density, and temperature distributions, as well as the dynamics. We also discuss the uncertainties on the derived values. The presence of an expanding component in AC Her is doubtful, but thanks to new maps and models, we estimate an upper limit to the mass of this outflow of < 3 10^-5 Mo, that is, the mass of the outflow is < 5 % of the total nebular mass. For 89 Her, we find a total nebular mass of 1.4 10^-2 Mo, of which ~ 50 % comes from an hourglass-shaped extended outflow. In the case of IRAS 19125+0343, the nebular mass is 1.1 10^-2 Mo, where the outflow contributes ~ 70 % of the total mass. The nebular mass of R Sct is 3.2 10^-2 Mo, of which ~ 75 % corresponds to a very extended outflow that surrounds the disk. Our results for IRAS 19125+0343 and R Sct lead us to introduce a new subclass of binary post-AGB stars, for which the outflow is the dominant component of the nebula. Moreover, the outflow mass fraction found in AC Her is smaller than those found in other disk-dominated binary post-AGB stars. 89 Her would represent an intermediate case between both subclasses.

S. W. Duchesne, M. Johnston-Hollitt, A. G. Wilber

17 pages, 31 individual figure files, submitted to PASA

We report on the detection of a giant radio halo in the cluster Abell 3404 as well as confirmation of the radio halo observed in Abell 141 (with linear extents $\sim 770$ kpc and $\sim 850$ kpc, respectively). We use the Murchison Widefield Array (MWA) in conjunction with the Australian Square Kilometre Array Pathfinder (ASKAP) and the Australia Telescope Compact Array (ATCA) to characterise the emission and intervening radio sources from $\sim100$-$1000$ MHz; power law models are fit to the spectral energy distributions with spectral indices $\alpha_{88}^{1110} = -1.80 \pm 0.09$ and $\alpha_{88}^{944} = -1.13 \pm 0.09$ for the radio halos in Abell 3404 and Abell 141, respectively. Each cluster is noted to have an atypical morphology for a radio-halo--hosting cluster, with Abell 141 having been previously reported to be in a pre-merging state, and Abell 3404 is largely relaxed with only minor evidence for a disturbed morphology. We find that the radio halo power is consistent with the current radio halo sample and $P_\nu$-$M$ scaling relations, but note that the radio halo in Abell 3404 is an ultra-steep-spectrum radio halo (USSRH) and, as with other USSRHs lies slightly below the best-fit $P_{1.4}$-$M$ relation. We find that an updated scaling relation is consistent with previous results and shifting the frequency to 150 MHz does not significantly alter the best-fit relations with a sample of 86 radio halos. We suggest that the USSRH halo in Abell 3404 represents the faint class of radio halos that will be found in clusters undergoing weak mergers.

Naoya Shimono, Tomonori Totani, Takahiro Sudoh

9 pages, 7 figures, submitted to MNRAS

Prospects of the Cherenkov Telescope Array (CTA) for the study of very high energy gamma-ray emission from nearby star-forming galaxies are investigated. In the previous work, we constructed a model to calculate luminosity and energy spectrum of pion-decay gamma-ray emission produced by cosmic-ray interaction with the interstellar medium (ISM), from four physical quantities of galaxies [star formation rate (SFR), gas mass, stellar mass, and effective radius]. The model is in good agreement with the observed GeV--TeV emission of several nearby galaxies. Applying this model to nearby galaxies that are not yet detected in TeV (mainly from the KINGFISH catalog), their hadronic gamma-ray luminosities and spectra are predicted. We identify galaxies of the highest chance of detection by CTA, including NGC 5236, M33, NGC 6946, and IC 342. Concerning gamma-ray spectra, NGC 1482 is particularly interesting because our model predicts that this galaxy is close to the calorimetric limit and its gamma-ray spectral index in GeV--TeV is close to that of cosmic-ray protons injected into ISM. Therefore this galaxy may be detectable by CTA even though its GeV flux is below the {\it Fermi} Large Area Telescope sensitivity limit. In the TeV regime, most galaxies are not in the calorimetric limit, and the predicted TeV flux is lower than that assuming a simple relation between the TeV luminosity and SFR of M82 and NGC 253, typically by a factor of 15. This means that a more sophisticated model beyond the calorimetric limit assumption is necessary to study TeV emission from star-forming galaxies.

Marco Molinaro, Mark Allen, Françoise Genova, André Schaaff, Margarida Castro Neves, Markus Demleitner, Sara Bertocco, Dave Morris, François Bonnarel, Stelios Voutsinas, Catherine Boisson, Giuliano Taffoni

14 pages, 9 figures Proceedings of SPIE 2020 Astronomical Telescopes + Instrumentation Conference 11449 Observatory Operations: Strategies, Processes, and Systems VIII

The International Virtual Observatory Alliance (IVOA) has developed and built, in the last two decades, an ecosystem of distributed resources, interoperable and based upon open shared technological standards. In doing so the IVOA has anticipated, putting into practice for the astrophysical domain, the ideas of FAIR-ness of data and service resources and the Open-ness of sharing scientific results, leveraging on the underlying open standards required to fill the above. In Europe, efforts in supporting and developing the ecosystem proposed by the IVOA specifications has been provided by a continuous set of EU funded projects up to current H2020 ESCAPE ESFRI cluster. In the meantime, in the last years, Europe has realised the importance of promoting the Open Science approach for the research communities and started the European Open Science Cloud (EOSC) project to create a distributed environment for research data, services and communities. In this framework the European VO community, had to face the move from the interoperability scenario in the astrophysics domain into a larger audience perspective that includes a cross-domain FAIR approach. Within the ESCAPE project the CEVO Work Package (Connecting ESFRI to EOSC through the VO) has one task to deal with this integration challenge: a challenge where an existing, mature, distributed e-infrastructure has to be matched to a forming, more general architecture. CEVO started its works in the first months of 2019 and has already worked on the integration of the VO Registry into the EOSC e-infrastructure. This contribution reports on the first year and a half of integration activities, that involve applications, services and resources being aware of the VO scenario and compatible with the EOSC architecture.

Daichi Tsuna, Kazumi Kashiyama, Toshikazu Shigeyama

22 pages, 12 figures. Fig 12 is a take-home figure for busy readers

Transients powered by interaction with the circumstellar medium (CSM) are often observed in wavelengths other than optical, and multi-wavelength modelling can be important when inferring the properties of the explosion and CSM, or for distinguishing from other powering mechanisms. We develop a model calculating time dependent emission spectrum of interaction-powered transients. We solve energy equations of electron-proton plasma in the shocked SN ejecta and CSM and a radiation transfer equation out to the outer edge of the CSM, incorporating the collisional relaxation and the comptonization of the bremsstrahlung radiation. We compare our model to observations of Type IIn supernovae covering frequency ranges from optical to X-rays. For SN 2010jl the observed optical and X-ray light curves can be consistently explained if clumpy or asymmetric structure in the CSM is assumed, in agreement with previous studies. For SN 2014C our model successfully reproduces the X-ray bremsstrahlung component and the emergence of H$\alpha$ emission at 400 days after explosion. Finally we find a parameter space where the supernova is extremely X-ray bright, reaching $10^{43}$-$10^{44}\ {\rm erg\ s^{-1}}$ for up to $100$ days. Such X-ray transients are likely detectable with all-sky surveys by e.g. eROSITA.

The thermal inertia of an asteroid is an indicator of the thermophysical properties of the regolith and is determined by the size of grains on the surface. Previous thermophysical modeling studies of asteroids have identified or suggested that object size, rotation period, and heliocentric distance (a proxy for temperature) as important factors that separately influence thermal inertia. In this work we present new thermal inertias for 239 asteroids and model all three factors in a multi-variate model of thermal inertia. Using multi-epoch infrared data of a large (239) set of objects observed by WISE, we derive the size, albedo, thermal inertia, surface roughness, and sense of spin using a thermophysical modelling approach that doesn't require a priori knowledge of an object's shape or spin axis direction. Our thermal inertia results are consistent with previous values from the literature for similarly sized asteroids, and we identify an excess of retrograde rotators among main-belt asteroids < 8 km. We then combine our results with thermal inertias from the literature to construct a multi-variate model and quantify the dependency on asteroid diameter, rotation period, and surface temperature. This multi-variate model, which accounts for co-dependencies between the three independent variables, identified asteroid diameter and surface temperature as strong controls on thermal inertia.

Lupin C.-C. Lin, H. H. Wang, C. Y. Hui, Jumpei Takata, Paul K. H. Yeung, Chin-Ping Hu, Albert K. H. Kong

11 pages, 5 figures and 1 table; accepted by MNRAS

PSR J1420-6048 is a young gamma-ray pulsar with recurrent glitches. Utilizing long-term monitoring data obtained from the Fermi Gamma-ray Space Telescope, we found that PSR J1420-6048 has shown gamma-ray flux variation and we also detected four glitches between 2008 and 2019. Two of the glitches are previously unknown, and their gamma-ray spectrum also shows variability between each glitch. Since the results might be contaminated by background sources, we discuss whether the observed changes in flux and spectra were caused by artificial misallocations of photons from a nearby pulsar wind nebula (HESS J1420-607) and a pulsar (PSR J1418-6058), or a change of the emission geometry from the target pulsar itself. We examine the correlation of the flux changes and the alternating pulse structure to investigate whether the emission geometry in the outer magnetosphere was changing. By assuming the observational features were not totally resulted from the background environment, we compare our results with similar phenomena observed in other gamma-ray pulsars and propose that a strong crust crack can cause timing anomaly of a neutron star, which can affect the particle accelerations or pair creation regions resulting in the changes of emission behaviors.

Kaustubh Rajwade, Benjamin Stappers, Christopher Williams, Ewan Barr, Mechiel Christiaan Bezuidenhout, Manisha Caleb, Laura Driessen, Fabian Jankowski, Mateusz Malenta, Vincent Morello, Sotirios Sanidas, Mayuresh Surnis

8 pages, 5 figures, published in the proceedings of SPIE Conference on Astronomy and Telescope Instrumentation, No. 114470J

Real-time detections of transients and rapid multi-wavelength follow-up are at the core of modern multi-messenger astrophysics. MeerTRAP is one such instrument that has been deployed on the MeerKAT radio telescope in South Africa to search for fast radio transients in real-time. This, coupled with the ability to rapidly localize the transient in combination with optical co-pointing by the MeerLICHT telescope gives the instrument the edge in finding and identifying the nature of the transient on short timescales. The commensal nature of the project means that MeerTRAP will keep looking for transients even if the telescope is not being used specifically for that purpose. Here, we present a brief overview of the MeerTRAP project. We describe the overall design, specifications and the software stack required to implement such an undertaking. We conclude with some science highlights that have been enabled by this venture over the last 10 months of operation.

Aleksandra Wołowska, Magdalena Kunert-Bajraszewska, Kunal P. Mooley, Aneta Siemiginowska, Preeti Kharb, C. H. Ishwara-Chandra, Gregg Hallinan, Mariusz Gromadzki, Dorota Kozieł-Wierzbowska

A recent multi-year Caltech-NRAO Stripe 82 Survey (CNSS) revealed a group of objects that appeared as new radio sources after $>$5--20 years of absence. They are transient phenomena with respect to the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) survey and constitute the first unbiased sample of renewed radio activity. Here we present the follow-up, radio, optical and X-ray study of them. The group consist of 12 sources, both quasars and galaxies with wide redshift ($\rm 0.04 < z < 1.7$) and luminosity ($\rm 22<log_{10}[L_{1.4GHz}/W~Hz^{-1}]>24.5$) distribution. Their radio properties in the first phase of activity, namely the convex spectra and compact morphology, allow them all to be classified as gigahertz-peaked spectrum (GPS) sources. We conclude that the spectral changes are a consequence of the evolution of newly-born radio jets. Our observations show that over the next few years of activity the GPS galaxies keep the convex shape of the spectrum, while GPS quasars rapidly transform into flat-spectrum sources, which may result in them not being recognized as young sources. The wide range of bolometric luminosities, black hole masses and jet powers among the transient sources indicates even greater population diversity in the group of young radio objects. We also suggest that small changes of the accretion disc luminosity (accretion rate) may be sufficient to ignite low-power radio activity that evolves on the scale of decades.

The Sun frequently accelerates near-relativistic electron beams that travel out through the solar corona and interplanetary space. Interacting with their plasma environment, these beams produce type III radio bursts, the brightest astrophysical radio sources seen from the Earth. The formation and motion of type III fine frequency structures is a puzzle but is commonly believed to be related to plasma turbulence in the solar corona and solar wind. Combining a theoretical framework with kinetic simulations and high-resolution radio type III observations using the Low Frequency Array, we quantitatively show that the fine structures are caused by the moving intense clumps of Langmuir waves in a turbulent medium. Our results show how type III fine structure can be used to remotely analyse the intensity and spectrum of compressive density fluctuations, and can infer ambient temperatures in astrophysical plasma, both significantly expanding the current diagnostic potential of solar radio emission.

Xiaodong Ren, Pablo Astudillo, Urs U. Graf, Richard E. Hills, Sebastian Jorquera, Bojan Nikolic, Stephen C. Parshley, Nicolás Reyes, Lars Weikert

Proceedings Volume 11445, Ground-based and Airborne Telescopes VIII SPIE Astronomical Telescopes + Instrumentation, 2020, Online Only Conference

We describe a system being developed for measuring the shapes of the mirrors of the Fred Young Submillimeter Telescope (FYST), now under construction for the CCAT Observatory. "Holographic" antenna-measuring techniques are an efficient and accurate way of measuring the surfaces of large millimeter-wave telescopes and they have the advantage of measuring the wave-front errors of the whole system under operational conditions, e.g. at night on an exposed site. Applying this to FYST, however, presents significant challenges because of the high accuracy needed, the fact that the telescope consists of two large off-axis mirrors, and a requirement that measurements can be made without personnel present. We use a high-frequency (~300GHz) source which is relatively close to the telescope aperture (<1/100th of the Fresnel distance) to minimize atmospheric effects. The main receiver is in the receiver cabin and can be moved under remote control to different positions, so that the wave-front errors in different parts of the focal plane can be measured. A second receiver placed on the yoke provides a phase reference. The signals are combined in a digital cross-correlation spectrometer. Scanning the telescope provides a map of the complex beam pattern. The surface errors are found by inference, i.e. we make models of the reflectors with errors and calculate the patterns expected, and then iterate to find the best match to the data. To do this we have developed a fast and accurate method for calculating the patterns using the Kirchhoff-Fresnel formulation. This paper presents details of the design and outlines the results from simulations of the measurement and inference process. These indicate that a measurement accuracy of ~3 microns rms is achievable.

T. M. Sitnova, L. I. Mashonkina, A. M. Tatarnikov, O. V. Voziakova, M. A. Burlak, Yu. V. Pakhomov, P. Jablonka, M. D. Neretina, A. Frebel

accepted for publication in MNRAS

We present atmospheric parameters and abundances for chemical elements from carbon to barium in metal-poor stars in Segue 1 (seven stars), Coma Berenices (three stars), and Triangulum II (one star) ultra-faint dwarf galaxies (UFDs). The effective temperatures rely on new photometric observations in the visible and infra-red bands, obtained with the 2.5 m telescope of the SAI MSU Caucasian observatory. Abundances of up to fourteen chemical elements were derived under the non-local thermodynamic equilibrium (NLTE) line formation, and LTE abundances were obtained for up to five more elements. For the first time we present abundance of oxygen in Seg 1 S1 and S4, silicon in ComaBer S2 and Tri II S40, potassium in Seg 1 S1-S6 and ComaBer S1-S3, and barium in Seg 1 S7. Three stars in Segue 1, two stars in Coma Berenices, and Triangulum II star have very low [Na/Mg] of -1.08 to -1.67 dex, which is usually attributed in the literature to an odd-even effect produced by nucleosynthesis in massive metal-free stars. We interpret this chemical property as a footprint of first stars, which is not blurred due to a small number of nucleosynthesis events that contributed to chemical abundance patterns of the sample stars. Our NLTE abundances of Sr and Ba in Coma Berenices, Segue 1, and Triangulum II report on lower [Sr/Ba] abundance ratio in the UFDs compared to that in classical dwarf spheroidal galaxies and the Milky Way halo. However, in UFDs, just as in massive galaxies, [Sr/Ba] is not constant and it can be higher than the pure r-process ratio. We suggest a hypothesis of Sr production in metal-poor binaries at the earliest epoch of galactic evolution.

S. Hunziker, H. M. Schmid, J. Ma, F. Menard, H. Avenhaus, A. Boccaletti, J. L. Beuzit, G. Chauvin, K. Dohlen, C. Dominik, N. Engler, C. Ginski, R. Gratton, T. Henning, M. Langlois, J. Milli, D. Mouillet, C. Tschudi, R. G. van Holstein, A. Vigan

20 pages, 14 figures

We present high-precision photometry and polarimetry for the protoplanetary disk around HD142527, with a focus on determining the light scattering parameters of the dust. We re-reduced polarimetric differential imaging data of HD142527 in the VBB (735 nm) and H-band (1625 nm) from the ZIMPOL and IRDIS subinstruments of SPHERE/VLT. With polarimetry and photometry based on reference star differential imaging, we were able to measure the linearly polarized intensity and the total intensity of the light scattered by the circumstellar disk with high precision. We used simple Monte Carlo simulations of multiple light scattering by the disk surface to derive constraints for three scattering parameters of the dust: the maximum polarization of $P_{\rm max}$, the asymmetry parameter $g$, and the single-scattering albedo $\omega$. We measure a reflected total intensity of $51.4\pm1.5$ mJy and $206\pm12$ mJy and a polarized intensity of $11.3\pm0.3$ mJy and $55.1\pm3.3$ mJy in the VBB and H-band, respectively. We also find in the visual range a degree of polarization that varies between $28\%$ on the far side of the disk and $17\%$ on the near side. The disk shows a red color for the scattered light intensity and the polarized intensity, which are about twice as high in the near-infrared when compared to the visual. We determine with model calculations the scattering properties of the dust particles and find evidence for strong forward scattering ($g\approx 0.5-0.75$), relatively low single-scattering albedo ($\omega \approx 0.2-0.5$), and high maximum polarization ($P_{\rm max} \approx 0.5-0.75$) at the surface on the far side of the disk for both observed wavelengths. The optical parameters indicate the presence of large aggregate dust particles, which are necessary to explain the high maximum polarization, the strong forward-scattering nature of the dust, and the observed red disk color.

Simon Guichandut, Andrew Cumming, Maurizio Falanga, Zhaosheng Li, Michael Zamfir

23 pages, 13 figures. Submitted to ApJ (March 2021)

We calculate steady-state models of radiation-driven super-Eddington winds and static expanded envelopes of neutron stars caused by high luminosities in type I X-ray bursts. We use flux-limited diffusion to model the transition from optically thick to optically thin, and include effects of general relativity, allowing us to study the photospheric radius close to the star as the hydrostatic atmosphere evolves into a wind. We find that the photospheric radius evolves monotonically from static envelopes ($r_{\rm ph}\lesssim 50-70$ km) to winds ($r_{\rm ph}\approx 100-1000$ km). Photospheric radii of less than $100$ km, as observed in most photospheric radius expansion bursts, can be explained by static envelopes, but only in a narrow range of luminosity. In most bursts, we would expect the luminosity to increase further, leading to a wind with photospheric radius $\gtrsim 100$ km. In the contraction phase, the expanded envelope solutions show that the photosphere is still $\approx 1$ km above the surface when the effective temperature is only $3\%$ away from its maximum value. This is a possible systematic uncertainty when interpreting the measured Eddington fluxes from bursts at touchdown. We also discuss the applicability of steady-state models to describe the dynamics of bursts. In particular, we show that the sub to super-Eddington transition during the burst rise is rapid enough that static models are not appropriate. Finally, we analyze the strength of spectral shifts in our models. Expected shifts at the photosphere are dominated by gravitational redshift, and are therefore predicted to be less than a few percent.

Fast radio bursts (FRBs) probe the total column density of free electrons in the intergalactic medium (IGM) along the path of propagation though the dispersion measures (DMs) which depend on the baryon mass fraction in the IGM, i.e., $f_{\rm IGM}$. In this letter, we investigate the large-scale clustering information of DMs to study the evolution of $f_{\rm IGM}$. When combining with the Planck 2018 measurements, we could give tight constraints on the evolution of $f_{\rm IGM}(z)$ from about $10^4$ FRBs with the intrinsic DM scatter of $30(1+z)~ \rm pc/cm^3$ spanning 80% of the sky and redshift range $z=0-3$. Firstly, we consider the Taylor expansion of $f_{\rm IGM}(z)$ up to second order, and find that the mean relative standard deviation $\sigma(f_{\rm IGM})\equiv\left\langle \sigma[f_{\rm IGM}(z)] /f_{\rm IGM}(z) \right\rangle$ is about 6.7%. In order to alleviate the dependence on fiducial model, we also adopt a non-parametric methods in this work, the local principle component analysis. We obtain the consistent, but weaker constraints on the evolution of $f_{\rm IGM}(z)$, namely the mean relative standard deviation $\sigma(f_{\rm IGM})$ is 21.4%. With the forthcoming surveys, this could be a complimentary method to investigate the baryon mass fraction in the IGM.

Ravit Helled, Stephanie Werner, Caroline Dorn, Tristan Guillot, Masahiro Ikoma, Yuichi Ito, Mihkel Kama, Tim Lichtenberg, Yamila Miguel, Oliver Shorttle, Paul J. Tackley, Diana Valencia, Allona Vazan

Accepted for publication in Experimental Astronomy

The recently adopted Ariel ESA mission will measure the atmospheric composition of a large number of exoplanets. This information will then be used to better constrain planetary bulk compositions. While the connection between the composition of a planetary atmosphere and the bulk interior is still being investigated, the combination of the atmospheric composition with the measured mass and radius of exoplanets will push the field of exoplanet characterisation to the next level, and provide new insights of the nature of planets in our galaxy. In this white paper, we outline the ongoing activities of the interior working group of the {\it Ariel} mission, and list the desirable theoretical developments as well as the challenges in linking planetary atmospheres, bulk composition and interior structure.

Sapna Mishra, M. Vivek, Hum Chand, Ravi Joshi

16 pages, 9 figures, Accepted for publication in MNRAS

We present a new set of 84 Broad absorption line (BAL) quasars ( 1.7 $<$ \zem $<$ 4.4) exhibiting an appearance of \civ BAL troughs over 0.3$-$4.8 rest-frame years by comparing the Sloan Digital Sky Survey Data Release (SDSSDR)-7, SDSSDR-12, and SDSSDR-14 quasar catalogs. We contrast the nature of BAL variability in this appearing BAL quasar sample with a disappearing BAL quasar sample studied in literature by comparing the quasar's intrinsic, BAL trough, and continuum parameters between the two samples. We find that appearing BAL quasars have relatively higher redshift and smaller probed timescales as compared to the disappearing BAL quasars. To mitigate the effect of any redshift bias, we created control samples of appearing and disappearing BAL quasars that have similar redshift distribution. We find that the appearing BAL quasars are relatively brighter and have shallower and wider BAL troughs compared to the disappearing BAL sample. The distribution of quasar continuum variability parameters between the two samples is clearly separated, with the appearance of the BAL troughs being accompanied by the dimming of the continuum and vice versa. Spectral index variations in the two samples also point to the anti-correlation between the BAL trough and continuum variations consistent with the "bluer when brighter" trend in quasars. We show that the intrinsic dust model is less likely to be a favorable scenario in explaining BAL appearance/disappearance. Our analysis suggests that the extreme variations of BAL troughs like BAL appearance/disappearance are mainly driven by changes in the ionization conditions of the absorbing gas.

David S. N. Rupke (1 and 2), Adam D. Thomas (2), Michael A. Dopita (2) ((1) Rhodes College, (2) Australian National University)

MNRAS, in press

Dusty, neutral outflows and inflows are a common feature of nearby star-forming galaxies. We characterize these flows in eight galaxies -- mostly AGN -- selected for their widespread NaI D signatures from the Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). This survey employs deep, wide field-of-view integral field spectroscopy at moderate spectral resolution (R=7000 at NaI D). We significantly expand the sample of sightlines in external galaxies in which the spatially-resolved relationship has been studied between cool, neutral gas properties -- N(NaI), Weq(NaI D) -- and dust -- E(B-V) from both stars and gas. Our sample shows strong, significant correlations of total Weq with E(B-V)_stars and g-i colour within individual galaxies; correlations with E(B-V)_gas are present but weaker. Regressions yield slope variations from galaxy to galaxy and intrinsic scatter ~1 Angstrom. The sample occupies regions in the space of N(NaI) and Weq^abs vs. E(B-V)_gas that are consistent with extrapolations from other studies to higher colour excess [E(B-V)_gas ~ 1]. For perhaps the first time in external galaxies, we detect inverse P Cygni profiles in the NaI D line, presumably due to inflowing gas. Via Doppler shifted NaI D absorption and emission lines, we find ubiquitous flows that differ from stellar rotation by >100 km/s or have |v,abs - v,em| > 100 km/s. Inflows and outflows extend toward the edge of the detected stellar disk/FOV, together subtend 10-40% of the projected disk, and have similar mean N(NaI) and Weq(NaI D). Outflows are consistent with minor-axis or jet-driven flows, while inflows tend toward the projected major axis. The inflows may result from non-axisymmetric potentials, tidal motions, or halo infall.

The possibility that primordial black holes (PBHs) form a part of dark matter has been considered for a long time but poorly constrained in the $1-100~M_{\odot}$ (or stellar mass range). However, a renewed special interest of PBHs in this mass window was triggered by the discovery at LIGO of the merger events of black-hole binaries. Fast radio bursts (FRBs) are bright radio transients with millisecond duration and high all-sky occurrence rate. Lensing effect of these bursts has been proposed as one of the cleanest probes for constraining the presence of PBHs in the stellar mass window. In this paper, we first investigate constraints on the abundance of PBHs from the latest FRB observations for both the monochromatic mass distribution and three other popular extended mass distributions (EMDs). We find that constraints from currently public FRB observations are relatively weaker than those from existing gravitational wave detections. Furthermore, we forecast constraining power of future FRB observations on the abundance of PBHs with different mass distributions of PBHs and different redshift distributions of FRBs taken into account. Finally, We find that constraints of parameter space on EMDs from $\sim10^5$ FRBs with $\overline{\Delta t}\leq1 ~\rm ms$ would be comparable with what can be constrained from gravitational wave events. It is foreseen that upcoming complementary multi-messenger observations will yield considerable constraints on the possibilities of PBHs in this intriguing mass window.

Hideaki Takemura, Fumitaka Nakamura, Shun Ishii, Yoshito Shimajiri, Patricio Sanhueza, Takashi Tsukagoshi, Ryohei Kawabe, Tomoya Hirota, Akimasa Kataoka

66 pages, 15 figures, accepted by PASJ

We have performed an unbiased dense core survey toward the Orion A Giant Molecular Cloud in the C$^{18}$O ($J$=1--0) emission line taken with the Nobeyama Radio Observatory (NRO) 45-m telescope. The effective angular resolution of the map is 26", which corresponds to $\sim$ 0.05 pc at a distance of 414 pc. By using the Herschel-Planck H$_2$ column density map, we calculate the C$^{18}$O fractional abundance and find that it is roughly constant over the column density range of $\lesssim$ 5 $\times$ 10$^{22}$ cm$^{-3}$, although a trend of C$^{18}$O depletion is determined toward higher column density. Therefore, C$^{18}$O intensity can follow the cloud structure reasonably well. The mean C$^{18}$O abundance in Orion A is estimated to be 5.7$\times$10$^{-7}$, which is about 3 times larger than the fiducial value. We identified 746 C$^{18}$O cores with astrodendro and classified 709 cores as starless cores. We compute the core masses by decomposing the Herschel-Planck dust column density using the relative proportions of the C$^{18}$O integrated intensities of line-of-sight components. Applying this procedure, we attempt to remove the contribution of the background emission, i.e., the ambient gas outside the cores. Then, we derived mass function for starless cores and found that it resembles the stellar initial mass function (IMF). The CMF for starless cores, $dN/dM$, is fitted with a power-law relation of $M^\alpha$ with a power index of $\alpha = -$2.25$\pm$ 0.16 at the high-mass slope ($\gtrsim$ 0.44 $M_\odot$). We also found that the ratio of each core mass to the total mass integrated along the line of sight is significantly large. Therefore, in the previous studies, the core masses derived from the dust image are likely to be overestimated at least by a factor of a few. Accordingly, such previous studies may underestimate the star formation efficiency of individual cores.

Hideaki Takemura, Fumitaka Nakamura, Shuo Kong, Héctor G. Arce, John M. Carpenter, Volker Ossenkopf-Okada, Ralf Klessen, Patricio Sanhueza, Yoshito Shimajiri, Takashi Tsukagoshi, Ryohei Kawabe, Shun Ishii, Kazuhito Dobashi, Tomomi Shimoikura, Paul F. Goldsmith, Álvaro Sánchez-Monge, Jens Kauffmann, Thushara Pillai, Paolo Padoan, Adam Ginsberg, Rowan J. Smith, John Bally, Steve Mairs, Jaime E. Pineda, Dariusz C. Lis, Blakesley Burkhart, Peter Schilke, Hope How-Huan Chen, Andrea Isella, Rachel K. Friesen, Alyssa A. Goodman, Doyal A. Harper

13 pages, 5 figures, accepted by ApJL

Applying dendrogram analysis to the CARMA-NRO C$^{18}$O ($J$=1--0) data having an angular resolution of $\sim$ 8", we identified 692 dense cores in the Orion Nebula Cluster (ONC) region. Using this core sample, we compare the core and initial stellar mass functions in the same area to quantify the step from cores to stars. About 22 \% of the identified cores are gravitationally bound. The derived core mass function (CMF) for starless cores has a slope similar to Salpeter's stellar initial mass function (IMF) for the mass range above 1 $M_\odot$, consistent with previous studies. Our CMF has a peak at a subsolar mass of $\sim$ 0.1 $M_\odot$, which is comparable to the peak mass of the IMF derived in the same area. We also find that the current star formation rate is consistent with the picture in which stars are born only from self-gravitating starless cores. However, the cores must gain additional gas from the surroundings to reproduce the current IMF (e.g., its slope and peak mass), because the core mass cannot be accreted onto the star with a 100\% efficiency. Thus, the mass accretion from the surroundings may play a crucial role in determining the final stellar masses of stars.

Daniel Schwartz, Cristiana Spingola, Anna Barnacka

9 pages, 5 figures, submitted to ApJ

We use a Chandra X-ray observation of the gravitationally lensed system MGB2016+112 at z=3.273 to elucidate presence of at least two X-ray sources. We find that these sources are consistent with the VLBI components measured by \citet{Spingola19}, which are separated by $\sim 200$ pc. Their intrinsic 0.5 -- 7 keV source frame luminosities are 2.6$\times$10$^{43}$ and 4.2$\times$10$^{44}$ erg s$^{-1}$. Most likely this system contains a dual active galactic nucleus (AGN), but we possibly are detecting an AGN plus a pc-scale X-ray jet, the latter lying in a region at very high magnification. The quadruply lensed X-ray source is within $\pm$40 pc (1$\sigma$) of its VLBI counterpart. Using a gravitational lens as a telescope, and a novel statistical application, we have achieved unprecedented accuracy for measuring metric distances at such large redshifts in X-ray astronomy, which is tens of mas if the source is located close to the caustics, while it is of hundreds of mas if the source is in a region at lower amplification. The present demonstration of this approach has implications for future X-ray investigations of large numbers of lensed systems.

J.M. Boxelaar, R.J. van Weeren, A. Botteon

Submitted to Astronomy & Computing

Here we present Halo-FDCA, a robust open source Python package for modeling and estimating total flux densities of radio (mini) halos in galaxy clusters. Radio halos are extended ( ~200 - 1500 kpc in size) synchrotron emitting sources found in galaxy clusters that trace the presence of cosmic rays and magnetic fields in the intracluster medium (ICM). These sources are centrally located and have a low surface brightness. Their exact origin is still unknown but they are likely related to cosmic rays being re-accelerated in-situ by merger or sloshing driven ICM turbulence. The presented algorithm combines the numerical power of the Markov Chain Monte Carlo routine and multiple theoretical models to estimate the total radio flux density of a radio halo from a radio image and its associated uncertainty. This method introduces a flexible analytic fitting procedure to replace existing simplistic manual measurements prone to biases and inaccuracies. It allows to easily determine the properties of the emission and is particularly suitable for future studies of large samples of clusters.

Fraser Kennedy, Philip Bull

20 pages, 14 figures. For the busy reader: see Figs. 4, 5, and 8

Despite being designed as an interferometer, the MeerKAT radio array (an SKA pathfinder) can also be used in autocorrelation (`single-dish') mode, where each dish scans the sky independently. Operating in this mode allows extremely high survey speeds to be achieved, albeit at significantly lower angular resolution. We investigate the recovery of the baryon acoustic oscillation (BAO) scale from multipoles of the redshift-space correlation function as measured by a low angular resolution 21cm intensity mapping survey of this kind. Our approach is to construct an analytic model of the multipoles of the correlation function and their covariance matrix that includes foreground contamination and beam resolution effects, which we then use to generate an ensemble of mock data vectors from which we attempt to recover the BAO scale. In line with previous studies, we find that recovery of the transverse BAO scale $\alpha_{\perp}$ is hampered by the strong smoothing effect of the instrumental beam with increasing redshift, while the radial scale $\alpha_\parallel$ is much more robust. The multipole formalism naturally incorporates transverse information when it is available however, and so there is no need to perform a radial-only analysis. In particular, the quadrupole of the correlation function preserves a distinctive BAO `bump' feature even for large smoothing scales. We also investigate the robustness of BAO scale recovery to beam model accuracy, severity of the foreground removal cuts, and accuracy of the covariance matrix model, finding in all cases that the radial BAO scale can be recovered in an accurate, unbiased manner.

J. F. Radcliffe (1, 2, 3), P. D. Barthel (2), A. P. Thomson (3), M. A. Garrett (3, 4), R. J. Beswick (3), T. W. B. Muxlow (3) ((1) University of Pretoria, SA, (2) University of Groningen, NL, (3) University of Manchester, UK, (4) University of Leiden, NL)

26 pages, 10 figures. Accepted for publication to A&A

Obtaining a census of active galactic nuclei (AGN) activity across cosmic time is critical to our understanding of galaxy evolution and formation. Many AGN classification techniques are compromised by dust obscuration. However, very long baseline interferometry (VLBI) can be used to identify compact emission that can only be attributed to AGN activity. This is the second in a series of papers dealing with the compact radio population in the GOODS-N field. We review 14 different AGN classification techniques in the context of a VLBI-detected sample, and use these to investigate the nature of the AGN as well as their host galaxies. We find that no single identification technique can identify all VLBI objects as AGN. Infrared colour-colour selection is most notably incomplete. However, the usage of multiple classification schemes can identify all VLBI-selected AGN, independently verifying similar approaches used in other deep field surveys. In the era of large area surveys with instruments such as the SKA and ngVLA, multi-wavelength coverage, which relies heavily upon observations from space, is often unavailable. Therefore, VLBI remains an integral component in detecting AGN of the jetted efficient and inefficient accretion types. A substantial fraction (46%) of the VLBI AGN have no X-ray counterpart, which is most likely due to lack of sensitivity in the X-ray band. A high fraction of the VLBI AGN reside in low or intermediate redshift dust-poor early-type galaxies. These most likely exhibit inefficient accretion. Finally, A significant fraction of the VLBI AGN reside in symbiotic dusty starburst - AGN systems. We present an extensive compilation of the multi-wavelength properties of all the VLBI-selected AGN in GOODS-N in the Appendix.

Adam S. Trotter, Daniel E. Reichart, Aaron P. LaCluyzé, Rachel Freed

15 pages, 2 figures, RTSRE Conference Proceedings, July 2018

Most students do not enroll in introductory astronomy as part of their major; for many, it is the last science course they will ever take. Thus, it has great potential to shape students' attitudes toward STEM fields for the rest of their life. We therefore argue that it is less important, when assessing the effectiveness of introductory astronomy courses, to explore traditional curricular learning gains than to explore the effects that various course components have on this attitude. We describe the results of our analysis of end-of-semester surveys returned by a total of 749 students in 2014-2015, at 10 institutions that employed at least part of the introductory astronomy lecture and lab curriculum we first implemented at UNC-Chapel Hill in 2009. Surveys were designed to measure each student's attitude, and to probe the correlation of attitude with their utilization of, and satisfaction with, various course components, along with other measures of their academic background and their self-assessed performance in the course. We find that students' attitudes are significantly positively correlated with the grade they expect to receive, and their rating of the course's overall effectiveness. To a lesser degree, we find that students' attitudes are positively correlated with their mathematical background, whether they intend to major or pursue a career in STEM, and their rating of the effectiveness of the instructor. We find that students' attitudes are negatively correlated with the amount of work they perceived the course to involve, and, surprisingly, the size and reputation of their home institution. We also find that, for the subsets of students who were exposed to them, students' attitudes are positively correlated with their perception of the helpfulness of the lecture component of the course, and of telescope-based labs that utilized UNC-CH's Skynet Robotic Telescope Network.

To understand the true nature of black holes, fundamental theoretical developments should be linked all the way to observational features of black holes in their natural astrophysical environments. Here, we take several steps to establish such a link. We construct a family of spinning, regular black-hole spacetimes based on a locality principle for new physics and analyze their shadow images. We identify characteristic image features associated to regularity (increased compactness and relative stretching) and to the locality principle (cusps and asymmetry) that persist in the presence of a simple analytical disk model. We conjecture that these occur as universal features of distinct classes of regular black holes based on different sets of construction principles for the corresponding spacetimes.

The impact of radiation dramatically increases at high altitudes in the Earth's atmosphere and in space. Therefore, monitoring and access to radiation environment measurements are critical for estimating the radiation exposure risks of aircraft and spacecraft crews and the impact of space weather disturbances on electronics. Addressing these needs requires convenient access to multi-source radiation environment data and enhancement of visualization and search capabilities. The Radiation Data Portal represents an interactive web-based application for search and visualization of in-flight radiation measurements. The Portal enhances the exploration capabilities of various properties of the radiation environment and provides measurements of dose rates along with information on space weather-related conditions. The Radiation Data Portal back-end is a MySQL relational database that contains the radiation measurements obtained from the Automated Radiation Measurements for Aerospace Safety (ARMAS) device and the soft X-ray and proton flux measurements from the Geostationary Operational Environmental Satellite (GOES). The implemented Application Programming Interface (API) and Python routines allow a user to retrieve the database records without interaction with the web interface. As a use case of the Radiation Data Portal, we examine ARMAS measurements during an enhancement of the Solar Proton (SP) fluxes, known as Solar Proton Events (SPEs), and compare them to measurements during SP-quiet periods.

Graciela B. Gelmini, Anna Simpson, Edoardo Vitagliano

7 pages, 3 figures

Axion-like particles (ALPs) are a compelling candidate for dark matter (DM). Their production is associated with the formation of a string-wall network. This system must annihilate, producing gravitational waves and non-relativistic ALPs. We show that these gravitational waves, if produced at temperatures below 100 eV, could be detected by future cosmological probes for ALPs with mass from $10^{-16}$ to $10^{6}$ eV that could constitute the entirety of the DM.

In hot dense plasmas of intermediate or high-Z elements in the state of local thermodynamic equilibrium, the number of electronic configurations contributing to key macroscopic quantities such as the spectral opacity and equation of state, can be enormous. In this work we present systematic methods for the analysis of the number of relativistic electronic configurations in a plasma. While the combinatoric number of configurations can be huge even for mid-Z elements, the number of configurations which have non negligible population is much lower and depends strongly and non-trivially on temperature and density. We discuss two useful methods for the estimation of the number of populated configurations: (i) using an exact calculation of the total combinatoric number of configurations within superconfigurations in a converged super-transition-array (STA) calculation, and (ii) by using an estimate for the multidimensional width of the probability distribution for electronic population over bound shells, which is binomial if electron exchange and correlation effects are neglected. These methods are analyzed, and the mechanism which leads to the huge number of populated configurations is discussed in detail. Comprehensive average atom finite temperature density functional theory (DFT) calculations are performed in a wide range of temperature and density for several low, mid and high Z plasmas. The effects of temperature and density on the number of populated configurations are discussed and explained.

In theories beyond the Standard Model, in particular in theories with extra spatial dimensions such as superstring theory, there are a large number of scalar fields which appear in the low energy effective action. These moduli fields must be stabilized. Often, moduli stabilization involves a stage during which the moduli fields oscillatte coherently about their ground state value. Here, we study moduli and graviton production during the period during which the background modulus field is oscillating, assuming that this period takes place in the radiation phase. We find a resonant production of moduli fluctuations, and a tachyonic instability to the generation of long wavelength gravitational waves. As a consequence, the period of moduli stabilization is short on a Hubble time scale.

Simone Bianchi

5 pages. Accepted by JAHH. Comments and references to the adoption of mean time in other cities/countries are welcomed

It is usually stated in the literature that Geneva was the first city to adopt mean solar time, in 1780, followed by London (or the whole of England) in 1792, Berlin in 1810 and Paris in 1816. In this short note I partially revise this statement, using primary references when available, and provide dates for a few other European cities. Though no exact date was found for the first public use of mean time, the primacy seems to belong to England, followed by Geneva in 1778-79 (for horologists), Berlin in 1810, Geneva in 1821 (for public clocks), Vienna in 1823, Paris in 1826, Rome in 1847, Turin in 1849, Milan, Bologna and Florence in 1860.

S.D. Odintsov, V.K. Oikonomou

PDU Accepted, Abstract reduced due to arXiv limitations

In this work we shall study the implications of a subclass of $E$-models cosmological attractors, namely of $a$-attractors, on hydrodynamically stable slowly rotating neutron stars. Specifically, we shall present the Jordan frame theory of the $a$-attractors, and by using a conformal transformation we shall derive the Einstein frame theory. We discuss the inflationary context of $a$-attractors in order to specify the allowed range of values for the free parameters of the model based on the latest cosmic-microwave-background-based Planck 2018 data. Accordingly, using the notation and physical units frequently used in theoretical astrophysics contexts, we shall derive the Tolman-Oppenheimer-Volkoff equations in the Einstein frame. Assuming a piecewise polytropic equation of state, the lowest density part of which shall be chosen to be the WFF1, or APR or the SLy EoS, we numerically solve the Tolman-Oppenheimer-Volkoff equations using a double shooting python-based "LSODA" numerical code. The resulting picture depends on the value of the parameter $a$ characterizing the $a$-attractors. As we show, for large values of $a$, which do not produce a viable inflationary era, the $M-R$ graphs are nearly identical to the general relativistic result, and these two are discriminated at large central densities values. Also, for large $a$-values, the WFF1 equation of state is excluded, due to the GW170817 constraints. In addition, the small $a$ cases produce larger masses and radii compared to the general relativistic case and are compatible with the GW170817 constraints on the radii of neutron stars. Our results indicate deep and not yet completely understood connections between non-minimal inflationary attractors and neutron stars phenomenology in scalar-tensor theory.

Jan Rafelski, Cheng Tao Yang

10 pages APPB, 3 figures, Prepared for the 60th anniversary of the Cracow School of Theoretical Physics

Muon abundance is required for the understanding of several fundamental questions regarding properties of the primordial Universe. In this paper we evaluate the production and decay rates of muons in the cosmic plasma as a function of temperature. This allows us to determine when exactly the muon abundance disappears. When the Universe cools below the temperature $kT_\mathrm{disappear}\approx 4.195$ MeV the muon decay rate overwhelms production rates and muons vanish quasi-instantaneously from the Universe. Interestingly, we show that at $T_\mathrm{disappear}$ the muon pair number is equal to baryon abundance within measurement error margin.

Luis A. Anchordoqui, Ignatios Antoniadis, Dieter Lust, Jorge F. Soriano

22 pages revtex, 6 figures

The Swampland de Sitter conjecture in combination with upper limits on the tensor-to-scalar ratio $r$ derived from observations of the cosmic microwave background endangers the paradigm of slow-roll single field inflation. This conjecture constrains the first and the second derivatives of the inflationary potential in terms of two ${\cal O} (1)$ constants $c$ and $c'$. In view of these restrictions we reexamine single-field inflationary potentials with $S$-duality symmetry, which ameliorate the unlikeliness problem of the initial condition. We compute $r$ at next-to-leading order in slow-roll parameters for the most general form of $S$-dual potentials and confront model predictions to constraints imposed by the de Sitter conjecture. We find that $c \sim {\cal O} (10^{-1})$ and $c' \sim {\cal O} (10^{-2})$ can accommodate the 95\% CL upper limit on $r$. By imposing at least 50 $e$-folds of inflation with the effective field theory description only valid over a field displacement ${\cal O} (1)$ when measured as a distance in the target space geometry, we further restrict $c \sim {\cal O} (10^{-2})$, while the constraint on $c'$ remains unchanged. We comment on how to accommodate the required small values of $c$ and $c'$.

R. A. Heinonen, P. H. Diamond, M. F. D. Katz, G. E. Ronimo

27 pages, 9 figures

Magnetohydrodynamic (MHD) turbulence on a $\beta$-plane with an in-plane mean field, a system which serves as a simple model for the solar tachocline, is investigated analytically and computationally. We show that conservation of squared magnetic potential in this system leads to a global increase in the cross-helicity, which is otherwise conserved in pure MHD (for which the Rossby parameter $\beta$ is zero). We perform a closure using weak turbulence theory and show that the cross-helicity spectrum entirely determines momentum transport. We also note that perturbation theory for small Rossby parameter is impossible in weak turbulence, since it changes the topology of surfaces on which resonant interactions occur. We support our results with numerical simulations, which clearly indicate that the cross-helicity is most important in the transitional regime between Rossby and Alfv\'enic turbulence.

Shota Nakagawa, Fuminobu Takahashi, Masaki Yamada

6 pages, 2 figures

Cosmic birefringence is predicted if an axion-like particle (ALP) moves after the recombination. We show that this naturally happens if the ALP is coupled to the dark matter density because it then acquires a large effective mass after the matter-radiation equality. Our scenario applies to a broad range of the ALP mass $m_\phi \lesssim 10^{-28}$ eV, even smaller than the present Hubble constant. We give a simple model to realize this scenario, where dark matter is made of hidden monopoles, which give the ALP such a large effective mass through the Witten effect. The mechanism works if the ALP decay constant is of order the GUT scale without a fine-tuning of the initial misalignment angle. For smaller decay constant, the hidden monopole can be a fraction of dark matter. We also study the implications for the QCD axion, and show that the domain wall problem can be solved by the effective mass.

Stefan Bender, Patrick J. Espy, Larry J. Paxton

15 pages, 5 figures, in review for Ann. Geophys

The coupling of the atmosphere to the space environment has become recognized as an important driver of atmospheric chemistry and dynamics. In order to quantify the effects of particle precipitation on the atmosphere, reliable global energy inputs on spatial scales commensurate with particle precipitation variations are required. To that end, we have validated the Special Sensor Ultraviolet Spectrographic Imagers (SSUSI) products for average electron energy and electron energy flux by comparing to EISCAT electron density profiles. This comparison shows that SSUSI FUV observations can be used to provide ionization rate profiles throughout the auroral region. The SSUSI on board the Defense Meteorological Satellite Program (DMSP) Block 5D3 satellites provide nearly hourly, high-resolution UV snapshots of auroral emissions. Here we use the SSUSI-derived energies and fluxes to drive standard parametrizations in order to obtain ionization-rate and electron-density profiles in the E-region (90--150 km), which are then compared to EISCAT ground-based electron density measurements. We compare the data from two satellites, DMSP F17 and F18, to the Troms{\o} UHF radar profiles. We find that differentiating between the magnetic local time (MLT) morning (3--11 h) and evening (15--23 h) provides the best fit to the ground-based data. The data agree well in the MLT morning sector using a Maxwellian electron spectrum, while in the evening sector using a Gaussian spectrum and accounting for bounce-electrons achieved optimum agreement with EISCAT. Depending on the satellite and MLT, the median of the differences varies between 0 and 20% above 105 km (F17) and $\pm$15% above 100 km (F18). Because of the large gradient below those altitudes, the relative differences get larger without substantially increasing absolute differences, with virtually no statistically significant differences at the 1-sigma level.

The TianQin space Gravitational Waves (GW) observatory will contain 3 geocentric and circularly orbiting spacecraft with an orbital radius of 10^5 km, to detect the GW in the milli-hertz frequency band. Each spacecraft pair will establish a 1.7*10^5 km-long laser interferometer immersed in the solar wind and the magnetospheric plasmas to measure the phase deviations induced by the GW. GW detection requires a high-precision measurement of the laser phase. The cumulative effects of the long distance and the periodic oscillations of the plasma density may induce an additional phase noise. This paper aims to model the plasma induced phase deviation of the inter-spacecraft laser signals, using a realistic orbit simulator and the Space Weather Modeling Framework (SWMF) model. Preliminary results show that the plasma density oscillation can induce the phase deviations close to 2*10^-6 rad/Hz^1/2 or 0.3pm/Hz^1/2 in the milli-hertz frequency band and it is within the error budget assigned to the displacement noise of the interferometry. The amplitude spectrum density of phases along three arms become more separated when the orbital plane is parallel to the Sun-Earth line or during a magnetic storm. Finally, the dependence of the phase deviations on the orbital radius is examined.

We study the observational predictions of the phenomenological AdS/QCD inspired model, in which the inflaton field emerges in a four dimensional strongly coupled gauge theory, in which the chiral symmetry breaking occurs through the formation of quark condensate. Based on a top-down approach of AdS/QCD, using a D7 brane in the background of $N_c$ D3 branes, it has already been shown that chiral symmetry breaking in a magnetic field through the generation of Higgs vacuum expectation value could be a second order phase transition, although it was doubted that this scenario could lead to enough amount of inflation. Using an iterative method, we consistently solve for the time-dependent parameters, including the embedding function of the D7 brane and the Hubble parameter of the expanding background. We show that with $N_c\sim 10^7$ and $g_{\rm{}_{UV}}\sim {\rm few}\times 0.1$, the predictions of the inflationary model is consistent with the most stringent constraints placed on the inflationary models by Planck 2018. Although the model is capable of producing a large amount of gravitational waves, $r\simeq 0.01$, the displacement of the canonical mass dimension one scalar field remains below the Planck mass, in violation of the Lyth bound.

$f(P)$ gravity is a novel extension of ECG in which the Ricci scalar in the action is replaced by a function of the curvature invariant $P$ which represents the contractions of the Riemann tensor at the cubic order \cite{p}. The present work is concentrated on bounding some $f(P)$ gravity models using the concept of energy conditions where the functional forms of $f(P)$ are represented as \textbf{a)} $f(P) = \alpha \sqrt{P}$, and \textbf{b)} $f(P) = \alpha \exp (P)$, where $\alpha$ is the sole model parameter. Energy conditions are interesting linear relationships between pressure and density and have been extensively employed to derive interesting results in Einstein's gravity, and are also an excellent tool to impose constraints on any cosmological model. To place the bounds, we ensured that the energy density must remain positive, the pressure must remain negative, and the EoS parameter must attain a value close to $-1$ to make sure that the bounds respect the accelerated expansion of the Universe and are also in harmony with the latest observational data. We report that for both the models, suitable parameter spaces exist which satisfy the aforementioned conditions and therefore posit the $f(P)$ theory of gravity to be a promising modified theory of gravitation.

L. Tan (Florida), N. C. Tsamis (Crete), R. P. Woodard (Florida)

42 pages, 2 figures, 10 tables, uses LaTeX 2e. This paper is dedicated to Stanley Deser on the occasion of his 90th birthday

Although matter contributions to the graviton self-energy $-i[\mbox{}^{\mu\nu} \Sigma^{\rho\sigma}](x;x')$ must be separately conserved on $x^{\mu}$ and ${x'}^{\mu}$, graviton contributions obey the weaker constraint of the Ward identity, which involves a divergence on both coordinates. On a general homogeneous and isotropic background this leads to just four structures functions for matter contributions but nine structure functions for graviton contributions. We propose a convenient parameterization for these nine structure functions. We also apply the formalism to explicit one loop computations of $-i[\mbox{}^{\mu\nu} \Sigma^{\rho\sigma}](x;x')$ on de Sitter background, one of the contributions from a massless, minimally coupled scalar and the other for the contribution from gravitons in the simplest gauge. We also specialize the linearized, quantum-corrected Einstein equation to the graviton mode function and to the gravitational response to a point mass.

James B. Hartle

23p, 4 figures, revtex4, history of ths article on title page. arXiv admin note: text overlap with arXiv:1805.12246

In a fundamental formulation of the quantum mechanics of a closed system such as the universe as a whole, three forms of information are needed to make predictions for the probabilities of alternative time histories of the closed system . These are the action functional of the elementary particles, the quantum istate of the universe, and the description of our specific history. We discuss the origin of the "quasiclassical realm" of familiar experience and Hamiltonian quantum mechanics with its preferred time in such a formulation of quantum cosmology. It is argued that these features of the universe are not general properties of quantum theory, but rather approximate features that are emergent after the Planck time as a consequence of theories of the closed system's quantum state and dynamics.