Search

The electrocardiogram (ECG) is an inexpensive and widely available tool for cardiovascular assessment. Despite its standardized format and small file size, the high complexity and inter-individual variability of ECG signals (typically a 60,000-size vector) make it challenging to use in deep learning models, especially when only small datasets are available. This study addresses these challenges by exploring feature generation methods from representative beat ECGs, focusing on Principal Component Analysis (PCA) and Autoencoders to reduce data complexity. We introduce three novel Variational Autoencoder (VAE) variants: Stochastic Autoencoder (SAE), Annealed beta-VAE (Abeta-VAE), and cyclical beta-VAE (Cbeta-VAE), and compare their effectiveness in maintaining signal fidelity and enhancing downstream prediction tasks. The Abeta-VAE achieved superior signal reconstruction, reducing the mean absolute error (MAE) to 15.7 plus-minus 3.2 microvolts, which is at the level of signal noise. Moreover, the SAE encodings, when combined with ECG summary features, improved the prediction of reduced Left Ventricular Ejection Fraction (LVEF), achieving an area under the receiver operating characteristic curve (AUROC) of 0.901. This performance nearly matches the 0.910 AUROC of state-of-the-art CNN models but requires significantly less data and computational resources. Our findings demonstrate that these VAE encodings are not only effective in simplifying ECG data but also provide a practical solution for applying deep learning in contexts with limited-scale labeled training data.

The prototypical nova super-remnant (NSR) was uncovered around the most rapidly recurring nova (RN), M31N 2008-12a. Simulations of the growth of NSRs revealed that these large structures should exist around all novae, whether classical or recurrent. NSRs consist of large shell-like structures surrounding excavated cavities. Predictions, informed by these simulations, led to the discovery of an extended cavity coincident with the Galactic RN, RS Ophiuchi, in far-infrared archival IRAS images. We propose that this cavity is associated with RS Oph and is therefore evidence of another NSR to be uncovered., Comment: 6 pages, 5 figures; Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters

There has been huge recent interest in the potential of making operational weather forecasts using machine learning techniques. As they become a part of the weather forecasting toolbox, there is a pressing need to understand how well current machine learning models can simulate high-impact weather events. We compare forecasts of Storm Ciar\'an, a European windstorm that caused sixteen deaths and extensive damage in Northern Europe, made by machine learning and numerical weather prediction models. The four machine learning models considered (FourCastNet, Pangu-Weather, GraphCast and FourCastNet-v2) produce forecasts that accurately capture the synoptic-scale structure of the cyclone including the position of the cloud head, shape of the warm sector and location of warm conveyor belt jet, and the large-scale dynamical drivers important for the rapid storm development such as the position of the storm relative to the upper-level jet exit. However, their ability to resolve the more detailed structures important for issuing weather warnings is more mixed. All of the machine learning models underestimate the peak amplitude of winds associated with the storm, only some machine learning models resolve the warm core seclusion and none of the machine learning models capture the sharp bent-back warm frontal gradient. Our study shows there is a great deal about the performance and properties of machine learning weather forecasts that can be derived from case studies of high-impact weather events such as Storm Ciar\'an.

We present early spectral observations of the very slow Galactic nova Gaia22alz, over its gradual rise to peak brightness that lasted 180 days. During the first 50 days, when the nova was only 3--4 magnitudes above its normal brightness, the spectra showed narrow (FWHM $\approx$ 400 km s$^{-1}$) emission lines of H Balmer, He I, He II, and C IV, but no P Cygni absorption. A few weeks later, the high-excitation He II and C IV lines disappeared, and P Cygni profiles of Balmer, He I, and eventually Fe II lines emerged, yielding a spectrum typical of classical novae before peak. We propose that the early spectra of Gaia22alz are produced in the white dwarf's envelope or accretion disk, reprocessing X-ray and ultraviolet emission from the white dwarf after a dramatic increase in the rate of thermonuclear reactions, during a phase known as the ``early X-ray/UV flash''. If true, this would be one of the rare times that the optical signature of the early X-ray/UV flash has been detected. While this phase might last only a few hours in other novae and thus be easily missed, it was possible to detect in Gaia22alz due to its very slow and gradual rise and thanks to the efficiency of new all-sky surveys in detecting transients on their rise. We also consider alternative scenarios that could explain the early spectral features of Gaia22alz and its unusually slow rise., Comment: 20 pages, 12 figures, 2 tables. Submitted to MNRAS

V906 Carinae was one of the best observed novae of recent times. It was a prolific dust producer and harboured shocks in the early evolving ejecta outflow. Here, we take a close look at the consequences of these early interactions through study of high-resolution UVES spectroscopy of the nebular stage and extrapolate backwards to investigate how the final structure may have formed. A study of ejecta geometry and shaping history of the structure of the shell is undertaken following a spectral line SHAPE model fit. A search for spectral tracers of shocks in the nova ejecta is undertaken and an analysis of the ionised environment. Temperature, density and abundance analyses of the evolving nova shell are presented.

The specific activity of the beta decay of $^{39}$Ar in atmospheric argon is measured using the DEAP-3600 detector. DEAP-3600, located 2 km underground at SNOLAB, uses a total of (3269 $\pm$ 24) kg of liquid argon distilled from the atmosphere to search for dark matter. This detector with very low background uses pulseshape discrimination to differentiate between nuclear recoils and electron recoils and is well-suited to measure the decay of $^{39}$Ar. With 167 live-days of data, the measured specific activity at the time of atmospheric extraction is [0.964 $\pm$ 0.001 (stat) $\pm$ 0.024 (sys)] Bq/kg$_{\rm atmAr}$ which is consistent with results from other experiments. A cross-check analysis using different event selection criteria provides a consistent result.

