Papers for Thursday, Mar 18 2021

In dense star clusters, such as globular and open clusters, dynamical interactions between stars and black holes (BHs) can be extremely frequent, leading to various astrophysical transients. Close encounters between a star and a stellar mass BH make it possible for the star to be tidally disrupted by the BH. Due to the relative low mass of the BH and the small cross section of the tidal disruption event (TDE) for cases with high penetration, disruptions caused by close encounters are usually partial disruptions. The existence of the remnant stellar core and its non-negligible mass compared to the stellar mass BH alters the accretion process significantly. We study this problem with SPH simulations using the code {\tt Phantom}, with the inclusion of radiation pressure, which is important for small mass BHs. Additionally, we develop a new, more general method of computing the fallback rate which does not rely on any approximation. Our study shows that the powerlaw slope of the fallback rate has a strong dependence on the mass of the BH in the stellar mass BH regime. Furthermore, in this regime, self-gravity of the fallback stream and local instabilities become more significant, and cause the disrupted material to collapse into small clumps before returning to the BH. This results in an abrupt increase of the fallback rate, which can significantly deviate from a powerlaw. Our results will help in the identification of TDEs by stellar mass BHs in dense clusters.

The Transit Timing Variation (TTV) technique is a powerful dynamical tool to measure exoplanetary masses by analysing transit light curves. We assessed the transit timing performances of the Ariel Fine Guidance Sensors (FGS1/2) based on the simulated light curve of a bright, 55 Cnc, and faint, K2-24, planet-hosting star. We estimated through a Markov-Chain Monte-Carlo analysis the transit time uncertainty at the nominal cadence of 1 second and, as a comparison, at a 30 and 60-s cadence. We found that at the nominal cadence Ariel will be able to measure the transit time with a precision of about 12s and 34s, for a star as bright as 55 Cnc and K2-24, respectively. We then ran dynamical simulations, also including the Ariel timing errors, and we found an improvement on the measurement of planetary masses of about $20-30\%$ in a K2-24-like planetary system through TTVs. We also simulated the conditions that allow us to detect the TTV signal induced by an hypothetical external perturber within the mass range between Earth and Neptune using 10 transit light curves by Ariel.

We present spatially-resolved two-dimensional maps and radial trends of the stellar populations and kinematics for a sample of six compact elliptical galaxies (cE) using spectroscopy from the Keck Cosmic Web Imager (KCWI). We recover their star formation histories, finding that all except one of our cEs are old and metal rich, with both age and metallicity decreasing toward their outer radii. We also use the integrated values within one effective radius to study different scaling relations. Comparing our cEs with others from the literature and from simulations we reveal the formation channel that these galaxies might have followed. All our cEs are fast rotators, with relatively high rotation values given their low ellipticites. In general, the properties of our cEs are very similar to those seen in the cores of more massive galaxies, and in particular, to massive compact galaxies. Five out of our six cEs are the result of stripping a more massive (compact or extended) galaxy, and only one cE is compatible with having been formed intrinsically as the low-mass, compact object that we see today. These results further confirm that cEs are a mixed-bag of galaxies that can be formed following different formation channels, reporting for the first time an evolutionary link within the realm of compact galaxies (at all stellar masses).

As the number of observed merging binary black holes (BHs) grows, accurate models are required to disentangle multiple formation channels. In models with isolated binaries, important uncertainties remain regarding the stability of mass transfer (MT) and common-envelope (CE) evolution. To study some of these uncertainties, we have computed simulations using MESA of a $30M_\odot$, low metallicity ($Z_\odot/10$) star with a BH companion. We developed a prescription to compute MT rates including possible outflows from outer Lagrangian points, and a method to self-consistently determine the core-envelope boundary in the case of CE evolution. We find that binaries survive a CE only if unstable MT happens after the formation of a deep convective envelope, resulting in a narrow range (0.2 dex) in period for envelope ejection. All cases where interaction is initiated with a radiative envelope have large binding energies ($\sim 10^{50}$ erg), and merge during CE even under the assumption that all the internal and recombination energy of the envelope, as well as the energy from an inspiral, is used for ejection. This is independent of core helium ignition for the donor, a condition under which various rapid-population synthesis calculations assume a successful ejection is possible. Moreover, we find that the critical mass ratio for instability is such that for periods between $\sim 1-1000$ days merging binary BHs can be formed via stable MT. A large fraction of these systems overflow their L$_2$ equipotential, in which case we find stable MT produces merging binary BHs even under extreme assumptions of mass and angular momentum outflows. Our conclusions are limited to the study of one donor star, but suggest that population synthesis calculations overestimate the formation rate of merging binary BHs produced by CE evolution, and that stable MT could dominate the rate from isolated binaries.

We analyze the rest-optical emission-line ratios of z~1.5 galaxies drawn from the MOSDEF survey. Based on composite spectra we investigate the mass-metallicity relation (MZR) at z~1.5 and measure its evolution to z=0. When using gas-phase metallicities based on the N2 line ratio, we find that the difference between the MZR at z~1.5 and z=0 depends on stellar mass, with a difference of $\Delta\rm log(\rm O/H)\sim0.25$ dex at $M_* <10^{9.75}M_{\odot}$ down to $\Delta\rm log(\rm O/H)\sim0.05$ for galaxies with $M_*>10^{10.5}M_{\odot}$. In contrast, the O3N2-based MZR shows a constant offset of $\Delta\rm log(\rm O/H)\sim0.30$ for all masses, consistent with previous MOSDEF results based on independent metallicity indicators, and suggesting that O3N2 provides a more robust metallicity calibration for our z~1.5 sample. We investigated the secondary dependence of the MZR on SFR by measuring correlated scatter about the mean $M_*$-sSFR and $M_*-\log(\rm O3N2)$ relations. We find an anti-correlation between $\log(\rm O/H)$ and sSFR offsets, indicating the presence of a $M_*$-SFR-Z relation, though with limited significance. Additionally, we find that our z~1.5 stacks lie along the z=0 metallicity sequence at fixed $\mu=\log(M_*/M_{\odot})-0.6\times\log(\rm SFR/M_{\odot}yr^{-1})$, suggesting that the z~1.5 stacks can be described by the z=0 fundamental metallicity relation (FMR). However, using different calibrations can shift the metallicities enough to be inconsistent with the local FMR, indicating that the choice of calibration is essential for understanding metallicity evolution with redshift. Finally, understanding how [NII]/H$\alpha$ scales with galaxy properties is crucial for describing the effects of blended [NII] and H$\alpha$ on redshift and H$\alpha$ flux measurements in future large surveys utilizing low-resolution spectra such as with Euclid and the Roman Space Telescope.

We present ASTRODEEP-GS43, a new multiwavelength photometric catalogue of the GOODS-South field, which builds and improves upon the previously released CANDELS catalogue. We provide photometric fluxes and corresponding uncertainties in 43 optical and infrared bands (25 wide and 18 medium filters), as well as photometric redshifts and physical properties of the 34930 CANDELS $H$-detected objects, plus an additional sample of 178 $H$-dropout sources, of which 173 are $Ks$-detected and 5 IRAC-detected. We keep the CANDELS photometry in 7 bands (CTIO $U$, Hubble Space Telescope WFC3 and ISAAC-$K$), and measure from scratch the fluxes in the other 36 (VIMOS, HST ACS, HAWK-I $Ks$, Spitzer IRAC, and 23 from Subaru SuprimeCAM and Magellan-Baade Fourstar) with state-of-the-art techniques of template-fitting. We then compute new photometric redshifts with three different software tools, and take the median value as best estimate. We finally evaluate new physical parameters from SED fitting, comparing them to previously published ones. Comparing to a sample of 3931 high quality spectroscopic redshifts, for the new photo-$z$'s we obtain a normalized median absolute deviation (NMAD) of 0.015 with 3.01$\%$ of outliers (0.011, 0.22$\%$ on the bright end at $I814$<22.5), similarly to the best available published samples of photometric redshifts, such as the COSMOS UltraVISTA catalogue. The ASTRODEEP-GS43 results are in qualitative agreement with previously published catalogues of the GOODS-South field, improving on them particularly in terms of SED sampling and photometric redshift estimates. The catalogue is available for download from the Astrodeep website.

As weak lensing surveys become deeper, they reveal more non-Gaussian aspects of the convergence field which require statistics beyond the power spectrum to extract. In Cheng et al. (2020), we showed that the scattering transform, a novel statistic borrowing concepts from convolutional neural networks, is a powerful tool to perform cosmological parameter estimation. Here, we extend this analysis to explore its sensitivity to dark energy and neutrino mass parameters with weak lensing surveys. We first use image synthesis to show visually that the scattering transform provides a better statistical vocabulary to characterize lensing mass maps compared to the power spectrum and bispectrum. We then show that it outperforms those two estimators in the low-noise regime. When the noise level increases and non-Gaussianity is diluted, though the constraint is not significantly tighter than that of the bispectrum, the scattering coefficients have much more Gaussian likelihood, which is essential for accurate cosmological inference. We argue that the scattering coefficients are preferred statistics considering both constraining power and likelihood properties.

We report the results of eROSITA and NICER observations of the June 2020 outburst of the Be/X-ray binary pulsar RX J0529.8-6556 in the Large Magellanic Cloud, along with the analysis of archival X-ray and optical data from this source. We find two anomalous features in the system's behavior. First, the pulse profile observed by NICER during maximum luminosity is similar to that observed by XMM-Newton in 2000, despite the fact that the X-ray luminosity was different by two orders of magnitude. By contrast, a modest decrease in luminosity in the 2020 observations generated a significant change in pulse profile. Second, we find that the historical optical outbursts are not strictly periodic, as would be expected if the outbursts were triggered by periastron passage, as is generally assumed. The optical peaks are also not coincident with the X-ray outbursts. We suggest that this behavior may result from a misalignment of the Be star disk and the orbital plane, which might cause changes in the timing of the passage of the neutron star through the disk as it precesses. We conclude that the orbital period of the source remains unclear.

