We present the first characterization of the statistical relationship between a large sample of novae in M31 and their progenitor stellar populations in the form of a delay time distribution. To this end, we leverage the spatially resolved stellar age distribution of the M31 disk derived from deep Hubble Space Telescope photometry by the Panchromatic Hubble Andromeda Treasury survey and a large catalog of novae in M31. Our delay time distribution has two statistically significant detections: one population of nova progenitors with ages between 2 and 3.2 Gyr and an unnormalized rate of (3.7 - 3.5 + 6.8 ± 2.1) × 10 - 9 events M ⊙ [-1], and another of ages between 7.9 Gyr and the age of the Universe with (4.8 - 0.9 + 1.0 ± 0.2) × 10 - 9 events M ⊙ [-1] (uncertainties are statistical and systematic, respectively). Together with the upper limits we derive at other time bins, these detections are consistent with either a constant production efficiency or a higher production efficiency of novae at earlier delay times.

Our interpretation of terrestrial exoplanet atmospheric spectra will always be limited by the accuracy of the data we use as input in our forward and retrieval models. Ultraviolet molecular absorption cross sections are one category of these essential model inputs; however, they are often poorly characterized at the longest wavelengths relevant to photodissociation. Photolysis reactions dominate the chemical kinetics of temperate terrestrial planet atmospheres. One molecule of particular importance is CO2, which is likely present in all terrestrial planet atmospheres. The photolysis of CO2 can introduce CO and O, as well as shield tropospheric water vapor from undergoing photolysis. This is important because H2O photolysis produces OH, which serves as a major reactive sink to many atmospheric trace gases. Here, we construct CO2 cross-section prescriptions at 195 K and 300 K extrapolated beyond 200 nm from measured cross sections. We compare results from the implementation of these new cross sections to the most commonly used CO2 prescriptions for temperate terrestrial planets with Archean-like atmospheres. We generally find that the observational consequences of CO2 dissociation beyond 200 nm are minimal so long as our least conservative (highest opacity) prescription can be ruled out. Moreover, implementing our recommended extended CO2 cross sections does not substantially alter previous results that show the consequential photochemical impact of extended H2O cross sections.

The structure of the low redshift Universe is dominated by a multiscale void distribution delineated by filaments and walls of galaxies. The characteristics of voids, such as morphology, average density profile, and correlation function, can be used as cosmological probes. However, their physical properties are difficult to infer due to shot noise and the general lack of tracer particles used to define them. In this work, we construct a robust, topology-based void-finding algorithm that utilizes Persistent Homology to detect persistent features in the data. We apply this approach to a volume-limited subsample of galaxies in the SDSS I/II Main Galaxy catalog with the r-band absolute magnitude brighter than M r = -20.19, and a set of mock catalogs constructed using the Horizon Run 4 cosmological N-body simulation. We measure the size distribution of voids, their averaged radial profile, sphericity, and the centroid nearest neighbor separation, using conservative values for the threshold and persistence. We find 32 topologically robust voids in the SDSS data over the redshift range 0.02 ≤ z ≤ 0.116, with effective radii in the range 21-56 h [-1] Mpc. The median nearest neighbor void separation is found to be ∼57 h [-1] Mpc, and the median radial void profile is consistent with the expected shape from the mock data.

M15 is a globular cluster with a known spread in neutron-capture elements. This paper presents abundances of neutron-capture elements for 62 stars in M15. Spectra were obtained with the Michigan/Magellan Fiber System spectrograph, covering a wavelength range from ∼4430 to 4630 Å. Spectral lines from Fe i, Fe ii, Sr i, Zr ii, Ba ii, La ii, Ce ii, Nd ii, Sm ii, Eu ii, and Dy ii were measured, enabling classifications and neutron-capture abundance patterns for the stars. Of the 62 targets, 44 are found to be highly Eu-enhanced r-II stars, another 17 are moderately Eu-enhanced r-I stars, and one star is found to have an s-process signature. The neutron-capture patterns indicate that the majority of the stars are consistent with enrichment by the r-process. The 62 target stars are found to show significant star-to-star spreads in Sr, Zr, Ba, La, Ce, Nd, Sm, Eu, and Dy, but no significant spread in Fe. The neutron-capture abundances are further found to have slight correlations with sodium abundances from the literature, unlike what has been previously found; follow-up studies are needed to verify this result. The findings in this paper suggest that the Eu-enhanced stars in M15 were enhanced by the same process, that the nucleosynthetic source of this Eu pollution was the r-process, and that the r-process source occurred as the first generation of cluster stars was forming.