The nova super-remnant (NSR) surrounding M31N 2008-12a (12a), the annually erupting recurrent nova (RN), is the only known example of this phenomenon. As this structure has grown as a result of frequent eruptions from 12a, we might expect to see NSRs around other RNe; this would confirm the RN--NSR association and strengthen the connection between novae and type Ia supernovae (SN Ia) as NSRs centered on SN Ia provide a lasting, unequivocal signpost to the single degenerate progenitor type of that explosion. The only previous NSR simulation used identical eruptions from a static white dwarf (WD). In this Paper, we simulate the growth of NSRs alongside the natural growth/erosion of the central WD, within a range of environments, accretion rates, WD temperatures, and initial WD masses. The subsequent evolving eruptions create dynamic NSRs tens of parsecs in radius comprising a low-density cavity, bordered by a hot ejecta pile-up region, and surrounded by a cool high-density, thin, shell. Higher density environments restrict NSR size, as do higher accretion rates, whereas the WD temperature and initial mass have less impact. NSRs form around growing or eroding WDs, indicating that NSRs also exist around old novae with low-mass WDs. Observables such as X-ray and H$\alpha$ emission from the modelled NSRs are derived to aid searches for more examples; only NSRs around high accretion rate novae will currently be observable. The observed properties of the 12a NSR can be reproduced when considering both the dynamically grown NSR and photoionisation by the nova system., Comment: 20 pages, 20 figures; Accepted for publication in Monthly Notices of the Royal Astronomical Society. For the animation in Figure 4, see https://www.dropbox.com/s/7a6zbu4eyv6kv8x/coolRun1.mp4?dl=0

A graph class $\mathcal{G}$ has linear growth if, for each graph $G \in \mathcal{G}$ and every positive integer $r$, every subgraph of $G$ with radius at most $r$ contains $O(r)$ vertices. In this paper, we show that every graph class with linear growth has bounded treewidth.

The recurrent nova RS Oph underwent a new outburst on August 8, 2021, reaching a visible brightness of V = 4.8 mag. Observations of the 2021 outburst made with the high resolution UVES spectrograph at the Kueyen-UT2 telescope of ESO-VLT in Paranal enabled detection of the possible presence of 7Be freshly made in the thermonuclear runaway reactions. The 7Be yields can be estimated in N(Be)/N(H) = 5.7 x 10^(-6), which are close to the lowest yields measured in classical novae so far. 7Be is short-lived and decays only into 7Li. By means of a spectrum taken during the nebular phase we estimated an ejected mass of about 1.1 x 10^(-5) Msun, providing an amount of about 4.4x 10^(-10) Msun of 7Li created in the 2021 event. Recurrent novae of the kind of RS Oph may synthesize slightly lower amount of 7Li per event as classical novae, but occur 10^3 times more frequently. The recurrent novae fraction is in the range of 10-30% and they could have contributed to the making of 7Li we observe today. The detection of 7Be in RS Oph provides further support to the recent suggestion that novae are the most effective source of 7Li in the Galaxy., Comment: Accepted 2022 September 17. Received 2022 September 16; in original form 2022 April 20. 12 pages, 13 figures

V392 Persei is a known dwarf nova (DN) that underwent a classical nova eruption in 2018. Here we report ground-based optical, Swift UV and X-ray, and Fermi-LAT \gamma-ray observations following the eruption for almost three years. V392 Per is one of the fastest evolving novae yet observed, with a $t_2$ decline time of 2 days. Early spectra present evidence for multiple and interacting mass ejections, with the associated shocks driving both the \gamma-ray and early optical luminosity. V392 Per entered Sun-constraint within days of eruption. Upon exit, the nova had evolved to the nebular phase, and we saw the tail of the super-soft X-ray phase. Subsequent optical emission captured the fading ejecta alongside a persistent narrow line emission spectrum from the accretion disk. Ongoing hard X-ray emission is characteristic of a standing accretion shock in an intermediate polar. Analysis of the optical data reveals an orbital period of 3.230 \pm 0.003 days, but we see no evidence for a white dwarf (WD) spin period. The optical and X-ray data suggest a high mass WD, the pre-nova spectral energy distribution (SED) indicates an evolved donor, and the post-nova SED points to a high mass accretion rate. Following eruption, the system has remained in a nova-like high mass transfer state, rather than returning to the pre-nova DN low mass transfer configuration. We suggest that this high state is driven by irradiation of the donor by the nova eruption. In many ways, V392 Per shows similarity to the well-studied nova and DN GK Persei., Comment: Accepted by MNRAS. 20 pages including references, with 37 pages of supplementary material

Prediction with experts' advice is one of the most fundamental problems in online learning and captures many of its technical challenges. A recent line of work has looked at online learning through the lens of differential equations and continuous-time analysis. This viewpoint has yielded optimal results for several problems in online learning. In this paper, we employ continuous-time stochastic calculus in order to study the discrete-time experts' problem. We use these tools to design a continuous-time, parameter-free algorithm with improved guarantees for the quantile regret. We then develop an analogous discrete-time algorithm with a very similar analysis and identical quantile regret bounds. Finally, we design an anytime continuous-time algorithm with regret matching the optimal fixed-time rate when the gains are independent Brownian Motions; in many settings, this is the most difficult case. This gives some evidence that, even with adversarial gains, the optimal anytime and fixed-time regrets may coincide., Comment: 30 pages, 1 figure. Version 2 diff: minor edits, reorganization for a journal submission, correct statement of Lemma 5.1 and a better formatted proof of the same lemma

Prediction with expert advice is a foundational problem in online learning. In instances with $T$ rounds and $n$ experts, the classical Multiplicative Weights Update method suffers at most $\sqrt{(T/2)\ln n}$ regret when $T$ is known beforehand. Moreover, this is asymptotically optimal when both $T$ and $n$ grow to infinity. However, when the number of experts $n$ is small/fixed, algorithms with better regret guarantees exist. Cover showed in 1967 a dynamic programming algorithm for the two-experts problem restricted to $\{0,1\}$ costs that suffers at most $\sqrt{T/2\pi} + O(1)$ regret with $O(T^2)$ pre-processing time. In this work, we propose an optimal algorithm for prediction with two experts' advice that works even for costs in $[0,1]$ and with $O(1)$ processing time per turn. Our algorithm builds up on recent work on the experts problem based on techniques and tools from stochastic calculus., Comment: 29 pages, 13 pages of main text, published in ALT 2022 (PMLR vol. 167)