The intensity of the Cosmic UV background (UVB), coming from all sources of ionising photons such as star-forming galaxies and quasars, determines the thermal evolution and ionization state of the intergalactic medium (IGM) and is, therefore, a critical ingredient for models of cosmic structure formation. Most of the previous estimates are based on the comparison between observed and simulated Lyman-$\alpha$ forest. We present the results of an independent method to constrain the product of the UVB photoionisation rate and the covering fraction of Lyman limit systems (LLSs) by searching for the fluorescent Lyman-$\alpha$ emission produced by self-shielded clouds. Because the expected surface brightness is well below current sensitivity limits for direct imaging, we developed a new method based on three-dimensional stacking of the IGM around Lyman-$\alpha$ emitting galaxies (LAEs) between 2.9<z<6.6 using deep MUSE observations. Combining our results with covering fractions of LLSs obtained from mock cubes extracted from the EAGLE simulation, we obtain new and independent constraints on the UVB at z>3 that are consistent with previous measurements, with a preference for relatively low UVB intensities at z=3, and which suggest a non-monotonic decrease of $\Gamma$HI with increasing redshift between 3<z<5. This could suggest a possible tension between some UVB models and current observations which however require deeper and wider observations in Lyman-$\alpha$ emission and absorption to be confirmed. Assuming instead a value of UVB from current models, our results constrain the covering fraction of LLSs at 3<z<4.5 to be less than 25% within 150kpc from LAEs.

We solve the equations of two-dimensional hydrodynamics describing a circumbinary disk accreting onto an eccentric, equal-mass binary. We compute the time rate of change of the binary semi-major axis $a$ and eccentricity $e$ over a continuous range of eccentricities spanning $e=0$ to $e=0.9$. We find that binaries with initial eccentricities $e_0\lesssim 0.1$ tend to $e=0$, where the binary semi-major axis expands. All others are attracted to $e \approx 0.4$, where the binary semi-major axis decays. The $e \approx 0.4$ attractor is caused by a rapid change in the disk response from an apsidally locked, nearly symmetric state to a precessing asymmetric state. The state change causes the time rates of change $\dot{a}$ and $\dot{e}$ to steeply change sign at the same critical eccentricity resulting in an attracting solution where $\dot{a} = \dot{e} = 0$. This does not, however, result in a stalled, eccentric binary. Due to a finite transition time between disk states, the binary eccentricity continually overshoots the attracting eccentricity resulting in a constant drift of the semi-major axis. For fiducial disk parameters, binaries with $e_0 \gtrsim 0.1$ evolve towards and then oscillate around $e \approx 0.4$ where they steadily shrink in semi-major axis. Because unequal mass binaries grow towards equal mass through preferential accretion, our results are applicable to a wide range of initial binary mass ratios. Hence, these findings merit further investigations of this disk transition; understanding its dependence on disk parameters is vital for determining the fate of binaries undergoing orbital evolution with a circumbinary disk.

Direct imaging of exoplanets is usually limited by quasi-static speckles. These uncorrected aberrations in a star's point spread function (PSF) obscure faint companions and limit the sensitivity of high-contrast imaging instruments. Most current approaches to processing differential imaging sequences like angular differential imaging (ADI) and spectral differential imaging (SDI) produce a self-calibrating dataset that are combined in a linear least squares solution to minimize the noise. Due to temporal and chromatic evolution of a telescope's PSF, the best correlated reference images are usually the most contaminated by the planet, leading to self-subtraction and reducing the planet throughput. In this paper, we present an algorithm that directly optimizes the non-linear equation for planet signal to noise ratio (SNR). This new algorithm does not require us to reject adjacent reference images and optimally balances noise reduction with self-subtraction. We then show how this algorithm can be applied to multiple images simultaneously for a further reduction in correlated noise, directly maximizing the SNR of the final combined image. Finally, we demonstrate the technique on an illustrative sequence of HR8799 using the new Julia-based Signal to Noise Analysis Pipeline (SNAP). We show that SNR optimization can provide up to a $5\times$ improvement in contrast close to the star. Applicable to both new and archival data, this technique will allow for the detection of lower mass, and closer in companions, or achieve the same sensitivity with less telescope time.

We provide a new method to derive heavy element abundances based on the unique suite of nebular lines in the mid- to far-infrared (IR) range. Using grids of photo-ionisation models that cover a wide range in O/H and N/O abundances, and ionisation parameter, our code HII-Chi-mistry-IR (HCm-IR) provides model-based abundances based on extinction free and temperature insensitive tracers, two significant advantages over optical diagnostics. The code is probed using a sample of 56 galaxies observed with $Spitzer$ and $Herschel$ covering a wide range in metallicity, $7.2 \lesssim 12+\log(O/H) \lesssim 8.9$. The IR model-based metallicities obtained are robust within a scatter of 0.03 dex when the hydrogen recombination lines, which are typically faint transitions in the IR range, are not available. When compared to the optical abundances obtained with the direct method, model-based methods, and strong-line calibrations, HCm-IR estimates show a typical dispersion of ~0.2 dex, in line with previous studies comparing IR and optical abundances, a do not introduce a noticeable systematic above $12+\log(O/H) \gtrsim 7.6$. This accuracy can be achieved using the lines [SIV]$_{10.5 \mu m}$, [SIII]$_{18.7,33.5 \mu m}$, [NeIII]$_{15.6 \mu m}$ and [NeII]$_{12.8 \mu m}$. Additionally, HCm-IR provides an independent N/O measurement when the [OIII]$_{52,88 \mu m}$ and [NIII]$_{57 \mu m}$ transitions are measured, and therefore the derived abundances in this case do not rely on particular assumptions in the N/O ratio. Large uncertainties (~0.4 dex) may affect the abundance determinations of galaxies at sub- or over-solar metallicities when a solar-like N/O ratio is adopted. Finally, the code has been applied to 8 galaxies located at $1.8 < z < 7.5$ with ground-based detections of far-IR lines redshifted in the submm range, revealing solar-like N/O and O/H abundances in agreement with recent studies.

Dust gaps and rings appear ubiquitous in bright protoplanetary disks. Disk-planet interaction with dust-trapping at the edges of planet-induced gaps is one plausible explanation. However, the sharpness of some observed dust rings indicate that sub-mm-sized dust grains have settled to a thin layer in some systems. We test whether or not such dust around gas gaps opened by planets can remain settled by performing three-dimensional, dust-plus-gas simulations of protoplanetary disks with an embedded planet. We find planets massive enough to open gas gaps stir small, sub-mm-sized dust grains to high disk elevations at the gap edges, where the dust scale-height can reach ~70% of the gas scale-height. We attribute this dust 'puff-up' to the planet-induced meridional gas flows previously identified by Fung & Chiang and others. We thus emphasize the importance of explicit 3D simulations to obtain the vertical distribution of sub-mm-sized grains around gas gaps opened by massive planets. We caution that the gas-gap-opening planet interpretation of well-defined dust rings is only self-consistent with large grains exceeding mm in size.

We present the second public data release (DR) of the VISTA EXtension to Auxiliary Surveys (VEXAS), where we classify objects into stars, galaxies and quasars based on an ensemble of machine learning algorithms. The aim of VEXAS is to build the widest multi-wavelength catalogue, providing reference magnitudes, colours and morphological information for a large number of scientific uses. We apply an ensemble of 32 different machine learning models, based on three different algorithms and on different magnitude sets, training samples and classification problems on the three VEXAS DR1 optical+infrared (IR) tables. The tables were created in DR1 cross-matching VISTA near-IR data with WISE far-IR data and with optical magnitudes from the Dark Energy Survey (VEXAS-DESW), the Sky Mapper Survey (VEXAS-SMW), and the PanSTARRS (VEXAS-PSW). We assemble a large table of spectroscopically confirmed objects (415 628 unique objects), based on the combination of 6 different spectroscopic surveys that we use for training. We develop feature imputation to classify also objects for which magnitudes in one or more bands are missing. We classify in total ~90 million objects in the Southern Hemisphere. Among these,~62.9M (~52.6M) are classified as 'high confidence' ('secure') stars, ~920k (~750k) as 'high confidence' ('secure') quasars and ~34.8M (~34.1M) as 'high confidence' ('secure') galaxies, with probabilities $p_{\rm class}\ge 0.7$ ($p_{\rm class}\ge 0.9$). The density of high-confidence extragalactic objects varies strongly with the survey depth: at $p_{\rm class}\ge 0.7$, there are 111/deg$^2$ quasars in the VEXAS-DESW footprint and 103/deg$^2$ in the VEXAS-PSW footprint, while only 10.7/deg$^2$ in the VEXAS-SM footprint. Improved depth in the midIR and coverage in the optical and nearIR are needed for the SM footprint that is not already covered by DESW and PSW.

The Earth's N2-dominated atmosphere is a very special feature. Firstly, N2 as main gas is unique on the terrestrial planets in the inner solar system and gives a hint for tectonic activity. Studying the origins of atmospheric nitrogen and its stability provides insights into the uniqueness of the Earth's habitat. Secondly, the coexistence of N2 and O2 within an atmosphere is unequaled in the entire solar system. Such a combination is strongly linked to the existence of aerobic lifeforms. The availability of nitrogen on the surface, in the ocean, and within the atmosphere can enable or prevent the habitability of a terrestrial planet, since nitrogen is vitally required by all known lifeforms. In the present work, the different origins of atmospheric nitrogen, the stability of nitrogen dominated atmospheres, and the development of early Earth's atmospheric N2 are discussed. We show why N2-O2-atmospheres constitute a biomarker not only for any lifeforms but for aerobic lifeforms, which was the first major step that led to higher developed life on Earth.

