The All-Sky Automated Survey for Supernovae (ASAS-SN) provides long baseline (${\sim}4$ yrs) $V-$band light curves for sources brighter than V$\lesssim17$ mag across the whole sky. We produced V-band light curves for a total of ${\sim}61.5$ million sources and systematically searched these sources for variability. We identified ${\sim} 426,000$ variables, including ${\sim} 219,000$ new discoveries. Most (${\sim}74\%$) of our discoveries are in the Southern hemisphere. Here we use spectroscopic information from LAMOST, GALAH, RAVE, and APOGEE to study the physical and chemical properties of these variables. We find that metal-poor eclipsing binaries have orbital periods that are shorter than metal-rich systems at fixed temperature. We identified rotational variables on the main-sequence, red giant branch and the red clump. A substantial fraction (${\gtrsim}80\%$) of the rotating giants have large $v_{\rm rot}$ or large NUV excesses also indicative of fast rotation. The rotational variables have unusual abundances suggestive of analysis problems. Semi-regular variables tend to be lower metallicity ($\rm [Fe/H]{\sim}-0.5$) than most giant stars. We find that the APOGEE DR16 temperatures of oxygen-rich semi-regular variables are strongly correlated with the $W_{RP}-W_{JK}$ color index for $\rm T_{eff}\lesssim3800$ K. Using abundance measurements from APOGEE DR16, we find evidence for Mg, O and N enrichment in the semi-regular variables. We find that the Aluminum abundances of the semi-regular variables are strongly correlated with the pulsation period, where the variables with $\rm P\gtrsim 60$ days are significantly depleted in Al.

AT 2019abn was discovered in the nearby M51 galaxy by the Zwicky Transient Facility at more than two magnitudes and around three weeks prior to its optical peak. We aim to conduct a detailed photometric and spectroscopic follow-up campaign for AT 2019abn, with early discovery allowing for significant pre-maximum observations of an intermediate luminosity red transient (ILRT) for the first time. This work is based on the analysis of u'BVi'z'H photometry and low-resolution spectroscopy using the Liverpool Telescope, medium-resolution spectroscopy with the Gran Telescopio Canarias (GTC), and near-infrared imaging with the GTC and the Nordic Optical Telescope. We present the most detailed optical light curve of an ILRT to date, with multi-band photometry starting around three weeks before peak brightness. The transient peaked at an observed absolute magnitude of M_r' = -13.1, although it is subject to significant reddening from dust in M51, implying an intrinsic M_r' ~ -15.2. The initial light curve showed a linear, achromatic rise in magnitude before becoming bluer at peak. After peak brightness, the transient gradually cooled. This is reflected in our spectra, which at later times show absorption from such species as Fe I, Ni I and Li I. A spectrum taken around peak brightness shows narrow, low-velocity absorption lines, which we interpret as likely to originate from pre-existing circumstellar material. We conclude that while there are some peculiarities, such as the radius evolution, AT 2019abn fits in well overall with the ILRT class of objects and is the most luminous member of the class seen to date.

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.

With an estimated diameter in the 320 to 350 km range, (704) Interamnia is the fifth largest main belt asteroid and one of the few bodies that fills the gap in size between the four largest bodies with $D$ > 400 km (Ceres, Vesta, Pallas and Hygiea) and the numerous smaller bodies with $D$ $\lesssim$ 200 km. However, despite its large size, little is known about the shape and spin state of Interamnia and, therefore, about its bulk composition and past collisional evolution. We aimed to test at what size and mass the shape of a small body departs from a nearly ellipsoidal equilibrium shape (as observed in the case of the four largest asteroids) to an irregular shape as routinely observed in the case of smaller ($D$ $\lesssim$ 200 km) bodies. We observed Interamnia as part of our ESO VLT/SPHERE large program (ID: 199.C-0074) at thirteen different epochs. In addition, several new optical lightcurves were recorded. These data, along with stellar occultation data from the literature, were fed to the All-Data Asteroid Modeling (ADAM) algorithm to reconstruct the 3D-shape model of Interamnia and to determine its spin state. Interamnia's volume-equivalent diameter of 332 $\pm$ 6 km implies a bulk density of $\rho$=1.98 $\pm$ 0.68 gcm$^{-3}$ , which suggests that Interamnia - like Ceres and Hygiea - contains a high fraction of water ice, consistent with the paucity of apparent craters. Our observations reveal a shape that can be well approximated by an ellipsoid, and that is compatible with a fluid hydrostatic equilibrium at the 2 $\sigma$ level. The rather regular shape of Interamnia implies that the size and mass limit, under which the shapes of minor bodies with a high amount of water ice in the subsurface become irregular, has to be searched among smaller ($D$ $\lesssim$ 300km) less massive ($m$ $\lesssim$ 3x10$^{19}$ kg) bodies.

We characterize ${\sim} 71,200$ W UMa type (EW) contact binaries, including ${\sim} 12,600$ new discoveries, using ASAS-SN $V$-band all-sky light curves along with archival data from Gaia, 2MASS, AllWISE, LAMOST, GALAH, RAVE, and APOGEE. There is a clean break in the EW period-luminosity relation at $\rm \log (\rm P/d){\simeq}-0.30$, separating the longer period early-type (A sub-type) EW binaries from the shorter period, late-type (W sub-type) systems. The two populations are even more cleanly separated in the space of period and effective temperature, by $\rm T_{eff}=6710\,K-1760\,K\,\log(P/0.5\,d)$. Early-type and late-type EW binaries follow opposite trends in $\rm T_{eff}$ with orbital period. For longer periods, early-type EW binaries are cooler, while late-type systems are hotter. We derive period-luminosity relationships (PLRs) in the $W_{JK}$, $V$, Gaia DR2 $G$, $J$, $H$, $K_s$ and $W_1$ bands for the late-type and early-type EW binaries separated both by period and effective temperature, and by period alone. The dichotomy of contact binaries is almost certainly related to the Kraft break and the related changes in envelope structure, winds and angular momentum loss.