We report the search for 7Be isotope in the outbursts of the classical nova V6595 Sgr by means of high resolution UVES observations taken at the ESO VLT in April 2021, about two weeks after discovery and under difficult circumstances due to the pandemic. Narrow absorption components with velocities at about -2620 and -2820 km/s, superposed on broader and shallow absorption, are observed in the outburst spectra for the 7BeII 313.0583, 313.1228 nm doublet resonance lines, as well as in several other elements such as CaII, FeI, MgI, NaI, HI but LiI. Using CaII K line as a reference element, we infer N(7Be)/N(H) ~ 7.4 x 10^{-6}, or ~ 9.8 x 10^{-6} when the 7Be decay is taken into account. The 7Be abundance is about half of the value most frequently measured in novae. The possible presence of over-ionization in the layers where 7Be is detected is also discussed. Observations taken at the Telescopio Nazionale Galileo (TNG) in La Palma 91 days after discovery showed prominent emission lines of Oxygen and Neon which allow to classify the nova as ONe type. Therefore, although 7Be is expected to be higher in CO novae, it is found at comparable levels in both nova types., Comment: 10 pages, 9 figures. Accepted MNRAS

The DEAP-3600 detector searches for the scintillation signal from dark matter particles scattering on a 3.3 tonne liquid argon target. The largest background comes from $^{39}$Ar beta decays and is suppressed using pulseshape discrimination (PSD). We use two types of PSD algorithm: the prompt-fraction, which considers the fraction of the scintillation signal in a narrow and a wide time window around the event peak, and the log-likelihood-ratio, which compares the observed photon arrival times to a signal and a background model. We furthermore use two algorithms to determine the number of photons detected at a given time: (1) simply dividing the charge of each PMT pulse by the charge of a single photoelectron, and (2) a likelihood analysis that considers the probability to detect a certain number of photons at a given time, based on a model for the scintillation pulseshape and for afterpulsing in the light detectors. The prompt-fraction performs approximately as well as the log-likelihood-ratio PSD algorithm if the photon detection times are not biased by detector effects. We explain this result using a model for the information carried by scintillation photons as a function of the time when they are detected., Comment: 14 pages, 9 figures

We consider the size of the smallest set of vertices required to intersect every longest path in a chordal graph. Such sets are known as longest path transversals. We show that if $\omega(G)$ is the clique number of a chordal graph $G$, then there is a transversal of order at most $4\lceil\frac{\omega(G)}{5}\rceil$. We also consider the analogous question for longest cycles, and show that if $G$ is a 2-connected chordal graph then there is a transversal intersecting all longest cycles of order at most $2\lceil\frac{\omega(G)}{3}\rceil$., Comment: 11 pages

In online convex optimization (OCO), Lipschitz continuity of the functions is commonly assumed in order to obtain sublinear regret. Moreover, many algorithms have only logarithmic regret when these functions are also strongly convex. Recently, researchers from convex optimization proposed the notions of "relative Lipschitz continuity" and "relative strong convexity". Both of the notions are generalizations of their classical counterparts. It has been shown that subgradient methods in the relative setting have performance analogous to their performance in the classical setting. In this work, we consider OCO for relative Lipschitz and relative strongly convex functions. We extend the known regret bounds for classical OCO algorithms to the relative setting. Specifically, we show regret bounds for the follow the regularized leader algorithms and a variant of online mirror descent. Due to the generality of these methods, these results yield regret bounds for a wide variety of OCO algorithms. Furthermore, we further extend the results to algorithms with extra regularization such as regularized dual averaging., Comment: 22 pages, camera-ready version, accepted in NeurIPS 2020 for poster presentation. Second version has a new reference (acknowledged in the acknowledgements section) and comments on the relationship of this work with the paper

The physical mechanism driving mass ejection during a nova eruption is still poorly understood. Possibilities include ejection in a single ballistic event, a common envelope interaction, a continuous wind, or some combination of these processes. Here we present a study of 12 Galactic novae, for which we have pre-maximum high-resolution spectroscopy. All 12 novae show the same spectral evolution. Before optical peak, they show a slow P Cygni component. After peak a fast component quickly arises, while the slow absorption remains superimposed on top of it, implying the presence of at least two physically distinct flows. For novae with high-cadence monitoring, a third, intermediate-velocity component is also observed. These observations are consistent with a scenario where the slow component is associated with the initial ejection of the accreted material and the fast component with a radiation-driven wind from the white dwarf. When these flows interact, the slow flow is swept up by the fast flow, producing the intermediate component. These colliding flows may produce the gamma-ray emission observed in some novae. Our spectra also show that the transient heavy element absorption lines seen in some novae have the same velocity structure and evolution as the other lines in the spectrum, implying an association with the nova ejecta rather than a pre-existing circumbinary reservoir of gas or material ablated from the secondary. While this basic scenario appears to qualitatively reproduce multi-wavelength observations of classical novae, substantial theoretical and observational work is still needed to untangle the rich diversity of nova properties., Comment: 39 pages, 35 figures, submitted to ApJ

A classical nova is an eruption on the surface of a white dwarf in an accreting binary system. The material ejected from the white dwarf surface generally forms an axisymmetric shell. The shaping mechanisms of nova shells are probes of the processes that take place at energy scales between planetary nebulae and supernova remnants. We report on the discovery of nova shells surrounding the post-nova systems V4362 Sagittarii (1994) and more limited observations of DO Aquilae (1925). Distance measurements of 0.5p/m1.4 kpc for V4362 Sgr and 6.7 p/m 3.5 kpc -0.2 for DO Aql are found based on the expansion parallax method. The growth rates are measured to be 0.07``/year for DO Aql and 0.32``/year for V4362 Sgr. A preliminary investigation into the ionisation structure of the nova shell associated with V4362 Sgr is presented. The observed ionisation structure of nova shells depends strongly on their morphology and the orientation of the central component towards the observer. X-ray, IR and UV observations as well as optical integral field unit spectroscopy are required to better understand these interesting objects., Comment: 18 pages, 14 Figures, 9 Tables, accepted MNRAS

