We present a differentiable, end-to-end Bayesian forward modeling framework for line intensity mapping cosmology experiments, with a specific focus on low-frequency radio telescopes targeting the redshifted 21 cm line from neutral hydrogen as a cosmological probe. Our framework is capable of posterior density estimation of the cosmological signal jointly with foreground and telescope parameters at the field level. Our key aim is to be able to optimize the model's high-dimensional, non-linear, and ill-conditioned parameter space, while also sampling from it to perform robust uncertainty quantification within a Bayesian framework. We show how a differentiable programming paradigm, accelerated by recent advances in machine learning software and hardware, can make this computationally-demanding, end-to-end Bayesian approach feasible. We demonstrate a proof-of-concept on a signal recovery problem for the Hydrogen Epoch of Reionization Array experiment, highlighting the framework's ability to build confidence in early 21 cm signal detections even in the presence of poorly understood foregrounds and instrumental systematics. We use a Hessian-preconditioned Hamiltonian Monte Carlo algorithm to efficiently sample our parameter space with a dimensionality approaching $N\sim10^5$, which enables joint, end-to-end nuisance parameter marginalization over foreground and instrumental terms. Lastly, we introduce a new spherical harmonic formalism that is a complete and orthogonal basis on the cut sky relevant to drift-scan radio surveys, which we call the spherical stripe harmonic formalism, and it's associated three-dimensional basis, the spherical stripe Fourier-Bessel formalism.

The Gamow Explorer will use Gamma Ray Bursts (GRBs) to: 1) probe the high redshift universe (z > 6) when the first stars were born, galaxies formed and Hydrogen was reionized; and 2) enable multi-messenger astrophysics by rapidly identifying Electro-Magnetic (IR/Optical/X-ray) counterparts to Gravitational Wave (GW) events. GRBs have been detected out to z ~ 9 and their afterglows are a bright beacon lasting a few days that can be used to observe the spectral fingerprints of the host galaxy and intergalactic medium to map the period of reionization and early metal enrichment. Gamow Explorer is optimized to quickly identify high-z events to trigger follow-up observations with JWST and large ground-based telescopes. A wide field of view Lobster Eye X-ray Telescope (LEXT) will search for GRBs and locate them with arc-minute precision. When a GRB is detected, the rapidly slewing spacecraft will point the 5 photometric channel Photo-z Infra-Red Telescope (PIRT) to identify high redshift (z > 6) long GRBs within 100s and send an alert within 1000s of the GRB trigger. An L2 orbit provides > 95% observing efficiency with pointing optimized for follow up by the James Webb Space Telescope (JWST) and ground observatories. The predicted Gamow Explorer high-z rate is >10 times that of the Neil Gehrels Swift Observatory. The instrument and mission capabilities also enable rapid identification of short GRBs and their afterglows associated with GW events. The Gamow Explorer will be proposed to the 2021 NASA MIDEX call and if approved, launched in 2028.

The so-called corner turning problem is a major bottleneck for radio telescopes with large numbers of antennas. The problem is essentially that of rapidly transposing a matrix that is too large to store on one single device; in radio interferometry, it occurs because data from each antenna needs to be routed to an array of processors that will each handle a limited portion of the data (a frequency range, say) but requires input from each antenna. We present a low-cost solution allowing the correlator to transpose its data in real time, without contending for bandwidth, via a butterfly network requiring neither additional RAM memory nor expensive general-purpose switching hardware. We discuss possible implementations of this using FPGA, CMOS, analog logic and optical technology, and conclude that the corner turner cost can be small even for upcoming massive radio arrays.