We characterize an all-sky catalog of ${\sim} 8,400$ $\delta$ Scuti variables in ASAS-SN, which includes ${\sim} 3,300$ new discoveries. Using distances from Gaia DR2, we derive period-luminosity relationships (PLRs) for both the fundamental mode and overtone pulsators in the $W_{JK}$, $V$, Gaia DR2 $G$, $J$, $H$, $K_s$ and $W_1$ bands. We find that the overtone pulsators have a dominant overtone mode, with many sources pulsating in the second overtone or higher order modes. The fundamental mode pulsators have metallicity dependent periods, with $\log_{10} (\rm P){\sim}-1.1$ for $\rm [Fe/H]<-0.3$ and $\log_{10} (\rm P){\sim}-0.9$ for $\rm [Fe/H]>0$ which leads to a period-dependant scale height. Stars with $\rm P>0.100\,d$ are predominantly located close to the Galactic disk ($\rm |Z|<0.5\,kpc$). The median period at a scale height of $\rm Z{\sim}0\, kpc$ also increases with the Galactocentric radius $\rm R$, from $\log_{10} (\rm P){\sim}-0.94$ for sources with $\rm R>9\, kpc$ to $\log_{10} (\rm P){\sim}-0.85$ for sources with $\rm R<7\, kpc$, which is indicative of a radial metallicity gradient. To illustrate potential applications of this all-sky catalog, we obtained 30 min cadence, image subtraction TESS light curves for a sample of 10 fundamental mode and 10 overtone $\delta$ Scuti stars discovered by ASAS-SN. From this sample, we identified two new $\delta$ Scuti eclipsing binaries, ASASSN-V J071855.62$-$434247.3 and ASASSN-V J170344.20$-$615941.2 with short orbital periods of $\rm P_{orb}=2.6096$ d and $\rm P_{orb}=2.5347$ d respectively.

Recent studies have indicated that triple star systems may play a role in the formation of an appreciable number of planetary nebulae, however only one triple central star is known to date (and that system is likely too wide to have had much influence on the evolution of its component stars). Here, we consider the possibility that Sh 2-71 was formed by a triple system which has since broken apart. We present the discovery of two regions of emission, seemingly aligned with the proposed tertiary orbit (i.e. in line with the axis formed by the two candidate central star systems previously considered in the literature). We also perform a few simple tests of the plausibility of the triple hypothesis based on the observed properties (coordinates, radial velocities, distances and proper motions) of the stars observed close to the projected centre of the nebula, adding further support through numerical integrations of binary orbits responding to mass loss. Although a number of open questions remain, we conclude that Sh 2-71 is currently one of the best candidates for planetary nebula formation influenced by triple-star interactions.

We investigate kinematics of mass loss from the vicinity of the second Lagrange point L2 with applications to merging binary stars, common envelope evolution and the associated transient brightenings. For ballistic particle trajectories, we characterize initial velocities and positional offsets from L2 which lead to unbound outflow, fall back followed by a formation of a decretion disk, collision with the binary surface, or a hydrodynamic shock close to the binary, where some particle trajectories loop and self-intersect. The latter two final states occur only when the trajectories are initiated with offset from L2 or with velocity vector different from corotation with the binary. We find that competition between the time-dependent and steeply radially decreasing tidal torques from the binary, Coriolis force and initial conditions lead to a non-trivial distribution of outcomes in the vicinity of L2. Specifically, even for initial velocities slower than corotation, we find that a set of initial position offsets lead to unbound outflows. Our results will aid in the interpretation of the morphology of mass loss streams in hydrodynamic simulations.

The All-Sky Automated Survey for Supernovae (ASAS-SN) provides long baseline (${\sim}4$ yrs) light curves for sources brighter than V$\lesssim17$ mag across the whole sky. As part of our effort to characterize the variability of all the stellar sources visible in ASAS-SN, we have produced ${\sim}30.1$ million V-band light curves for sources in the southern hemisphere using the APASS DR9 catalog as our input source list. We have systematically searched these sources for variability using a pipeline based on random forest classifiers. We have identified ${\sim} 220,000$ variables, including ${\sim} 88,300$ new discoveries. In particular, we have discovered ${\sim}48,000$ red pulsating variables, ${\sim}23,000$ eclipsing binaries, ${\sim}2,200$ $\delta$-Scuti variables and ${\sim}10,200$ rotational variables. The light curves and characteristics of the variables are all available through the ASAS-SN variable stars database (https://asas-sn.osu.edu/variables). The pre-computed ASAS-SN V-band light curves for all the ${\sim}30.1$ million sources are available through the ASAS-SN photometry database (https://asas-sn.osu.edu/photometry). This effort will be extended to provide ASAS-SN light curves for sources in the northern hemisphere and for V$\lesssim17$ mag sources across the whole sky that are not included in APASS DR9.

We explore the synergy between photometric and spectroscopic surveys by searching for periodic variable stars among the targets observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) using photometry from the All-Sky Automated Survey for Supernovae (ASAS-SN). We identified 1924 periodic variables among more than 258000 APOGEE targets; 465 are new discoveries. We homogeneously classified 430 eclipsing and ellipsoidal binaries, 139 classical pulsators (Cepheids, RR Lyrae and delta Scuti), 719 long period variables (pulsating red giants) and 636 rotational variables. The search was performed using both visual inspection and machine learning techniques. The light curves were also modeled with the damped random walk stochastic process. We find that the median [Fe/H] of variable objects is lower by 0.3 dex than that of the overall APOGEE sample. Eclipsing binaries and ellipsoidal variables are shifted to a lower median [Fe/H] by 0.2 dex. Eclipsing binaries and rotational variables exhibit significantly broader spectral lines than the rest of the sample. We make ASAS-SN light curves for all the APOGEE stars publicly available and provide parameters for the variable objects.