We present a study of the mass transfer and wind outflows of SS433, focusing on the so-called stationary lines based on archival high and low resolution optical spectra, and new optical multifilter polarimetry and low resolution optical spectra spanning an interval of a decade and a broad range of precessional and orbital phases. We derive $\text{E(B-V)}=0.86\pm0.10$ and revised UV and U band polarizations and polarization angles that yield the same position angle as the optical. The polarization wavelength dependence is consistent with optical-dominating electron scattering with a Rayleigh component in U and the UV filters; no polarization changes were observed during a flare event. Using profile orbital and precessional modulation of multiple lines we derive properties for the accretion disk, present evidence for a strong disk wind, determine its velocity structure, and demonstrate its variability on timescales unrelated to the orbit. We derive a mass ratio $q=0.37\pm0.04$, and masses $\text{M}_X=4.2\pm0.4\ \text{M}_\odot$, $\text{M}_A=11.3\pm 0.6\ \text{M}_\odot$, and show that the A star fills its Roche surface. The O I 7772 \r{A} and 8446 \r{A} lines show different but related orbital modulation and no evidence for a circumbinary disk component. Instead, the spectral line profile variability can be understood with an ionization stratified outflow predicted by thermal wind modeling, which also accounts for an extended equatorial structure detected at long wavelength., Comment: Typos corrected. Hyperref added. Accepted for publication in A&A (2020 May 27). Acknowledgements modified

A search has been performed for neutrinos from two sources, the $hep$ reaction in the solar $pp$ fusion chain and the $\nu_e$ component of the diffuse supernova neutrino background (DSNB), using the full dataset of the Sudbury Neutrino Observatory with a total exposure of 2.47 kton-years after fiducialization. The $hep$ search is performed using both a single-bin counting analysis and a likelihood fit. We find a best-fit flux that is compatible with solar model predictions while remaining consistent with zero flux, and set a one-sided upper limit of $\Phi_{hep} < 30\times10^{3}~\mathrm{cm}^{-2}~\mathrm{s}^{-1}$ [90% credible interval (CI)]. No events are observed in the DSNB search region, and we set an improved upper bound on the $\nu_e$ component of the DSNB flux of $\Phi^\mathrm{DSNB}_{\nu_e} < 19~\textrm{cm}^{-2}~\textrm{s}^{-1}$ (90% CI) in the energy range $22.9 < E_\nu < 36.9$~MeV., Comment: 11 pages, 6 figures

Online mirror descent (OMD) and dual averaging (DA) -- two fundamental algorithms for online convex optimization -- are known to have very similar (and sometimes identical) performance guarantees when used with a fixed learning rate. Under dynamic learning rates, however, OMD is provably inferior to DA and suffers a linear regret, even in common settings such as prediction with expert advice. We modify the OMD algorithm through a simple technique that we call stabilization. We give essentially the same abstract regret bound for OMD with stabilization and for DA by modifying the classical OMD convergence analysis in a careful and modular way that allows for straightforward and flexible proofs. Simple corollaries of these bounds show that OMD with stabilization and DA enjoy the same performance guarantees in many applications -- even under dynamic learning rates. We also shed light on the similarities between OMD and DA and show simple conditions under which stabilized-OMD and DA generate the same iterates., Comment: 27 pages main text, 37 pages in total, 1 figure. Version 2: Revision for camera-ready version of ICML 2020, with a new abstract, new discussion and acknowledgements sections, and some other minor modifications. Version 3: Technical report version of JMLR submission, with minor revisions, full proofs, and more details on the setting with composite functions