It is well known that the power spectrum is not able to fully characterize the statistical properties of non-Gaussian density fields. Recently, many different statistics have been proposed to extract information from non-Gaussian cosmological fields that perform better than the power spectrum. The Fisher matrix formalism is commonly used to quantify the accuracy with which a given statistic can constrain the value of the cosmological parameters. However, these calculations typically rely on the assumption that the sampling distribution of the considered statistic follows a multivariate Gaussian distribution. In this work, we follow Sellentin & Heavens and use two different statistical tests to identify non-Gaussianities in different statistics such as the power spectrum, bispectrum, marked power spectrum, and wavelet scattering transform (WST). We remove the non-Gaussian components of the different statistics and perform Fisher matrix calculations with the Gaussianized statistics using Quijote simulations. We show that constraints on the parameters can change by a factor of ∼2 in some cases. We show with simple examples how statistics that do not follow a multivariate Gaussian distribution can achieve artificially tight bounds on the cosmological parameters when using the Fisher matrix formalism. We think that the non-Gaussian tests used in this work represent a powerful tool to quantify the robustness of Fisher matrix calculations and their underlying assumptions. We release the code used to compute the power spectra, bispectra, and WST that can be run on both CPUs and GPUs.

Overshooting of turbulent motions from convective regions into adjacent stably stratified zones plays a significant role in stellar interior dynamics, as this process may lead to mixing of chemical species and contribute to the transport of angular momentum and magnetic fields. We present a series of fully nonlinear, three-dimensional (3D) anelastic simulations of overshooting convection in a spherical shell that are focused on the dependence of the overshooting dynamics on the density stratification and the rotation, both key ingredients in stars that however have not been studied systematically together via global simulations. We demonstrate that the overshoot lengthscale is not simply a monotonic function of the density stratification in the convective region, but instead it depends on the ratio of the density stratifications in the two zones. Additionally, we find that the overshoot lengthscale decreases with decreasing Rossby number Ro and scales as Ro0.23 while it also depends on latitude with higher Rossby cases leading to a weaker latitudinal variation. We examine the mean flows arising due to rotation and find that they extend beyond the base of the convection zone into the stable region. Our findings may provide a better understanding of the dynamical interaction between stellar convective and radiative regions, and motivate future studies particularly related to the solar tachocline and the implications of its overlapping with the overshoot region.

We study the nature of energy release and transfer for two sub-A class solar microflares observed during the second Focusing Optics X-ray Solar Imager (FOXSI-2) sounding rocket flight on 2014 December 11. FOXSI is the first solar-dedicated instrument to utilize focusing optics to image the Sun in the hard X-ray (HXR) regime, sensitive to energies of 4-20 keV. Through spectral analysis of the microflares using an optically thin isothermal plasma model, we find evidence for plasma heated to ~10 MK and emission measures down to ~10[44] cm[-3]. Though nonthermal emission was not detected for the FOXSI-2 microflares, a study of the parameter space for possible hidden nonthermal components shows that there could be enough energy in nonthermal electrons to account for the thermal energy in microflare 1, indicating that this flare is plausibly consistent with the standard thick-target model. With a solar-optimized design and improvements in HXR focusing optics, FOXSI-2 offers approximately five times greater sensitivity at 10 keV than the Nuclear Spectroscopic Telescope Array for typical microflare observations and allows for the first direct imaging spectroscopy of solar HXRs with an angular resolution at scales relevant for microflares. Harnessing these improved capabilities to study small-scale events, we find evidence for spatial and temporal complexity during a sub-A class flare. This analysis, combined with contemporaneous observations by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory, indicates that these microflares are more similar to large flares in their evolution than to the single burst of energy expected for a nanoflare.

Solar flares are explosive releases of magnetic energy. Hard X-ray (HXR) flare emission originates from both hot (millions of Kelvin) plasma and nonthermal accelerated particles, giving insight into flare energy release. The Nuclear Spectroscopic Telescope ARray (NuSTAR) utilizes direct-focusing optics to attain much higher sensitivity in the HXR range than that of previous indirect imagers. This paper presents 11 NuSTAR microflares from two active regions (AR 12671 on 2017 August 21 and AR 12712 on 2018 May 29). The temporal, spatial, and energetic properties of each are discussed in context with previously published HXR brightenings. They are seen to display several "large flare" properties, such as impulsive time profiles and earlier peak times in higher-energy HXRs. For two events where the active region background could be removed, microflare emission did not display spatial complexity; differing NuSTAR energy ranges had equivalent emission centroids. Finally, spectral fitting showed a high-energy excess over a single thermal model in all events. This excess was consistent with additional higher-temperature plasma volumes in 10/11 microflares and only with an accelerated particle distribution in the last. Previous NuSTAR studies focused on one or a few microflares at a time, making this the first to collectively examine a sizable number of events. Additionally, this paper introduces an observed variation in the NuSTAR gain unique to the extremely low livetime (<1%) regime and establishes a correction method to be used in future NuSTAR solar spectral analysis.