In intermediate polars (IPs), the intrinsic thermal emissions from white dwarfs (WDs) have typically been studied. Few reports have analyzed X-ray reflections from WDs. We recently developed an elaborate IP-reflection spectral model. Herein, we report the first application of a reflection model with an IP thermal model to the spectra of the brightest typical IP V1223 Sagittarii observed by the Suzaku and NuSTAR satellites. The model reasonably reproduces the spectra within the range of 5-78 keV and estimates the WD mass as 0.92$\pm$0.02 $M_\odot$. The WD mass estimated by the proposed model is consistent with that measured using an active galactic nuclei reflection model and a partial covering absorption model. However, the choice of incorrect parameter values, such as an unsuitable fitting energy band and an incorrect metal abundance, was found to introduce systematic errors (e.g., $<\sim$ 0.2 $M_\odot$ in the WD mass) in the WD mass measurement. Our spin phase-resolved analysis resulted in discoveries regarding the modulations of the equivalent width of the fluorescent iron K$_{\alpha}$ line and the angle between the post-shock accretion column and the line-of-sight (viewing angle). The viewing angle anti-correlates approximately with the X-ray flux and has average and semi-amplitude values of 55$^\circ$ and 7$^\circ$, respectively, which points toward two WD spin axis angles from the line-of-sight of 55$^\circ$ and 7$^\circ$, respectively. Both estimated spin axis angles are different from the reported system inclination of 24$^\circ$.

The accurate measurement of stellar masses over a wide range of galaxy properties is essential for better constraining models of galaxy evolution. Emission line galaxies (ELGs) tend to have better redshift estimates than continuum-selected objects and have been shown to span a large range of physical properties, including stellar mass. Using data from the 3D-HST Treasury program, we construct a carefully vetted sample of 4350 ELGs at redshifts 1.16<z<1.90. We combine the 3D-HST emission line fluxes with far-UV through near-IR photometry and use the MCSED spectral energy distribution fitting code to constrain the galaxies' physical parameters, such as their star formation rate (SFRs) and stellar masses. Our sample is consistent with the z~2 mass-metallicity relation. More importantly, we show there is a simple but tight correlation between stellar mass and absolute magnitude in a near-IR filter that will be particularly useful in quickly calculating accurate stellar masses for millions of galaxies in upcoming missions such as Euclid and the Nancy Grace Roman Space Telescope.

A sunspot catalogue was published by the Coimbra Astronomical Observatory (Portugal), now named Geophysical and Astronomical Observatory of the University of Coimbra, for the period 1929-1941. We digitalized data included in that catalogue and provide a machine-readable version. We show the reconstructions for the (total and hemispheric) sunspot number index and sunspot area according to this catalogue, comparing it with the sunspot number index (version 2) and Balmaceda sunspot area series (Balmaceda et al., J. Geophys. Res. 114, A07104, 2009). Moreover, we also compared the Coimbra catalogue with records made at the Royal Greenwich Observatory. The results demonstrate that the historical catalogue compiled by the Coimbra Astronomical Observatory contain reliable sunspot data and therefore can be considered for studies about solar activity.

Across a large range of scales, accreting sources show remarkably similar patterns of variability, most notably the log-normality of the luminosity distribution and the linear root-mean square (rms)-flux relationship. These results are often explained using the theory of propagating fluctuations in which fluctuations in the viscosity create perturbations in the accretion rate at all radii, propagate inwards and combine multiplicatively. While this idea has been extensively studied analytically in a linear regime, there has been relatively little numerical work investigating the non-linear behaviour. In this paper, we present a suite of stochastically driven 1-d $\alpha$-disc simulations, exploring the behaviour of these discs. We find that the eponymous propagating fluctuations are present in all simulations across a wide range of model parameters, in contradiction to previous work. Of the model parameters, we find by far the most important to be the timescale on which the viscosity fluctuations occur. Physically, this timescale will depend on the underlying physical mechanism, thought to be the magnetorotational instability (MRI). We find a close relationship between this fluctuation timescale and the break frequency in the power spectral density (PSD) of the luminosity, a fact which could allow observational probes of the behaviour of the MRI dynamo. We report a fitting formula for the break frequency as a function of the fluctuation timescale, the disc thickness and the mass of the central object.

The precise characterization of terrestrial atmospheres with the James Webb Space Telescope (JWST) is one of the utmost goals of exoplanet astronomy in the next decade. With JWST's impending launch, it is crucial we are well prepared to understand the subtleties of terrestrial atmospheres - particularly ones we may have not needed to consider before due to instrumentation limitations. In this work we show that patchy ice cloud variability is present in the upper atmospheres of M-dwarf terrestrial planets, particularly along the limbs. Here we test whether these variable clouds will introduce unexpected biases in the multi-epoch observations necessary to constrain atmospheric abundances. Using 3D ExoCAM general circulation models (GCMs) of TRAPPIST-1e, we simulate five different climates with varying pCO$_2$ to explore the strength of this variability. These models are post-processed using NASA Goddard's Planetary Spectrum Generator (PSG) and PandExo to generate simulated observations with JWST's NIRSpec PRISM mode at 365 different temporal outputs from each climate. Assuming the need for 10 transits of TRAPPIST-1e to detect molecular features at great confidence, we then use CHIMERA to retrieve on several randomly selected weighted averages of our simulated observations to explore the effect of multi-epoch observations with variable cloud cover along the limb on retrieved abundances. We find that the variable spectra do not affect retrieved abundances at detectable levels for our sample of TRAPPIST-1e models.

We present multi-wavelength fast timing observations of the black hole X-ray binary MAXI J1820+070 (ASASSN-18ey), taken with the Karl G. Jansky Very Large Array (VLA), Atacama Large Millimeter/Sub-Millimeter Array (ALMA), Very Large Telescope (VLT), New Technology Telescope (NTT), Neutron Star Interior Composition Explorer (NICER), and XMM-Newton. Our data set simultaneously samples ten different electromagnetic bands (radio - X-ray) over a 7-hour period during the hard state of the 2018-2019 outburst. The emission we observe is highly variable, displaying multiple rapid flaring episodes. To characterize the variability properties in our data, we implemented a combination of cross-correlation and Fourier analyses. We find that the emission is highly correlated between different bands, measuring time-lags ranging from hundreds of milliseconds between the X-ray/optical bands to minutes between the radio/sub-mm bands. Our Fourier analysis also revealed, for the first time in a black hole X-ray binary, an evolving power spectral shape with electromagnetic frequency. Through modelling these variability properties, we find that MAXI J1820+070 launches a highly relativistic ($\Gamma=6.81^{+1.06}_{-1.15}$) and confined ($\phi=0.45^{+0.13}_{-0.11}$ deg) jet, which is carrying a significant amount of power away from the system (equivalent to $\sim0.6 \, L_{1-100{\rm keV}}$). We additionally place constraints on the jet composition and magnetic field strength in the innermost jet base region. Overall, this work demonstrates that time-domain analysis is a powerful diagnostic tool for probing jet physics, where we can accurately measure jet properties with time-domain measurements alone.

The spectra arising from the disks of nova-like variables show many of the features seen in stellar atmospheres. They are typically modelled either from an appropriated weighted set of stellar atmospheres or a disk atmosphere with energy is dissipated near the disk plane, with the effective temperature distribution expected from a steady state accretion disk. However these models generally over-predict the depth of the Balmer jump and the slope of the spectrum in the ultraviolet. The problem is likely due to energy dissipation in the disk atmosphere, which produces a flatter vertical temperature profile than is observed in stars. Here, we provide validation for this hypothesis in the form of spectra generated using the stellar atmosphere code TLUSTY using a parametric prescription for energy dissipation as a function of depth and closely match the spectrum of the nova-like IX Vel over the wavelength range 1150-6000 \AA.

Starshade in formation flight with a space telescope is a rapidly maturing technology that would enable imaging and spectral characterization of small planets orbiting nearby stars in the not-too-distant future. While performance models of the starshade-assisted exoplanet imaging have been developed and used to design future missions, their results have not been verified from the analyses of synthetic images. Following a rich history of using community data challenges to evaluate image-processing capabilities in astronomy and exoplanet fields, the Starshade Technology Development to TRL5 (S5), a focused technology development activity managed by the NASA Exoplanet Exploration Program, is organizing and implementing a starshade exoplanet data challenge. The purpose of the data challenge is to validate the flow down of requirements from science to key instrument performance parameters and to quantify the required accuracy of noisy background calibration with synthetic images. This data challenge distinguishes itself from past efforts in the exoplanet field in that (1) it focuses on the detection and spectral characterization of small planets in the habitable zones of nearby stars, and (2) it develops synthetic images that simultaneously include multiple background noise terms -- some specific to starshade observations -- including residual starlight, solar glint, exozodiacal light, detector noise, as well as variability resulting from starshade's motion and telescope jitter. In this paper, we provide an overview of the design and rationale of the data challenge. Working with data challenge participants, we expect to achieve improved understanding of the noise budget and background calibration in starshade-assisted exoplanet observations in the context of both Starshade Rendezvous with Roman and HabEx.

