We report results from searches for new physics with low-energy electronic recoil data recorded with the XENON1T detector, with an exposure of 0.65 t-y and an unprecedentedly low background rate of $76 \pm 2_{stat}$ events/(t y keV) between 1-30 keV. An excess over known backgrounds is observed below 7 keV, rising towards lower energies and prominent between 2-3 keV. The solar axion model has a 3.5$\sigma$ significance, and a three-dimensional 90% confidence surface is reported for axion couplings to electrons, photons, and nucleons. This surface is inscribed in the cuboid defined by $g_{ae} < 3.7 \times 10^{-12}$, $g_{ae}g_{an}^{eff} < 4.6 \times 10^{-18}$, and $g_{ae}g_{a\gamma} < 7.6\times10^{-22}~{GeV}^{-1}$, and excludes either $g_{ae}=0$ or $g_{ae}g_{a\gamma}=g_{ae}g_{an}^{eff}=0$. The neutrino magnetic moment signal is similarly favored over background at 3.2$\sigma$ and a confidence interval of $\mu_{\nu} \in (1.4,~2.9)\times10^{-11}\mu_B$ (90% C.L.) is reported. Both results are in strong tension with stellar constraints. The excess can also be explained by $\beta$ decays of tritium at 3.2$\sigma$ significance with a corresponding tritium concentration in xenon of $(6.2 \pm 2.0) \times 10^{-25}$ mol/mol. Such a trace amount can be neither confirmed nor excluded with current knowledge of production and reduction mechanisms. The significances of the solar axion and neutrino magnetic moment hypotheses are decreased to 2.1$\sigma$ and 0.9$\sigma$, respectively, if an unconstrained tritium component is included in the fitting. With respect to bosonic dark matter, the excess favors a monoenergetic peak at ($2.3 \pm 0.2$) keV (68% C.L.) with a 3.0$\sigma$ global (4.0$\sigma$ local) significance over background. This analysis sets the most restrictive direct constraints to date on pseudoscalar and vector bosonic dark matter for most masses between 1 and 210 keV/c$^2$.

The hard X-ray transient source Swift J1845.7-0037 was discovered in 2012 by Swift/BAT. However, at that time no dedicated observations of the source were performed. On Oct 2019 the source became active again, and X-ray pulsations with a period of ~199s were detected with Swift/XRT. This triggered follow-up observations with NuSTAR. Here we report on the timing and spectral analysis of the source properties using NuSTAR and Swift/XRT. The main goal was to confirm pulsations and search for possible cyclotron lines in the broadband spectrum of the source to probe its magnetic field. Despite highly significant pulsations with period of 207.379(2) were detected, no evidence for a cyclotron line was found in the spectrum of the source. We therefore discuss the strength of the magnetic field based on the source flux and the detection of the transition to the "cold-disc" accretion regime during the 2012 outburst. Our conclusion is that, most likely, the source is a highly magnetized neutron star with B 1e13G at a large distance of d~10 kpc. The latter one consistent with the non-detection of a cyclotron line in the NuSTAR energy band.

We analytically compute time domain gravitational waveforms produced in the final stages of extreme mass ratio inspirals of non-spinning compact objects into supermassive nearly extremal Kerr black holes. Conformal symmetry relates all corotating equatorial orbits in the geodesic approximation to circular orbits through complex conformal transformations. We use this to obtain the time domain Teukolsky perturbations for generic equatorial corotating plunges in closed form. The resulting gravitational waveforms consist of an intermediate polynomial ringdown phase in which the decay rate depends on the impact parameters, followed by an exponential quasi-normal mode decay. The waveform amplitude exhibits critical behavior when the orbital angular momentum tends to a minimal value determined by the innermost stable circular orbit. We show that either near-critical or large angular momentum leads to a significant extension of the LISA observable volume of gravitational wave sources of this kind.