Cosmic-ray (CR) sources temporarily enhance the relativistic particle density in their vicinity over the background distribution accumulated from the Galaxy-wide past injection activity and propagation. If individual sources are close enough to the solar system, their localized enhancements may present as features in the measured spectra of the CRs and in the associated secondary electromagnetic emissions. Large-scale loop-like structures visible in the radio sky are possible signatures of such nearby CR sources. If so, these loops may also have counterparts in the high-latitude γ-ray sky. Using ~10 yr of data from the Fermi Large Area Telescope, applying Bayesian analysis including Gaussian Processes, we search for extended enhanced emission associated with putative nearby CR sources in the energy range from 1 GeV to 1 TeV for the sky region |b| > 30°. We carefully control the systematic uncertainty due to imperfect knowledge of the interstellar gas distribution. Radio Loop IV is identified for the first time as a γ-ray emitter, and we also find significant emission from Loop I. Strong evidence is found for asymmetric features about the Galactic l = 0° meridian that may be associated with parts of the so-called "Fermi Bubbles," and some evidence is also found for γ-ray emission from other radio loops. Implications for the CRs producing the features and possible locations of the sources of the emissions are discussed.

Herein we present a laboratory rotational study of cyanoacetic acid (CH2(CN)C(O)OH), an organic acid as well as a -CN bearing molecule, that is a candidate molecular system to be detected in the interstellar medium (ISM). Our investigation aims to provide direct experimental frequencies of cyanoacetic acid to guide its eventual astronomical search in low-frequency surveys. Using different jet-cooled rotational spectroscopic techniques in the time domain, we have determined a precise set of the relevant rotational spectroscopic constants, including the [14]N nuclear quadrupole coupling constants for the two distinct structures, cis- and gauche- cyanoacetic acid. We believe this work will potentially allow the detection of cyanoacetic acid in the interstellar medium, whose rotational features have remained unknown until now.

Hot stars with hot Jupiters have a wide range of obliquities, while cool stars with hot Jupiters tend to have low obliquities. An enticing explanation for this pattern is tidal realignment of the cool host stars, although this explanation assumes that obliquity damping occurs faster than orbital decay, an assumption that needs further exploration. Here we revisit this tidal realignment problem, building on previous work identifying a low-frequency component of the time-variable tidal potential that affects the obliquity but not the orbital separation. We adopt a recent empirically based model for the stellar tidal quality factor and its sharp increase with forcing frequency. This leads to enhanced dissipation at low frequencies, and efficient obliquity damping. We model the tidal evolution of 46 observed hot Jupiters orbiting cool stars. A key parameter is the stellar age, which we determine in a homogeneous manner for the sample, taking advantage of Gaia DR2 data. We explore a variety of tidal histories and futures for each system, finding in most cases that the stellar obliquity is successfully damped before the planet is destroyed. A testable prediction of our model is that hot Jupiter hosts with orbital periods shorter than 2-3 days should have obliquities much smaller than 1°. With the possible exception of WASP-19b, the predicted future lifetimes of the planets range from 108 yr to more than 1010 yr. Thus, our model implies that these hot Jupiters are probably not in immediate danger of being devoured by their host stars while they are on the main sequence.

Since its launch, the Alpha Magnetic Spectrometer-02 (AMS-02) has delivered outstanding quality measurements of the spectra of cosmic-ray (CR) species (p ¯ , e [±], and nuclei, 1H-8O, 10Ne, 12Mg, 14Si) which resulted in a number of breakthroughs. One of the latest long-awaited surprises is the spectrum of 26Fe just published by AMS-02. Because of the large fragmentation cross section and large ionization energy losses, most of CR iron at low energies is local and may harbor some features associated with relatively recent supernova (SN) activity in the solar neighborhood. Our analysis of the new AMS-02 results, together with Voyager 1 and ACE-CRIS data, reveals an unexpected bump in the iron spectrum and in the Fe/He, Fe/O, and Fe/Si ratios at 1-2 GV, while a similar feature in the spectra of He, O, and Si and in their ratios is absent, hinting at a local source of low-energy CRs. The found excess extends the recent discoveries of radioactive [60]Fe deposits in terrestrial and lunar samples and in CRs. We provide an updated local interstellar spectrum (LIS) of iron in the energy range from 1 MeV nucleon[-1] to ~10 TeV nucleon[-1]. Our calculations employ the GALPROP-HELMOD framework, which has proved to be a reliable tool in deriving the LIS of CR p ¯ , e [-], and nuclei Z ⩽ 28.