Research into the processes of photoionised nebulae plays a significant part in our understanding of stellar evolution. It is extremely difficult to visually represent or model ionised nebula, requiring astronomers to employ sophisticated modelling code to derive temperature, density and chemical composition. Existing codes are available that often require steep learning curves and produce models derived from mathematical functions. In this article we will introduce PyCross: PyCloudy Rendering Of Shape Software. This is a pseudo 3D modelling application that generates photoionisation models of optically thin nebulae, created using the Shape software. Currently PyCross has been used for novae and planetary nebulae, and it can be extended to Active Galactic Nuclei or any other type of photoionised axisymmetric nebulae. Functionality, an operational overview, and a scientific pipeline will be described with scenarios where PyCross has been adopted for novae (V5668 Sagittarii (2015) & V4362 Sagittarii (1994)) and a planetary nebula (LoTr1). Unlike the aforementioned photoionised codes this application does not require any coding experience, nor the need to derive complex mathematical models, instead utilising the select features from Cloudy/PyCloudy and Shape. The software was developed using a formal software development lifecycle, written in Python and will work without the need to install any development environments or additional python packages. This application, Shape models and PyCross archive examples are freely available to students, academics and research community on GitHub for download (https://github.com/karolfitzgerald/PyCross_OSX_App)., Comment: 15 pages, 12 figures

Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $\gamma$-rays from classical novae have hinted that shocks internal to the nova ejecta may dominate the nova emission. Shocks have also been suggested to power the luminosity of events as diverse as stellar mergers (Metzger & Pejcha 2017), supernovae (Moriya et al. 2018), and tidal disruption events (Roth et al. 2016), but observational confirmation has been lacking. Here we report simultaneous space-based optical and $\gamma$-ray observations of the 2018 nova V906 Carinae (ASASSN-18fv), revealing a remarkable series of distinct correlated flares in both bands. The optical and $\gamma$-ray flares occur simultaneously, implying a common origin in shocks. During the flares, the nova luminosity doubles, implying that the bulk of the luminosity is shock-powered. Furthermore, we detect concurrent but weak X-ray emission from deeply embedded shocks, confirming that the shock power does not appear in the X-ray band and supporting its emergence at longer wavelengths. Our data, spanning the spectrum from radio to $\gamma$-ray, provide direct evidence that shocks can power substantial luminosity in classical novae and other optical transients., Comment: 49 pages including supplementary information. Accepted and published in Nature Astronomy

The spectrum of cosmogenic neutrons at Earth's surface covers a wide energy range, from thermal to several GeV. The flux of secondary neutrons varies with latitude, elevation, solar activity, and nearby material, including ground moisture. We report the results from a campaign to measure count rates in neutron detectors responding to three different energy ranges conducted near the geomagnetic North Pole at CFS Alert, Nunavut, Canada (82.5 degrees , 62.5 degrees W; vertical geomagnetic cutoff rigidity, RC = 0 GV) in June of 2016. In November 2016, we performed a follow-on measurement campaign in southern Canada at similar RC (1.5 GV) and elevations. We conducted these measurements, at varying elevation and ground moisture content, with unmoderated and moderated 3He detectors for thermal and epithermal-to-MeV sensitivity, and with EJ-299-33 pulse shape discrimination plastic scintillator detectors for fast neutrons. Background gamma rays were monitored with NaI(Tl) detectors. Using these data sets, we compared the measured count rates to a predictive model. This is the first ever data set taken from this location on Earth. We find that for the thermal and epithermal-to-MeV neutron measurements the predictive model and data are in good agreement, except at one location on rock-covered ground near 1 km elevation. The discrepancy at that location may be attributable to ground moisture variability. Other measurements, during this campaign and prior, support the assertion that ground moisture plays a critical role in determining neutron flux., Comment: 28 pages, 9 figures

Ultra-cool white dwarfs are among the oldest stellar remnants in the Universe. Their efficient gravitational settling and low effective temperatures are responsible for the smooth spectra they exhibit. For that reason, it is not possible to derive their radial velocities or to find the chemistry of the progenitors. The best that can be done is to infer such properties from associated sources, which are coeval. The simplest form of such a system is a common proper motion pair where one star is an evolved degenerate and the other a main-sequence star. In this work, we present the discovery of the first of such a system, the M dwarf LHS 6328 and the ultra-cool white dwarf PSO J1801+625, from the Pan-STARRS 1 3$\pi$ survey and the $Gaia$ Data Release 2. Follow-up spectra were collected covering a usable wavelength range of 3750-24500 {\AA}. Their spectra show that the white dwarf has an effective temperature of 3550K and surface gravity of $\log g = 7.45 \pm 0.13$ or $\log g = 7.49 \pm 0.13$ for a CO or He core, respectively, when compared against synthetic spectra of ultra-cool white dwarf atmosphere models. The system has slightly subsolar metallicity with -0.25 < [Fe/H] < 0.0, and a spatial velocity of (U, V, W) = ($-114.26 \pm 0.24$, $222.94 \pm 0.60$, $10.25 \pm 0.34$) km s $^{-1}$ the first radial velocity and metallicity measurements of an ultra-cool white dwarf. This makes it the first and only benchmark of its kind to date., Comment: 11 pages, 4 figures, 5 tables. Some typographical changes

We consider the classical problem of prediction with expert advice. In the fixed-time setting, where the time horizon is known in advance, algorithms that achieve the optimal regret are known when there are two, three, or four experts or when the number of experts is large. Much less is known about the problem in the anytime setting, where the time horizon is not known in advance. No minimax optimal algorithm was previously known in the anytime setting, regardless of the number of experts. Even for the case of two experts, Luo and Schapire have left open the problem of determining the optimal algorithm. We design the first minimax optimal algorithm for minimizing regret in the anytime setting. We consider the case of two experts, and prove that the optimal regret is $\gamma \sqrt{t} / 2$ at all time steps $t$, where $\gamma$ is a natural constant that arose 35 years ago in studying fundamental properties of Brownian motion. The algorithm is designed by considering a continuous analogue of the regret problem, which is solved using ideas from stochastic calculus., Comment: 47 pages, 1 figure

DEAP-3600 is a liquid-argon scintillation detector looking for dark matter. Scintillation events in the liquid argon (LAr) are registered by 255 photomultiplier tubes (PMTs), and pulseshape discrimination (PSD) is used to suppress electromagnetic background events. The excellent PSD performance of LAr makes it a viable target for dark matter searches, and the LAr scintillation pulseshape discussed here is the basis of PSD. The observed pulseshape is a combination of LAr scintillation physics with detector effects. We present a model for the pulseshape of electromagnetic background events in the energy region of interest for dark matter searches. The model is composed of a) LAr scintillation physics, including the so-called intermediate component, b) the time response of the TPB wavelength shifter, including delayed TPB emission at $\mathcal O$(ms) time-scales, and c) PMT response. TPB is the wavelength shifter of choice in most LAr detectors. We find that approximately 10\% of the intensity of the wavelength-shifted light is in a long-lived state of TPB. This causes light from an event to spill into subsequent events to an extent not usually accounted for in the design and data analysis of LAr-based detectors.

Neutrons produced in nuclear interactions initiated by cosmic-ray muons present an irreducible background to many rare-event searches, even in detectors located deep underground. Models for the production of these neutrons have been tested against previous experimental data, but the extrapolation to deeper sites is not well understood. Here we report results from an analysis of cosmogenically produced neutrons at the Sudbury Neutrino Observatory. A specific set of observables are presented, which can be used to benchmark the validity of GEANT4 physics models. In addition, the cosmogenic neutron yield, in units of $10^{-4}\;\text{cm}^{2}/\left(\text{g}\cdot\mu\right)$, is measured to be $7.28 \pm 0.09\;\text{stat.} ^{+1.59}_{-1.12}\;\text{syst.}$ in pure heavy water and $7.30 \pm 0.07\;\text{stat.} ^{+1.40}_{-1.02}\;\text{syst.}$ in NaCl-loaded heavy water. These results provide unique insights into this potential background source for experiments at SNOLAB.

We consider stochastic gradient descent algorithms for minimizing a non-smooth, strongly-convex function. Several forms of this algorithm, including suffix averaging, are known to achieve the optimal $O(1/T)$ convergence rate in expectation. We consider a simple, non-uniform averaging strategy of Lacoste-Julien et al. (2011) and prove that it achieves the optimal $O(1/T)$ convergence rate with high probability. Our proof uses a recently developed generalization of Freedman's inequality. Finally, we compare several of these algorithms experimentally and show that this non-uniform averaging strategy outperforms many standard techniques, and with smaller variance.

Fast and accurate molecular force field (FF) parameterization is still an unsolved problem. Accurate FFs are not generally available for all molecules, like novel druglike molecules. While methods based on quantum mechanics (QM) exist to parameterize them with better accuracy, they are computationally expensive and slow, which limits applicability to a small number of molecules. Here, we present an automated FF parameterization method which can utilize either DFT calculations or approximate QM energies produced by different neural network potentials (NNPs), to obtain improved parameters for molecules. We demonstrate that for the case of torchani-ANI-1x NNP, we can parameterize small molecules in a fraction of time compared with an equivalent parameterization using DFT QM calculations while producing more accurate parameters than FF (GAFF2). We expect our method to be of critical importance in computational structure-based drug discovery. The current version is available at PlayMolecule (www.playmolecule.org) and implemented in HTMD, allowing to parameterize molecules with different QM and NNP options.