Many studies have shown that RR Lyrae variable stars (RRL) are powerful stellar tracers of Galactic halo structure and satellite galaxies. The Dark Energy Survey (DES), with its deep and wide coverage (g ~ 23.5 mag) in a single exposure; over 5000 deg$^{2}$) provides a rich opportunity to search for substructures out to the edge of the Milky Way halo. However, the sparse and unevenly sampled multiband light curves from the DES wide-field survey (median 4 observations in each of grizY over the first three years) pose a challenge for traditional techniques used to detect RRL. We present an empirically motivated and computationally efficient template fitting method to identify these variable stars using three years of DES data. When tested on DES light curves of previously classified objects in SDSS stripe 82, our algorithm recovers 89% of RRL periods to within 1% of their true value with 85% purity and 76% completeness. Using this method, we identify 5783 RRL candidates, ~31% of which are previously undiscovered. This method will be useful for identifying RRL in other sparse multiband data sets.

We use deep \textit{Hubble Space Telescope} spectroscopy to constrain the metallicities and (\editone{light-weighted}) ages of massive ($\log M_\ast/M_\odot\gtrsim10$) galaxies selected to have quiescent stellar populations at $1.0<z<1.8$. The data include 12--orbit depth coverage with the WFC3/G102 grism covering $\sim$ $8,000<\lambda<11,500$~\AA\, at a spectral resolution of $R\sim 210$ taken as part of the CANDELS Lyman-$\alpha$ Emission at Reionization (CLEAR) survey. At $1.0<z<1.8$, the spectra cover important stellar population features in the rest-frame optical. We simulate a suite of stellar population models at the grism resolution, fit these to the data for each galaxy, and derive posterior likelihood distributions for metallicity and age. We stack the posteriors for subgroups of galaxies in different redshift ranges that include different combinations of stellar absorption features. Our results give \editone{light-weighted ages of $t_{z \sim 1.1}= 3.2\pm 0.7$~Gyr, $t_{z \sim 1.2}= 2.2\pm 0.6$~Gyr, $t_{z\sim1.3}= 3.1\pm 0.6$~Gyr, and $t_{z\sim1.6}= 2.0 \pm 0.6$~Gyr, \editone{for galaxies at $z\sim 1.1$, 1.2, 1.3, and 1.6. This} implies that most of the massive quiescent galaxies at $1<z<1.8$ had formed $>68$\% of their stellar mass by a redshift of $z>2$}. The posteriors give metallicities of \editone{$Z_{z\sim1.1}=1.16 \pm 0.29$~$Z_\odot$, $Z_{z\sim1.2}=1.05 \pm 0.34$~$Z_\odot$, $Z_{z\sim1.3}=1.00 \pm 0.31$~$Z_\odot$, and $Z_{z\sim1.6}=0.95 \pm 0.39$~$Z_\odot$}. This is evidence that massive galaxies had enriched rapidly to approximately Solar metallicities as early as $z\sim3$.