A growing body of evidence suggests that the solar wind is powered to a large extent by an Alfv\'en-wave (AW) energy flux. AWs energize the solar wind via two mechanisms: heating and work. We use high-resolution direct numerical simulations of reflection-driven AW turbulence (RDAWT) in a fast-solar-wind stream emanating from a coronal hole to investigate both mechanisms. In particular, we compute the fraction of the AW power at the coronal base ($P_{\rm AWb}$) that is transferred to solar-wind particles via heating between the coronal base and heliocentric distance $r$, which we denote $\chi_{\rm H}(r)$, and the fraction that is transferred via work, which we denote $\chi_{\rm W}(r)$. We find that $\chi_{\rm W}(r_{\rm A})$ ranges from 0.15 to 0.3, where $r_{\rm A}$ is the Alfv\'en critical point. This value is small compared to~one because the Alfv\'en speed $v_{\rm A} $ exceeds the outflow velocity $U$ at $r<r_{\rm A}$, so the AWs race through the plasma without doing much work. At $r>r_{\rm A}$, where $v_{\rm A} < U$, the AWs are in an approximate sense "stuck to the plasma", which helps them do pressure work as the plasma expands. However, much of the AW power has dissipated by the time the AWs reach $r=r_{\rm A}$, so the total rate at which AWs do work on the plasma at $r>r_{\rm A}$ is a modest fraction of $P_{\rm AWb}$. We find that heating is more effective than work at $r<r_{\rm A}$, with $\chi_{\rm H}(r_{\rm A})$ ranging from 0.5 to 0.7. The reason that $\chi_{\rm H} \geq 0.5$ in our simulations is that an appreciable fraction of the local AW power dissipates within each Alfv\'en-speed scale height in RDAWT, and there are a few Alfv\'en-speed scale heights between the coronal base and $r_{\rm A}$.

The Canada-France Imaging Survey (CFIS) will consist of deep, high-resolution r-band imaging over ~5000 square degrees of the sky, representing a first-rate opportunity to identify recently-merged galaxies. Due to the large number of galaxies in CFIS, we investigate the use of a convolutional neural network (CNN) for automated merger classification. Training samples of post-merger and isolated galaxy images are generated from the IllustrisTNG simulation processed with the observational realism code RealSim. The CNN's overall classification accuracy is 88 percent, remaining stable over a wide range of intrinsic and environmental parameters. We generate a mock galaxy survey from IllustrisTNG in order to explore the expected purity of post-merger samples identified by the CNN. Despite the CNN's good performance in training, the intrinsic rarity of post-mergers leads to a sample that is only ~6 percent pure when the default decision threshold is used. We investigate trade-offs in purity and completeness with a variable decision threshold and find that we recover the statistical distribution of merger-induced star formation rate enhancements. Finally, the performance of the CNN is compared with both traditional automated methods and human classifiers. The CNN is shown to outperform Gini-M20 and asymmetry methods by an order of magnitude in post-merger sample purity on the mock survey data. Although the CNN outperforms the human classifiers on sample completeness, the purity of the post-merger sample identified by humans is frequently higher, indicating that a hybrid approach to classifications may be an effective solution to merger classifications in large surveys.

We present newly-calibrated period-$\phi_{31}$-[Fe/H] relations for fundamental mode RR Lyrae stars in the optical and, for the first time, mid-infrared. This work's calibration dataset provides the largest and most comprehensive span of parameter space to date with homogeneous metallicities from $-3<\textrm{[Fe/H]}<0.4$ and accurate Fourier parameters derived from 1980 ASAS-SN ($V$-band) and 1083 WISE (NEOWISE extension, $W1$ and $W2$ bands) RR Lyrae stars with well-sampled light curves. We compare our optical period-$\phi_{31}$-[Fe/H] with those available in the literature and demonstrate that our relation minimizes systematic trends in the lower and higher metallicity range. Moreover, a direct comparison shows that our optical photometric metallicities are consistent with both those from high-resolution spectroscopy and globular clusters, supporting the good performance of our relation. We found an intrinsic scatter in the photometric metallicities (0.41 dex in the $V$-band and 0.50 dex in the infrared) by utilizing large calibration datasets covering a broad metallicity range. This scatter becomes smaller when optical and infrared bands are used together (0.37 dex). Overall, the relations derived in this work have many potential applications, including large-area photometric surveys with JWST in the infrared and LSST in the optical.

We report the time variability of the late-time radio emission in a Type-I superluminous supernova (SLSN), PTF10hgi, at z = 0.0987. The Karl G. Jansky Very Large Array 3 GHz observations at 8.6 and 10 years after the explosion both detected radio emission with a ~40% decrease in flux density in the second epoch. This is the first report of a significant variability of the late-time radio light curve in a SLSN. Through combination with previous measurements in two other epochs, we constrained both the rise and decay phases of the radio light curve over three years, peaking at approximately 8-9 years after the explosion with a peak luminosity of L(3GHz) = 2 x 10^21 W/Hz. Possible scenarios for the origin of the variability are an active galactic nucleus (AGN) in the host galaxy, an afterglow caused by the interaction between an off-axis jet and circumstellar medium, and a wind nebula powered by a newly-born magnetar. Comparisons with models show that the radio light curve can be reproduced by both the afterglow model and magnetar wind nebula model. Considering the flat radio spectrum at 1-15 GHz and an upper limit at 0.6 GHz obtained in previous studies, plausible scenarios are a low-luminosity flat-spectrum AGN or a magnetar wind nebula with a shallow injection spectral index.

A starshade suppresses starlight by a factor of 1E11 in the image plane of a telescope, which is crucial for directly imaging Earth-like exoplanets. The state of the art in high contrast post-processing and signal detection methods were developed specifically for images taken with an internal coronagraph system and focus on the removal of quasi-static speckles. These methods are less useful for starshade images where such speckles are not present. This paper is dedicated to investigating signal processing methods tailored to work efficiently on starshade images. We describe a signal detection method, the generalized likelihood ratio test (GLRT), for starshade missions and look into three important problems. First, even with the light suppression provided by the starshade, rocky exoplanets are still difficult to detect in reflected light due to their absolute faintness. GLRT can successfully flag these dim planets. Moreover, GLRT provides estimates of the planets' positions and intensities and the theoretical false alarm rate of the detection. Second, small starshade shape errors, such as a truncated petal tip, can cause artifacts that are hard to distinguish from real planet signals; the detection method can help distinguish planet signals from such artifacts. The third direct imaging problem is that exozodiacal dust degrades detection performance. We develop an iterative generalized likelihood ratio test to mitigate the effect of dust on the image. In addition, we provide guidance on how to choose the number of photon counting images to combine into one co-added image before doing detection, which will help utilize the observation time efficiently. All the methods are demonstrated on realistic simulated images.

Nuclear isomers are populated in the rapid neutron capture process (r process) of nucleosynthesis. The r process may cover a wide range of temperatures, potentially starting from several tens of GK (several MeV) and then cooling as material is ejected from the event. As the r-process environment cools, isomers can freeze out of thermal equilibrium or be directly populated as astrophysical isomers (astromers). Two key behaviors of astromers -- ground state<->isomer transition rates and thermalization temperatures -- are determined by direct transition rates between pairs of nuclear states. We perform a sensitivity study to constrain the effects of unknown transitions on astromer behavior. We also introduce a categorization of astromers that describes their potential effects in hot environments. We provide a table of neutron-rich isomers that includes the astromer type, thermalization temperature, and key unmeasured transition rates.

Heating of neutral gas by energetic sources is crucial for the prediction of the 21 cm signal during the epoch of reionization (EoR). To investigate differences induced on statistics of the 21 cm signal by various source types, we use five radiative transfer simulations which have the same stellar UV emission model and varying combinations of more energetic sources, such as X-ray binaries (XRBs), accreting nuclear black holes (BHs) and hot interstellar medium emission (ISM). We find that the efficient heating from the ISM increases the average global 21~cm signal, while reducing its fluctuations and thus power spectrum. A clear impact is also observed in the bispectrum in terms of scale and timing of the transition between a positive and a negative value. The impact of XRBs is similar to that of the ISM, although it is delayed in time and reduced in intensity because of the less efficient heating. Due to the paucity of nuclear BHs, the behaviour of the 21~cm statistics in their presence is very similar to that of a case when only stars are considered, with the exception of the latest stages of reionization, when the effect of BHs is clearly visible. We find that differences between the source scenarios investigated here are larger than the instrumental noise of SKA1-low at $z \gtrsim 7-8$, suggesting that in the future it might be possible to constrain the spectral energy distribution of the sources contributing to the reionization process.

We present results from a spectroscopically blind search for HI 21-cm and OH 18-cm absorption lines at 2 < z < 5 towards 88 AGNs with 1.4GHz spectral luminosity of 27 - 29.3 W/Hz. One 21-cm absorber, which is associated with M1540-1453 at $z_{abs}$ = 2.1139, is detected. This is only the fourth known associated absorption at $z>2$. The detection rate ($1.6^{+3.8}_{-1.4}$%) suggests low covering factor of cold neutral medium (CNM; T~100K) associated with these powerful AGNs. The intervening HI 21-cm and OH 18-cm absorption searches, with a sensitivity to detect CNM in damped Ly-alpha systems (DLAs), have comoving absorption path lengths of $\Delta$X = 130.1 and 167.7, respectively. Using these we estimate the number of HI and OH absorber per unit comoving path lengths to be <0.014 and <0.011, respectively. The former is at least 4.5 times lower than that of DLAs and consistent with the CNM cross-section estimated using H$_2$ and CI absorbers at z>2. The AGNs in our sample selected using mid-infrared colors are optically fainter compared to the optical- and radio-selected quasars used to search for DLAs. In our optical spectra obtained using SALT and NOT, we detect 5 DLAs (redshift path ~9.3) and 2 proximate DLAs (within 3000 km/s of the AGN redshift). This is slightly excessive compared to the statistics based on optically selected quasars. The non-detection of HI 21-cm absorption from these DLAs suggests small CNM covering fraction around galaxies at z>2.