The DEAP-3600 experiment is searching for WIMP dark matter with a 3.3 tonne single phase liquid argon (LAr) target, located 2.1 km underground at SNOLAB. The experimental signature of dark matter interactions is keV-scale $^{40}$Ar nuclear recoils (NR) producing 128 nm LAr scintillation photons observed by PMTs. The largest backgrounds in DEAP-3600 are electronic recoils (ER) induced by $\beta$ and $\gamma$-rays originating from internal and external radioactivity in the detector material. A background model of the ER interactions in DEAP-3600 was developed and is described in this work. The model is based on several components which are expected from radioisotopes in the LAr, from ex-situ material assay measurements, and from dedicated independent in-situ analyses. This prior information is used in a Bayesian fit of the ER components to a 247.2 d dataset to model the radioactivity in the surrounding detector materials. While excellent discrimination between ERs and NRs is reached with pulse shape discrimination, utilizing the large difference between fast and slow components of LAr scintillation light, detailed knowledge of the ER background and activity of detector components, sets valuable constraints on other key types of backgrounds in the detector: neutrons and alphas. In addition, the activity of $^{42}$Ar in LAr in DEAP-3600 is determined by measuring the daughter decay of $^{42}$K. This cosmogenically activated trace isotope is a relevant background at higher energies for other rare event searches using atmospheric argon e.g. DarkSide-20k, GERDA or LEGEND. The specific activity of $^{42}$Ar in the atmosphere is found to be $40.4 \pm 5.9$ $\mu$Bq/kg of argon.

Neutron production in GeV-scale neutrino interactions is a poorly studied process. We have measured the neutron multiplicities in atmospheric neutrino interactions in the Sudbury Neutrino Observatory experiment and compared them to the prediction of a Monte Carlo simulation using GENIE and a minimally modified version of GEANT4. We analyzed 837 days of exposure corresponding to Phase I, using pure heavy water, and Phase II, using a mixture of Cl in heavy water. Neutrons produced in atmospheric neutrino interactions were identified with an efficiency of $15.3\%$ and $44.3\%$, for Phase I and II respectively. The neutron production is measured as a function of the visible energy of the neutrino interaction and, for charged current quasi-elastic interaction candidates, also as a function of the neutrino energy. This study is also performed classifying the complete sample into two pairs of event categories: charged current quasi-elastic and non charged current quasi-elastic, and $\nu_{\mu}$ and $\nu_e$. Results show good overall agreement between data and Monte Carlo for both phases, with some small tension with a statistical significance below $2\sigma$ for some intermediate energies.

We present a detailed study of the 2017 eruption of the classical nova ASASSN-17pf (LMCN 2017-11a), which is located in the Large Magellanic Cloud, including data from AAVSO, ASAS-SN, SALT, SMARTS, SOAR, and the Neil Gehrels \textit{Swift} Observatory. The optical light-curve is characterized by multiple maxima (flares) on top of a slowly evolving light-curve (with a decline time, $t_2>$ 100 d). The maxima correlate with the appearance of new absorption line systems in the optical spectra characterized by increasing radial velocities. We suggest that this is evidence of multiple episodes of mass-ejection with increasing expansion velocities. The line profiles in the optical spectra indicate very low expansion velocities (FWHM $\sim$ 190 km s$^{-1}$), making this nova one of the slowest expanding ever observed, consistent with the slowly evolving light-curve. The evolution of the colors and spectral energy distribution show evidence of decreasing temperatures and increasing effective radii for the pseudo-photosphere during each maximum. The optical and infrared light-curves are consistent with dust formation 125 days post-discovery. We speculate that novae showing several optical maxima have multiple mass-ejection episodes leading to shocks that may drive $\gamma$-ray emission and dust formation., Comment: Submitted to ApJ, 31 pages, 17 figures, and 7 tables

DEAP-3600 is a single-phase liquid argon (LAr) direct-detection dark matter experiment, operating 2 km underground at SNOLAB (Sudbury, Canada). The detector consists of 3279 kg of LAr contained in a spherical acrylic vessel. This paper reports on the analysis of a 758 tonne\cdot day exposure taken over a period of 231 live-days during the first year of operation. No candidate signal events are observed in the WIMP-search region of interest, which results in the leading limit on the WIMP-nucleon spin-independent cross section on a LAr target of $3.9\times10^{-45}$ cm$^{2}$ ($1.5\times10^{-44}$ cm$^{2}$) for a 100 GeV/c$^{2}$ (1 TeV/c$^{2}$) WIMP mass at 90\% C. L. In addition to a detailed background model, this analysis demonstrates the best pulse-shape discrimination in LAr at threshold, employs a Bayesian photoelectron-counting technique to improve the energy resolution and discrimination efficiency, and utilizes two position reconstruction algorithms based on PMT charge and photon arrival times.

Consider the problem of minimizing functions that are Lipschitz and strongly convex, but not necessarily differentiable. We prove that after $T$ steps of stochastic gradient descent, the error of the final iterate is $O(\log(T)/T)$ with high probability. We also construct a function from this class for which the error of the final iterate of deterministic gradient descent is $\Omega(\log(T)/T)$. This shows that the upper bound is tight and that, in this setting, the last iterate of stochastic gradient descent has the same general error rate (with high probability) as deterministic gradient descent. This resolves both open questions posed by Shamir (2012). An intermediate step of our analysis proves that the suffix averaging method achieves error $O(1/T)$ with high probability, which is optimal (for any first-order optimization method). This improves results of Rakhlin (2012) and Hazan and Kale (2014), both of which achieved error $O(1/T)$, but only in expectation, and achieved a high probability error bound of $O(\log \log(T)/T)$, which is suboptimal. We prove analogous results for functions that are Lipschitz and convex, but not necessarily strongly convex or differentiable. After $T$ steps of stochastic gradient descent, the error of the final iterate is $O(\log(T)/\sqrt{T})$ with high probability, and there exists a function for which the error of the final iterate of deterministic gradient descent is $\Omega(\log(T)/\sqrt{T})$.