PDS 70b is a recently discovered and directly imaged exoplanet within the wide ($\gtrsim$40 au) cavity around PDS 70 (Keppler et al. 2018, M\"uller et al. 2018). Ongoing accretion onto the central star suggests that accretion onto PDS 70b may also be ongoing. We present the first high contrast images at H$\alpha$ (656 nm) and nearby continuum (643 nm) of PDS 70 utilizing the MagAO system. The combination of these filters allows for the accretion rate of the young planet to be inferred, as hot infalling hydrogen gas will emit strongly at H$\alpha$ over the optical continuum. We detected a source in H$\alpha$ at the position of PDS 70b on two sequential nights in May 2018, for which we establish a false positive probability of $<$0.1%. We conclude that PDS 70b is a young, actively accreting planet. We utilize the H$\alpha$ line luminosity to derive a mass accretion rate of $\dot M= 10^{-8\pm1}$ M$_{Jup}/yr$, where the large uncertainty is primarily due to the unknown amount of optical extinction from the circumstellar and circumplanetary disks. PDS 70b represents the second case of an accreting planet interior to a disk gap, and is among the early examples of a planet observed during its formation.

Despite promising astrometric signals, to date there has been no success in direct imaging of a hypothesized third member of the Sirius system. Using the Clio instrument and MagAO adaptive optics system on the Magellan Clay 6.5 m telescope, we have obtained extensive imagery of Sirius through a vector apodizing phase plate (vAPP) coronagraph in a narrowband filter at 3.9 microns. The vAPP coronagraph and MagAO allow us to be sensitive to planets much less massive than the limits set by previous non-detections. However, analysis of these data presents challenges due to the target's brightness and unique characteristics of the instrument. We present a comparison of dimensionality reduction techniques to construct background illumination maps for the whole detector using the areas of the detector that are not dominated by starlight. Additionally, we describe a procedure for sub-pixel alignment of vAPP data using a physical-optics-based model of the coronagraphic PSF.

MagAO-X is an entirely new "extreme" adaptive optics system for the Magellan Clay 6.5 m telescope, funded by the NSF MRI program starting in Sep 2016. The key science goal of MagAO-X is high-contrast imaging of accreting protoplanets at H$\alpha$. With 2040 actuators operating at up to 3630 Hz, MagAO-X will deliver high Strehls (>70%), high resolution (19 mas), and high contrast ($< 1\times10^{-4}$) at H$\alpha$ (656 nm). We present an overview of the MagAO-X system, review the system design, and discuss the current project status.

Adaptive optics systems correct atmospheric turbulence in real time. Most adaptive optics systems used routinely correct in the near infrared, at wavelengths greater than 1 micron. MagAO- X is a new extreme adaptive optics (ExAO) instrument that will offer corrections at visible-to- near-IR wavelengths. MagAO-X will achieve Strehl ratios greater than 70% at H-alpha when running the 2040 actuator deformable mirror at 3.6 kHz. A visible pyramid wavefront sensor (PWFS) optimized for sensing at 600-1000 nm wavelengths will provide the high-order wavefront sensing on MagAO- X. We present the optical design and predicted performance of the MagAO-X pyramid wavefront sensor.

We present a review of data types and statistical methods often encountered in astronomy. The aim is to provide an introduction to statistical applications in astronomy for statisticians and computer scientists. We highlight the complex, often hierarchical, nature of many astronomy inference problems and advocate for cross-disciplinary collaborations to address these challenges.

The Axion Resonant InterAction Detection Experiment (ARIADNE) is a collaborative effort to search for the QCD axion using techniques based on nuclear magnetic resonance. In the experiment, axions or axion-like particles would mediate short-range spin-dependent interactions between a laser-polarized 3He gas and a rotating (unpolarized) tungsten source mass, acting as a tiny, fictitious "magnetic field". The experiment has the potential to probe deep within the theoretically interesting regime for the QCD axion in the mass range of 0.1-10 meV, independently of cosmological assumptions. The experiment relies on a stable rotary mechanism and superconducting magnetic shielding, required to screen the 3He sample from ordinary magnetic noise. Progress on testing the stability of the rotary mechanism is reported, and the design for the superconducting shielding is discussed.

We present the discovery of 1847 Mira candidates in the Local Group galaxy M33 using a novel semi-parametric periodogram technique coupled with a Random Forest classifier. The algorithms were applied to ~2.4x10^5 I-band light curves previously obtained by the M33 Synoptic Stellar Survey. We derive preliminary Period-Luminosity relations at optical, near- & mid-infrared wavelengths and compare them to the corresponding relations in the Large Magellanic Cloud.

We develop a non-linear semi-parametric Gaussian process model to estimate periods of Miras with sparsely-sampled light curves. The model uses a sinusoidal basis for the periodic variation and a Gaussian process for the stochastic changes. We use maximum likelihood to estimate the period and the parameters of the Gaussian process, while integrating out the effects of other nuisance parameters in the model with respect to a suitable prior distribution obtained from earlier studies. Since the likelihood is highly multimodal for period, we implement a hybrid method that applies the quasi-Newton algorithm for Gaussian process parameters and search the period/frequency parameter over a dense grid. A large-scale, high-fidelity simulation is conducted to mimic the sampling quality of Mira light curves obtained by the M33 Synoptic Stellar Survey. The simulated data set is publicly available and can serve as a testbed for future evaluation of different period estimation methods. The semi-parametric model outperforms an existing algorithm on this simulated test data set as measured by period recovery rate and quality of the resulting Period-Luminosity relations.

Six planetary nebulae (PN) are known in the Kepler space telescope field of view, three newly identified. Of the 5 central stars of PN with useful Kepler data, one, J193110888+4324577, is a short-period, post common envelope binary exhibiting relativistic beaming effects. A second, the central star of the newly identified PN Pa5, has a rare O(He) spectral type and a periodic variability consistent with an evolved companion, where the orbital axis is almost aligned with the line of sight. The third PN, NGC~6826 has a fast rotating central star, something that can only be achieved in a merger. Fourth, the central star of the newly identified PN Kn61, has a PG1159 spectral type and a mysterious semi-periodic light variability which we conjecture to be related to the interplay of binarity with a stellar wind. Finally, the central star of the circular PN A61 does not appear to have a photometric variability above 2 mmag. With the possible exception of the variability of Kn61, all other variability behaviour, whether due to binarity or not, would not easily have been detected from the ground. We conclude, based on very low numbers, that there may be many more close binary or close binary products to be discovered with ultra-high precision photometry. With a larger number of high precision photometric observations we will be able to determine how much higher than the currently known 15 per cent, the short period binary fraction for central stars of PN is likely to be.

Dust attenuation affects nearly all observational aspects of galaxy evolution, yet very little is known about the form of the dust-attenuation law in the distant Universe. Here, we model the spectral energy distributions (SEDs) of galaxies at z = 1.5--3 from CANDELS with rest-frame UV to near-IR imaging under different assumptions about the dust law, and compare the amount of inferred attenuated light with the observed infrared (IR) luminosities. Some individual galaxies show strong Bayesian evidence in preference of one dust law over another, and this preference agrees with their observed location on the plane of infrared excess (IRX, $L_{\text{TIR}}/L_{\text{UV}}$) and UV slope ($\beta$). We generalize the shape of the dust law with an empirical model, $A_{\lambda,\delta}=E(B-V)\ k_\lambda\ (\lambda/\lambda_V)^\delta$ where $k_\lambda$ is the dust law of Calzetti et al. (2000), and show that there exists a correlation between the color excess ${E(B-V)}$ and tilt $\delta$ with ${\delta=(0.62\pm0.05)\log(E(B-V))}$+ ${(0.26~\pm~0.02)}$. Galaxies with high color excess have a shallower, starburst-like law, and those with low color excess have a steeper, SMC-like law. Surprisingly, the galaxies in our sample show no correlation between the shape of the dust law and stellar mass, star-formation rate, or $\beta$. The change in the dust law with color excess is consistent with a model where attenuation is caused by by scattering, a mixed star-dust geometry, and/or trends with stellar population age, metallicity, and dust grain size. This rest-frame UV-to-near-IR method shows potential to constrain the dust law at even higher ($z>3$) redshifts.

We present a new method of extending the single band Analysis of Variance period estimation algorithm to multiple bands. We use SDSS Stripe 82 RR Lyrae to show that in the case of low number of observations per band and non-simultaneous observations, improvements in period recovery rates of up to $\approx$60\% are observed. We also investigate the effect of inter-band observing cadence on period recovery rates. We find that using non-simultaneous observation times between bands is ideal for the multiband method, and using simultaneous multiband data is only marginally better than using single band data. These results will be particularly useful in planning observing cadences for wide-field astronomical imaging surveys such as LSST. They also have the potential to improve the extraction of transient data from surveys with few ($\lesssim 30$) observations per band across several bands, such as the Dark Energy Survey.

Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Spitzer Extended Deep Survey. We describe an updated, Bayesian spectral-energy distribution fitting method that incorporates effects of nebular line emission, star formation histories that are constant or rising with time, and different dust attenuation prescriptions (starburst and Small Magellanic Cloud). From $z$=6.5 to $z$=3.5 star-forming galaxies in CANDELS follow a nearly unevolving correlation between stellar mass and SFR that follows SFR $\sim$ $M_\star^a$ with $a = 0.54 \pm 0.16$ at $z\sim 6$ and $0.70 \pm 0.21$ at $z\sim 4$. This evolution requires a star formation history that increases with decreasing redshift (on average, the SFRs of individual galaxies rise with time). The observed scatter in the SFR-stellar mass relation is tight, $\sigma(\log \mathrm{SFR}/\mathrm{M}_\odot$ yr$^{-1})< 0.3\ - $ 0.4 dex, for galaxies with $\log M_\star/\mathrm{M}_\odot > 9$ dex. Assuming that the SFR is tied to the net gas inflow rate (SFR $\sim$ $\dot{M}_\mathrm{gas}$), then the scatter in the gas inflow rate is also smaller than 0.3$-$0.4 dex for star-forming galaxies in these stellar mass and redshift ranges, at least when averaged over the timescale of star formation. We further show that the implied star formation history of objects selected on the basis of their co-moving number densities is consistent with the evolution in the SFR-stellar mass relation.

We consider the estimation of a common period for a set of functions sampled at irregular intervals. The problem arises in astronomy, where the functions represent a star's brightness observed over time through different photometric filters. While current methods can estimate periods accurately provided that the brightness is well--sampled in at least one filter, there are no existing methods that can provide accurate estimates when no brightness function is well--sampled. In this paper we introduce two new methods for period estimation when brightnesses are poorly--sampled in all filters. The first, multiband generalized Lomb-Scargle (MGLS), extends the frequently used Lomb-Scargle method in a way that na\"{i}vely combines information across filters. The second, penalized generalized Lomb-Scargle (PGLS), builds on the first by more intelligently borrowing strength across filters. Specifically, we incorporate constraints on the phases and amplitudes across the different functions using a non--convex penalized likelihood function. We develop a fast algorithm to optimize the penalized likelihood by combining block coordinate descent with the majorization-minimization (MM) principle. We illustrate our methods on synthetic and real astronomy data. Both advance the state-of-the-art in period estimation; however, PGLS significantly outperforms MGLS when all functions are extremely poorly--sampled.

The observed deficit of strongly Lyman-alpha emitting galaxies at z>6.5 is attributed to either increasing neutral hydrogen in the intergalactic medium (IGM) and/or to the evolving galaxy properties. To investigate this, we have performed very deep near-IR spectroscopy of z>7 galaxies using MOSFIRE on the Keck-I Telescope. We measure the Lyman-alpha fraction at z~8 (combined photometric redshift peak at z=7.7) using two methods. First, we derived NLy{\alpha}/Ntot directly using extensive simulations to correct for incompleteness. Second, we used a Bayesian formalism (introduced by Treu et al. 2012) that compares the z>7 galaxy spectra to models of the Lyman-alpha equivalent width (WLy{\alpha}) distribution at z~6. We explored two simple evolutionary scenarios: smooth evolution where Lyman-alpha is attenuated in all galaxies by a constant factor (perhaps owing to processes from galaxy evolution or a slowly increasing IGM opacity), and patchy evolution where Lyman-alpha is blocked in some fraction of galaxies (perhaps due to the IGM being opaque along only some fraction of sightlines). The Bayesian formalism places stronger constraints compared with the direct method. Combining our data with that in the literature we find that at z~8 the Lyman-alpha fraction has dropped by a factor >3(84% confidence interval) using both the smooth and patchy scenarios compared to the z~6 values. Furthermore, we find a tentative evidence that the data favor the patchy scenario over smooth (with "positive" Bayesian evidence), extending trends observed at z~7 to higher redshift. If this decrease is a result of reionization as predicted by theory, then our data imply the volume averaged neutral hydrogen fraction in the IGM to be >0.3 suggesting that the reionization of the universe is in progress at z~8.

We report on a study of GS 1826--238 using all available {\it RXTE} observations, concentrating on the behavior of the hard X-rays during type-I bursts. We find a hard X-ray shortage at 30--50 keV promoted by the shower of soft X-rays coming from type-I bursts. This shortage happens with a time delay after the peak of the soft flux of 3.6 $\pm$ 1.2 seconds.The behavior of hard X-rays during bursts indicates cooling and reheating of the corona, during which a large amount of energy is required. We speculate that this energy originates from the feedback of the type-I bursts to the accretion process, resulting in a rapid temporary increase of the accretion rate.

Efficient and automated classification of periodic variable stars is becoming increasingly important as the scale of astronomical surveys grows. Several recent papers have used methods from machine learning and statistics to construct classifiers on databases of labeled, multi--epoch sources with the intention of using these classifiers to automatically infer the classes of unlabeled sources from new surveys. However, the same source observed with two different synoptic surveys will generally yield different derived metrics (features) from the light curve. Since such features are used in classifiers, this survey-dependent mismatch in feature space will typically lead to degraded classifier performance. In this paper we show how and why feature distributions change using OGLE and \textit{Hipparcos} light curves. To overcome survey systematics, we apply a method, \textit{noisification}, which attempts to empirically match distributions of features between the labeled sources used to construct the classifier and the unlabeled sources we wish to classify. Results from simulated and real--world light curves show that noisification can significantly improve classifier performance. In a three--class problem using light curves from \textit{Hipparcos} and OGLE, noisification reduces the classifier error rate from 27.0% to 7.0%. We recommend that noisification be used for upcoming surveys such as Gaia and LSST and describe some of the promises and challenges of applying noisification to these surveys.

As the number of observed Gamma-Ray Bursts (GRBs) continues to grow, follow-up resources need to be used more efficiently in order to maximize science output from limited telescope time. As such, it is becoming increasingly important to rapidly identify bursts of interest as soon as possible after the event, before the afterglows fade beyond detectability. Studying the most distant (highest redshift) events, for instance, remains a primary goal for many in the field. Here we present our Random forest Automated Triage Estimator for GRB redshifts (RATE GRB-z) for rapid identification of high-redshift candidates using early-time metrics from the three telescopes onboard Swift. While the basic RATE methodology is generalizable to a number of resource allocation problems, here we demonstrate its utility for telescope-constrained follow-up efforts with the primary goal to identify and study high-z GRBs. For each new GRB, RATE GRB-z provides a recommendation - based on the available telescope time - of whether the event warrants additional follow-up resources. We train RATE GRB-z using a set consisting of 135 Swift bursts with known redshifts, only 18 of which are z > 4. Cross-validated performance metrics on this training data suggest that ~56% of high-z bursts can be captured from following up the top 20% of the ranked candidates, and ~84% of high-z bursts are identified after following up the top ~40% of candidates. We further use the method to rank 200+ Swift bursts with unknown redshifts according to their likelihood of being high-z.

Despite the great promise of machine-learning algorithms to classify and predict astrophysical parameters for the vast numbers of astrophysical sources and transients observed in large-scale surveys, the peculiarities of the training data often manifest as strongly biased predictions on the data of interest. Typically, training sets are derived from historical surveys of brighter, more nearby objects than those from more extensive, deeper surveys (testing data). This sample selection bias can cause catastrophic errors in predictions on the testing data because a) standard assumptions for machine-learned model selection procedures break down and b) dense regions of testing space might be completely devoid of training data. We explore possible remedies to sample selection bias, including importance weighting (IW), co-training (CT), and active learning (AL). We argue that AL---where the data whose inclusion in the training set would most improve predictions on the testing set are queried for manual follow-up---is an effective approach and is appropriate for many astronomical applications. For a variable star classification problem on a well-studied set of stars from Hipparcos and OGLE, AL is the optimal method in terms of error rate on the testing data, beating the off-the-shelf classifier by 3.4% and the other proposed methods by at least 3.0%. To aid with manual labeling of variable stars, we developed a web interface which allows for easy light curve visualization and querying of external databases. Finally, we apply active learning to classify variable stars in the ASAS survey, finding dramatic improvement in our agreement with the ACVS catalog, from 65.5% to 79.5%, and a significant increase in the classifier's average confidence for the testing set, from 14.6% to 42.9%, after a few AL iterations.