The Gamow Explorer will use Gamma Ray Bursts (GRBs) to: 1) probe the high redshift universe (z > 6) when the first stars were born, galaxies formed and Hydrogen was reionized; and 2) enable multi-messenger astrophysics by rapidly identifying Electro-Magnetic (IR/Optical/X-ray) counterparts to Gravitational Wave (GW) events. GRBs have been detected out to z ~ 9 and their afterglows are a bright beacon lasting a few days that can be used to observe the spectral fingerprints of the host galaxy and intergalactic medium to map the period of reionization and early metal enrichment. Gamow Explorer is optimized to quickly identify high-z events to trigger follow-up observations with JWST and large ground-based telescopes. A wide field of view Lobster Eye X-ray Telescope (LEXT) will search for GRBs and locate them with arc-minute precision. When a GRB is detected, the rapidly slewing spacecraft will point the 5 photometric channel Photo-z Infra-Red Telescope (PIRT) to identify high redshift (z > 6) long GRBs within 100s and send an alert within 1000s of the GRB trigger. An L2 orbit provides > 95% observing efficiency with pointing optimized for follow up by the James Webb Space Telescope (JWST) and ground observatories. The predicted Gamow Explorer high-z rate is >10 times that of the Neil Gehrels Swift Observatory. The instrument and mission capabilities also enable rapid identification of short GRBs and their afterglows associated with GW events. The Gamow Explorer will be proposed to the 2021 NASA MIDEX call and if approved, launched in 2028.

We used TESS observations to search for pulsations in six known Of?p stars in the Magellanic Clouds. We find evidence for pulsational variability in three Of?p stars: UCAC4 115-008604, OGLE SMC-SC6 237339, and AzV 220. Two of them, UCAC4 115-008604 and OGLE SMC-SC6 237339, have been reported to possess kG-order magnetic fields. The obtained results are important to constrain and improve stellar evolution models.

Exoplanet catalogs produced by surveys suffer from a lack of completeness (not every planet is detected) and less than perfect reliability (not every planet in the catalog is a true planet), particularly near the survey's detection limit. Exoplanet occurrence rate studies based on such a catalog must be corrected for completeness and reliability. The final Kepler data release, DR25, features a uniformly vetted planet candidate catalog and data products that facilitate corrections. We present a new probabilistic approach to the characterization of Kepler completeness and reliability, making full use of the Kepler DR25 products. We illustrate the impact of completeness and reliability corrections with a Poisson-likelihood occurrence rate method, using a recent stellar properties catalog that incorporates Gaia stellar radii and essentially uniform treatment of the stellar population. Correcting for reliability has a significant impact: the exoplanet occurrence rate for orbital period and radius within 20% of Earth's around GK dwarf stars, corrected for reliability, is 0.015+0.011-0.007, whereas not correcting results in 0.034+0.018-0.012 - correcting for reliability reduces this occurrence rate by more than a factor of two. We further show that using Gaia-based vs. DR25 stellar properties impacts the same occurrence rate by a factor of two. We critically examine the the DR25 catalog and the assumptions behind our occurrence rate method. We propose several ways in which confidence in both the Kepler catalog and occurrence rate calculations can be improved. This work provides an example of how the community can use the DR25 completeness and reliability products.

We derive new terms in the post-Newtonian (PN) expansion of the generalized redshift invariant $\langle u^t \rangle_\tau$ for a small body in eccentric, equatorial orbit about a massive Kerr black hole. The series is computed analytically using the Teukolsky formalism for first-order black hole perturbation theory (BHPT), along with the CCK method for metric reconstruction using the Hertz potential in ingoing radiation gauge. Modal contributions with small values of $l$ are derived via the semi-analytic solution of Mano-Suzuki-Takasugi (MST), while the remaining values of $l$ to infinity are determined via direct expansion of the Teukolsky equation. Each PN order is calculated as a series in eccentricity $e$ but kept exact in the primary black hole's spin parameter $a$. In total, the PN terms are expanded to $e^{16}$ through 6PN relative order, and separately to $e^{10}$ through 8PN relative order. Upon grouping eccentricity coefficients by spin dependence, we find that many resulting component terms can be simplified to closed-form functions of eccentricity, in close analogy to corresponding terms derived previously in the Schwarzschild limit. We use numerical calculations to compare convergence of the full series to its Schwarzschild counterpart and discuss implications for gravitational wave analysis.