We numerically investigate the dynamics of a supernova fallback accretion confronting with a relativistic wind from a newborn neutron star (NS). The time evolution of the accretion shock in the radial direction is basically characterized by the encounter radius of the flows $r_\mathrm{enc}$ and a dimensionless parameter $\zeta \equiv L/\dot M_\mathrm{fb}c^2$, where $L$ is the NS wind luminosity and $\dot M_\mathrm{fb}$ is the fallback mass accretion rate. We find that the critical condition for the fallback matter to reach the near NS surface can be simply described as $\zeta < \zeta_\mathrm{min} \approx GM_*/c^2r_\mathrm{enc}$ independent of the wind Lorentz factor, where $M_*$ is the NS mass. By combining the condition for the fallback matter to bury the surface magnetic field under the NS crust, we discuss the possibility that the trifurcation of NSs into rotation-powered pulsars, central compact objects (CCOs), and magnetars can be induced by supernova fallback.

Eric H. Strach (1914-2011) studied medicine at University of Prague and graduated in 1938. Strach dedicated a great part of his life to astronomy becoming a consistent and meticulous observer. He joined the Liverpool Astronomical Association and British Astronomical Association during the 1960s and obtained two recognitions as proof of his great work in solar physics: the BAA's Merlin Medal and Gift in 1999 and Walter Goodacre Medal and Gift, ten years later. Strach recorded four decades (1969-2008) of systematic solar records in his observation notebooks although he started his observations from the late 1950s. In this work, we document the valuable effort made by Strach in getting four decades of solar records and the importance of this kind of long observation series for studies of space weather and climate. We present the sunspot group number series according to Strach's data and a long observation series of prominences recorded by Strach.

A revision is presented of the sunspot observations made by Charles Malapert from 1618 to 1626, studying several documentary sources that include those observations. The revised accounting of the group numbers recorded by Malapert for that period shows new information unavailable in the current sunspot group database. The average solar activity level calculated from these revised records of Malapert is by almost one third greater than that calculated from his records included in the current group database. Comparison of the sunspot observations made by Malapert and by other astronomers of that time with regard to the number of recorded groups and sunspot positions on the solar disk shows good agreement. Malapert reported that he only recorded one sunspot group in each sunspot drawing presented in Austriaca Sidera Heliocyclia (the documentary source which includes most of the sunspot records made by Malapert), although he sometimes observed several groups. Therefore, the sunspot counts obtained in this present work on Malapert's sunspot observations represents the lower limit of the solar activity level corresponding to those records.

An analysis of the sunspot observations made by Hevelius during 1642-1645 is presented. These records are the only systematic sunspot observations just before the Maunder Minimum. We have studied different phenomena meticulously recorded by Hevelius after translating the original Latin texts. We re-evaluate the observations of sunspot groups by Hevelius during this period and obtain an average value 7% greater than that calculated from his observations given in the current group database. Furthermore, the average of the active day fraction obtained in this work from Hevelius' records previous to the Maunder Minimum is significantly greater than the solar activity level obtained from Hevelius' sunspot observations made during the Maunder Minimum (70% vs. 30%). We also present the butterfly diagram obtained from the sunspot positions recorded by Hevelius for the period 1642-1645. It can be seen that no hemispheric asymmetry exists during this interval, in contrast with the Maunder Minimum. Hevelius noted a ~3-month period that appeared to lack sunspots in early 1645 that gave the first hint of the impending Maunder Minimum. Recent studies claim that the Maunder Minimum was not a grand minimum period speculating that astronomers of that time, due to the Aristotelian ideas, did not record all sunspots that they observed, producing thus an underestimation of the solar activity level. However, we show the good quality of the sunspot records made by Hevelius indicates that his reports of sunspots were true to the observations.

Nova outbursts play an important role in the chemical evolution of galaxies, especially they are the main source of synthetic $^{13}\rm C$, $^{15}\rm N$, $^{17}\rm O$ and some radioactive isotopes like $^{22}\rm Na$ and $^{26}\rm Al$. The enrichment of He in nova ejecta indicates that the accreted material may mix with the He-shell (He-mixing). The purpose of this work is to investigate how the He-mixing affects the nova outbursts in a systematic way. We evolved a series of accreting WD models, and found that the mass fraction of H and He in nova ejecta can be influenced by different He-mixing fractions significantly. We also found that both the nova cycle duration and ejected mass increase with the He-mixing fractions. Meanwhile, the nuclear energy production from $p$-$p$ chains decreases with the He-mixing fraction during the nova outbursts, whereas the CNO-cycle increases. The present work can reproduce the chemical abundances in the ejecta of some novae, such as GQ Mus, ASASSN-18fv, HR Del, T Aur and V443 Sct. This implies that the He-mixing process cannot be neglected when studying nova outbursts. This study also develops a He-mixing meter (i.e. $\rm He/H$) that can be used to estimate the He-mixing fraction in classical nova systems.

Context: The long-term carbonate-silicate cycle plays an important role in the evolution of Earth's climate and, therefore, may also be an important mechanism in the evolution of the climates of Earth-like exoplanets. Aims: We investigate the effects of radiogenic mantle heating, core size, and planetary mass on the evolution of the atmospheric partial $CO_2$ pressure, and the ability of a long-term carbon cycle driven by plate tectonics to control the atmospheric $CO_2$ pressure. Methods: We developed a box-model which connects carbon cycling to parametrized mantle convection. The carbon cycle was coupled to the thermal evolution via the plate speed, which depends on the global Rayleigh number. Results: We find decreasing atmospheric $CO_2$ pressure with time, up to an order of magnitude over 10 Gyr. Higher abundances of radioactive isotopes result in higher $CO_2$ pressures. We find a decreasing Rayleigh number and plate speed toward planets with larger core mass fractions $f_c$, which leads to lower atmospheric $CO_2$ pressure. More massive planets may favor the development of more $CO_2$ rich atmospheres due to hotter interiors. Conclusions: The dependence of plate tectonics on mantle cooling has a significant effect on the long-term evolution of the atmospheric $CO_2$ pressure. Carbon cycling mediated by plate tectonics is efficient in regulating planetary climates for a wide range of mantle radioactive isotope abundances, planet masses and core sizes. More efficient carbon cycling on planets with a high mantle abundance of thorium or uranium highlights the importance of mapping the abundances of these elements in host stars of potentially habitable exoplanets. Inefficient carbon recycling on planets with a large core mass fraction ($f_c\gtrsim 0.8$) emphasizes the importance of precise mass-radius measurements of Earth-sized exoplanets.

Aims: Interstellar molecules form early in the evolutionary sequence of interstellar material that eventually forms stars and planets. To understand this evolutionary sequence, it is important to characterize the chemical composition of its first steps. Methods: In this paper, we present the result of a 2 and 3 mm survey of five cold clumps identified by the Planck mission. We carried out a radiative transfer analysis on the detected lines in order to put some constraints on the physical conditions within the cores and on the molecular column densities. We also performed chemical models to reproduce the observed abundances in each source using the gas-grain model Nautilus. Results: Twelve molecules were detected: H2CO, CS, SO, NO, HNO, HCO+, HCN, HNC, CN, CCH, CH3OH, and CO. Here, CCH is the only carbon chain we detected in two sources. Radiative transfer analyses of HCN, SO, CS, and CO were performed to constrain the physical conditions of each cloud with limited success. The sources have a density larger than $10^4$ cm$^{-3}$ and a temperature lower than 15 K. The derived species column densities are not very sensitive to the uncertainties in the physical conditions, within a factor of 2. The different sources seem to present significant chemical differences with species abundances spreading over one order of magnitude. The chemical composition of these clumps is poorer than the one of Taurus Molecular Cloud 1 Cyanopolyyne Peak (TMC-1 CP) cold core. Our chemical model reproduces the observational abundances and upper limits for 79 to 83\% of the species in our sources. The "best" times for our sources seem to be smaller than those of TMC-1, indicating that our sources may be less evolved and explaining the smaller abundances and the numerous non-detections. Also, CS and HCN are always overestimated by our models.

The Fermi-LAT observational data of the diffuse $\gamma$ ray emission from the Large Magellanic Cloud (LMC), were examined to test for the existence of underlying long range correlations. A statistical test applied to the data indicated that the probability that data are random is $\sim 10^{-99}$. Thus we proceeded and have used the counts-number data to compute 2D spatial autocorrelation, power spectrum, and structure function. The most important result of the present study is a clear indication for large scale spatial underlying correlations. This is evident in {\bf all} the functions mentioned above. The 2D power spectrum has a logarithmic slope of -3 on large spatial scales and a logarithmic slope of -4 on small spatial scales. The structure function has logarithmic slopes equaling 1 and 2 for the large and small scales respectively. The logarithmic slopes of the structure function and the power spectrum are consistent. A plausible interpretation of these results is the existence of a large scale {\it compressible turbulence} with a 3D logarithmic slope of -4 extending over the entire extent of the LMC. This may reflect the fact that the $\gamma$ Ray emission is in star forming regions, where jets and shocks are abundant. Both the power spectrum and structure function exhibit steeper logarithmic slopes for smaller spatial scales. This is interpreted as an indication that the turbulent region has an effective depth of about 1.5 kpc.

We report on a multi-wavelength analysis of the 26 January 2014 solar eruption involving a coronal mass ejection (CME) and a Type-II radio burst, performed by combining data from various space-and ground-based instruments. An increasing standoff distance with height shows the presence of a strong shock, which further manifests itself in the continuation of the metric Type-II burst into the decameter-hectometric (DH) domain. A plot of speed versus position angle (PA) shows different points on the CME leading edge travelled with different speeds. From the starting frequency of the Type-II burst and white-light data, we find that the shock signature producing the Type-II burst might be coming from the flanks of the CME. Measuring the speeds of the CME flanks, we find the southern flank to be at a higher speed than the northern flank; further the radio contours from Type-II imaging data showed that the burst source was coming from the southern flank of the CME. From the standoff distance at the CME nose, we find that the local Alfven speed is close to the white-light shock speed, thus causing the Mach number to be small there. Also, the presence of a streamer near the southern flank appears to have provided additional favorable conditions for the generation of shock-associated radio emission. These results provide conclusive evidence that the Type-II emission could originate from the flanks of the CME, which in our study is from the the southern flank of the CME.