The long baseline between the Earth and the Sun makes solar neutrinos an excellent test beam for exploring possible neutrino decay. The signature of such decay would be an energy-dependent distortion of the traditional survival probability which can be fit for using well-developed and high precision analysis methods. Here a model including neutrino decay is fit to all three phases of $^8$B solar neutrino data taken by the Sudbury Neutrino Observatory. This fit constrains the lifetime of neutrino mass state $\nu_2$ to be ${>8.08\times10^{-5}}$ s/eV at $90\%$ confidence. An analysis combining this SNO result with those from other solar neutrino experiments results in a combined limit for the lifetime of mass state $\nu_2$ of ${>1.04\times10^{-3}}$ s/eV at $99\%$ confidence.

Experimental tests of Lorentz symmetry in systems of all types are critical for ensuring that the basic assumptions of physics are well-founded. Data from all phases of the Sudbury Neutrino Observatory, a kiloton-scale heavy water Cherenkov detector, are analyzed for possible violations of Lorentz symmetry in the neutrino sector. Such violations would appear as one of eight possible signal types in the detector: six seasonal variations in the solar electron neutrino survival probability differing in energy and time dependence, and two shape changes to the oscillated solar neutrino energy spectrum. No evidence for such signals is observed, and limits on the size of such effects are established in the framework of the Standard Model Extension, including 40 limits on perviously unconstrained operators and improved limits on 15 additional operators. This makes limits on all minimal, Dirac-type Lorentz violating operators in the neutrino sector available for the first time.

By the construction of suitable graphs and the determination of their chromatic polynomials, we resolve two open questions concerning real chromatic roots. First we exhibit graphs for which the Beraha number $B_{10} = (5 + \sqrt{5})/2$ is a chromatic root. As it was previously known that no other non-integer Beraha number is a chromatic root, this completes the determination of precisely which Beraha numbers can be chromatic roots. Next we construct an infinite family of $3$-connected graphs such that for any $k \geqslant 1$, there is a member of the family with $q=2$ as a chromatic root of multiplicity at least $k$. The former resolves a question of Salas and Sokal [J. Statist. Pys. 104 (2001) pp. 609--699] and the latter a question of Dong and Koh [J. Graph Theory 70 (2012) pp. 262--283]., Comment: This first section of this paper was previously uploaded as arxiv:1611.01941

MapReduce has become the de facto standard model for designing distributed algorithms to process big data on a cluster. There has been considerable research on designing efficient MapReduce algorithms for clustering, graph optimization, and submodular optimization problems. We develop new techniques for designing greedy and local ratio algorithms in this setting. Our randomized local ratio technique gives $2$-approximations for weighted vertex cover and weighted matching, and an $f$-approximation for weighted set cover, all in a constant number of MapReduce rounds. Our randomized greedy technique gives algorithms for maximal independent set, maximal clique, and a $(1+\epsilon)\ln \Delta$-approximation for weighted set cover. We also give greedy algorithms for vertex colouring with $(1+o(1))\Delta$ colours and edge colouring with $(1+o(1))\Delta$ colours., Comment: 16 pages