A recently observed bump in the cosmic-ray (CR) spectrum from 0.3 to 30 TV is likely caused by a stellar bow shock that reaccelerates preexisting CRs, which further propagate to the Sun along the magnetic field lines. Along their way, these particles generate an Iroshnikov-Kraichnan (I-K) turbulence that controls their propagation and sustains the bump. Ad hoc fitting of the bump shape requires six adjustable parameters. Our model requires none, merely depending on three physical unknowns that we constrain using the fit. These are the shock Mach number, M, its size, l ⊥, and the distance to it, ζ obs. Altogether, they define the bump rigidity R 0. With M ≈ 1.5-1.6 and R 0 ≈ 4.4 TV, the model fits the data with ≈0.08% accuracy. The fit critically requires the I-K spectrum predicted by the model and rules out the alternatives. These attributes of the fit make an accidental agreement highly unlikely. In turn, the R 0 and M derived from the fit impose the distance-size relation on the shock: ζ obs (pc) ~ 10 2 l ⊥ (pc) . For sufficiently large bow shocks, l ⊥ = 10[-3]-10[-2] pc, we find the distance of ζ obs = 3-10 pc. Three promising stars in this range are the Scholz's Star at 6.8 pc, Epsilon Indi at 3.6 pc, and Epsilon Eridani at 3.2 pc. Based on their current positions and velocities, we propose that Epsilon Indi and Epsilon Eridani can produce the observed spectral bump. Moreover, Epsilon Eridani's position is only ~6°.7 off of the magnetic field direction in the solar neighborhood, which also changes the CR arrival direction distribution. Given the proximity of these stars, the bump appearance may change in a relatively short time.

We present elemental abundance data of C, N, O, Na, Mg, Al, Ca, and Cr in Genesis silicon targets. For Na, Mg, Al, and Ca, data from three different SW regimes are also presented. Data were obtained by backside depth profiling using Secondary Ion Mass Spectrometry. The accuracy of these measurements exceeds those obtained by in-situ observations; therefore the Genesis data provide new insights into elemental fractionation between Sun and solar wind, including differences between solar wind regimes. We integrate previously published noble gas and hydrogen elemental abundances from Genesis targets, as well as preliminary values for K and Fe. The abundances of the solar wind elements measured display the well-known fractionation pattern that correlates with each elements' First Ionization Potential (FIP). When normalized either to spectroscopic photospheric solar abundances or to those derived from CI-chondritic meteorites, the fractionation factors of low-FIP elements (K, Na, Al, Ca, Cr, Mg, Fe) are essentially identical within uncertainties, but the data are equally consistent with an increasing fractionation with decreasing FIP. The elements with higher FIPs between ~11 and ~16 eV (C, N, O, H, Ar, Kr, Xe) display a relatively well-defined trend of increasing fractionation with decreasing FIP, if normalized to modern 3D photospheric model abundances. Among the three Genesis regimes, the Fast SW displays the least elemental fractionation for almost all elements (including the noble gases) but differences are modest: for low-FIP elements the precisely measured Fast-Slow SW variations are less than 3%.

We apply the helioseismic methodology of Fourier Legendre decomposition to 88 months of Dopplergrams obtained by the Helioseismic and Magnetic Imager (HMI) as the basis of inferring the depth variation of the mean meridional flow, as averaged between 20° and 60° latitude and in time, in both the northern and southern hemispheres. We develop and apply control procedures designed to assess and remove center-to-limb artifacts using measurements obtained by performing the analysis with respect to artificial poles at the east and west limbs. Forward modeling is carried out using sensitivity functions proportional to the mode kinetic energy density to evaluate the consistency of the corrected frequency shifts with models of the depth variation of the meridional circulation in the top half of the convection zone. The results, taken at face value, imply substantial differences between the meridional circulation in the northern and southern hemispheres. The inferred presence of a return (equatorward propagating) flow at a depth of approximately 40 Mm below the photosphere in the northern hemisphere is surprising and appears to be inconsistent with many other helioseismic analyses. This discrepancy may be the result of the inadequacy of our methodology to remove systematic errors in HMI data. Our results appear to be at least qualitatively similar to those by Gizon et al., which point to an anomaly in HMI data that is not present in MDI or GONG data.