Optical observations of the quasar J0015+1842 at z=2.6 have recently revealed the presence of multi-phase outflowing winds intercepted by the line of sight to the active nucleus. Here, we present 3-mm observations of this system with the NOrthern Extended Millimeter Array (NOEMA). Our data reveals molecular gas, traced via a Gaussian CO(3-2) line, with a remarkably large velocity dispersion (FWHM=1014+/-120 km/s). The integrated line flux indicates a total molecular mass in the range M_H2 ~ (3.4-17)x10^10 Msun, with the spread dominated by the assumed CO-to-H2 conversion factor. Assuming the 3-mm continuum emission is thermal, we derive a dust mass of Mdust ~ 2x10^9 Msun. J0015+1842 is located in the molecular gas-rich region in the IR vs CO line luminosity diagram, in-between the main locus of main-sequence and sub-millimeter galaxies and that of most other AGNs targeted so far for CO measurements. While the large velocity dispersion of the CO line suggests a post merger system, J0015+1842 is observed to be a regular, only very moderately dust-reddened (Av~0.3-0.4) type-I quasar from its UV-optical spectrum, from which we infer a mass of the super-massive black hole be around M_BH~6x10^8 Msun. We suggest that J0015+1842 is observed at a galaxy evolutionary stage where a massive merger has brought significant amounts of gas towards an actively accreting super-massive black hole (quasar) that already cleared the way in the circum-nuclear region towards the observer through powerful winds, while the host still contains a large amount of dust and molecular gas with high velocity dispersion. Observations of a sample of similar systems should help determining better the respective importance of evolution and orientation in the appearance of quasars and their host galaxies and have the potential to investigate early feedback and star-formation processes in galaxies in their quasar phases.

Active regions often show S-shaped structures in the corona called sigmoids. These are highly sheared and twisted loops formed along the polarity inversion line. They are considered to be one of the best pre-eruption signatures for CMEs. Here, we investigate the thermodynamic evolution of an on-disk sigmoid observed during December 24-28, 2015. For this purpose, we have employed Emission Measure (EM) and filter-ratio techniques on the observations recorded by the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) and X-ray Telescope (XRT) onboard Hinode. The EM analysis showed multi-thermal plasma along the sigmoid and provided a peak temperature of 10-12.5 MK for all observed flares. The sigmoidal structure showed emission from Fe XVIII (93.93 {\AA}) and Fe XXI 128.75 {\AA}) lines in the AIA 94 and 131 {\AA} channels, respectively. Our results show that the hot plasma is often confined to very hot strands. The temperature obtained from the EM analysis was found to be in good agreement with that obtained using the XRT, AIA, and GOES filter-ratio methods. These results provide important constraints for the thermodynamic modeling of sigmoidal structures in the core of active regions. Moreover, this study also benchmarks different techniques available for temperature estimation in solar coronal structures.

To investigate the feasibility of ancillary target observations with ESA's Ariel mission, we compiled a list of potentially interesting young stars: FUors, systems harbouring extreme debris discs and a larger sample of young stellar objects showing strong near/mid-infrared excess. These objects can be observed as additional targets in the waiting times between the scheduled exoplanet transit and occultation observations. After analyzing the schedule for Ariel an algorithm was constructed to find the optimal target to be observed in each gap. The selection was mainly based on the slew and stabilization time needed to observe the selected YSO, but it also incorporated the scientific importance of the targets and whether they have already been sufficiently measured. After acquiring an adequately large sample of simulation data, it was concluded that approximately 99.2% of the available -- at least one hour long -- gaps could be used effectively. With an average slewing and stabilization time of about 16.7 minutes between scheduled exoplanet transits and ancillary targets, this corresponds to an additional $2881 \pm 56$ hours of active data gathering. When this additional time is used to observe our selected 200 ancillary targets, a typical signal-to-noise ratio of $\sim$10$^4$ can be achieved along the whole spectral window covered by Ariel.

Astrochemistry lies at the nexus of astronomy, chemistry, and molecular physics. On the basis of precise laboratory data, a rich collection of more than 200 familiar and exotic molecules have been identified in the interstellar medium, the vast majority by their unique rotational fingerprint. Despite this large body of work, there is scant evidence in the radio band for the basic building blocks of chemistry on earth -- five and six-membered rings -- despite long standing and sustained efforts during the past 50 years. In contrast, a peculiar structural motif, highly unsaturated carbon in a chain-like arrangement, is instead quite common in space. The recent astronomical detection of cyanobenzene, the simplest aromatic nitrile, in the dark molecular cloud TMC-1, and soon afterwards in additional pre-stellar, and possibly protostellar sources, establishes that aromatic chemistry is likely widespread in the earliest stages of star formation. The subsequent discovery of cyanocyclopentadienes and even cyanonapthlenes in TMC-1 provides further evidence that organic molecules of considerable complexity are readily synthesized in regions with high visual extinction but where the low temperature and pressure are remarkably low. This review focuses on laboratory efforts now underway to understand the rich transition region between linear and planar carbon structures using microwave spectroscopy. We present key features, advantages, and disadvantages of current detection methods, a discussion of the types of molecules found in space and in the laboratory, and approaches under development to identify entirely new species in complex mixtures. Studies focusing on the cyanation of hydrocarbons and the formation of benzene from acyclic precursors are highlighted, as is the role that isotopic studies might play in elucidating the chemical pathways to ring formation.

We use Gaia DR2 data to survey the classic Monoceros OB1 region and look for the existence of a dispersed young population, co-moving with the cloud complex. An analysis of the distribution of proper motions reveals a 20-30 Myr association of young stars, about 300-400 pc away from the far side of the Mon OB1 complex, along the same general line-of-sight. We characterize the new association, Monoceros OB4, and estimate it contains between 1400 and 2500 stars, assuming a standard IMF, putting it on par in size with NGC\,2264. We find from the internal proper motions that Mon OB4 is unbound and expanding. Our results seem to unveil a larger and more complex Monoceros star formation region, suggesting an elongated arrangement that seems to be at least 300 x 60 pc.

The Vela supernova remnant (SNR) shows several ejecta fragments protruding beyond the forward shock (shrapnel). Recent studies have revealed high Si abundance in two shrapnel (A and G), located in opposite directions with respect to the SNR center. This suggests the possible existence of a Si-rich jet-counterjet structure. We analyzed an XMM-Newton observation of a bright clump, behind shrapnel G, which lies along the direction connecting A and G. The aim is to study the physical and chemical properties of this clump to ascertain whether it is part of this putative jet-like structure. We produced background-corrected and adaptively-smoothed count-rate images and median photon energy maps, and performed a spatially resolved spectral analysis. We identified two structures with different physical properties. The first one is remarkably elongated along the direction connecting A and G. Its X-ray spectrum is much softer than that of the other two shrapnel, to the point of hindering the determination of the Si abundance, however its physical and chemical properties are consistent with those of shrapnel A and G. The second structure, running along the southeast-northwest direction, has a higher temperature and appears like a thin filament. By analyzing the ROSAT data, we have found that this filament is part of a very large and coherent structure that we identified in the western rim of the shell. We obtained a thorough description of the tail of Shrapnel G. In addition we discovered a coherent and very extended feature that we interpret as a signature of an earlier interaction of the remnant with the stellar wind of its progenitor star. The peculiar Ne/O ratio we found in the wind residual may be suggestive of a Wolf-Rayet progenitor for Vela SNR, though further analysis is required to address this point.

Context. The high number of super-Earth and Earth-like planets in the habitable zone (HZ) detected around M-dwarf stars in the last years has revealed these stellar objects to be the key for planetary radial velocity (RV) searches. Aims. Using the HARPS-N spectrograph within The HArps-n red Dwarf Exoplanet Survey (HADES) we reach the precision needed to detect small planets with a few Earth masses using the RV technique. Methods. We obtained 138 HARPS-N RV measurements between 2013 May and 2020 September of GJ 720 A, classified as an M0.5V star located at a distance of 15.56 pc. To characterize the stellar variability and to discern the periodic variation due to the Keplerian signals from those related to stellar activity, the HARPS-N spectroscopic activity indicators and the simultaneous photometric observations were analyzed. The combined analysis of HARPS-N RVs and activity indicators let us to address the nature of the periodic signals. The final model and the orbital planetary parameters were obtained by fitting simultaneously the stellar variability and the Keplerian signal using a Gaussian process regression and following a Bayesian criterion. Results. The HARPS-N RV periodic signals around 40 d and 100 d have counterparts at the same frequencies in HARPS-N activity indicators and photometric light curves. Then we attribute these periodicities to stellar activity the former period being likely associated with the stellar rotation. GJ 720 A shows the most significant signal at 19.466$\pm$0.005 d with no counterparts in any stellar activity indices. We hence ascribe this RV signal, having a semiamplitude of 4.72$\pm$0.27 m/s , to the presence of a sub-Neptune mass planet. The planet GJ 720 Ab has a minimum mass of 13.64$\pm$0.79 M$_{\oplus}$, it is in circular orbit at 0.119$\pm$0.002 AU from its parent star, and lies inside the inner boundary of the HZ around its parent star.