Studies of the magnetic characteristics of massive stars have recently received significant attention because they are progenitors of highly magnetised compact objects. Stars initially more massive than about 8M_sun leave behind neutron stars and black holes by the end of their evolution. The merging of binary compact remnant systems produces astrophysical transients detectable by gravitational wave observatories. Studies of magnetic fields in massive stars with low metallicities are of particular interest because they provide important information on the role of magnetic fields in the star formation of the early Universe. While several detections of massive Galactic magnetic stars have been reported in the last few decades, the impact of a low-metallicity environment on the occurrence and strength of stellar magnetic fields has not yet been explored. Because of the similarity between Of?p stars in the Magellanic Clouds (MCs) and Galactic magnetic Of?p stars, which possess globally organised magnetic fields, we searched for magnetic fields in Of?p stars in the MCs. Additionally, we observed the massive contact binary Cl* NGC 346 SSN7 in the Small Magellanic Cloud to test the theoretical scenario that the origin of magnetic fields involves a merger event or a common envelope evolution. We obtained and analysed measurements of the magnetic field in four massive Of?p stars in the MCs and the binary Cl* NGC 346 SSN7 using the ESO/VLT FORS2 spectrograph in spectropolarimetric mode. We detected kilogauss-scale magnetic fields in two Of?p-type stars and in the contact binary Cl* NGC 346 SSN7. These results suggest that the impact of low metallicity on the occurrence and strength of magnetic fields in massive stars is low. However, because the explored stellar sample is very small, additional observations of massive stars in the MCs are necessary.

Exoplanet catalogs produced by surveys suffer from a lack of completeness (not every planet is detected) and less than perfect reliability (not every planet in the catalog is a true planet), particularly near the survey's detection limit. Exoplanet occurrence rate studies based on such a catalog must be corrected for completeness and reliability. The final Kepler data release, DR25, features a uniformly vetted planet candidate catalog and data products that facilitate corrections. We present a new probabilistic approach to the characterization of Kepler completeness and reliability, making full use of the Kepler DR25 products. We illustrate the impact of completeness and reliability corrections with a Poisson-likelihood occurrence rate method, using a recent stellar properties catalog that incorporates Gaia stellar radii and essentially uniform treatment of the stellar population. Correcting for reliability has a significant impact: the exoplanet occurrence rate for orbital period and radius within 20% of Earth's around GK dwarf stars, corrected for reliability, is 0.015+0.011-0.007, whereas not correcting results in 0.034+0.018-0.012 - correcting for reliability reduces this occurrence rate by more than a factor of two. We further show that using Gaia-based vs. DR25 stellar properties impacts the same occurrence rate by a factor of two. We critically examine the the DR25 catalog and the assumptions behind our occurrence rate method. We propose several ways in which confidence in both the Kepler catalog and occurrence rate calculations can be improved. This work provides an example of how the community can use the DR25 completeness and reliability products.

Complex periodic variables (CPVs) are stars that exhibit highly structured and periodic optical light curves. Previous studies have indicated that these stars are typically disk-free pre-main-sequence M dwarfs with rotation periods ranging from 0.2 to 2 days. To advance our understanding of these enigmatic objects, we conducted a blind search using TESS 2-minute data of 65,760 K and M dwarfs with $T$<16 mag and $d$<150 pc. We found 50 high-quality CPVs, and subsequently determined that most are members of stellar associations. Among the new discoveries are the brightest ($T$$\approx$9.5 mag), closest ($d$$\approx$20 pc), and oldest ($\approx$200 Myr) CPVs known. One exceptional object, LP 12-502, exhibited up to eight flux dips per cycle. Some of these dips coexisted with slightly different periods, and the shortest-duration dips precisely matched the expected timescale for transiting small bodies at the corotation radius. Broadly, our search confirms that CPVs are mostly young ($\lesssim$150 Myr) and low-mass ($\lesssim$0.4 $M_\odot$). The flux dips characteristic of the class have lifetimes of $\approx$100 cycles, although stellar flares seem to induce sudden dip collapse once every few months. The most plausible explanation for these phenomena remains corotating concentrations of gas or dust. The gas or dust is probably entrained by the star's magnetic field, and the sharp features could result from a multipolar field topology, a hypothesis supported by correspondences between the light curves of CPVs and of rapidly rotating B stars known to have multipolar magnetic fields.