Classical novae are eruptions on the surface of a white dwarf in a binary system. The material ejected from the white dwarf surface generally forms an axisymmetric shell of gas and dust around the system. The three-dimensional structure of these shells is difficult to untangle when viewed on the plane of the sky. In this work a geometrical model is developed to explain new observations of the 2015 nova V5668 Sagittarii. To understand the ionisation structure in terms of the nova shell morphology and estimate the emission distribution directly following the light-curve's dust-dip. High-cadence optical polarimetry and spectroscopy observations of a nova are presented. The ejecta is modelled in terms of morpho-kinematics and photoionisation structure. Initially observational results are presented, including broadband polarimetry and spectroscopy of V5668 Sgr nova during eruption. Variability over these observations provides clues towards the evolving structure of the nova shell. The position angle of the shell is derived from polarimetry, which is attributed to scattering from small dust grains. Shocks in the nova outflow are suggested in the photometry and the effect of these on the nova shell are illustrated with various physical diagnostics. Changes in density and temperature as the super soft source phase of the nova began are discussed. Gas densities are found to be of the order of 10$^{9}$ cm$^{-3}$ for the nova in its auroral phase. The blackbody temperature of the central stellar system is estimated to be around $2.2\times10^{5}$ K at times coincident with the super soft source turn-on. It was found that the blend around 4640 $\rm{\AA}$ commonly called `nitrogen flaring' is more naturally explained as flaring of the O~{\sc ii} multiplet (V1) from 4638 - 4696 $\rm{\AA}$, i.e. `oxygen flaring'., Comment: 17 pages, 13 figures, 3 tables, Accepted A&A 22 Nov 2017

The Dark matter Experiment using Argon Pulse-shape discrimination (DEAP) has been designed for a direct detection search for particle dark matter using a single-phase liquid argon target. The projected cross section sensitivity for DEAP-3600 to the spin-independent scattering of Weakly Interacting Massive Particles (WIMPs) on nucleons is $10^{-46}~\rm{cm}^{2}$ for a 100 GeV/$c^2$ WIMP mass with a fiducial exposure of 3 tonne-years. This paper describes the physical properties and construction of the DEAP-3600 detector., Comment: minor corrections in version 2. 70 pages, 24 figures, to be submitted to Astroparticle Physics Journal

The Lov\'asz Local Lemma is a seminal result in probabilistic combinatorics. It gives a sufficient condition on a probability space and a collection of events for the existence of an outcome that simultaneously avoids all of those events. Over the years, more general conditions have been discovered under which the conclusion of the lemma continues to hold. In this note we provide short proofs of two of those more general results: Shearer's lemma and the cluster expansion lemma, in their "lopsided" form. We conclude by using the cluster expansion lemma to prove that the symmetric form of the local lemma holds with probabilities bounded by $1/ed$, rather than the bound $1/e(d+1)$ required by the traditional proofs.

Molecular simulations produce very high-dimensional data-sets with millions of data points. As analysis methods are often unable to cope with so many dimensions, it is common to use dimensionality reduction and clustering methods to reach a reduced representation of the data. Yet these methods often fail to capture the most important features necessary for the construction of a Markov model. Here we demonstrate the results of various dimensionality reduction methods on two simulation data-sets, one of protein folding and another of protein-ligand binding. The methods tested include a k-means clustering variant, a non-linear auto encoder, principal component analysis and tICA. The dimension-reduced data is then used to estimate the implied timescales of the slowest process by a Markov state model analysis to assess the quality of the projection. The projected dimensions learned from the data are visualized to demonstrate which conformations the various methods choose to represent the molecular process., Comment: 11 pages, 10 figures

This paper reports the first results of a direct dark matter search with the DEAP-3600 single-phase liquid argon (LAr) detector. The experiment was performed 2 km underground at SNOLAB (Sudbury, Canada) utilizing a large target mass, with the LAr target contained in a spherical acrylic vessel of 3600 kg capacity. The LAr is viewed by an array of PMTs, which would register scintillation light produced by rare nuclear recoil signals induced by dark matter particle scattering. An analysis of 4.44 live days (fiducial exposure of 9.87 tonne-days) of data taken with the nearly full detector during the initial filling phase demonstrates the detector performance and the best electronic recoil rejection using pulse-shape discrimination in argon, with leakage $<1.2\times 10^{-7}$ (90% C.L.) between 16 and 33 keV$_{ee}$. No candidate signal events are observed, which results in the leading limit on WIMP-nucleon spin-independent cross section on argon, $<1.2\times 10^{-44}$ cm$^2$ for a 100 GeV/c$^2$ WIMP mass (90% C.L.)., Comment: Accepted for publication in Physical Review Letters

The Hamamatsu R5912-HQE photomultiplier-tube (PMT) is a novel high-quantum efficiency PMT. It is currently used in the DEAP-3600 dark matter detector and is of significant interest for future dark matter and neutrino experiments where high signal yields are needed. We report on the methods developed for in-situ characterization and monitoring of DEAP's 255 R5912-HQE PMTs. This includes a detailed discussion of typical measured single-photoelectron charge distributions, correlated noise (afterpulsing), dark noise, double, and late pulsing characteristics. The characterization is performed during the detector commissioning phase using laser light injected through a light diffusing sphere and during normal detector operation using LED light injected through optical fibres.

Tests on $B-L$ symmetry breaking models are important probes to search for new physics. One proposed model with $\Delta(B-L)=2$ involves the oscillations of a neutron to an antineutron. In this paper a new limit on this process is derived for the data acquired from all three operational phases of the Sudbury Neutrino Observatory experiment. The search was concentrated in oscillations occurring within the deuteron, and 23 events are observed against a background expectation of 30.5 events. These translate to a lower limit on the nuclear lifetime of $1.48\times 10^{31}$ years at 90% confidence level (CL) when no restriction is placed on the signal likelihood space (unbounded). Alternatively, a lower limit on the nuclear lifetime was found to be $1.18\times 10^{31}$ years at 90% CL when the signal was forced into a positive likelihood space (bounded). Values for the free oscillation time derived from various models are also provided in this article. This is the first search for neutron-antineutron oscillation with the deuteron as a target., Comment: 14 pages, 8 figures

We study an algorithm for approximating the multivariate independence polynomial $Z(\mathbf{z})$, with negative and complex arguments, an object that has strong connections to combinatorics and to statistical physics. In particular, the independence polynomial with negative arguments, $Z(-\mathbf{p})$, determines the Shearer region, the maximal region of probabilities to which the Lovasz Local Lemma (LLL) can be extended (Shearer 1985). In statistical physics, complex zeros of the independence polynomial relate to existence of phase transitions. Our main result is a deterministic algorithm to compute approximately the independence polynomial in any root-free complex polydisc centered at the origin. Our algorithm is essentially the same as Weitz's algorithm for positive parameters up to the tree uniqueness threshold, and the core of our analysis is a novel multivariate form of the correlation decay technique, which can handle non-uniform complex parameters. In particular, in the univariate real setting our work implies that Weitz's algorithm works in an interval between two critical points $(\lambda'_c(d), \lambda_c(d))$, and outside of this interval an approximation of $Z(\mathbf{z})$ is known to be NP-hard. As an application, we give a sub-exponential time algorithm for testing approximate membership in the Shearer region. We also give a new rounding based deterministic algorithm for Shearer's lemma (an extension of the LLL), which, however, runs in sub-exponential time. On the hardness side, we prove that evaluating $Z(\mathbf{z})$ at an arbitrary point in Shearer's region, and testing membership in Shearer's region, are #P-hard problems. We also establish the best possible dependence of the exponent of the run time of Weitz's correlation decay technique in the negative regime on the distance to the boundary of the Shearer region., Comment: 35 pages. Extended abstract to appear in Proceedings of ACM-SIAM SODA, 2018

Mader first proved that high average degree forces a given graph as a minor. Often motivated by Hadwiger's Conjecture, much research has focused on the average degree required to force a complete graph as a minor. Subsequently, various authors have consider the average degree required to force an arbitrary graph $H$ as a minor. Here, we strengthen (under certain conditions) a recent result by Reed and Wood, giving better bounds on the average degree required to force an $H$-minor when $H$ is a sparse graph with many high degree vertices. This solves an open problem of Reed and Wood, and also generalises (to within a constant factor) known results when $H$ is an unbalanced complete bipartite graph.

We generalise a result of Corr\'{a}di and Hajnal and show that every graph with average degree at least $\tfrac{4}{3}kr$ contains $k$ vertex disjoint cycles, each of order at least $r$, as long as $k \geq 6$. This bound is sharp when $r=3$., Comment: 15 pages

The treewidth of a graph is an important invariant in structural and algorithmic graph theory. This paper studies the treewidth of line graphs. We show that determining the treewidth of the line graph of a graph $G$ is equivalent to determining the minimum vertex congestion of an embedding of $G$ into a tree. Using this result, we prove sharp lower bounds in terms of both the minimum degree and average degree of $G$. These results are precise enough to exactly determine the treewidth of the line graph of a complete graph and other interesting examples. We also improve the best known upper bound on the treewidth of a line graph. Analogous results are proved for pathwidth., Comment: 18 pages (including appendices)