In this work, we discuss the application of convolutional neural networks (CNNs) as a tool to advantageously initialize Stokes profile inversions. To demonstrate the usefulness of CNNs, we concentrate in this paper on the inversion of LTE Stokes profiles. We use observations taken with the spectropolarimeter onboard the Hinode spacecraft as a test benchmark. First, we carefully analyze the data with the SIR inversion code using a given initial atmospheric model. The code provides a set of atmospheric models that reproduce the observations. These models are then used to train a CNN. Afterwards, the same data are again inverted with SIR but using the trained CNN to provide the initial guess atmospheric models for SIR. The CNNs allow us to significantly reduce the number of inversion cycles when used to compute initial guess model atmospheres, decreasing the computational time for LTE inversions by a factor of two to four. CNN's alone are much faster than assisted inversions, but the latter are more robust and accurate. The advantages and limitations of machine learning techniques for estimating optimum initial atmospheric models for spectral line inversions are discussed. Finally, we describe a python wrapper for the SIR and DeSIRe codes that allows for the easy setup of parallel inversions. The assisted inversions can speed up the inversion process, but the efficiency and accuracy of the inversion results depend strongly on the solar scene and the data used for the CNN training. This method (assisted inversions) will not obviate the need for analyzing individual events with the utmost care but will provide solar scientists with a much better opportunity to sample large amounts of inverted data, which will undoubtedly broaden the physical discovery space.

W 601 (NGC 6611 601) is one of the handful of known magnetic Herbig Ae/Be stars. We report the analysis of a large dataset of high-resolution spectropolarimetry. The star is a previously unreported spectroscopic binary, consisting of 2 B2 stars with a mass ratio of 1.8, masses of 12 M$_\odot$ and 6.2 $M_\odot$, in an eccentric 110-day orbit. The magnetic field belongs to the secondary, W 601 B. The H$\alpha$ emission is consistent with an origin in W 601 B's centrifugal magnetosphere; the star is therefore not a classical Herbig Be star in the sense that its emission is not formed in an accretion disk. However, the low value of $\log{g} = 3.8$ determined via spectroscopic analysis, and the star's membership in the young NGC 6611 cluster, are most consistent with it being on the pre-main sequence. The rotational period inferred from the variability of the H$\alpha$ line and the longitudinal magnetic field $\langle B_z \rangle$ is 1.13 d. Modelling of Stokes $V$ and $\langle B_z \rangle$ indicates a surface dipolar magnetic field $B_{\rm d}$ between 6 and $11$ kG. With its strong emission, rapid rotation, and strong surface magnetic field, W 601 B is likely a precursor to H$\alpha$-bright magnetic B-type stars such as $\sigma$ Ori E. By contrast, the primary is an apparently non-magnetic ($B_{\rm d} < 300$ G) pre-main sequence early B-type star. In accordance with expectations from magnetic braking, the non-magnetic primary is apparently more rapidly rotating than the magnetic star.

NASA's Stardust mission utilized a sample collector composed of aerogel and aluminum foil to return cometary and interstellar particles to Earth. Analysis of the aluminum foil begins with locating craters produced by hypervelocity impacts of cometary and interstellar dust. Interstellar dust craters are typically less than one micrometer in size and are sparsely distributed, making them difficult to find. In this paper, we describe a convolutional neural network based on the VGG16 architecture that achieves high specificity and sensitivity in locating impact craters in the Stardust interstellar collector foils. We evaluate its implications for current and future analyses of Stardust samples.

Most models of the low frequency quasi periodic oscillations (QPOs) in low-mass X-ray binaries (LMXBs) explain the dynamical properties of those QPOs. On the other hand, in recent years reverberation models that assume a lamp-post geometry have been successfull in explaining the energy-dependent time lags of the broad-band noise component in stellar mass black-holes and active galactic nuclei. We have recently shown that Comptonisation can explain the spectral-timing properties of the kilo-hertz (kHz) QPOs observed in neutron star (NS) LMXBs. It is therefore worth exploring whether the same family of models would be as successful in explaining the low-frequency QPOs. In this work, we use a Comptonisation model to study the frequency dependence of the phase lags of the type-C QPO in the BH LMXB GRS 1915+105. The phase lags of the QPO in GRS 1915+105 make a transition from hard to soft at a QPO frequency of around 1.8 Hz. Our model shows that at high QPO frequencies a large corona of ~ 100-150 R_g covers most of the accretion disc and makes it 100\% feedback dominated, thus producing soft lags. As the observed QPO frequency decreases, the corona gradually shrinks down to around 3-17 R_g, and at 1.8 Hz feedback onto the disc becomes inefficient leading to hard lags. We discuss how changes in the accretion geometry affect the timing properties of the type-C QPO.

Classification of galactic morphologies is a crucial task in galactic astronomy, and identifying fine structures of galaxies (e.g., spiral arms, bars, and clumps) is an essential ingredient in such a classification task. However, seeing effects can cause images we obtain to appear blurry, making it difficult for astronomers to derive galaxies' physical properties and, in particular, distant galaxies. Here, we present a method that converts blurred images obtained by the ground-based Subaru Telescope into quasi Hubble Space Telescope (HST) images via machine learning. Using an existing deep learning method called generative adversarial networks (GANs), we can eliminate seeing effects, effectively resulting in an image similar to an image taken by the HST. Using multiple Subaru telescope image and HST telescope image pairs, we demonstrate that our model can augment fine structures present in the blurred images in aid for better and more precise galactic classification. Using our first of its kind machine learning-based deblurring technique on space images, we can obtain up to 18% improvement in terms of CW-SSIM (Complex Wavelet Structural Similarity Index) score when comparing the Subaru-HST pair versus SeeingGAN-HST pair. With this model, we can generate HST-like images from relatively less capable telescopes, making space exploration more accessible to the broader astronomy community. Furthermore, this model can be used not only in professional morphological classification studies of galaxies but in all citizen science for galaxy classifications.

GG Carinae is a binary whose primary component is a \sgBe. Using photometric data from TESS, ASAS, OMC, and ASAS-SN, and spectroscopic data from the Global Jet Watch to study visible He\,I, Fe\,II and Si\,II emission lines, we investigate the short-period variations which are exhibited in GG Car. We find a hitherto neglected periodicity of $1.583156\pm0.0002$\,days that is present in both its photometry and the radial velocities of its emission lines, alongside variability at the well-established $\sim$31-day orbital period. We find that the amplitudes of the shorter-period variations in both photometry and some of the emission lines are modulated by the orbital phase of the binary, such that the short-period variations have largest amplitudes when the binary is at periastron. There are no significant changes in the phases of the short-period variations over the orbital period. We investigate potential causes of the 1.583-day variability, and find that the observed period agrees well with the expected period of the $l=2$ f-mode of the primary given its mass and radius. We propose that the primary is periodically pulled out of hydrostatic equilibrium by the quadrupolar tidal forces when the components are near periastron in the binary's eccentric orbit ($e=0.5$) and the primary almost fills its Roche lobe. This causes an oscillation at the $l=2$ f-mode frequency which is damped as the distance between the components increases.

In pre-recombination resolutions of the Hubble tension, such as early dark energy, new physics before recombination shifts the values of relevant cosmological parameters so that the models can fit with CMB and BAO observations as well as $\Lambda$CDM does. In this paper, we clarify how the parameter shifts are related with $\delta H_0$, particularly we find the shift of primordial scalar spectral index scales as ${\delta n_s}\simeq 0.4{\delta H_0\over H_0}$ by performing the Monte Carlo Markov chain analysis with Planck2018+BAO+Pantheon+R19+Keck Array/BICEP dataset. A novel point of our result is that if the current $H_0$ measured locally is correct, complete resolution of the Hubble tension seems to be pointing to a scale invariant Harrison-Zeldovich spectrum, i.e. $n_s= 1$ for $H_0\sim 73$km/s/Mpc.

This paper proposes a new approach for the selection of low-energy neutrino bursts, such as the ones detected after a supernova. It exploits the temporal structure of the expected signal with respect to the more diffuse background by defining a "Real-time Test Statistic" (RTS) that would allow identifying very weak signals, hard to select using standard clustering methods. For a given background rate, the new method (RTS$^2$: RTS for Supernovae) increases signal efficiency while keeping the same false alarm rate for Poisson-distributed background. By adding a spatial penalty term to the definition of RTS, one can also reject spatially-correlated backgrounds such as the ones due to spallation events. Furthermore, the method is easy to implement in a real-time monitoring system as RTS can be computed recursively for successive events, and it can be easily adapted for detectors of all scales that may want to send prompt alerts e.g. through SNEWS 2.0 network.

We present a new analysis of the light curve of the young planet-hosting star TOI 451 in the light of new observations from TESS Cycle 3. Our joint analysis of the transits of all three planets, using all available TESS data, results in an improved ephemeris for TOI 451 b and TOI 451 c, which will help to plan follow-up observations. The updated mid-transit times are $\textrm{BJD}-2,457\,000=$ $1410.9896_{ - 0.0029 }^{ + 0.0032 }$, $1411.7982_{-0.0020}^{+0.0022}$, and $1416.63407_{-0.00100}^{+0.00096}$ for TOI 451 b, c, and d, respectively, and the periods are $1.8587028_{-10e-06}^{+08e-06}$, $9.192453_{-3.3e-05}^{+4.1e-05}$, and $16.364932_{-3.5e-05}^{+3.6e-05 }$ days. We also model the out-of-transit light curve using a Gaussian Process with a quasi-periodic kernel and infer a change in the properties of the active regions on the surface of TOI 451 between TESS Cycles 1 and 3.