We study the effects of a thin gaseous accretion disk on the inspiral of a stellar--mass black hole into a supermassive black hole. We construct a phenomenological angular momentum transport equation that reproduces known disk effects. Disk torques modify the gravitational wave phase evolution to detectable levels with LISA for reasonable disk parameters. The Fourier transform of disk-modified waveforms acquires a correction with a different frequency trend than post-Newtonian vacuum terms. Such inspirals could be used to detect accretion disks with LISA and to probe their physical parameters.

We calculate the eccentricity dependence of the high-order post-Newtonian (PN) series for the generalized redshift invariant $\langle u^t \rangle_\tau$ for eccentric-orbit extreme-mass-ratio inspirals on a Schwarzschild background. These results are calculated within first-order black hole perturbation theory (BHPT) using Regge-Wheeler-Zerilli (RWZ) gauge. Our \textsc{Mathematica} code is based on a familiar procedure, using PN expansion of the Mano-Suzuki-Takasugi (MST) analytic function formalism for $l$ modes up to a certain maximum and then using a direct general-$l$ PN expansion of the RWZ equation for arbitrarily high $l$. We calculate dual expansions in PN order and in powers of eccentricity, reaching 10PN relative order and $e^{20}$. Detailed knowledge of the eccentricity expansion at each PN order allows us to find within the eccentricity dependence numerous closed-form expressions and multiple infinite series with known coefficients. We find leading logarithm sequences in the PN expansion of the redshift invariant that reflect a similar behavior in the PN expansion of the energy flux to infinity. A set of flux terms and special functions that appear in the energy flux, like the Peters-Mathews flux itself, are shown to reappear in the redshift PN expansion.

Malignant melanoma is a common skin cancer that is mostly curable before metastasis -when growths spawn in organs away from the original site. Melanoma is the most dangerous type of skin cancer if left untreated due to the high risk of metastasis. This paper presents Melatect, a machine learning (ML) model embedded in an iOS app that identifies potential malignant melanoma. Melatect accurately classifies lesions as malignant or benign over 96.6% of the time with no apparent bias or overfitting. Using the Melatect app, users have the ability to take pictures of skin lesions (moles) and subsequently receive a mole classification. The Melatect app provides a convenient way to get free advice on lesions and track these lesions over time. A recursive computer image analysis algorithm and modified MLOps pipeline was developed to create a model that performs at a higher accuracy than existing models. Our training dataset included 18,400 images of benign and malignant lesions, including 18,000 from the International Skin Imaging Collaboration (ISIC) archive, as well as 400 images gathered from local dermatologists; these images were augmented using DeepAugment, an AutoML tool, to 54,054 images.

We present an IR-monitoring survey with the $Spitzer$ Space Telescope of the star forming region GGD 12-15. Over 1000 objects were monitored including about 350 objects within the central 5 arcminutes which is found to be especially dense in cluster members. The monitoring took place over 38 days and is part of the Young Stellar Object VARiability (YSOVAR) project. The region was also the subject of a contemporaneous 67ks $Chandra$ observation. The field includes 119 previously identified pre-main sequence star candidates. X-rays are detected from 164 objects, 90 of which are identified with cluster members. Overall, we find that about half the objects in the central 5 arcminutes are young stellar objects based on a combination of their spectral energy distribution, IR variability and X-ray emission. Most of the stars with IR excess relative to a photosphere show large amplitude (>0.1 mag) mid-IR variability. There are 39 periodic sources, all but one of these is found to be a cluster member. Almost half of the periodic sources do not show IR excesses. Overall, more than 85% of the Class I, flat spectrum, and Class II sources are found to vary. The amplitude of the variability is larger in more embedded young stellar objects. Most of the Class~I/II objects exhibit redder colors in a fainter state, compatible with time-variable extinction. A few become bluer when fainter, which can be explained with significant changes in the structure of the inner disk. A search for changes in the IR due to X-ray events is carried out, but the low number of flares prevented an analysis of the direct impact of X-ray flares on the IR lightcurves. However, we find that X-ray detected Class II sources have longer timescales for change in the mid-IR than a similar set of non-X-ray detected Class IIs.