The nearby SN 2017eaw is a Type II-P (``plateau') supernova showing early-time, moderate CSM interaction. We present a comprehensive study of this SN including the analysis of high-quality optical photometry and spectroscopy covering the very early epochs up to the nebular phase, as well as near-UV and near-infrared spectra, and early-time X-ray and radio data. The combined data of SNe 2017eaw and 2004et allow us to get an improved distance to the host galaxy, NGC 6946, as $D \sim 6.85$ $\pm 0.63$ Mpc; this fits in recent independent results on the distance of the host and disfavors the previously derived (30% shorter) distances based on SN 2004et. From modeling the nebular spectra and the quasi-bolometric light curve, we estimate the progenitor mass and some basic physical parameters for the explosion and the ejecta. Our results agree well with previous reports on a RSG progenitor star with a mass of $\sim15-16$ M$_\odot$. Our estimation on the pre-explosion mass-loss rate ($\dot{M} \sim3 \times 10^{-7} -$ $1\times 10^{-6} M_{\odot}$ yr$^{-1}$) agrees well with previous results based on the opacity of the dust shell enshrouding the progenitor, but it is orders of magnitude lower than previous estimates based on general light-curve modeling of Type II-P SNe. Combining late-time optical and mid-infrared data, a clear excess at 4.5 $\mu$m can be seen, supporting the previous statements on the (moderate) dust formation in the vicinity of SN 2017eaw.

The mid-infrared (mid-IR) wavelength regime offers several advantages for following the late-time evolution of supernovae (SNe). First, the peaks of the SN spectral energy distributions shift toward longer wavelengths following the photospheric phase. Second, mid-IR observations suffer less from effects of interstellar extinction. Third, and perhaps most important, the mid-IR traces dust formation and circumstellar interaction at late-times (>100 days) after the radioactive ejecta component fades. The Spitzer Space Telescope has provided substantial mid-IR observations of SNe since its launch in 2003. More than 200 SNe have been targeted, but there are even more SNe that have been observed serendipitously. Here we present the results of a comprehensive study based on archival Spitzer/IRAC images of more than 1100 SN positions; from this sample, 119 SNe of various subclasses have been detected, including 45 SNe with previously unpublished mid-IR photometry. The photometry reveal significant amounts of warm dust in some cases. We perform an in-depth analysis to constrain the origin and heating mechanism of the dust, and present the resulting statistics.

Shock revival in core-collapse supernovae (CCSNe) may be due to the neutrino mechanism. While it is known that in a neutrino-powered CCSN, explosion begins when the neutrino luminosity of the proto-neutron star exceeds a critical value, the physics of this condition in time-dependent, multidimensional simulations are not fully understood. \citet{Pejcha2012} found that an `antesonic condition' exists for time-steady spherically symmetric models, potentially giving a physical explanation for the critical curve observed in simulations. In this paper, we extend that analysis to time-dependent, spherically symmetric polytropic models. We verify the critical antesonic condition in our simulations, showing that models exceeding it drive transonic winds whereas models below it exhibit steady accretion. In addition, we find that (1) high spatial resolution is needed for accurate determination of the antesonic ratio and shock radius at the critical curve, and that low resolution simulations systematically underpredict these quantities, making explosion more difficult at lower resolution; (2) there is an important physical connection between the critical mass accretion rate at explosion and the mass loss rate of the post-explosion wind: the two are directly proportional at criticality, implying that, at criticality, the wind kinetic power is tied directly to the accretion power; (3) the value of the post-shock adiabatic index $\Gamma$ has a large effect on the length and time scales of the post-bounce evolution of the explosion larger values of $\Gamma$ result in a longer transition from the accretion to wind phases.

The All-Sky Automated Survey for Supernovae (ASAS-SN) provides long baseline (${\sim}4$ yrs) light curves for sources brighter than V$\lesssim17$ mag across the whole sky. The Transiting Exoplanet Survey Satellite (TESS) has started to produce high-quality light curves with a baseline of at least 27 days, eventually for most of the sky. The combination of ASAS-SN and TESS light curves probes both long and short term variability in great detail, especially towards the TESS continuous viewing zones (CVZ) at the ecliptic poles. We have produced ${\sim}1.3$ million V-band light curves covering a total of ${\sim}1000 \, \rm deg^2$ towards the southern TESS CVZ and have systematically searched these sources for variability. We have identified ${\sim} 11,700$ variables, including ${\sim} 7,000$ new discoveries. The light curves and characteristics of the variables are all available through the ASAS-SN variable stars database (https://asas-sn.osu.edu/variables). We also introduce an online resource to obtain pre-computed ASAS-SN V-band light curves (https://asas-sn.osu.edu/photometry) starting with the light curves of the ${\sim}1.3$ million sources studied in this work. This effort will be extended to provide ASAS-SN light curves for ${\sim}50\;$million sources over the entire sky.

We report the identification of the bright (V${\sim}13.3$ mag) star FY Sct as a long period detached eclipsing binary using a combined ASAS-SN and ASAS light curve spanning 2000-2018. The orbital period is P${\sim}2.57$ years and the primary eclipse lasts ${\sim}73$ days. The eclipse profile is suggestive of a disk eclipsing binary rather than a stellar component. We also detect ${\sim}0.4$ mag pulsations with a period of P${\sim78}$ d. The next eclipse begins on September, 28, 2018. Further photometric and spectroscopic observations are encouraged, particularly when the system is in eclipse.

The variable stars in the VSX catalog are derived from a multitude of inhomogeneous data sources and classification tools. This inhomogeneity complicates our understanding of variable star types, statistics, and properties, and it directly affects attempts to build training sets for current (and next) generation all-sky, time-domain surveys. We homogeneously analyze the ASAS-SN V-band light curves of ${\sim}412,000$ variables from the VSX catalog. The variables are classified using an updated random forest classifier with an $F_1$ score of 99.4% and refinement criteria for individual classifications. We have derived periods for ${\sim}44,000$ variables in the VSX catalog that lack a period, and have reclassified ${\sim} 17,000$ sources into new broad variability groups with high confidence. We have also reclassified ${\sim} 94,000$ known variables with miscellaneous/generic classifications. The light curves, classifications, and a range of properties of the variables are all available through the ASAS-SN variable stars database (https://asas-sn.osu.edu/variables). We also include the V-band light curves for a set of ${\sim}4,000$ rare variables and transient sources, including cataclysmic variables, symbiotic binaries and flare stars.

Binary stars commonly pass through phases of direct interaction which result in the rapid loss of mass, energy, and angular momentum. Though crucial to understanding the fates of these systems, including their potential as gravitational wave sources, this short-lived phase is poorly understood and has thus far been unambiguously observed in only a single event, V1309 Sco. Here we show that the complex and previously-unexplained photometric behavior of V1309 Sco prior to its main outburst results naturally from the runaway loss of mass and angular momentum from the outer Lagrange point, which lasts for thousands of orbits prior to the final dynamical coalescence, much longer than predicted by contemporary models. This process enshrouds the binary in a "death spiral" outflow, which affects the amplitude and phase modulation of its light curve, and contributes to driving the system together. The total amount of mass lost during this gradual phase ($\sim 0.05 M_\odot$) rivals the mass lost during the subsequent dynamical interaction phase, which has been the main focus of "common envelope" modeling so far. Analogous features in related transients suggest that this behavior is ubiquitous.

We present the first release of the Czech Variable star catalogue that currently contains 1228 stars whose variability was discovered by 60 Czech observers. The catalogue contains confirmed variable stars of various types, but also candidates. We give precise coordinates, cross identification with other catalogues, information about constellation, variability type, brightness, light elements, name of the discoverer and year of discovery. In eighty-eight percent of stars the variability type is estimated, for more than 60 % of the stars the light ephemerides are given.

Core-collapse supernova explosions expose the structure and environment of massive stars at the moment of their death. We use the global fitting technique of Pejcha & Prieto (2015a,b) to estimate a set of physical parameters of 19 normal Type II SNe, such as their distance moduli, reddenings, $^{56}$Ni masses $M_{\rm Ni}$, and explosion energies $E_{\rm exp}$ from multicolor light curves and photospheric velocity curves. We confirm and characterize known correlations between $M_{\rm Ni}$ and bolometric luminosity at 50 days after the explosion, and between $M_{\rm Ni}$ and $E_{\rm exp}$. We pay special attention to the observed distribution of $M_{\rm Ni}$ coming from a joint sample of 38 Type~II SNe, which can be described as a skewed-Gaussian-like distribution between $0.005~M_{\odot}$ and $0.280~M_{\odot}$, with a median of $0.031~M_{\odot}$, mean of $0.046~M_{\odot}$, standard deviation of $0.048~M_{\odot}$ and skewness of $3.050$. We use two-sample Kolmogorov-Smirnov test and two-sample Anderson-Darling test to compare the observed distribution of $M_{\rm Ni}$ to results from theoretical hydrodynamical codes of core-collapse explosions with the neutrino mechanism presented in the literature. Our results show that the theoretical distributions obtained from the codes tested in this work, KEPLER and Prometheus Hot Bubble, are compatible with the observations irrespective of different pre-supernova calibrations and different maximum mass of the progenitors.

Luminous red novae (LRN) are a class of optical transients believed to originate from the mergers of binary stars, or "common envelope" events. Their light curves often show secondary maxima, which cannot be explained in the previous models of thermal energy diffusion or hydrogen recombination without invoking multiple independent shell ejections. We propose that double-peaked light curves are a natural consequence of a collision between dynamically-ejected fast shell and pre-existing equatorially-focused material, which was shed from the binary over many orbits preceding the dynamical event. The fast shell expands freely in the polar directions, powering the initial optical peak through cooling envelope emission. Radiative shocks from the collision in the equatorial plane power the secondary light curve peak on the radiative diffusion timescale of the deeper layers, similar to luminous Type IIn supernovae and some classical novae. Using a detailed 1D analytic model, informed by complementary 3D hydrodynamical simulations, we show that shock-powered emission can explain the observed range of peak timescales and luminosities of the secondary peaks in LRN for realistic variations in the binary parameters and fraction of the binary mass ejected. The dense shell created by the radiative shocks in the equatorial plane provides an ideal location for dust nucleation consistent with the the inferred aspherical geometry of dust in LRN. For giant stars, the ejecta forms dust when the shock-powered luminosity is still high, which could explain the infrared transients recently discovered by Spitzer. Our results suggest that pre-dynamical mass loss is common if not ubiquitous in stellar mergers, providing insight into the instabilities responsible for driving the binary merger.

We present photometric and spectroscopic observations of ASASSN-13co, an unusually luminous Type II supernova and the first core-collapse supernova discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN). First detection of the supernova was on UT 2013 August 29 and the data presented span roughly 3.5 months after discovery. We use the recently developed model from Pejcha & Prieto (2015) to model the multi-band light curves of ASASSN-13co and derive the bolometric luminosity curve. We compare ASASSN-13co to other Type II supernovae to show that it was unusually luminous for a Type II supernova and that it exhibited an atypical light curve shape that does not cleanly match that of either a standard Type II-L or Type II-P supernova.

We study mass loss from the outer Lagrange point (L2) in binary stellar mergers and their luminous transients by means of radiative hydrodynamical simulations. Previously, we showed that for binary mass ratios 0.06 < q < 0.8, synchronous L2 mass loss results in a radiatively inefficient, dust-forming unbound equatorial outflow. A similar outflow exists irrespective of q if the ratio of the sound speed to the orbital speed at the injection point is sufficiently large, \epsilon = c_T/v_orb > 0.15. By contrast, for cold L2 mass-loss (\epsilon < 0.15) from binaries with q < 0.06 or q > 0.8, the equatorial outflow instead remains marginally-bound and falls back to the binary over tens to hundreds of binary orbits, where it experiences additional tidal torqueing and shocking. As the bound gas becomes virialized with the binary, the luminosity of the system increases slowly at approximately constant photosphere radius, causing the temperature to rise. Subsequent evolution depends on the efficiency of radiative cooling. If the bound atmosphere is able to cool efficiently, as quantified by radiative diffusion time being shorter than the advection time (t_diff/t_adv < 1), then the virialized gas collapses to an excretion disk, while for t_diff/t_adv > 1 an isotropic wind is formed. Between these two extremes, an inflated envelope transports the heat generated near the binary to the surface by meridional flows. In all cases, the radiated luminosity reaches a fraction ~0.01 to 0.1 of Mdot v_orb^2/2, where Mdot is the mass outflow rate. We discuss the implications of our results for transients in the luminosity gap between classical novae and supernovae, such as V1309 Sco and V838 Mon.

We study transients produced by equatorial disk-like outflows from catastrophically mass-losing binary stars with an asymptotic velocity and energy deposition rate near the inner edge which are proportional to the binary escape velocity v_esc. As a test case, we present the first smoothed-particle radiation-hydrodynamics calculations of the mass loss from the outer Lagrange point with realistic equation of state and opacities. The resulting spiral stream becomes unbound for binary mass ratios 0.06 < q < 0.8. For synchronous binaries with non-degenerate components, the spiral-stream arms merge at a radius of ~10a, where a is the binary semi-major axis, and the accompanying shock thermalizes about 10% of the kinetic power of the outflow. The mass-losing binary outflows produce luminosities reaching up to ~10^6 L_Sun and effective temperatures spanning 500 < T_eff < 6000 K, which is compatible with many of the class of recently-discovered red transients such as V838 Mon and V1309 Sco. Dust readily forms in the outflow, potentially in a catastrophic global cooling transition. The appearance of the transient is viewing angle-dependent due to vastly different optical depths parallel and perpendicular to the binary plane. We predict a correlation between the peak luminosity and the outflow velocity, which is roughly obeyed by the known red transients. Outflows from mass-losing binaries can produce luminous (10^5 L_Sun) and cool (T_eff < 1500 K) transients lasting a year or longer, as has potentially been detected by Spitzer surveys of nearby galaxies.

Hydrogen-rich Type II-Plateau supernovae exhibit correlations between the plateau luminosity $L_{\rm pl}$, the nickel mass $M_{\rm Ni}$, the explosion energy $E_{\rm exp}$, and the ejecta mass $M_{\rm ej}$. Using our global, self-consistent, multi-band model of nearby well-observed supernovae, we find that the covariances of these quantities are strong and that the confidence ellipsoids are oriented in the direction of the correlations, which reduces their significance. By proper treatment of the covariance matrix of the model, we discover a significant intrinsic width to the correlations between $L_{\rm pl}$, $E_{\rm exp}$, and $M_{\rm Ni}$, where the uncertainties due to the distance and the extinction dominate. For fixed $E_{\rm exp}$, the spread in $M_{\rm Ni}$ is about 0.25 dex, which we attribute to the differences in the progenitor internal structure. We argue that the effects of incomplete $\gamma$-ray trapping are not important in our sample. Similarly, the physics of the Type II-Plateau supernova light curves leads to inherently degenerate estimates of $E_{\rm exp}$ and $M_{\rm ej}$, which makes their observed correlation weak. Ignoring the covariances of supernova parameters or the intrinsic width of the correlations causes significant biases in the slopes of the fitted relations. Our results imply that Type II-Plateau supernova explosions are not described by a single physical parameter or a simple one-dimensional trajectory through the parameter space, but instead reflect the diversity of the core and surface properties of their progenitors. We discuss the implications for the physics of the explosion mechanism and possible future observational constraints.

If the neutrino luminosity from the proto-neutron star formed during a massive star core collapse exceeds a critical threshold, a supernova (SN) results. Using spherical quasi-static evolutionary sequences for hundreds of progenitors over a range of metallicities, we study how the explosion threshold maps onto observables, including the fraction of successful explosions, the neutron star (NS) and black hole (BH) mass functions, the explosion energies (E_SN) and nickel yields (M_Ni), and their mutual correlations. Successful explosions are intertwined with failures in a complex pattern that is not simply related to initial progenitor mass or compactness. We predict that progenitors with initial masses of 15 +/- 1, 19 +/- 1, and 21-26 M_Sun are most likely to form BHs, that the BH formation probability is non-zero at solar-metallicity and increases significantly at low metallicity, and that low luminosity, low Ni-yield SNe come from progenitors close to success/failure interfaces. We qualitatively reproduce the observed E_SN-M_Ni correlation, we predict a correlation between the mean and width of the NS mass and E_SN distributions, and that the means of the NS and BH mass distributions are correlated. We show that the observed mean NS mass of ~1.33 M_Sun implies that the successful explosion fraction is higher than 0.35. Overall, we show that the neutrino mechanism can in principle explain the observed properties of SNe and their compact objects. We argue that the rugged landscape of progenitors and outcomes mandates that SN theory should focus on reproducing the wide ranging distributions of observed SN properties.

We present a new self-consistent and versatile method that derives photospheric radius and temperature variations of Type II-Plateau supernovae based on their expansion velocities and photometric measurements. We apply the method to a sample of 26 well-observed, nearby supernovae with published light curves and velocities. We simultaneously fit ~230 velocity and ~6800 magnitude measurements distributed over 21 photometric passbands spanning wavelengths from 0.19 to 2.2 microns. The light curve differences among the Type II-Plateau supernovae are well-modeled by assuming different rates of photospheric radius expansion, which we explain as different density profiles of the ejecta and we argue that steeper density profiles result in flatter plateaus, if everything else remains unchanged. The steep luminosity decline of Type II-Linear supernovae is due to fast evolution of the photospheric temperature, which we verify with a successful fit of SN1980K. Eliminating the need for theoretical supernova atmosphere models, we obtain self-consistent relative distances, reddenings, and nickel masses fully accounting for all internal model uncertainties and covariances. We use our global fit to estimate the time evolution of any missing band tailored specifically for each supernova and we construct spectral energy distributions and bolometric light curves. We produce bolometric corrections for all filter combinations in our sample. We compare our model to the theoretical dilution factors and find good agreement for the B and V filters. Our results differ from the theory when the I, J, H, or K bands are included. We investigate the reddening law towards our supernovae and find reasonable agreement with standard R_V ~ 3.1 reddening law in UBVRI bands. Results for other bands are inconclusive. We make our fitting code publicly available.

A large reverberation mapping study of the Seyfert 1 galaxy NGC 7469 has yielded emission-line lags for Hbeta 4861 and He II 4686 and a central black hole mass measurement of about 10 million solar masses, consistent with previous measurements. A very low level of variability during the monitoring campaign precluded meeting our original goal of recovering velocity-delay maps from the data, but with the new Hbeta measurement, NGC 7469 is no longer an outlier in the relationship between the size of the Hbeta-emitting broad-line region and the AGN luminosity. It was necessary to detrend the continuum and Hbeta and He II 4686 line light curves and those from archival UV data for different time-series analysis methods to yield consistent results.

V1309 Sco was proposed to be a stellar merger and a common envelope transient based on the pre-outburst light curve of a contact eclipsing binary with a rapidly decaying orbital period. Using published data, I show that the period decay timescale P/Pdot of V1309 Sco decreased from ~1000 to ~170 years in less than about 6 years, which implies a very high value of second period derivative. I argue that V1309 Sco experienced an onset of dynamical mass loss through the outer Lagrange point, which eventually obscured the binary. The photosphere of the resulting continuous optically-thick outflow expands as the mass-loss rate increases, explaining the ~200 day rise to optical maximum. The model yields the mass-loss rate of the binary star as a function of time and fits the observed light curve remarkably well. It is also possible to observationally constrain the properties of the surface layers undergoing the dynamical mass loss. V1309 Sco is thus a prototype of a new class of stellar transients distinguished by a slow rise to optical maximum that are driven by dynamical mass loss from a binary. I discuss implications of these findings for stellar transients and other suggested common envelope events.

We study the orbital evolution of hierarchical quadruple systems composed of two binaries on a long mutual orbit, where each binary acts as a Kozai-Lidov (KL) perturber on the other. We find that the coupling between the two binaries qualitatively changes the behavior of their KL cycles. The binaries can experience coherent eccentricity oscillations as well as excursions to very high eccentricity that occur over a much larger fraction of the parameter space than in triple systems. For a ratio of outer to inner semi-major axes of 10 to 20, about 30 to 50% of equal-mass quadruples reach eccentricity 1-e < 10^{-3} in one of the binaries. This is about 4 to 12 times more than for triples with equivalent parameters. Orbital "flips" and collisions without previous tidal interaction are similarly enhanced in quadruples relative to triples. We argue that the frequency of evolutionary paths influenced by KL cycles is comparable in the triple and quadruple populations even though field quadruples are a factor of ~5 less frequent than triples. Additionally, quadruples might be a non-negligible source of triples and provide fundamentally new evolutionary outcomes involving close binaries, mergers, collisions, and associated transients, which occur without any fine tuning of parameters. Finally, we study the perturbations to a planetary orbit due to a distant binary and we find that the fraction of orbital flips is a factor of 3 to 4 higher than for the corresponding triple system given our fiducial parameters with implications for hot Jupiters and star-planet collisions.

We present velocity-resolved reverberation results for five active galactic nuclei. We recovered velocity-delay maps using the maximum-entropy method for four objects: Mrk 335, Mrk 1501, 3C120, and PG2130+099. For the fifth, Mrk 6, we were only able to measure mean time delays in different velocity bins of the Hbeta emission line. The four velocity-delay maps show unique dynamical signatures for each object. For 3C120, the Balmer lines show kinematic signatures consistent with both an inclined disk and infalling gas, but the HeII 4686 emission line is suggestive only of inflow. The Balmer lines in Mrk 335, Mrk 1501, and PG 2130+099 show signs of infalling gas, but the HeII emission in Mrk 335 is consistent with an inclined disk. We also see tentative evidence of combined virial motion and infalling gas from the velocity-binned analysis of Mrk 6. The maps for 3C120 and Mrk 335 are two of the most clearly defined velocity-delay maps to date. These maps constitute a large increase in the number of objects for which we have resolved velocity-delay maps and provide evidence supporting the reliability of reverberation-based black hole mass measurements.

We present the results from a detailed analysis of photometric and spectrophotometric data on five Seyfert 1 galaxies observed as a part of a recent reverberation mapping program. The data were collected at several observatories over a 140-day span beginning in 2010 August and ending in 2011 January. We obtained high sampling-rate light curves for Mrk 335, Mrk 1501, 3C120, Mrk 6, and PG2130+099, from which we have measured the time lag between variations in the 5100 Angstrom continuum and the H-beta broad emission line. We then used these measurements to calculate the mass of the supermassive black hole at the center of each of these galaxies. Our new measurements substantially improve previous measurements of MBH and the size of the broad line-emitting region for four sources and add a measurement for one new object. Our new measurements are consistent with photoionization physics regulating the location of the broad line region in active galactic nuclei.

The evolution of multiple stellar systems can be driven by Kozai cycles and tidal friction (KCTF), which shrink the orbit of the inner binary. There is an interesting possibility that two close binaries on a common long-period orbit experience mutually-induced KCTF. We present the discovery of a possible new quadruple system composed of two unresolved eclipsing binaries (EBs), CzeV343 (V~13.5 mag). We obtained photometric observations of CzeV343 that completely cover the two orbital periods and we successfully model the light curves as the sum of two detached EBs. We provide confidence intervals for the model parameters and minima timings by bootstrap resampling of our data. One of the EBs shows a distinctly eccentric orbit with a total eccentricity of about 0.18. The two orbital periods, 1.20937 and 0.80693 days, are within 0.1% of a 3:2 ratio. We speculate that this might be the result of KCTF-driven evolution of a quadruple system and we discuss this hypothesis in the context of other quadruple systems composed of two EBs. We make our double EB fitting code publicly available to provide a tool for long-term monitoring of the mutual orbit in such systems.

In the seconds after collapse of a massive star, the newborn proto-neutron star (PNS) radiates neutrinos of all flavors. The absorption of electron-type neutrinos below the radius of the stalled shockwave may drive explosions (the "neutrino mechanism"). Because the heating rate is proportional to the square of neutrino energy, flavor conversion of mu and tau neutrinos to electron-type neutrinos via collective neutrino oscillations (CnuO) may in principle increase the heating rate and drive explosions. In order to assess the potential importance of CnuO for the shock revival, we solve the steady-state boundary value problem of spherically-symmetric accretion between the PNS surface (r_nu) and the shock (r_S), including a scheme for flavor conversion via CnuO. For a given r_nu, PNS mass (M), accretion rate (Mdot), and assumed values of the neutrino energies from the PNS, we calculate the critical neutrino luminosity above which accretion is impossible and explosion results. We show that CnuO can decrease the critical luminosity by a factor of at most ~1.5, but only if the flavor conversion is fully completed inside r_S and if there is no matter suppression. The magnitude of the effect depends on the model parameters (M, Mdot, and r_nu) through the shock radius and the physical scale for flavor conversion. We quantify these dependencies and find that CnuO could lower the critical luminosity only for small M and Mdot, and large r_nu. However, for these parameter values CnuO are suppressed due to matter effects. By quantifying the importance of CnuO and matter suppression at the critical neutrino luminosity for explosion, we show in agreement with previous studies that CnuO are unlikely to affect the neutrino mechanism of core-collapse supernovae significantly.

The observed distribution of neutron star (NS) masses reflects the physics of core-collapse supernova explosions and the structure of the massive stars that produce them at the end of their evolution. We present a Bayesian analysis that directly compares the NS mass distribution observed in double NS systems to theoretical models of NS formation. We find that models with standard binary mass ratio distributions are strongly preferred over independently picking the masses from the initial mass function, although the strength of the inference depends on whether current assumptions for identifying the remnants of the primary and secondary stars are correct. Second, NS formation models with no mass fallback are favored because they reduce the dispersion in NS masses. The double NS system masses thus directly point to the mass coordinate where the supernova explosion was initiated, making them an excellent probe of the supernova explosion mechanism. If we assume no fallback and simply vary the mass coordinate separating the remnant and the supernova ejecta, we find that for solar metallicity stars the explosion most likely develops at the edge of the iron core at a specific entropy of about 2.8 k_B. The primary limitations of our study are the poor knowledge of the supernova explosion mechanism and the lack of broad range of SN model explosions of LMC to solar metallicity.

We perform a global fit to ~5,000 radial velocity and ~177,000 magnitude measurements in 29 photometric bands covering 0.3 to 8.0 microns distributed among 287 Galactic, LMC, and SMC Cepheids with P > 10 days. We assume that the Cepheid light curves and radial velocities are fully characterized by distance, reddening, and time-dependent radius and temperature variations. We construct phase curves of radius and temperature for periods between 10 and 100 days, which yield light curve templates for all our photometric bands and can be easily generalized to any additional band. With only 4 to 6 parameters per Cepheid, depending on the existence of velocity data and the amount of freedom in the distance, the models have typical rms light and velocity curve residuals of 0.05 mag and 3.5 km/s. The model derives the mean Cepheid spectral energy distribution and its derivative with respect to temperature, which deviate from a black body in agreement with metal-line and molecular opacity effects. We determine a mean reddening law towards the Cepheids in our sample, which is not consistent with standard assumptions in either the optical or near-IR. Based on stellar atmosphere models we predict the biases in distance, reddening, and temperature determinations due to the metallicity and we quantify the metallicity signature expected for our fit residuals. The observed residuals as a function of wavelength show clear differences between the individual galaxies, which are compatible with these predictions. In particular, we find that metal-poor Cepheids appear hotter. Finally, we provide a framework for optimally selecting filters that yield the smallest overall errors in Cepheid parameter determination, or filter combinations for suppressing or enhancing the metallicity effects on distance determinations. We make our templates publicly available.

(Abridged) Neutrino heating may drive core-collapse supernova explosions. Although it is known that the stalled accretion shock turns into explosion when the neutrino luminosity from the collapsed core exceeds a critical value (L_crit) (the "neutrino mechanism"), the physics of L_crit, as well as its dependence on the properties of the proto-neutron star (PNS) and changes to the microphysics has never been systematically explored. We solve the one-dimensional steady-state accretion problem between the PNS surface and the accretion shock. We quantify the deep connection between the solution space of steady-state accretion flows with bounding shocks and the neutrino mechanism. We show that there is a maximum, critical sound speed above which it is impossible to maintain accretion with a standoff shock, because the shock jump conditions cannot be satisfied. The physics of this critical sound speed is general and does not depend on a specific heating mechanism. For the simple model of pressure-less free-fall onto a shock bounding an isothermal accretion flow with sound speed c_T, we show that if c_T^2/v_escape^2 > 3/16 explosion results. We generalize this result to the more complete supernova problem, showing explicitly that the same physics determines L_crit. We find that the critical condition for explosion can be written as c_S^2/v_escape^2 = 0.19, where c_S is the adiabatic sound speed. This "antesonic" condition describes L_crit over a broad range in accretion rate and microphysics. We show that the addition of the accretion luminosity (L_acc) reduces L_crit non-trivially. As in previous work, we find that L_crit is always significantly higher than the maximum possible value of L_acc. Finally, we provide evidence that the reduction in L_crit seen in recent multi-dimensional simulations results from a reduction in the efficiency of cooling, rather than an increase in the heating rate.

We present the first results from a detailed analysis of photometric and spectrophotometric data on the narrow-line Seyfert 1 galaxy Mrk 335, collected over a 120-day span in the fall of 2010. From these data we measure the lag in the He II 4686 broad emission line relative to the optical continuum to be 2.7 \pm 0.6 days and the lag in the H\beta 4861 broad emission line to be 13.9 \pm 0.9 days. Combined with the line width, the He II lag yields a black hole mass, MBH = (2.6 \pm 0.8)\times 10^7 Msun. This measurement is consistent with measurements made using the H\beta 4861 line, suggesting that the He II emission originates in the same region as H\beta, but at a much smaller radius. This constitutes the first robust lag measurement for a high-ionization line in a narrow-line Seyfert 1 galaxy.

We present a catalog of photometrically variable stars discovered within two 21'.3 X 21'.3 fields centered on the Cygnus OB2 association. There have hitherto been no deep optical variability studies of Cyg OB2 despite it being replete with early-type massive stars, perhaps due to the high and variable extinction (up to A_V ~ 20) that permeates much of the region. Here we provide results of the first variability study with this combination of spatial coverage (~ 0.5 deg) and photometric depth (R ~ 21 mag). We find 121 stars to be variable in both R- and I-band, 116 of them newly discovered. Of the 121 variables, we identify 27 eclipsing binaries (EBs) and eclipsing binary candidates, 20 potential Herbig Ae/Be stars, and 52 pulsating variables. Confirming both the status and the cluster membership of the Herbig Ae/Be stars would address the uncertainty regarding the age and star formation history of Cyg OB2. We match our catalog to known variables and binaries in the region, 2MASS near-IR (NIR) data, and Chandra X-ray observations to find counterparts to new variables in other wavelengths.

We present the discovery of a Neptune-mass planet OGLE-2007-BLG-368Lb with a planet-star mass ratio of q=[9.5 +/- 2.1] x 10^{-5} via gravitational microlensing. The planetary deviation was detected in real-time thanks to the high cadence of the MOA survey, real-time light curve monitoring and intensive follow-up observations. A Bayesian analysis returns the stellar mass and distance at M_l = 0.64_{-0.26}^{+0.21} M_\sun and D_l = 5.9_{-1.4}^{+0.9} kpc, respectively, so the mass and separation of the planet are M_p = 20_{-8}^{+7} M_\oplus and a = 3.3_{-0.8}^{+1.4} AU, respectively. This discovery adds another cold Neptune-mass planet to the planetary sample discovered by microlensing, which now comprise four cold Neptune/Super-Earths, five gas giant planets, and another sub-Saturn mass planet whose nature is unclear. The discovery of these ten cold exoplanets by the microlensing method implies that the mass ratio function of cold exoplanets scales as dN_{\rm pl}/d\log q \propto q^{-0.7 +/- 0.2} with a 95% confidence level upper limit of n < -0.35 (where dN_{\rm pl}/d\log q \propto q^n). As microlensing is most sensitive to planets beyond the snow-line, this implies that Neptune-mass planets are at least three times more common than Jupiters in this region at the 95% confidence level.

We use the recently released OGLE-III catalog of 17692 fundamental mode RR Lyr (RRab) stars in the Large Magellanic Cloud (LMC) to investigate the structure of its stellar halo. We apply conservative cuts in period, amplitude and magnitude to remove blends and other contamination. We use period--luminosity and period--color relations to determine distance and extinction of every star in our final sample of 9393 stars. In order to determine the scatter of our method, we compare the distributions of distances in two regions at the edges of the covered area with a central region. We determine the intrinsic line-of-sight dispersion in the center to be 0.135 mag or 3.21 kpc (FWHM of 0.318 mag or 7.56 kpc), assuming zero depth in one of the edge regions. The conservative cuts we apply reduce the derived depth significantly. Furthermore, we find that the distribution of RRab stars is deformed in the sense that stars on the Eastern side are closer than on the Western side. We model the RRab distribution as a triaxial ellipsoid and determine its axes ratios to be 1:2.00:3.50 with the longest axis inclined by 6 degrees from the line of sight. Another result of our analysis is an extinction map of the LMC and a map of internal reddening, which we make publicly available.

We demonstrate the newly developed resource for exoplanet researchers - The Exoplanet Transit Database. This database is designed to be a web application and it is open for any exoplanet observer. It came on-line in September 2008. The ETD consists of three individual sections. One serves for predictions of the transits, the second one for processing and uploading new data from the observers. We use a simple analytical model of the transit to calculate the central time of transit, its duration and the depth of the transit. These values are then plotted into the observed - computed diagrams (O-C), that represent the last part of the application.

A significant fraction of novae exhibit a series of rebrightenings on the decline branch of their light curves. We use visual observations to study this phase in several well-observed novae. We find that these rebrightenings are isolated flare-like events on otherwise smooth light curves and we show that in most novae in our sample the time intervals between consecutive flares gradually increase as a geometric series; rebrightenings are equally spaced in logarithmic time. We also find a correlation between the rate of increase in the time between rebrightenings with the speed class of the nova in the sense that slower novae tend to increase the time intervals faster. We attribute these rebrightenings to instabilities in the envelope hydrogen burning and we mention other cases of such timing pattern in astronomy and natural sciences in general.

We explore the sensitivity of two-point-mass gravitational microlensing to the extended nature of the source star, as well as the related sensitivity to its limb darkening. We demonstrate that the sensitive region, usually considered to be limited to a source-diameter-wide band along the caustic, is strongly expanded near cusps, most prominently along their outer axis. In the case of multi-component caustics, facing cusps may form a region with a non-negligible extended-source effect spanning the gap between them. We demonstrate that for smaller sources the size of the sensitive region extending from a cusp measured in units of source radii increases, scaling as the inverse cube root of the radius. We study the extent of different sensitivity contours and show that for a microlensed Galactic bulge giant the probability of encountering at least a 1% extended-source effect is higher than the probability of caustic crossing by 40-60% when averaged over a typical range of lens-component separations, with the actual value depending on the mass ratio of the components. We derive analytical expressions for the extended-source effect and chromaticity for a source positioned off the caustic. These formulae are more generally applicable to any gravitational lens with a sufficiently small source. Using exactly computed amplifications we test the often used linear-fold caustic approximation and show that it may lead to errors on the level of a few percent even in near-ideal caustic-crossing events. Finally, we discuss several interesting cases of observed binary and planetary microlensing events and point out the importance of our results for the measurement of stellar limb darkening from microlensing light curves.

We present an analysis of light variations in UBVRI of the anomalous object in the center of planetary nebula Sh 2-71. We refined the linear ephemeris of the light maxima to JD{max} = 2449862.0 + 68.101 (E-96), but also identified long-term, obviously non-periodic variations. The latter manifest themselves in large O-C shifts, a variable profile of light curves (hereafter LC) and changes in the mean brightness of the object. Our spectroscopic observations suggested the presence of a superdense nebula in the center of Sh 2-71.

We confirmed the presence of regular UBVRIc light variations of the object in the centre of planetary nebula Sh 2-71, the improved period of them being P =(68.132+/-0.005) days. The shapes and amplitudes of light curves in particular colours are briefly discussed.

We present visual and CCD VRI observations of a suspected symbiotic binary V335 Vul. Long term photometry suggests usual Mira or SRa pulsations. High speed V band CCD photometry didn't reveal anything unusual above 0.02 mag.