The aim of this work is to determine abundances of neutron-capture elements for thin- and thick-disc F, G, and K stars in several sky fields near the north ecliptic pole and to compare the results with the Galactic chemical evolution models, to explore elemental gradients according to stellar ages, mean galactocentric distances, and maximum heights above the Galactic plane. The observational data were obtained with the 1.65m telescope at the Moletai Astronomical Observatory and a fibre-fed high-resolution spectrograph. Elemental abundances were determined using a differential spectrum synthesis with the MARCS stellar model atmospheres and accounting for the hyperfine-structure effects. We determined abundances of Sr, Y, Zr, Ba, La, Ce, Pr, Nd, Sm, and Eu for 424 thin- and 82 thick-disc stars. The sample of thick-disc stars shows a clearly visible decrease in [Eu/Mg] with increasing [Fe/H] compared to the thin-disc stars, bringing more evidence of a different chemical evolution in these two Galactic components. Abundance correlation with age slopes for the investigated thin-disc stars are slightly negative for the majority of s-process dominated elements, while r-process dominated elements have positive correlations. Our sample of thin-disc stars with ages spanning from 0.1 to 9 Gyrs give the [Y/Mg]=0.022 ($\pm$0.015)-0.027 ($\pm$0.003)*age [Gyr] relation. For the thick-disc stars, when we also took data from other studies into account, we found that [Y/Mg] cannot serve as an age indicator. The radial [El/Fe] gradients in the thin disc are negligible for the s-process dominated elements and become positive for the r-process dominated elements. The vertical gradients are negative for the light s-process dominated elements and become positive for the r-process dominated elements. In the thick disc, the radial [El/Fe] slopes are negligible, and the vertical slopes are predominantly negative.

Accretion disks and coronae around massive black holes have been studied extensively, and they are known to be coupled. Over a period of 30 years, however, the X-ray (coronal) flux of Mrk 817 increased by a factor of 40 while its UV (disk) flux remained relatively steady. Recent high-cadence monitoring finds that the X-ray and UV continua in Mrk 817 are also decoupled on time scales of weeks and months. These findings could require mechanical beaming of the innermost accretion flow, and/or an absorber that shields the disk and/or broad line region (BLR) from the X-ray corona. Herein, we report on a 135 ks observation of Mrk 817 obtained with NuSTAR, complemented by simultaneous X-ray coverage via the Neil Gehrels Swift Observatory. The X-ray data strongly prefer a standard relativistic disk reflection model over plausible alternatives. Comparable fits with related models constrain the spin to lie in the range 0.5 < a < 1, and the viewing angle to lie between 10 deg. < theta < 22 deg. (including 1-sigma statistical errors and small systematic errors related to differences between the models). The spectra also reveal strong evidence of moderately ionized absorption, similar to but likely less extreme than obscuring events in NGC 5548 and NGC 3783. Archival Swift data suggest that the absorption may be variable. Particularly if the column density of this absorber is higher along the plane of the disk, it may intermittently mask or prevent coupling between the central engine, disk, and BLR in Mrk 817.

The first exploration of Pluto was motivated by (i) the many intriguing aspects of this body, its atmosphere, and its giant impact binary-planet formation; as well as (ii) the scientific desire to initiate the reconnaissance of the newly-discovered population of dwarf planets in the Kuiper Belt. That exploration took place in the form of a single spacecraft flyby that yielded an impressive array of exciting results that have transformed our understanding of this world and its satellites, and which opened our eyes to the exciting nature of the dwarf planet population of the Kuiper Belt. From Pluto's five-object satellite system, to its hydrocarbon haze-laden atmosphere, to its variegated distribution of surface volatiles, to its wide array of geologic expressions that include extensive glaciation and suspected cryovolcanoes, plus the tantalizing possibility of an interior ocean, the Pluto system has proven to be as complex as larger terrestrial bodies like Mars, and it begs for future exploration. Owing to Pluto's high obliquity (and consequently, current-epoch southern hemisphere polar winter darkness) and the single spacecraft nature of the New Horizons flyby, only about 40% of Pluto and its binary satellite, Charon, could be mapped at high pixel scales (better than 10 km/pix). Additionally, due to their distances from New Horizons at closest approach, none of Pluto's small moons could be studied at high resolution during the flyby. Furthermore, studies of the time variability of atmospheric, geologic, and surface-atmosphere interactions cannot be practically made by additional flybys, and they cannot be made from Earth-based observations. We find that these limitations, combined with Pluto's many important, open scientific questions, strongly motivate a Pluto System follow on orbiter mission.

In the standard approach to deriving inflationary predictions, we evolve a vacuum state in time according to the rules of a given model. Since the only observables are the future values of correlators and not their time evolution, this brings about a large degeneracy: a vast number of different models are mapped to the same minute number of observables. Furthermore, due to the lack of time-translation invariance, even tree-level calculations require an increasing number of nested integrals that quickly become intractable. Here we ask how much of the final observables can be "bootstrapped" directly from locality, unitarity and symmetries. To this end, we introduce two new bootstrap tools to efficiently compute cosmological correlators/wavefunctions. The first is a Manifestly Local Test (MLT) that any $n$-point (wave)function of massless scalars or gravitons must satisfy if it is to arise from a manifestly local theory. When combined with a sub-set of the recently proposed Bootstrap Rules, this allows us to compute explicitly all bispectra to all orders in derivatives for a single scalar. Since we don't invoke soft theorems, this can also be extended to multi-field inflation. The second is a partial energy recursion relation that allows us to compute exchange correlators. Combining a bespoke complex shift of the partial energies with Cauchy's integral theorem and the Cosmological Optical Theorem, we fix exchange correlators up to a boundary term. The latter can be determined up to contact interactions using unitarity and manifest locality. As an illustration, we use these tools to bootstrap scalar inflationary trispectra due to graviton exchange and inflaton self-interactions.

A distinct signature of compact extra dimensions would be a Kaluza-Klein tower of gravitational waves. Motivated by this prospect, we compute the corresponding spectrum on a warped toroidal background. We evaluate in particular the impact of the warp factor on the spectrum. To that end, we use the complete warp factor H of standard string compactifications, generated by D-branes and orientifolds, thus connecting to recent works on stringy de Sitter constructions. The problematic region close to an orientifold where H < 0 leads to unphysical tachyonic modes in the spectrum. We develop tools that overcome this difficulty and lead to a tachyon-free spectrum. We show, in particular, that the warp factor can lower the first Kaluza-Klein mass by at least 69%.

Inflationary perturbations are approximately Gaussian and deviations from Gaussianity are usually calculated using in-in perturbation theory. This method, however, fails for unlikely events on the tail of the probability distribution: in this regime non-Gaussianities are important and perturbation theory breaks down for $|\zeta| \gtrsim |f_{\rm \scriptscriptstyle NL}|^{-1}$. In this paper we show that this regime is amenable to a semiclassical treatment, $\hbar \to 0$. In this limit the wavefunction of the Universe can be calculated in saddle-point, corresponding to a resummation of all the tree-level Witten diagrams. The saddle can be found by solving numerically the classical (Euclidean) non-linear equations of motion, with prescribed boundary conditions. We apply these ideas to a model with an inflaton self-interaction $\propto \lambda \dot\zeta^4$. Numerical and analytical methods show that the tail of the probability distribution of $\zeta$ goes as $\exp(-\lambda^{-1/4}\zeta^{3/2})$, with a clear non-perturbative dependence on the coupling. Our results are relevant for the calculation of the abundance of primordial black holes.

The $^{23}$Na($\alpha,p$)$^{26}$Mg reaction has been identified as having a significant impact on the nucleosynthesis of several nuclei between Ne and Ti in type-Ia supernovae, and of $^{23}$Na and $^{26}$Al in massive stars. The reaction has been subjected to renewed experimental interest recently, motivated by high uncertainties in early experimental data and in the statistical Hauser-Feshbach models used in reaction rate compilations. Early experiments were affected by target deterioration issues and unquantifiable uncertainties. Three new independent measurements instead are utilizing inverse kinematics and Rutherford scattering monitoring to resolve this. In this work we present directly measured angular distributions of the emitted protons to eliminate a discrepancy in the assumptions made in the recent reaction rate measurements, which results in cross sections differing by a factor of 3. We derive a new combined experimental reaction rate for the $^{23}$Na($\alpha,p$)$^{26}$Mg reaction with a total uncertainty of 30% at relevant temperatures. Using our new $^{23}$Na($\alpha,p$)$^{26}$Mg rate, the $^{26}$Al and $^{23}$Na production uncertainty is reduced to within 8%. In comparison, using the factor of 10 uncertainty previously recommended by the rate compilation STARLIB, $^{26}$Al and $^{23}$Na production was changing by more than a factor of 2. In type-Ia supernova conditions, the impact on production of $^{23}$Na is constrained to within 15%.

We explore the benefits of adapted gauges to small mass ratio binary black hole evolutions in the moving puncture formulation. We find expressions that approximate the late time behavior of the lapse and shift, $(\alpha_0,\beta_0)$, and use them as initial values for their evolutions. We also use a position and black hole mass dependent damping term, $\eta[\vec{x}_1(t),\vec{x}_2(t),m_1,m_2]$, in the shift evolution, rather than a constant or conformal-factor dependent choice. We have found that this substantially reduces noise generation at the start of the numerical integration and keeps the numerical grid stable around both black holes, allowing for more accuracy with lower resolutions. We test our choices for this gauge in detail in a case study of a binary with a 7:1 mass ratio, and then use 15:1 and 32:1 binaries for a convergence study. Finally, we apply our new gauge to a 64:1 binary and a 128:1 binary to well cover the comparable and small mass ratio regimes.

The space-based gravitational wave (GW) detector LISA is expected to observe signals from a large population of compact object binaries, comprised predominantly of white dwarfs, in the Milky Way. Resolving individual sources from this population against its self-generated confusion noise poses a major data analysis problem. We present an iterative source estimation and subtraction method to address this problem based on the use of particle swarm optimization (PSO). In addition to PSO, a novel feature of the method is the cross-validation of sources estimated from the same data using different signal parameter search ranges. This is found to greatly reduce contamination by spurious sources and may prove to be a useful addition to any multi-source resolution method. Applied to a recent mock data challenge, the method is able to find $O(10^4)$ Galactic binaries across a signal frequency range of $[0.1,15]$ mHz, and, for frequency $\gtrsim 4$ mHz, reduces the residual data after subtracting out estimated signals to the instrumental noise level.