Recent observations support the hypothesis that a large fraction of "short-hard" gamma-ray bursts (SHBs) are associated with compact binary inspiral. Since gravitational-wave (GW) measurements of well-localized inspiraling binaries can measure absolute source distances, simultaneous observation of a binary's GWs and SHB would allow us to independently determine both its luminosity distance and redshift. Such a "standard siren" (the GW analog of a standard candle) would provide an excellent probe of the relatively nearby universe's expansion, complementing other standard candles. In this paper, we examine binary measurement using a Markov Chain Monte Carlo technique to build the probability distributions describing measured parameters. We assume that each SHB observation gives both sky position and the time of coalescence, and we take both binary neutron stars and black hole-neutron star coalescences as plausible SHB progenitors. We examine how well parameters particularly distance) can be measured from GW observations of SHBs by a range of ground-based detector networks. We find that earlier estimates overstate how well distances can be measured, even at fairly large signal-to-noise ratio. The fundamental limitation to determining distance proves to be a degeneracy between distance and source inclination. Overcoming this limitation requires that we either break this degeneracy, or measure enough sources to broadly sample the inclination distribution. (Abridged)

Next generation radio observatories such as the MWA, LWA, LOFAR, CARMA and SKA provide a number of challenges for interferometric data analysis. These challenges include heterogeneous arrays, direction-dependent instrumental gain, and refractive and scintillating atmospheric conditions. From the analysis perspective, this means that calibration solutions can not be described using a single complex gain per antenna. In this paper we use the optimal map-making formalism developed for CMB analyses to extend traditional interferometric radio analysis techniques--removing the assumption of a single complex gain per antenna and allowing more complete descriptions of the instrumental and atmospheric conditions. Due to the similarity with holographic mapping of radio antenna surfaces, we call this extended analysis approach software holography. The resulting analysis algorithms are computationally efficient, unbiased, and optimally sensitive. We show how software holography can be used to solve some of the challenges of next generation observations, and how more familiar analysis techniques can be derived as limiting cases.

The key to detecting neutral hydrogen during the epoch of reionization (EoR) is to separate the cosmological signal from the dominating foreground radiation. We developed direct optimal mapping (DOM) to map interferometric visibilities; it contains only linear operations, with full knowledge of point spread functions from visibilities to images. Here, we demonstrate a fast Fourier transform-based image power spectrum and its window functions computed from the DOM images. We use noiseless simulation, based on the Hydrogen Epoch of Reionization Array Phase I configuration, to study the image power spectrum properties. The window functions show &lt;10^{-11} of the integrated power leaks from the foreground-dominated region into the EoR window; the 2D and 1D power spectra also verify the separation between the foregrounds and the EoR.

We argue that the stellar velocity dispersion observed in an elliptical galaxy is a good proxy for the halo velocity dispersion. As dark matter halos are almost completely characterized by a single scale parameter, the stellar velocity dispersion tells us the virial radius of the halo and the mass contained within. This permits non-dimensionalizing of the stellar mass and effective radius axes of the stellar mass fundamental plane by the virial radius and halo mass, respectively.

We study the X-ray luminosity function (XLF) of low mass X-ray binaries (LMXB) in the nearby early-type galaxy Centaurus A, concentrating primarily on two aspects of binary populations: the XLF behavior at the low luminosity limit and comparison between globular cluster and field sources. The 800 ksec exposure of the deep Chandra VLP program allows us to reach a limiting luminosity of 8e35 erg/s, about 2-3 times deeper than previous investigations. We confirm the presence of the low luminosity break in the overall LMXB XLF at log(L_X)=37.2-37.6 below which the luminosity distribution follows a constant dN/d(ln L). Separating globular cluster and field sources, we find a statistically significant difference between the two luminosity distributions with a relative underabundance of faint sources in the globular cluster population. This demonstrates that the samples are drawn from distinct parent populations and may disprove the hypothesis that the entire LMXB population in early type galaxies is created dynamically in globular clusters. As a plausible explanation for this difference in the XLFs, we suggest that there is an enhanced fraction of helium accreting systems in globular clusters, which are created in collisions between red giants and neutron stars. Due to the 4 times higher ionization temperature of He, such systems are subject to accretion disk instabilities at approximately 20 times higher mass accretion rate, and therefore are not observed as persistent sources at low luminosities.