Showing 201–250 of 590 results for author: Drake, J

Globally Distributed Energetic Neutral Atom Maps for the "Croissant" Heliosphere

Authors: Marc Kornbleuth, Merav Opher, Adam T. Michael, James F. Drake

Abstract: A recent study by Opher et al. (2015) suggested the heliosphere has a "croissant" shape, where the heliosheath plasma is confined by the toroidal solar magnetic field. The "croissant" heliosphere is in contrast to the classically accepted view of a comet-like tail. We investigate the effect of the "croissant" heliosphere model on energetic neutral atom (ENA) maps. Regardless of the existence of a… ▽ More A recent study by Opher et al. (2015) suggested the heliosphere has a "croissant" shape, where the heliosheath plasma is confined by the toroidal solar magnetic field. The "croissant" heliosphere is in contrast to the classically accepted view of a comet-like tail. We investigate the effect of the "croissant" heliosphere model on energetic neutral atom (ENA) maps. Regardless of the existence of a split tail, the confinement of the heliosheath plasma should appear in ENA maps. ENA maps from the Interstellar Boundary Explorer (IBEX) have shown two high latitude lobes with excess ENA flux at higher energies in the tail of the heliosphere. These lobes could be a signature of the confinement of the heliosheath plasma, while some have argued they are caused by the fast/slow solar wind profile. Here we present ENA maps of the "croissant" heliosphere, focusing on understanding the effect of the heliosheath plasma collimation by the solar magnetic field while using a uniform solar wind. We incorporate pick-up ions (PUIs) into our model based on Malama et al. (2006) and Zank et al. (2010). We use the neutral solution from our MHD model to determine the angular variation of the PUIs, and include the extinction of PUIs in the heliosheath. In the presence of a uniform solar wind, we find that the collimation in the "croissant" heliosphere does manifest itself into two high latitude lobes of increased ENA flux in the downwind direction. △ Less

Submitted 17 August, 2018; originally announced August 2018.

Comments: 14 pages, 1 table, 7 figures, Accepted for publication in ApJ

Spectral curvature of shock-accelerated particles in solar cycle 23

Authors: Connie Zhou, Federico Fraschetti, Jeremy J. Drake, Martin Pohl

Abstract: We have fitted the protons momentum distribution of the 16 GLEs of solar cycle 23 with a log-parabola and have found a correlation between the previously determined energy break and the global spectral curvature of the log-parabola. We have fitted the protons momentum distribution of the 16 GLEs of solar cycle 23 with a log-parabola and have found a correlation between the previously determined energy break and the global spectral curvature of the log-parabola. △ Less

Submitted 14 August, 2018; originally announced August 2018.

Comments: 2 pages, 1 figure, accepted in RNAAS

The Fast, Luminous Ultraviolet Transient AT2018cow: Extreme Supernova, or Disruption of a Star by an Intermediate-Mass Black Hole?

Authors: Daniel A. Perley, Paolo A. Mazzali, Lin Yan, S. Bradley Cenko, Suvi Gezari, Kirsty Taggart, Nadia Blagorodnova, Christoffer Fremling, Brenna Mockler, Avinash Singh, Nozomu Tominaga, Masaomi Tanaka, Alan M. Watson, Tomás Ahumada, G. C. Anupama, Chris Ashall, Rosa L. Becerra, David Bersier, Varun Bhalerao, Joshua S. Bloom, Nathaniel R. Butler, Chris Copperwheat, Michael W. Coughlin, Kishalay De, Andrew J. Drake , et al. (38 additional authors not shown)

Abstract: Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwi… ▽ More Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwide campaign to observe AT2018cow (ATLAS18qqn), the first fast-luminous optical transient to be found in real time at low redshift. Our first spectra (<2 days after discovery) are entirely featureless. A very broad absorption feature suggestive of near-relativistic velocities develops between 3-8 days, then disappears. Broad emission features of H and He develop after >10 days. The spectrum remains extremely hot throughout its evolution, and the photospheric radius contracts with time (receding below R<10^14 cm after 1 month). This behaviour does not match that of any known supernova, although a relativistic jet within a fallback supernova could explain some of the observed features. Alternatively, the transient could originate from the disruption of a star by an intermediate-mass black hole, although this would require long-lasting emission of highly super-Eddington thermal radiation. In either case, AT2018cow suggests that the population of fast luminous transients represents a new class of astrophysical event. Intensive follow-up of this event in its late phases, and of any future events found at comparable distance, will be essential to better constrain their origins. △ Less

Submitted 23 November, 2018; v1 submitted 2 August, 2018; originally announced August 2018.

Comments: Corrected Figure 8 / Table 4 to use final fits. Includes machine-readable photometry table (hopefully for real this time)

Wave generation and heat flux suppression in astrophysical plasma systems

Authors: G. T. Roberg-Clark, J. F. Drake, M. Swisdak, C. S. Reynolds

Abstract: Heat flux suppression in collisionless plasmas for a large range of plasma $β$ is explored using two-dimensional particle-in-cell simulations with a strong, sustained thermal gradient. We find that a transition takes place between whistler-dominated (high-$β$) and double-layer-dominated (low-$β$) heat flux suppression. Whistlers saturate at small amplitude in the low beta limit and are unable to e… ▽ More Heat flux suppression in collisionless plasmas for a large range of plasma $β$ is explored using two-dimensional particle-in-cell simulations with a strong, sustained thermal gradient. We find that a transition takes place between whistler-dominated (high-$β$) and double-layer-dominated (low-$β$) heat flux suppression. Whistlers saturate at small amplitude in the low beta limit and are unable to effectively suppress the heat flux. Electrostatic double layers suppress the heat flux to a mostly constant factor of the free streaming value once this transition happens. The double layer physics is an example of ion-electron coupling and occurs on a scale of roughly the electron Debye length. The scaling of ion heating associated with the various heat flux driven instabilities is explored over the full range of $β$ explored. The range of plasma-$β$s studied in this work makes it relevant to the dynamics of a large variety of astrophysical plasmas, including the intracluster medium of galaxy clusters, hot accretion flows, stellar and accretion disk coronae, and the solar wind. △ Less

Submitted 24 September, 2018; v1 submitted 12 July, 2018; originally announced July 2018.

Comments: 16 pages, 6 figures. Accepted for publication in Astrophys. J

Journal ref: G. T. Roberg-Clark et al 2018 ApJ 867 154

The stellar rotation--activity relationship in fully convective M dwarfs

Authors: Nicholas J. Wright, Elisabeth R. Newton, Peter K. G. Williams, Jeremy J. Drake, Rakesh K. Yadav

Abstract: The coronal activity-rotation relationship is considered to be a proxy for the underlying stellar dynamo responsible for magnetic activity in solar and late-type stars. While this has been studied in considerable detail for partly-convective stars that are believed to operate an interface dynamo, it is poorly unconstrained in fully-convective stars that lack the necessary shear layer between radia… ▽ More The coronal activity-rotation relationship is considered to be a proxy for the underlying stellar dynamo responsible for magnetic activity in solar and late-type stars. While this has been studied in considerable detail for partly-convective stars that are believed to operate an interface dynamo, it is poorly unconstrained in fully-convective stars that lack the necessary shear layer between radiative core and the convective envelope. We present new X-ray observations of 19 slowly-rotating fully-convective stars with rotation periods from the MEarth Project. We use these to calculate X-ray luminosities (or upper limits for undetected sources) and combine these with existing measurements from \citet{wrig16b}. We confirm the existence of fully-convective stars in the X-ray unsaturated regime and find that these objects follow the same rotation-activity relationship seen for partly-convective stars. We measure the power-law slope of the relationship between Rossby number (the ratio of the rotation period to the convective turnover time) and the fractional X-ray luminosity for X-ray unsaturated fully-convective stars for the first time, and find it to be consistent with that found for partly-convective stars. We discuss this implications of this result for our understanding of stellar magnetic dynamos in fully- and partly-convective stars. Finally, we also use this data to improve empirical estimates of the convective turnover time for fully-convective stars. △ Less

Submitted 9 July, 2018; originally announced July 2018.

Comments: 12 pages, 5 figures, accepted for publication in MNRAS

Suppression of Coronal Mass Ejections in active stars by an overlying large-scale magnetic field: A numerical study

Authors: Julián D. Alvarado-Gómez, Jeremy J. Drake, Ofer Cohen, Sofia P. Moschou, Cecilia Garraffo

Abstract: We present results from a set of numerical simulations aimed at exploring the mechanism of coronal mass ejection (CME) suppression in active stars by an overlying large-scale magnetic field. We use a state-of-the-art 3D magnetohydrodynamic (MHD) code which considers a self-consistent coupling between an Alfvén wave-driven stellar wind solution, and a first-principles CME model based on the eruptio… ▽ More We present results from a set of numerical simulations aimed at exploring the mechanism of coronal mass ejection (CME) suppression in active stars by an overlying large-scale magnetic field. We use a state-of-the-art 3D magnetohydrodynamic (MHD) code which considers a self-consistent coupling between an Alfvén wave-driven stellar wind solution, and a first-principles CME model based on the eruption of a flux-rope anchored to a mixed polarity region. By replicating the driving conditions used in simulations of strong solar CMEs, we show that a large-scale dipolar magnetic field of $75$ G is able to fully confine eruptions within the stellar corona. Our simulations also consider CMEs exceeding the magnetic energy used in solar studies, which are able to escape the large-scale magnetic field confinement. The analysis includes a qualitative and quantitative description of the simulated CMEs and their dynamics, which reveals a drastic reduction of the radial speed caused by the overlying magnetic field. With the aid of recent observational studies, we place our numerical results in the context of solar and stellar flaring events. In this way, we find that this particular large-scale magnetic field configuration establishes a suppression threshold around $\sim$$3 \times 10^{32}$~erg in the CME kinetic energy. Extending the solar flare-CME relations to other stars, such CME kinetic energies could be typically achieved during erupting flaring events with total energies larger than $6 \times 10^{32}$ erg (GOES class $\sim$X70). △ Less

Submitted 7 June, 2018; originally announced June 2018.

Comments: 16 Pages, 8 figures, Accepted for publication in The Astrophysical Journal

Diffuse X-ray emission in the Cygnus OB2 association

Authors: J. F. Albacete Colombo, J. J. Drake, E. Flaccomio, N. J. Wright, V. Kashyap, M. G. Guarcello, K. Briggs, J. E. Drew, D. M. Fenech, G. Micela, M. McCollough, R. K. Prinja, N. Schneider, S. Sciortino, J. S. Vink

Abstract: We present a large-scale study of diffuse X-ray emission in the nearby massive stellar association Cygnus OB2 as part of the Chandra Cygnus OB2 Legacy Program. We used 40 Chandra X-ray ACIS-I observations covering $\sim$1.0 deg$^2$. After removing 7924 point-like sources detected in our survey, background-corrected X-ray emission, the adaptive smoothing reveals large-scale diffuse X-ray emission.… ▽ More We present a large-scale study of diffuse X-ray emission in the nearby massive stellar association Cygnus OB2 as part of the Chandra Cygnus OB2 Legacy Program. We used 40 Chandra X-ray ACIS-I observations covering $\sim$1.0 deg$^2$. After removing 7924 point-like sources detected in our survey, background-corrected X-ray emission, the adaptive smoothing reveals large-scale diffuse X-ray emission. Diffuse emission was detected in the sub-bands Soft [0.5 : 1.2] and Medium [1.2 : 2.5], and marginally in the Hard [2.5 : 7.0] keV band. From X-ray spectral analysis of stacked spectra we compute a total [0.5 : 7.0 keV] diffuse X-ray luminosity of L$_{\rm x}^{\rm diff}\approx$4.2$\times$10$^{\rm 34}$ erg s$^{-1}$, characterized with plasma temperature components at kT$\approx$ 0.11, 0.40 and 1.18 keV, respectively. The HI absorption column density corresponding to these temperatures has a distribution consistent with N$_{\rm H}$ = 0.43, 0.80 and 1.39 $\times$10$^{22}$ cm$^{-2}$. The extended medium band energy emission likely arises from O-type stellar winds thermalized by wind-wind collisions in the most populated regions of the association, while the soft band emission probably arises from less energetic termination shocks against the surrounding Interstellar-Medium. Super-soft and Soft diffuse emission appears more widely dispersed and intense than the medium band emission. The diffuse X-ray emission is generally spatially coincident with low-extinction regions that we attribute to the ubiquitous influence of powerful stellar winds from massive stars and their interaction with the local Interstellar-Medium. Diffuse X-ray emission is volume-filling, rather than edge-brightened, oppositely to other star-forming regions. We reveal the first observational evidence of X-ray haloes around some evolved massive stars. △ Less

Submitted 4 June, 2018; originally announced June 2018.

Comments: Accepted for publication on ApJSS, Cygnus OB2 Special Issue (2018). 20 pages, 19 figures

A new physical interpretation of optical and infrared variability in quasars

Authors: Nicholas P. Ross, K. E. Saavik Ford, Matthew Graham, Barry McKernan, Daniel Stern, Aaron M. Meisner, Roberto J. Assef, Arjun Dey, Andrew J. Drake, Hyunsung D. Jun, Dustin Lang

Abstract: Changing-look quasars are a recently identified class of active galaxies in which the strong UV continuum and/or broad optical hydrogen emission lines associated with unobscured quasars either appear or disappear on timescales of months to years. The physical processes responsible for this behaviour are still debated, but changes in the black hole accretion rate or accretion disk structure appear… ▽ More Changing-look quasars are a recently identified class of active galaxies in which the strong UV continuum and/or broad optical hydrogen emission lines associated with unobscured quasars either appear or disappear on timescales of months to years. The physical processes responsible for this behaviour are still debated, but changes in the black hole accretion rate or accretion disk structure appear more likely than changes in obscuration. Here we report on four epochs of spectroscopy of SDSS J110057.70-005304.5, a quasar at a redshift of $z=0.378$ whose UV continuum and broad hydrogen emission lines have faded, and then returned over the past $\approx$20 years. The change in this quasar was initially identified in the infrared, and an archival spectrum from 2010 shows an intermediate phase of the transition during which the flux below rest-frame $\approx$3400Å has decreased by close to an order of magnitude. This combination is unique compared to previously published examples of changing-look quasars, and is best explained by dramatic changes in the innermost regions of the accretion disk. The optical continuum has been rising since mid-2016, leading to a prediction of a rise in hydrogen emission line flux in the next year. Increases in the infrared flux are beginning to follow, delayed by a $\sim$3 year observed timescale. If our model is confirmed, the physics of changing-look quasars are governed by processes at the innermost stable circular orbit (ISCO) around the black hole, and the structure of the innermost disk. The easily identifiable and monitored changing-look quasars would then provide a new probe and laboratory of the nuclear central engine. △ Less

Submitted 27 September, 2018; v1 submitted 17 May, 2018; originally announced May 2018.

Comments: 13 pages, 4 figures, 3 tables. Published in MNRAS. All code and data links on GitHub, https://github.com/d80b2t/WISE_LC

Journal ref: 2018MNRAS.480.4468R

Anti-correlation between X-ray luminosity and pulsed fraction in the Small Magellanic Cloud pulsar SXP 1323

Authors: Jun Yang, Andreas Zezas, Malcolm J. Coe, Jeremy J. Drake, JaeSub Hong, Silas G. T. Laycock, Daniel R. Wik

Abstract: We report the evidence for the anti-correlation between pulsed fraction (PF) and luminosity of the X-ray pulsar SXP 1323, found for the first time in a luminosity range $10^{35}$--$10^{37}$ erg s$^{-1}$ from observations spanning 15 years. The phenomenon of a decrease in X-ray PF when the source flux increases has been observed in our pipeline analysis of other X-ray pulsars in the Small Magellani… ▽ More We report the evidence for the anti-correlation between pulsed fraction (PF) and luminosity of the X-ray pulsar SXP 1323, found for the first time in a luminosity range $10^{35}$--$10^{37}$ erg s$^{-1}$ from observations spanning 15 years. The phenomenon of a decrease in X-ray PF when the source flux increases has been observed in our pipeline analysis of other X-ray pulsars in the Small Magellanic Cloud (SMC). It is expected that the luminosity under a certain value decreases as the PF decreases due to the propeller effect. Above the propeller region, an anti-correlation between the PF and flux might occur either as a result of an increase in the un-pulsed component of the total emission or a decrease of the pulsed component. Additional modes of accretion may also be possible, such as spherical accretion and a change in emission geometry. At higher mass accretion rates, the accretion disk could also extend closer to the neutron star (NS) surface, where a reduced inner radius leads to hotter inner disk emission. These modes of plasma accretion may affect the change in the beam configuration to fan-beam dominant emission. △ Less

Submitted 11 May, 2018; originally announced May 2018.

Comments: It has been accepted for publication in Monthly Notices of the Royal Astronomical Society Letters

The Reduction of Magnetic Reconnection Outflow Jets to Sub-Alfvénic Speeds

Authors: Colby C. Haggerty, Michael A. Shay, Alexandros Chasapis, Tai D. Phan, James F. Drake, Kittipat Malakit, Paul A. Cassak, Rungployphan Kieokaew

Abstract: The outflow velocity of jets produced by collisionless magnetic reconnection is shown to be reduced by the ion exhaust temperature in simulations and observations. We derive a scaling relationship for the outflow velocity based on the upstream Alfvén speed and the parallel ion exhaust temperature, which is verified in kinetic simulations and observations. The outflow speed reduction is shown to be… ▽ More The outflow velocity of jets produced by collisionless magnetic reconnection is shown to be reduced by the ion exhaust temperature in simulations and observations. We derive a scaling relationship for the outflow velocity based on the upstream Alfvén speed and the parallel ion exhaust temperature, which is verified in kinetic simulations and observations. The outflow speed reduction is shown to be due to the firehose instability criterion, and so for large enough guide fields this effect is suppressed and the outflow speed reaches the upstream Alfvén speed based on the reconnecting component of the magnetic field. △ Less

Submitted 25 April, 2018; originally announced April 2018.

The Revolution Revolution: magnetic morphology driven spin-down

Authors: C. Garraffo, J. J. Drake, A. Dotter, J. Choi, D. J. Burke, S. P. Moschou, J. D. Alvarado-Gomez, V. L. Kashyap, O. Cohen

Abstract: Observations of young open clusters show a bimodal distribution of rotation periods that has been difficult to explain with existing stellar spin-down models. Detailed MHD stellar wind simulations have demonstrated that surface magnetic field morphology has a strong influence on wind-driven angular momentum loss. Observations suggest that faster rotating stars store a larger fraction of their magn… ▽ More Observations of young open clusters show a bimodal distribution of rotation periods that has been difficult to explain with existing stellar spin-down models. Detailed MHD stellar wind simulations have demonstrated that surface magnetic field morphology has a strong influence on wind-driven angular momentum loss. Observations suggest that faster rotating stars store a larger fraction of their magnetic flux in higher-order multipolar components of the magnetic field. In this work, we present a new model for stellar spin-down that, for the first time, accounts for the stellar surface magnetic field configuration. We show how a magnetic complexity that evolves from complex toward simple configurations as a star spins down can explain the salient features of stellar rotation evolution, including the bimodal distribution of both slow and fast rotators seen in young open clusters. △ Less

Submitted 15 June, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Comments: 8 pages, 4 figures, to appear in ApJ

Characterizing Ion Flows Across a Dipolarization Jet

Authors: Harry Arnold, Marc Swisdak, James Drake

Abstract: The structure of dipolarization jets with finite width in the dawn-dusk direction relevant to magnetic reconnection in the Earth's magnetotail is explored with particle-in-cell simulations. We carry out Riemann simulations of the evolution of the jet in the dawn-dusk, north-south plane to investigate the dependence of the jet structure on the jet width in the dawn-dusk direction. We find that the… ▽ More The structure of dipolarization jets with finite width in the dawn-dusk direction relevant to magnetic reconnection in the Earth's magnetotail is explored with particle-in-cell simulations. We carry out Riemann simulations of the evolution of the jet in the dawn-dusk, north-south plane to investigate the dependence of the jet structure on the jet width in the dawn-dusk direction. We find that the magnetic field and Earth-directed ion flow structure depend on the dawn-dusk width. A reversal in the usual Hall magnetic field near the center of the current sheet on the dusk side of larger jets is observed. For small widths, the maximum velocity of the Earthward flow is significantly reduced below the theoretical limit of the upstream Alfvén speed. However, the ion flow speed approaches this limit once the width exceeds the ion Larmor radius based on the normal magnetic field, $B_z$. △ Less

Submitted 16 March, 2018; originally announced March 2018.

Energy Dissipation in the Upper Atmospheres of Trappist-1 Planets

Authors: O. Cohen, A. Glocer, C. Garraffo, J. D. Drake, J. M. Bell

Abstract: We present a method to quantify the upper-limit of the energy transmitted from the intense stellar wind to the upper atmospheres of three of the Trappist-1 planets (e, f, and g). We use a formalism that treats the system as two electromagnetic regions, where the efficiency of the energy transmission between one region (the stellar wind at the planetary orbits) to the other (the planetary ionospher… ▽ More We present a method to quantify the upper-limit of the energy transmitted from the intense stellar wind to the upper atmospheres of three of the Trappist-1 planets (e, f, and g). We use a formalism that treats the system as two electromagnetic regions, where the efficiency of the energy transmission between one region (the stellar wind at the planetary orbits) to the other (the planetary ionospheres) depends on the relation between the conductances and impedances of the two regions. Since the energy flux of the stellar wind is very high at these planetary orbits, we find that for the case of high transmission efficiency (when the conductances and impedances are close in magnitude), the energy dissipation in the upper planetary atmospheres is also very large. On average, the Ohmic energy can reach $0.5-1~W/m^2$, about 1\% of the stellar irradiance and 5-15 times the EUV irradiance. Here, using constant values for the ionospheric conductance, we demonstrate that the stellar wind energy could potentially drive large atmospheric heating in terrestrial planets, as well as in hot jupiters. More detailed calculations are needed to assess the ionospheric conductance and to determine more accurately the amount of heating the stellar wind can drive in close-orbit planets. △ Less

Submitted 13 March, 2018; originally announced March 2018.

Comments: 6 pages, 1 table, 2 figures, accepted to ApJ Letters

Exploring Extreme Space Weather Factors of Exoplanetary Habitability

Authors: V. S. Airapetian, V. Adibekyan, M. Ansdell, O. Cohen, M. Cuntz, W. Danchi, C. F. Dong, J. J. Drake, A. Fahrenbach, K. France, K. Garcia-Sage, A. Glocer, J. L. Grenfell, G. Gronoff, H. Hartnett, W. Henning, N. R. Hinkel, A. G. Jensen, M. Jin, P. Kalas, S. R. Kane, K. Kobayashi, R. Kopparapu, J. Leake, M. López-Puertas , et al. (24 additional authors not shown)

Abstract: It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems.… ▽ More It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years. △ Less

Submitted 9 March, 2018; originally announced March 2018.

Comments: 6 pages, the white paper submitted to the US National Academy of Sciences call on Exoplanet Science Strategy

The Solar Wind Environment in Time

Authors: Quentin Pognan, Cecilia Garraffo, Ofer Cohen, Jeremy J. Drake

Abstract: We use magnetograms of 8 solar analogues of ages 30~Myr to 3.6~Gyr obtained from Zeeman Doppler Imaging (ZDI) and taken from the literature, together with two solar magnetograms, to drive magnetohydrodynamical (MHD) wind simulations and construct an evolutionary scenario of the solar wind environment and its angular momentum loss rate. With observed magnetograms of the radial field strength as the… ▽ More We use magnetograms of 8 solar analogues of ages 30~Myr to 3.6~Gyr obtained from Zeeman Doppler Imaging (ZDI) and taken from the literature, together with two solar magnetograms, to drive magnetohydrodynamical (MHD) wind simulations and construct an evolutionary scenario of the solar wind environment and its angular momentum loss rate. With observed magnetograms of the radial field strength as the only variant in the wind model, we find that power law model fitted to the derived angular momentum loss rate against time, $t$, results in a spin down relation $Ω\propto t^{-0.51}$, for angular speed $Ω$, which is remarkably consistent with the well-established Skumanich law $Ω\propto t^{-0.5}$. We use the model wind conditions to estimate the magnetospheric standoff distances for an Earth-like test planet situated at 1~AU for each of the stellar cases, and to obtain trends of minimum and maximum wind ram pressure and average ram pressure in the solar system through time. The wind ram pressure declines with time as $\overline{P_{ram}}\propto t^{2/3}$, amounting to a factor of 50 or so over the present lifetime of the solar system. △ Less

Submitted 16 February, 2018; v1 submitted 14 February, 2018; originally announced February 2018.

Comments: 21 pages, 11 figures, accepted for publication in ApJ, 8th February 2018

SNhunt151: an explosive event inside a dense cocoon

Authors: N. Elias-Rosa, S. Benetti, E. Cappellaro, A. Pastorello, G. Terreran, A. Morales-Garoffolo, S. C. Howerton, S. Valenti, E. Kankare, A. J. Drake, S. G. Djorgovski, L. Tomasella, L. Tartaglia, T. Kangas, P. Ochner, A. V. Filippenko, F. Ciabattari, S. Geier, D. A. Howell, J. Isern, S. Leonini, G. Pignata, M. Turatto

Abstract: SNhunt151 was initially classified as a supernova (SN) impostor (nonterminal outburst of a massive star). It exhibited a slow increase in luminosity, lasting about 450 d, followed by a major brightening that reaches M_V ~ -18 mag. No source is detected to M_V > -13 mag in archival images at the position of SNhunt151 before the slow rise. Low-to-mid-resolution optical spectra obtained during the pr… ▽ More SNhunt151 was initially classified as a supernova (SN) impostor (nonterminal outburst of a massive star). It exhibited a slow increase in luminosity, lasting about 450 d, followed by a major brightening that reaches M_V ~ -18 mag. No source is detected to M_V > -13 mag in archival images at the position of SNhunt151 before the slow rise. Low-to-mid-resolution optical spectra obtained during the pronounced brightening show very little evolution, being dominated at all times by multicomponent Balmer emission lines, a signature of interaction between the material ejected in the new outburst and the pre-existing circumstellar medium. We also analyzed mid-infrared images from the Spitzer Space Telescope, detecting a source at the transient position in 2014 and 2015. Overall, SNhunt151 is spectroscopically a Type IIn SN, somewhat similar to SN2009ip. However, there are also some differences, such as a slow pre-discovery rise, a relatively broad light-curve peak showing a longer rise time (~ 50 d) and a slower decline, along with a negligible change in the temperature around the peak (T < 10^4 K). We suggest that SNhunt151 is the result of an outburst, or a SN explosion, within a dense circumstellar nebula, similar to those embedding some luminous blue variables like Eta Carinae and originating from past mass-loss events. △ Less

Submitted 9 January, 2018; originally announced January 2018.

Comments: Accepted for publication in MNRAS. 19 pages with 10 tables and 11 figures

A and F stars as probes of outer Galactic disc kinematics

Authors: A. Harris, J. E. Drew, H. J. Farnhill, M. Monguió, M. Gebran, N. J. Wright, J. J. Drake, S. E. Sale

Abstract: Previous studies of the rotation law in the outer Galactic disc have mainly used gas tracers or clump giants. Here, we explore A and F stars as alternatives: these provide a much denser sampling in the outer disc than gas tracers and have experienced significantly less velocity scattering than older clump giants. This first investigation confirms the suitability of A stars in this role. Our work i… ▽ More Previous studies of the rotation law in the outer Galactic disc have mainly used gas tracers or clump giants. Here, we explore A and F stars as alternatives: these provide a much denser sampling in the outer disc than gas tracers and have experienced significantly less velocity scattering than older clump giants. This first investigation confirms the suitability of A stars in this role. Our work is based on spectroscopy of $\sim$ 1300 photometrically-selected stars in the red calcium-triplet region, chosen to mitigate against the effects of interstellar extinction. The stars are located in two low Galactic latitude sightlines, at longitudes $\ell = 118^{\circ}$, sampling strong Galactic rotation shear, and $\ell = 178^{\circ}$, near the Anticentre. With the use of Markov Chain Monte Carlo parameter fitting, stellar parameters and radial velocities are measured, and distances computed. The obtained trend of radial velocity with distance is inconsistent with existing flat or slowly rising rotation laws from gas tracers (Brand & Blitz 1993; Reid et al. 2014). Instead, our results fit in with those obtained by Huang et al. (2016) from disc clump giants that favoured rising circular speeds. An alternative interpretation in terms of spiral arm perturbation is not straight forward. We assess the role that undetected binaries in the sample and distance error may have in introducing bias, and show that the former is a minor factor. The random errors in our trend of circular velocity are within $\pm 5$ km s$^{-1}$. △ Less

Submitted 3 January, 2018; originally announced January 2018.

Localized Oscillatory Dissipation in Magnetopause Reconnection

Authors: J. L. Burch, R. E. Ergun, P. A. Cassak, J. M. Webster, R. B. Torbert, B. L. Giles, J. C. Dorelli, A. C. Rager, K. -J. Hwang, T. D. Phan, K. J. Genestreti, R. C. Allen, L. -J. Chen, S. Wang, D. Gershman, O. Le Contel, C. T. Russell, R. J. Strangeway, F. D. Wilder, D. B. Graham, M. Hesse, J. F. Drake, M. Swisdak, L. M. Price, M. A. Shay , et al. (4 additional authors not shown)

Abstract: Data from the NASA Magnetospheric Multiscale (MMS) mission are used to investigate asymmetric magnetic reconnection at the dayside boundary between the Earth's magnetosphere and the solar wind (the magnetopause). High-resolution measurements of plasmas, electric and magnetic fields, and waves are used to identify highly localized (~15 electron Debye lengths) standing wave structures with large ele… ▽ More Data from the NASA Magnetospheric Multiscale (MMS) mission are used to investigate asymmetric magnetic reconnection at the dayside boundary between the Earth's magnetosphere and the solar wind (the magnetopause). High-resolution measurements of plasmas, electric and magnetic fields, and waves are used to identify highly localized (~15 electron Debye lengths) standing wave structures with large electric-field amplitudes (up to 100 mV/m). These wave structures are associated with spatially oscillatory dissipation, which appears as alternatingly positive and negative values of J dot E (dissipation). For small guide magnetic fields the wave structures occur in the electron stagnation region at the magnetosphere edge of the EDR. For larger guide fields the structures also occur near the reconnection x-line. This difference is explained in terms of channels for the out-of-plane current (agyrotropic electrons at the stagnation point and guide-field-aligned electrons at the x-line). △ Less

Submitted 13 December, 2017; originally announced December 2017.

Unmixing the Galactic Halo with RR Lyrae tagging

Authors: V. Belokurov, A. J. Deason, S. E. Koposov, M. Catelan, D. Erkal, A. J. Drake, N. W. Evans

Abstract: We show that tagging RR Lyrae stars according to their location in the period-amplitude diagram can be used to shed light on the genesis of the Galactic stellar halo. The mixture of RR Lyrae of ab type, separated into classes along the lines suggested by Oosterhoff, displays a strong and coherent evolution with Galactocentric radius. The change in the RR Lyrae composition appears to coincide with… ▽ More We show that tagging RR Lyrae stars according to their location in the period-amplitude diagram can be used to shed light on the genesis of the Galactic stellar halo. The mixture of RR Lyrae of ab type, separated into classes along the lines suggested by Oosterhoff, displays a strong and coherent evolution with Galactocentric radius. The change in the RR Lyrae composition appears to coincide with the break in the halo's radial density profile at ~25 kpc. Using simple models of the stellar halo, we establish that at least three different types of accretion events are necessary to explain the observed RRab behavior. Given that there exists a correlation between the RRab class fraction and the total stellar content of a dwarf satellite, we hypothesize that the field halo RRab composition is controlled by the mass of the progenitor contributing the bulk of the stellar debris at the given radius. This idea is tested against a suite of cosmological zoom-in simulations of Milky Way-like stellar halo formation. Finally, we study some of the most prominent stellar streams in the Milky Way halo and demonstrate that their RRab class fractions follow the trends established previously. △ Less

Submitted 20 November, 2017; originally announced November 2017.

Comments: submitted to MNRAS, comments welcome. 12 pages, 9 figures

A Monster CME Obscuring A Demon Star Flare

Authors: Sofia-Paraskevi Moschou, Jeremy J. Drake, Ofer Cohen, Julian D. Alvarado-Gomez, Cecilia Garraffo

Abstract: We explore the scenario of a Coronal Mass Ejection (CME) being the cause of the observed continuous X-ray absorption of the August 30 1997 superflare on the eclipsing binary Algol (the Demon Star). The temporal decay of the absorption is consistent with absorption by a CME undergoing self-similar evolution with uniform expansion velocity. We investigate the kinematic and energetic properties of th… ▽ More We explore the scenario of a Coronal Mass Ejection (CME) being the cause of the observed continuous X-ray absorption of the August 30 1997 superflare on the eclipsing binary Algol (the Demon Star). The temporal decay of the absorption is consistent with absorption by a CME undergoing self-similar evolution with uniform expansion velocity. We investigate the kinematic and energetic properties of the CME using the ice-cream cone model for its three-dimensional structure in combination with the observed profile of the hydrogen column density decline with time. Different physically justified length scales were used that allowed us to estimate lower and upper limits of the possible CME characteristics. Further consideration of the maximum available magnetic energy in starspots leads us to quantify its mass as likely lying in the range $2\times 10^{21}$-$2\times 10^{22}$ g and kinetic energy in the range $7\times 10^{35}$-$3\times 10^{38}$ erg. The results are in reasonable agreement with extrapolated relations between flare X-ray fluence and CME mass and kinetic energy derived for solar CMEs. △ Less

Submitted 19 October, 2017; originally announced October 2017.

Comments: 10 pages, 5 figures, accepted for publication in ApJ

Localized and intense energy conversion in the diffusion region of asymmetric magnetic reconnection

Authors: M. Swisdak, J. F. Drake, L. Price, J. L. Burch, P. A. Cassak, T. -D. Phan

Abstract: We analyze a high-resolution simulation of magnetopause reconnection observed by the Magnetospheric Multiscale (MMS) mission and explain the occurrence of strongly localized dissipation with an amplitude more than an order of magnitude larger than expected. Unlike symmetric reconnection, wherein reconnection of the ambient reversed magnetic field drives the dissipation, we find the annihilation of… ▽ More We analyze a high-resolution simulation of magnetopause reconnection observed by the Magnetospheric Multiscale (MMS) mission and explain the occurrence of strongly localized dissipation with an amplitude more than an order of magnitude larger than expected. Unlike symmetric reconnection, wherein reconnection of the ambient reversed magnetic field drives the dissipation, we find the annihilation of the self-generated, out-of-plane (Hall) magnetic field plays the dominant role. Electrons flow along the magnetosheath separatrices, converge in the diffusion region, and jet past the X-point into the magnetosphere. The resulting accumulation of negative charge generates intense parallel electric fields that eject electrons along the magnetospheric separatrices and produce field-aligned beams. Many of these features match MMS observations. △ Less

Submitted 14 May, 2018; v1 submitted 12 October, 2017; originally announced October 2017.

Comments: Accepted version

Far beyond the Sun: I. The beating magnetic heart in Horologium

Authors: Julián D. Alvarado-Gómez, Gaitee A. J. Hussain, Jeremy J. Drake, Jean-François Donati, Jorge Sanz-Forcada, Beate Stelzer, Ofer Cohen, Eliana M. Amazo-Gómez, Jason H. Grunhut, Cecilia Garraffo, Sofia P. Moschou, James Silvester, Mary E. Oksala

Abstract: A former member of the Hyades cluster, $ι$ Horologii ($ι$ Hor) is a planet-hosting Sun-like star which displays the shortest coronal activity cycle known to date ($P_{\rm cyc}\sim$ 1.6 yr). With an age of $\sim$$625$ Myr, $ι$ Hor is also the youngest star with a detected activity cycle. The study of its magnetic properties holds the potential to provide fundamental information to understand the or… ▽ More A former member of the Hyades cluster, $ι$ Horologii ($ι$ Hor) is a planet-hosting Sun-like star which displays the shortest coronal activity cycle known to date ($P_{\rm cyc}\sim$ 1.6 yr). With an age of $\sim$$625$ Myr, $ι$ Hor is also the youngest star with a detected activity cycle. The study of its magnetic properties holds the potential to provide fundamental information to understand the origin of cyclic activity and stellar magnetism in late-type stars. In this series of articles, we present the results of a comprehensive project aimed at studying the evolving magnetic field in this star and how this evolution influences its circumstellar environment. This paper summarizes the first stage of this investigation, with results from a long-term observing campaign of $ι$ Hor using ground-based high-resolution spectropolarimetry. The analysis includes precise measurements of the magnetic activity and radial velocity of the star, and their multiple time-scales of variability. In combination with values reported in the literature, we show that the long-term chromospheric activity evolution of $ι$ Hor follows a beating pattern, caused by the superposition of two periodic signals of similar amplitude at $P_{1} \simeq 1.97 \pm 0.02$ yr and $P_{2} \simeq 1.41 \pm 0.01$ yr. Additionally, using the most recent parameters for $ι$ Hor b in combination with our activity and radial velocity measurements, we find that stellar activity dominates the radial velocity residuals, making the detection of additional planets in this system challenging. Finally, we report here the first measurements of the surface longitudinal magnetic field strength of $ι$ Hor, which displays varying amplitudes within $\pm4$ G and served to estimate the rotation period of the star ($P_{\rm rot} =~7.70^{+0.18}_{-0.67}$ d). △ Less

Submitted 6 October, 2017; originally announced October 2017.

Comments: 14 pages, 7 Figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

Mottled protoplanetary disk ionization by magnetically-channeled T Tauri star energetic particles

Authors: Federico Fraschetti, Jeremy J. Drake, Ofer Cohen, Cecilia Garraffo

Abstract: The evolution of protoplanetary disks is believed to be driven largely by angular momentum transport resulting from magnetized disk winds and turbulent viscosity. The ionization of the disk that is essential for these processes has been thought due to host star coronal X-rays but could also arise from energetic particles produced by coronal flares or by travelling shock waves and advected by the s… ▽ More The evolution of protoplanetary disks is believed to be driven largely by angular momentum transport resulting from magnetized disk winds and turbulent viscosity. The ionization of the disk that is essential for these processes has been thought due to host star coronal X-rays but could also arise from energetic particles produced by coronal flares or by travelling shock waves and advected by the stellar wind. We have performed test-particle numerical simulations of energetic protons propagating into a realistic T~Tauri stellar wind, including a superposed small-scale magnetostatic turbulence. The isotropic (Kolmogorov power spectrum) turbulent component is synthesised along the individual particle trajectories. We have investigated the energy range $[0.1 - 10]$ GeV, consistent with expectations from {\it Chandra} X-ray observations of large flares on T~Tauri stars and with recent indications by the {\it Herschel} Space Observatory of a significant contribution of energetic particles to the disk ionization of young stars. In contrast with a previous theoretical study finding dominance of energetic particles over X-ray in the ionization throughout the disk, we find that the disk ionization is likely dominated by X-rays over much of its area except within narrow regions where particles are channeled onto the disk by the strongly-tangled and turbulent magnetic field. The radial thickness of such regions is $\sim 5$ stellar radii close to the star and broadens with increasing radial distance. This likely continues out to large distances from the star ($10$ AU or greater) where particles can be copiously advected and diffused by the turbulent wind. △ Less

Submitted 29 December, 2017; v1 submitted 3 October, 2017; originally announced October 2017.

Comments: 18 pages, 14 figures, ApJ, in press

Suppression of electron thermal conduction by whistler turbulence in a sustained thermal gradient

Authors: G. T. Roberg-Clark, J. F. Drake, C. S. Reynolds, M. Swisdak

Abstract: The dynamics of weakly magnetized collisionless plasmas in the presence of an imposed temperature gradient along an ambient magnetic field is explored with particle-in-cell simulations and modeling. Two thermal reservoirs at different temperatures drive an electron heat flux that destabilizes off-angle whistler-type modes. The whistlers grow to large amplitude, $δB / B_{0} \simeq 1$, and resonantl… ▽ More The dynamics of weakly magnetized collisionless plasmas in the presence of an imposed temperature gradient along an ambient magnetic field is explored with particle-in-cell simulations and modeling. Two thermal reservoirs at different temperatures drive an electron heat flux that destabilizes off-angle whistler-type modes. The whistlers grow to large amplitude, $δB / B_{0} \simeq 1$, and resonantly scatter the electrons, significantly reducing the heat flux. A surprise is that the resulting steady state heat flux is largely independent of the thermal gradient. The rate of thermal conduction is instead controlled by the finite propagation speed of the whistlers, which act as mobile scattering centers that convect the thermal energy of the hot reservoir. The results are relevant to thermal transport in high $β$ astrophysical plasmas such as hot accretion flows and the intracluster medium of galaxy clusters. △ Less

Submitted 19 December, 2017; v1 submitted 31 August, 2017; originally announced September 2017.

Comments: Accepted for publication in Physical Review Letters. 5 pages, 5 figures. Supplementary material: 1 page, 1 figure

Journal ref: Phys. Rev. Lett. 120, 035101 (2018)

Redshifted X-rays from the material accreting onto TW Hya: evidence of a low-latitude accretion spot

Authors: C. Argiroffi, J. J. Drake, R. Bonito, S. Orlando, G. Peres, M. Miceli

Abstract: High resolution spectroscopy, providing constraints on plasma motions and temperatures, is a powerful means to investigate the structure of accretion streams in CTTS. In particular, the accretion shock region, where the accreting material is heated to temperatures of a few MK as it continues its inward bulk motion, can be probed by X-ray spectroscopy. To attempt to detect for the first time the mo… ▽ More High resolution spectroscopy, providing constraints on plasma motions and temperatures, is a powerful means to investigate the structure of accretion streams in CTTS. In particular, the accretion shock region, where the accreting material is heated to temperatures of a few MK as it continues its inward bulk motion, can be probed by X-ray spectroscopy. To attempt to detect for the first time the motion of this X-ray-emitting post-shock material, we searched for a Doppler shift in the deep Chandra/HETGS observation of the CTTS TW Hya. This test should unveil the nature of this X-ray emitting plasma component in CTTS, and constrain the accretion stream geometry. We searched for a Doppler shift in the X-ray emission from TW Hya with two different methods, by measuring the position of a selected sample of emission lines, and by fitting the whole TW Hya X-ray spectrum, allowing the line-of-sight velocity to vary. We found that the plasma at T~2-4 MK has a line-of-sight velocity of 38.3+/-5.1 km/s with respect to the stellar photosphere. This result definitively confirms that this X-ray-emitting material originates in the post-shock region, at the base of the accretion stream, and not in coronal structures. The comparison of the observed velocity along the line of sight, 38.3+/-5.1 km/s, with the inferred intrinsic velocity of the post shock of TW Hya, v_post~110-120 km/s, indicates that the footpoints of the accretion streams on TW Hya are located at low latitudes on the stellar surface. Our results indicate that complex magnetic field geometries, such as that of TW Hya, permit low-latitude accretion spots. Moreover, since on TW Hya the redshift of the soft X-ray emission is very similar to that of the narrow component of the CIV resonance doublet at 1550 Ang, as found by Ardila et al. (2013), then the plasma at 2-4 MK and that at 0.1 MK likely originate in the same post-shock regions. △ Less

Submitted 4 September, 2017; v1 submitted 10 August, 2017; originally announced August 2017.

Comments: Accepted for publication in Astronomy & Astrophysics; 2nd version after language editor corrections; 16 pages, 8 figures, 6 tables

Supernovae 2016bdu and 2005gl, and their link with SN 2009ip-like transients: another piece of the puzzle

Authors: A. Pastorello, C. S. Kochanek, M. Fraser, S. Dong, N. Elias-Rosa, S. Benetti, E. Cappellaro, L. Tomasella, A. J. Drake, J. Hermanen, T. Reynolds, B. J. Shappee, S. J. Smartt, K. C. Chambers, M. E. Huber, K. Smith, K. Z. Stanek, A. V. Filippenko, E. J. Christensen, L. Denneau, S. G. Djorgovski, H. Flewelling, C. Gall, A. Gal-Yam, S. Geier , et al. (30 additional authors not shown)

Abstract: Supernova (SN) 2016bdu is an unusual transient resembling SN 2009ip. SN 2009ip-like events are characterized by a long-lasting phase of erratic variability which ends with two luminous outbursts a few weeks apart. The second outburst is significantly more luminous (about 3 mag) than the first. In the case of SN 2016bdu, the first outburst (Event A) reached an absolute magnitude M(r) ~ -15.3 mag, w… ▽ More Supernova (SN) 2016bdu is an unusual transient resembling SN 2009ip. SN 2009ip-like events are characterized by a long-lasting phase of erratic variability which ends with two luminous outbursts a few weeks apart. The second outburst is significantly more luminous (about 3 mag) than the first. In the case of SN 2016bdu, the first outburst (Event A) reached an absolute magnitude M(r) ~ -15.3 mag, while the second one (Event B) occurred over one month later and reached M(r) ~ -18 mag. By inspecting archival data, a faint source at the position of SN 2016bdu is detectable several times in the past few years. We interpret these detections as signatures of a phase of erratic variability, similar to that experienced by SN 2009ip between 2008 and mid-2012, and resembling the currently observed variability of the luminous blue variable SN 2000ch in NGC 3432. Spectroscopic monitoring of SN 2016bdu during the second peak initially shows features typical of a SN IIn. One month after the Event B maximum, the spectra develop broad Balmer lines with P Cygni profiles and broad metal features. At these late phases, the spectra resemble those of a typical Type II SN. All members of this SN 2009ip-like group are remarkably similar to the Type IIn SN 2005gl. For this object, the claim of a terminal SN explosion is supported by the disappearance of the progenitor star. The similarity with SN 2005gl suggests that all members of this family may finally explode as genuine SNe, although the unequivocal detection of nucleosynthesised elements in their nebular spectra is still missing. △ Less

Submitted 3 July, 2017; originally announced July 2017.

Comments: Submitted to MNRAS on April 10, 2017; re-submitted on June 23 including suggestions from the referee. 24 pages, 12 figures, 5 tables

Long-term Periodicities of Cataclysmic Variables with Synoptic Surveys

Authors: Michael Ting-Chang Yang, Yi Chou, Chow-Choong Ngeow, Chin-Ping Hu, Yi-Hao Su, Thomas A. Prince, Shrinivas R. Kulkarni, David Levitan, Russ Laher, Jason Surace, Andrew J. Drake, Stanislav G. Djorgovski, Ashish A. Mahabal, Matthew J. Graham, Ciro Donalek

Abstract: A systematic study on the long-term periodicities of known Galactic cataclysmic variables (CVs) was conducted. Among 1580 known CVs, 344 sources were matched and extracted from the Palomar Transient Factory (PTF) data repository. The PTF light curves were combined with the Catalina Real-Time Transient Survey (CRTS) light curves and analyzed. Ten targets were found to exhibit long-term periodic var… ▽ More A systematic study on the long-term periodicities of known Galactic cataclysmic variables (CVs) was conducted. Among 1580 known CVs, 344 sources were matched and extracted from the Palomar Transient Factory (PTF) data repository. The PTF light curves were combined with the Catalina Real-Time Transient Survey (CRTS) light curves and analyzed. Ten targets were found to exhibit long-term periodic variability, which is not frequently observed in the CV systems. These long-term variations are possibly caused by various mechanisms, such as the precession of the accretion disk, hierarchical triple star system, magnetic field change of the companion star, and other possible mechanisms. We discuss the possible mechanisms in this study. If the long-term period is less than several tens of days, the disk precession period scenario is favored. However, the hierarchical triple star system or the variations in magnetic field strengths are most likely the predominant mechanisms for longer periods. △ Less

Submitted 20 June, 2017; originally announced June 2017.

Comments: 33 pages, 9 figures (manuscript form), Accepted for publication in PASP

The Threatening Environment of the TRAPPIST-1 Planets

Authors: Cecilia Garraffo, Jeremy J. Drake, Ofer Cohen, Julian D. Alvarado-Gomez, Sofia P. Moschou

Abstract: Recently, four additional Earth-mass planets were discovered orbiting the nearby ultracool M8 dwarf TRAPPIST-1, making a remarkable total of seven planets with equilibrium temperatures compatible with the presence of liquid water on their surface. Temperate terrestrial planets around an M-dwarf orbit close to their parent star, rendering their atmospheres vulnerable to erosion by the stellar wind… ▽ More Recently, four additional Earth-mass planets were discovered orbiting the nearby ultracool M8 dwarf TRAPPIST-1, making a remarkable total of seven planets with equilibrium temperatures compatible with the presence of liquid water on their surface. Temperate terrestrial planets around an M-dwarf orbit close to their parent star, rendering their atmospheres vulnerable to erosion by the stellar wind and energetic electromagnetic and particle radiation. Here, we use state-of-the-art 3D magnetohydrodynamic models to simulate the wind around TRAPPIST-1 and study the conditions at each planetary orbit. All planets experience a stellar wind pressure between $10^3$ and $10^5$ times the solar wind pressure on Earth. All orbits pass through wind pressure changes of an order of magnitude and most planets spend a large fraction of their orbital period in the sub-Alfvénic regime. For plausible planetary magnetic field strengths, all magnetospheres are greatly compressed and undergo much more dynamic change than that of the Earth. The planetary magnetic fields connect with the stellar radial field over much of the planetary surface, allowing direct flow of stellar wind particles onto the planetary atmosphere. These conditions could result in strong atmospheric stripping and evaporation and should be taken into account for any realistic assessment of the evolution and habitability of the TRAPPIST-1 planets. △ Less

Submitted 14 June, 2017; originally announced June 2017.

Comments: 12 pages, 5 figures, to appear in ApJ Letters

Understanding extreme quasar optical variability with CRTS: I. Major AGN flares

Authors: Matthew J. Graham, S. G. Djorgovski, Andrew J. Drake, Daniel Stern, Ashish A. Mahabal, Eilat Glikman, Steve Larson, Eric Christensen

Abstract: There is a large degree of variety in the optical variability of quasars and it is unclear whether this is all attributable to a single (set of) physical mechanism(s). We present the results of a systematic search for major flares in AGN in the Catalina Real-time Transient Survey as part of a broader study into extreme quasar variability. Such flares are defined in a quantitative manner as being a… ▽ More There is a large degree of variety in the optical variability of quasars and it is unclear whether this is all attributable to a single (set of) physical mechanism(s). We present the results of a systematic search for major flares in AGN in the Catalina Real-time Transient Survey as part of a broader study into extreme quasar variability. Such flares are defined in a quantitative manner as being atop of the normal, stochastic variability of quasars. We have identified 51 events from over 900,000 known quasars and high probability quasar candidates, typically lasting 900 days and with a median peak amplitude of $Δm = 1.25$ mag. Characterizing the flare profile with a Weibull distribution, we find that nine of the sources are well described by a single-point single-lens model. This supports the proposal by Lawrence et al. (2016) that microlensing is a plausible physical mechanism for extreme variability. However, we attribute the majority of our events to explosive stellar-related activity in the accretion disk: superluminous supernovae, tidal disruption events, and mergers of stellar mass black holes. △ Less

Submitted 9 June, 2017; originally announced June 2017.

Comments: 25 pages, 18 figures, accepted for publication by MNRAS

The role of three-dimensional transport in driving enhanced electron acceleration during magnetic reconnection

Authors: J. T. Dahlin, J. F. Drake, M. Swisdak

Abstract: Magnetic reconnection is an important driver of energetic particles in many astrophysical phenomena. Using kinetic particle-in-cell (PIC) simulations, we explore the impact of three-dimensional reconnection dynamics on the efficiency of particle acceleration. In two-dimensional systems, Alfvénic outflows expel energetic electrons into flux ropes where they become trapped and disconnected from acce… ▽ More Magnetic reconnection is an important driver of energetic particles in many astrophysical phenomena. Using kinetic particle-in-cell (PIC) simulations, we explore the impact of three-dimensional reconnection dynamics on the efficiency of particle acceleration. In two-dimensional systems, Alfvénic outflows expel energetic electrons into flux ropes where they become trapped and disconnected from acceleration regions. However, in three-dimensional systems these flux ropes develop axial structure that enables particles to leak out and return to acceleration regions. This requires a finite guide field so that particles may move quickly along the flux rope axis. We show that greatest energetic electron production occurs when the guide field is of the same order as the reconnecting component: large enough to facilitate strong transport, but not so large as to throttle the dominant Fermi mechanism responsible for efficient electron acceleration. This suggests a natural explanation for the envelope of electron acceleration during the impulsive phase of eruptive flares. △ Less

Submitted 1 June, 2017; originally announced June 2017.

Deep Chandra Survey of the Small Magellanic Cloud. II. Timing Analysis of X-ray Pulsars

Authors: JaeSub Hong, Vallia Antoniou, Andreas Zezas, Frank Haberl, Manami Sasaki, Jeremy Drake, Paul Plucinsky, Silas Laycock

Abstract: We report the timing analysis results of X-ray pulsars from a recent deep Chandra survey of the Small Magellanic Cloud (SMC). We have analyzed a total exposure of 1.4 Ms from 31 observations over a 1.2 deg$^2$ region in the SMC under a Chandra X-ray Visionary Program. Using the Lomb-Scargle and epoch folding techniques, we have detected periodic modulations from 20 pulsars and a new candidate puls… ▽ More We report the timing analysis results of X-ray pulsars from a recent deep Chandra survey of the Small Magellanic Cloud (SMC). We have analyzed a total exposure of 1.4 Ms from 31 observations over a 1.2 deg$^2$ region in the SMC under a Chandra X-ray Visionary Program. Using the Lomb-Scargle and epoch folding techniques, we have detected periodic modulations from 20 pulsars and a new candidate pulsar. The survey also covers 11 other pulsars with no clear sign of periodic modulation. The 0.5-8 keV X-ray luminosity ($L_X$) of the pulsars ranges from $10^{34}$ to $10^{37}$ erg s$^{-1}$ at 60 kpc. All the Chandra sources with $L_X$ $\gtrsim 4 \times 10^{35}$ erg s$^{-1}$ exhibit X-ray pulsations. The X-ray spectra of the SMC pulsars (and high mass X-ray binaries) are in general harder than those of the SMC field population. All but SXP~8.02 can be fitted by an absorbed power-law model with a photon index of $Γ$ $\lesssim$ 1.5. The X-ray spectrum of the known magnetar SXP~8.02 is better fitted with a two-temperature blackbody model. Newly measured pulsation periods of SXP~51.0, SXP~214 and SXP~701 are significantly different from the previous XMM-Newton and RXTE measurements. This survey provides a rich data set for energy-dependent pulse profile modeling. Six pulsars show an almost eclipse-like dip in the pulse profile. Phase-resolved spectral analysis reveals diverse spectral variation during pulsation cycle: e.g., for an absorbed power-law model, some exhibit an (anti)-correlation between absorption and X-ray flux, while others show more intrinsic spectral variation. △ Less

Submitted 26 May, 2017; originally announced May 2017.

Comments: 24 pages, 19 figures, 11 tables, submitted to ApJ

Extreme Variability in a Broad Absorption Line Quasar

Authors: Daniel Stern, Matthew J. Graham, Nahum Arav, S. G. Djorgovski, Carter Chamberlain, Aaron J. Barth, Ciro Donalek, Andrew J. Drake, Eilat Glikman, Hyunsung D. Jun, Ashish A. Mahabal, Charles C. Steidel

Abstract: CRTS J084133.15+200525.8 is an optically bright quasar at z=2.345 that has shown extreme spectral variability over the past decade. Photometrically, the source had a visual magnitude of V~17.3 between 2002 and 2008. Then, over the following five years, the source slowly brightened by approximately one magnitude, to V~16.2. Only ~1 in 10,000 quasars show such extreme variability, as quantified by t… ▽ More CRTS J084133.15+200525.8 is an optically bright quasar at z=2.345 that has shown extreme spectral variability over the past decade. Photometrically, the source had a visual magnitude of V~17.3 between 2002 and 2008. Then, over the following five years, the source slowly brightened by approximately one magnitude, to V~16.2. Only ~1 in 10,000 quasars show such extreme variability, as quantified by the extreme parameters derived for this quasar assuming a damped random walk model. A combination of archival and newly acquired spectra reveal the source to be an iron low-ionization broad absorption line (FeLoBAL) quasar with extreme changes in its absorption spectrum. Some absorption features completely disappear over the 9 years of optical spectra, while other features remain essentially unchanged. We report the first definitive redshift for this source, based on the detection of broad H-alpha in a Keck/MOSFIRE spectrum. Absorption systems separated by several 1000 km/s in velocity show coordinated weakening in the depths of their troughs as the continuum flux increases. We interpret the broad absorption line variability to be due to changes in photoionization, rather than due to motion of material along our line of sight. This source highlights one sort of rare transition object that astronomy will now be finding through dedicated time-domain surveys. △ Less

Submitted 12 April, 2017; originally announced April 2017.

Comments: 6 pages, 4 figures; accepted for publication in ApJ

Turbulence in Three-Dimensional Simulations of Magnetopause Reconnection

Authors: L. Price, M. Swisdak, J. F. Drake, J. L. Burch, P. A. Cassak, R. E. Ergun

Abstract: We present detailed analysis of the turbulence observed in three-dimensional particle-in-cell simulations of magnetic reconnection at the magnetopause. The parameters are representative of an electron diffusion region encounter of the Magnetospheric Multiscale (MMS) mission. The turbulence is found to develop around both the magnetic X line and separatrices, is electromagnetic in nature, is charac… ▽ More We present detailed analysis of the turbulence observed in three-dimensional particle-in-cell simulations of magnetic reconnection at the magnetopause. The parameters are representative of an electron diffusion region encounter of the Magnetospheric Multiscale (MMS) mission. The turbulence is found to develop around both the magnetic X line and separatrices, is electromagnetic in nature, is characterized by a wave vector $k$ given by $kρ_e\sim(m_eT_e/m_iT_i)^{0.25}$ with $ρ_e$ the electron Larmor radius, and appears to have the ion pressure gradient as its source of free energy. Taken together, these results suggest the instability is a variant of the lower hybrid drift instability. The turbulence produces electric field fluctuations in the out-of-plane direction (the direction of the reconnection electric field) with an amplitude of around $\pm 10$~mV/m, which is much greater than the reconnection electric field of around $0.1$~mV/m. Such large values of the out-of-plane electric field have been identified in the MMS data. The turbulence in the simulations controls the scale lengths of the density profile and current layers in asymmetric reconnection, driving them closer to $\sqrt{ρ_eρ_i}$ than the $ρ_e$ or $d_e$ scalings seen in 2-D reconnection simulations, and produces significant anomalous resistivity and viscosity in the electron diffusion region. △ Less

Submitted 3 May, 2018; v1 submitted 5 April, 2017; originally announced April 2017.

Comments: 11 pages, 10 figures

Journal ref: Journal of Geophysical Research: Space Physics, 122, 11086--11099 (2017)

A Comprehensive Library of X-ray Pulsars in the Small Magellanic Cloud: Time Evolution of their Luminosities and Spin Periods

Authors: J. Yang, S. G. T. Laycock, D. M. Christodoulou, S. Fingerman, M. J. Coe, J. J. Drake

Abstract: We have collected and analyzed the complete archive of {\itshape XMM-Newton\} (116), {\itshape Chandra\} (151), and {\itshape RXTE\} (952) observations of the Small Magellanic Cloud (SMC), spanning 1997-2014. The resulting observational library provides a comprehensive view of the physical, temporal and statistical properties of the SMC pulsar population across the luminosity range of… ▽ More We have collected and analyzed the complete archive of {\itshape XMM-Newton\} (116), {\itshape Chandra\} (151), and {\itshape RXTE\} (952) observations of the Small Magellanic Cloud (SMC), spanning 1997-2014. The resulting observational library provides a comprehensive view of the physical, temporal and statistical properties of the SMC pulsar population across the luminosity range of $L_X= 10^{31.2}$--$10^{38}$~erg~s$^{-1}$. From a sample of 67 pulsars we report $\sim$1654 individual pulsar detections, yielding $\sim$1260 pulse period measurements. Our pipeline generates a suite of products for each pulsar detection: spin period, flux, event list, high time-resolution light-curve, pulse-profile, periodogram, and spectrum. Combining all three satellites, we generated complete histories of the spin periods, pulse amplitudes, pulsed fractions and X-ray luminosities. Some pulsars show variations in pulse period due to the combination of orbital motion and accretion torques. Long-term spin-up/down trends are seen in 12/11 pulsars respectively, pointing to sustained transfer of mass and angular momentum to the neutron star on decadal timescales. Of the sample 30 pulsars have relatively very small spin period derivative and may be close to equilibrium spin. The distributions of pulse-detection and flux as functions of spin-period provide interesting findings: mapping boundaries of accretion-driven X-ray luminosity, and showing that fast pulsars ($P<$10 s) are rarely detected, which yet are more prone to giant outbursts. Accompanying this paper is an initial public release of the library so that it can be used by other researchers. We intend the library to be useful in driving improved models of neutron star magnetospheres and accretion physics. △ Less

Submitted 15 March, 2017; originally announced March 2017.

Comments: 17 pages, 11 + 58 (appendix) figures. To appear in the Astrophysical Journal Supplement

A multi-observatory database of X-ray pulsars in the Magellanic Clouds

Authors: J. Yang, S. G. T. Laycock, J. J. Drake, M. J. Coe, S. Fingerman, J. Hong, V. Antoniou, A. Zezas

Abstract: Using hundreds of XMM-Newton and Chandra archival observations and nearly a thousand RXTE observations, we have generated a comprehensive library of the known pulsars in the Small and Large Magellanic Clouds (SMC, LMC). The pulsars are detected multiple times across the full parameter spaces of X-ray luminosity ($L_X= 10^{31-38}$~erg/s) and spin period ( P$<$1s -- P$>$1000s) and the library enable… ▽ More Using hundreds of XMM-Newton and Chandra archival observations and nearly a thousand RXTE observations, we have generated a comprehensive library of the known pulsars in the Small and Large Magellanic Clouds (SMC, LMC). The pulsars are detected multiple times across the full parameter spaces of X-ray luminosity ($L_X= 10^{31-38}$~erg/s) and spin period ( P$<$1s -- P$>$1000s) and the library enables time-domain studies at a range of energy scales. The high time-resolution and sensitivity of the EPIC cameras are complemented by the angular resolution of Chandra and the regular monitoring of RXTE. Our processing %$\sim$15 year pipeline uses the latest calibration files and software to generate a suite of useful products for each pulsar detection: event lists, high time-resolution light curves, periodograms, spectra, and complete histories of $\dot{P}$, the pulsed fraction, etc., in the broad (0.2-12 keV), soft (0.2-2 keV), and hard (2-12 keV) energy bands. After combining the observations from these telescopes, we found that 28 pulsars show long-term spin up and 25 long-term spin down. We also used the faintest and brightest sources to map out the lower and upper boundaries of accretion-powered X-ray emission: the propeller line and the Eddington line, respectively. We are in the process of comparing the observed pulse profiles to geometric models of X-ray emission in order to constrain the physical parameters of the pulsars. Finally we are preparing a public release of the library so that it can be used by others in the astronomical community. △ Less

Submitted 14 March, 2017; originally announced March 2017.

Comments: 6 pages, 6 figures

Enabling New ALMA Science with Improved Support for Time-Domain Observations

Authors: Kate D. Alexander, Edo Berger, Geoff Bower, Sarah Casewell, S. Brad Cenko, Shami Chatterjee, Ilse Cleeves, Jim Cordes, Jeremy Drake, Maria Drout, Trent Dupuy, Tarraneh Eftekhari, Giovanni Fazio, Wen-fai Fong, James Guillochon, Mark Gurwell, Michael Johnson, Tomasz Kaminski, Albert Kong, Tanmoy Laskar, Casey Law, Stuart P. Littlefair, Meredith MacGregor, W. Peter Maksym, Lynn Matthews , et al. (12 additional authors not shown)

Abstract: While the Atacama Large Millimeter/submillimeter Array (ALMA) is a uniquely powerful telescope, its impact in certain fields of astrophysics has been limited by observatory policies rather than the telescope's innate technical capabilities. In particular, several observatory policies present challenges for observations of variable, mobile, and/or transient sources --- collectively referred to here… ▽ More While the Atacama Large Millimeter/submillimeter Array (ALMA) is a uniquely powerful telescope, its impact in certain fields of astrophysics has been limited by observatory policies rather than the telescope's innate technical capabilities. In particular, several observatory policies present challenges for observations of variable, mobile, and/or transient sources --- collectively referred to here as "time-domain" observations. In this whitepaper we identify some of these policies, describe the scientific applications they impair, and suggest changes that would increase ALMA's science impact in Cycle 6 and beyond. Parties interested in time-domain science with ALMA are encouraged to join the ALMA Time-domain Special Interest Group (ATSIG) by signing up for the ATSIG mailing list at https://groups.google.com/group/alma-td-sig . △ Less

Submitted 14 March, 2017; originally announced March 2017.

Comments: 9 pages; whitepaper submitted to the ALMA Science Advisory Council; corresponding author P. K. G. Williams (pwilliams@cfa.harvard.edu)

The Twist of the Draped Interstellar Magnetic Field Ahead of the Heliopause: A Magnetic Reconnection Driven Rotational Discontinuity

Authors: M. Opher, J. F. Drake, M. Swisdak, B. Zieger, G. Toth

Abstract: Based on the difference between the orientation of the interstellar $B_{ISM}$ and the solar magnetic fields, there was an expectation that the magnetic field direction would rotate dramatically across the heliopause (HP). However, the Voyager 1 spacecraft measured very little rotation across the HP. Previously we showed that the $B_{ISM}$ twists as it approaches the HP and acquires a strong T comp… ▽ More Based on the difference between the orientation of the interstellar $B_{ISM}$ and the solar magnetic fields, there was an expectation that the magnetic field direction would rotate dramatically across the heliopause (HP). However, the Voyager 1 spacecraft measured very little rotation across the HP. Previously we showed that the $B_{ISM}$ twists as it approaches the HP and acquires a strong T component (East-West). Here we establish that reconnection in the eastern flank of the heliosphere is responsible for the twist. On the eastern flank the solar magnetic field has twisted into the positive N direction and reconnects with the Southward pointing component of the $B_{ISM}$. Reconnection drives a rotational discontinuity (RD) that twists the $B_{ISM}$ into the -T direction and propagates upstream in the interstellar medium towards the nose. The consequence is that the N component of $B_{ISM}$ is reduced in a finite width band upstream of the HP. Voyager 1 currently measures angles ($δ=sin^{-1}(B_{N}/B)$) close to solar values. We present MHD simulations to support this scenario, suppressing reconnection in the nose region while allowing it in the flanks, consistent with recent ideas about reconnection suppression from diamagnetic drifts. The jump in plasma $β$ (the plasma to magnetic pressure) across the nose of HP is much greater than in the flanks because the heliosheath $β$ is greater there than in the flanks. Large-scale reconnection is therefore suppressed in the nose but not at the flanks. Simulation data suggest that $B_{ISM}$ will return to its pristine value $10-15~AU$ past the HP. △ Less

Submitted 20 February, 2017; originally announced February 2017.

Comments: 19 pages, 5 figures, submitted

Journal ref: Astrophysical Journal Letters 839, Number 1, 2017

The formation of magnetic depletions and flux annihilation due to reconnection in the heliosheath

Authors: J. F. Drake, M. Swisdak, M. Opher, J. D. Richardson

Abstract: The misalignment of the solar rotation axis and the magnetic axis of the Sun produces a periodic reversal of the Parker spiral magnetic field and the sectored solar wind. The compression of the sectors is expected to lead to reconnection in the heliosheath (HS). We present particle-in-cell simulations of the sectored HS that reflect the plasma environment along the Voyager 1 and 2 trajectories, sp… ▽ More The misalignment of the solar rotation axis and the magnetic axis of the Sun produces a periodic reversal of the Parker spiral magnetic field and the sectored solar wind. The compression of the sectors is expected to lead to reconnection in the heliosheath (HS). We present particle-in-cell simulations of the sectored HS that reflect the plasma environment along the Voyager 1 and 2 trajectories, specifically including unequal positive and negative azimuthal magnetic flux as seen in the Voyager data \citep{Burlaga03}. Reconnection proceeds on individual current sheets until islands on adjacent current layers merge. At late time bands of the dominant flux survive, separated by bands of deep magnetic field depletion. The ambient plasma pressure supports the strong magnetic pressure variation so that pressure is anti-correlated with magnetic field strength. There is little variation in the magnetic field direction across the boundaries of the magnetic depressions. At irregular intervals within the magnetic depressions are long-lived pairs of magnetic islands where the magnetic field direction reverses so that spacecraft data would reveal sharp magnetic field depressions with only occasional crossings with jumps in magnetic field direction. This is typical of the magnetic field data from the Voyager spacecraft \citep{Burlaga11,Burlaga16}. Voyager 2 data reveals that fluctuations in the density and magnetic field strength are anti-correlated in the sector zone as expected from reconnection but not in unipolar regions. The consequence of the annihilation of subdominant flux is a sharp reduction in the "number of sectors" and a loss in magnetic flux as documented from the Voyager 1 magnetic field and flow data \citep{Richardson13}. △ Less

Submitted 6 February, 2017; originally announced February 2017.

Exploring impulsive solar magnetic energy release and particle acceleration with focused hard X-ray imaging spectroscopy

Authors: Steven Christe, Samuel Krucker, Lindsay Glesener, Albert Shih, Pascal Saint-Hilaire, Amir Caspi, Joel Allred, Marina Battaglia, Bin Chen, James Drake, Brian Dennis, Dale Gary, Szymon Gburek, Keith Goetz, Brian Grefenstette, Mikhail Gubarev, Iain Hannah, Gordon Holman, Hugh Hudson, Andrew Inglis, Jack Ireland, Shinosuke Ishikawa, James Klimchuk, Eduard Kontar, Adam Kowalski , et al. (15 additional authors not shown)

Abstract: How impulsive magnetic energy release leads to solar eruptions and how those eruptions are energized and evolve are vital unsolved problems in Heliophysics. The standard model for solar eruptions summarizes our current understanding of these events. Magnetic energy in the corona is released through drastic restructuring of the magnetic field via reconnection. Electrons and ions are then accelerate… ▽ More How impulsive magnetic energy release leads to solar eruptions and how those eruptions are energized and evolve are vital unsolved problems in Heliophysics. The standard model for solar eruptions summarizes our current understanding of these events. Magnetic energy in the corona is released through drastic restructuring of the magnetic field via reconnection. Electrons and ions are then accelerated by poorly understood processes. Theories include contracting loops, merging magnetic islands, stochastic acceleration, and turbulence at shocks, among others. Although this basic model is well established, the fundamental physics is poorly understood. HXR observations using grazing-incidence focusing optics can now probe all of the key regions of the standard model. These include two above-the-looptop (ALT) sources which bookend the reconnection region and are likely the sites of particle acceleration and direct heating. The science achievable by a direct HXR imaging instrument can be summarized by the following science questions and objectives which are some of the most outstanding issues in solar physics (1) How are particles accelerated at the Sun? (1a) Where are electrons accelerated and on what time scales? (1b) What fraction of electrons is accelerated out of the ambient medium? (2) How does magnetic energy release on the Sun lead to flares and eruptions? A Focusing Optics X-ray Solar Imager (FOXSI) instrument, which can be built now using proven technology and at modest cost, would enable revolutionary advancements in our understanding of impulsive magnetic energy release and particle acceleration, a process which is known to occur at the Sun but also throughout the Universe. △ Less

Submitted 3 January, 2017; originally announced January 2017.

Comments: Next Generation Solar Physics Mission white paper, 2 figures

Detection of quasars in the time domain

Authors: Matthew J. Graham, S. G. Djorgovski, Daniel J. Stern, Andrew J. Drake, Ashish Mahabal

Abstract: The time domain is the emerging forefront of astronomical research with new facilities and instruments providing unprecedented amounts of data on the temporal behavior of astrophysical populations. Dealing with the size and complexity of this requires new techniques and methodologies. Quasars are an ideal work set for developing and applying these: they vary in a detectable but not easily quantifi… ▽ More The time domain is the emerging forefront of astronomical research with new facilities and instruments providing unprecedented amounts of data on the temporal behavior of astrophysical populations. Dealing with the size and complexity of this requires new techniques and methodologies. Quasars are an ideal work set for developing and applying these: they vary in a detectable but not easily quantifiable manner whose physical origins are poorly understood. In this paper, we will review how quasars are identified by their variability and how these techniques can be improved, what physical insights into their variability can be gained from studying extreme examples of variability, and what approaches can be taken to increase the number of quasars known. These will demonstrate how astroinformatics is essential to discovering and understanding this important population. △ Less

Submitted 21 December, 2016; originally announced December 2016.

Comments: 10 pages, 6 figures, IAU Symposium 325, "Astroinformatics"

A fork in the Sagittarius trailing debris

Authors: Camila Navarrete, Vasily Belokurov, Sergey E. Koposov, Mike Irwin, Márcio Catelan, Sonia Duffau, Andrew J. Drake

Abstract: We take advantage of the deep and wide coverage of the VST ATLAS survey to study the line-of-sight structure of the Sagittarius stellar stream in the Southern hemisphere, only ~40° away from the progenitor. We use photometrically selected Sub-Giant Branch (SGB) stars to reveal a complex debris morphology of the trailing arm and detect at least two clear peaks in the SGB distance modulus distributi… ▽ More We take advantage of the deep and wide coverage of the VST ATLAS survey to study the line-of-sight structure of the Sagittarius stellar stream in the Southern hemisphere, only ~40° away from the progenitor. We use photometrically selected Sub-Giant Branch (SGB) stars to reveal a complex debris morphology of the trailing arm and detect at least two clear peaks in the SGB distance modulus distribution. The separation between the two line-of-sight components is at least 5 kpc at the edge of the VST ATLAS footprint, but appears to change along the stream, which allows us to conclude that these detections correspond to two physically independent stellar structures, rather than a mix of co-distant stellar populations within a single stream. Our discovery of a fork in the Sgr trailing arm is verified using Blue Horizontal Branch stars and our distance measurements are calibrated using RR Lyrae stars from the Catalina Real-Time Transient Survey. Comparing with numerical simulations of the Sgr dwarf disruption, the more distant of the two components in the fork matches perfectly with the track of the trailing debris. However, no obvious counterpart exists in the simulation for the closer line-of-sight component. △ Less

Submitted 20 December, 2016; originally announced December 2016.

Comments: 15 pages, 12 figures, accepted for publication in MNRAS

Giant coronal loops dominate the quiescent X-ray emission in rapidly rotating M stars

Authors: O. Cohen, R. Yadav, C. Garraffo, S. H. Saar, S. J. Wolk, V. L. Kashyap, J. J. Drake, I. Pillitteri

Abstract: Observations indicate that magnetic fields in rapidly rotating stars are very strong, on both small and large scales. What is the nature of the resulting corona? Here we seek to shed some light on this question. We use the results of an anelastic dynamo simulation of a rapidly rotating fully-convective M-star to drive a physics-based model for the stellar corona. We find that due to the several ki… ▽ More Observations indicate that magnetic fields in rapidly rotating stars are very strong, on both small and large scales. What is the nature of the resulting corona? Here we seek to shed some light on this question. We use the results of an anelastic dynamo simulation of a rapidly rotating fully-convective M-star to drive a physics-based model for the stellar corona. We find that due to the several kilo Gauss large-scale magnetic fields at high latitudes, the corona and its X-ray emission are dominated by star-size large hot loops, while the smaller, underlying colder loops are not visible much in the X-ray. Based on this result we propose that, in rapidly rotating stars, emission from such coronal structures dominates the quiescent, cooler but saturated X-ray emission. △ Less

Submitted 7 November, 2016; originally announced November 2016.

Comments: 9 pages, 4 figures, accepted to ApJ

Origin of asymmetries in X-ray emission lines from the blast wave of the 2014 outburst of nova V745 Sco

Authors: S. Orlando, J. J. Drake, M. Miceli

Abstract: The symbiotic nova V745 Sco was observed in outburst on 2014 February 6. Its observations by the Chandra X-ray Observatory at days 16 and 17 have revealed a spectrum characterized by asymmetric and blue-shifted emission lines. Here we investigate the origin of these asymmetries through three-dimensional hydrodynamic simulations describing the outburst during the first 20 days of evolution. The mod… ▽ More The symbiotic nova V745 Sco was observed in outburst on 2014 February 6. Its observations by the Chandra X-ray Observatory at days 16 and 17 have revealed a spectrum characterized by asymmetric and blue-shifted emission lines. Here we investigate the origin of these asymmetries through three-dimensional hydrodynamic simulations describing the outburst during the first 20 days of evolution. The model takes into account thermal conduction and radiative cooling and assumes a blast wave propagates through an equatorial density enhancement. From the simulations, we synthesize the X-ray emission and derive the spectra as they would be observed with Chandra. We find that both the blast wave and the ejecta distribution are efficiently collimated in polar directions due to the presence of the equatorial density enhancement. The majority of the X-ray emission originates from the interaction of the blast with the equatorial density enhancement and is concentrated on the equatorial plane as a ring-like structure. Our "best-fit" model requires a mass of ejecta in the outburst $M_{ej} \approx 3\times 10^{-7}\,M_{\odot}$ and an explosion energy $E_b \approx 3 \times 10^{43}$ erg and reproduces the distribution of emission measure vs temperature and the evolution of shock velocity and temperature inferred from the observations. The model predicts asymmetric and blue-shifted line profiles similar to those observed and explains their origin as due to substantial X-ray absorption of red-shifted emission by ejecta material. The comparison of predicted and observed Ne and O spectral line ratios reveals no signs of strong Ne enhancement and suggests the progenitor is a CO white dwarf. △ Less

Submitted 18 October, 2016; originally announced October 2016.

Comments: 16 pages, 17 Figures; accepted for publication on MNRAS

Stellar flares and the dark energy of CMEs

Authors: Jeremy J. Drake, Ofer Cohen, Cecilia Garraffo, Vinay Kashyap

Abstract: Flares we observe on stars in white light, UV or soft X-rays are probably harbingers of coronal mass ejections (CMEs). If we use the Sun as a guide, large stellar flares will dissipate two orders of magnitude less X-ray radiative energy than the kinetic energy in the associated CME. Since coronal emission on active stars appears to be dominated by flare activity, CMEs pose a quandary for understan… ▽ More Flares we observe on stars in white light, UV or soft X-rays are probably harbingers of coronal mass ejections (CMEs). If we use the Sun as a guide, large stellar flares will dissipate two orders of magnitude less X-ray radiative energy than the kinetic energy in the associated CME. Since coronal emission on active stars appears to be dominated by flare activity, CMEs pose a quandary for understanding the fraction of their energy budget stars can spend on magnetic activity. One answer is magnetic suppression of CMEs, in which the strong large-scale fields of active stars entrap and prevent CMEs unless their free energy exceeds a critical value. The CME-less flaring active region NOAA 2192 presents a possible solar analogue of this. Monster CMEs will still exist, and have the potential to ravage planetary atmospheres. △ Less

Submitted 17 October, 2016; originally announced October 2016.

Comments: 6 pages, to appear in Proceedings of IAU Symposium No. 320

Optical, UV, and X-Ray Evidence for a 7-Year Stellar Cycle in Proxima Centauri

Authors: B. J. Wargelin, S. H. Saar, G. Pojmański, J. J. Drake, V. L. Kashyap

Abstract: Stars of stellar type later than about M3.5 are believed to be fully convective and therefore unable to support magnetic dynamos like the one that produces the 11-year solar cycle. Because of their intrinsic faintness, very few late M stars have undergone long-term monitoring to test this prediction, which is critical to our understanding of magnetic field generation in such stars. Magnetic activi… ▽ More Stars of stellar type later than about M3.5 are believed to be fully convective and therefore unable to support magnetic dynamos like the one that produces the 11-year solar cycle. Because of their intrinsic faintness, very few late M stars have undergone long-term monitoring to test this prediction, which is critical to our understanding of magnetic field generation in such stars. Magnetic activity is also of interest as the driver of UV and X-ray radiation, as well as energetic particles and stellar winds, that affect the atmospheres of close-in planets that lie within habitable zones, such as the recently discovered Proxima b. We report here on several years of optical, UV, and X-ray observations of Proxima Centauri (GJ 551; dM5.5e): 15 years of ASAS photometry in the V band (1085 nights) and 3 years in the I band (196 nights), 4 years of Swift XRT and UVOT observations (more than 120 exposures), and 9 sets of X-ray observations from other X-ray missions (ASCA, XMM-Newton, and three Chandra instruments) spanning 22 years. We confirm previous reports of an 83-day rotational period and find strong evidence for a 7-year stellar cycle, along with indications of differential rotation at about the solar level. X-ray/UV intensity is anti-correlated with optical V-band brightness for both rotational and cyclical variations. From comparison with other stars observed to have X-ray cycles we deduce a simple empirical relationship between X-ray cyclic modulation and Rossby number, and we also present Swift UV grism spectra covering 2300-6000 Å. △ Less

Submitted 11 October, 2016; originally announced October 2016.

Comments: 19 pages with 15 figures. Submitted to MNRAS 2016 Jul 1; revised version accepted 2016 Oct 3 and published online 2016 Oct 8. Available at http://dx.doi.org/10.1093/mnras/stw2570

The Space Weather of Proxima Centauri b

Authors: Cecilia Garraffo, Jeremy J. Drake, Ofer Cohen

Abstract: A planet orbiting in the "habitable zone" of our closest neighboring star, Proxima Centauri, has recently been discovered, and the next natural question is whether or not Proxima b is "habitable". Stellar winds are likely a source of atmospheric erosion that could be particularly severe in the case of M dwarf habitable zone planets that reside close to their parent star. Here we study the stellar… ▽ More A planet orbiting in the "habitable zone" of our closest neighboring star, Proxima Centauri, has recently been discovered, and the next natural question is whether or not Proxima b is "habitable". Stellar winds are likely a source of atmospheric erosion that could be particularly severe in the case of M dwarf habitable zone planets that reside close to their parent star. Here we study the stellar wind conditions that Proxima b experiences over its orbit. We construct 3-D MHD models of the wind and magnetic field around Proxima Centauri using a surface magnetic field map for a star of the same spectral type and scaled to match the observed ~ 600 G surface magnetic field strength of Proxima. We examine the wind conditions and dynamic pressure over different plausible orbits that sample the constrained parameters of the orbit of Proxima b. For all the parameter space explored, the planet is subject to stellar wind pressures of more than 2000 times those experienced by Earth from the solar wind. During an orbit, Proxima b is also subject to pressure changes of 1 to 3 orders of magnitude on timescales of a day. Its magnetopause standoff distance consequently undergoes sudden and periodic changes by a factor of 2 to 5. Proxima b will traverse the interplanetary current sheet twice each orbit, and likely crosses into regions of subsonic wind quite frequently. These effects should be taken into account in any physically realistic assessment or prediction of its atmospheric reservoir, characteristics and loss. △ Less

Submitted 9 December, 2016; v1 submitted 28 September, 2016; originally announced September 2016.

Comments: 14 pages, 5 figures, ApJ Letters, 833, 1

Simulating the environment around planet-hosting stars - II. Stellar winds and inner astrospheres

Authors: J. D. Alvarado-Gómez, G. A. J. Hussain, O. Cohen, J. J. Drake, C. Garraffo, J. Grunhut, T. I. Gombosi

Abstract: We present the results of a comprehensive numerical simulation of the environment around three exoplanet-host stars (HD 1237, HD 22049, and HD 147513). Our simulations consider one of the latest models currently used for space weather studies in the Heliosphere. Large-scale magnetic field maps, recovered with two implementations of the tomographic technique of Zeeman-Doppler imaging, serve to driv… ▽ More We present the results of a comprehensive numerical simulation of the environment around three exoplanet-host stars (HD 1237, HD 22049, and HD 147513). Our simulations consider one of the latest models currently used for space weather studies in the Heliosphere. Large-scale magnetic field maps, recovered with two implementations of the tomographic technique of Zeeman-Doppler imaging, serve to drive steady-state solutions in each system. This paper contains the description of the stellar wind and inner astrosphere, while the coronal structure was previously discussed in Alvarado-Gómez et al. (2016). The analysis includes the magneto-hydrodynamical properties of the stellar wind, the associated mass and angular momentum loss rates, as well as the topology of the astrospheric current sheet in each system. A systematic comparison among the considered cases is performed, including two reference solar simulations covering activity minimum and maximum. For HD 1237, we investigate the interactions between the structure of the developed stellar wind, and a possible magnetosphere around the Jupiter-mass planet in this system. We find that the process of particle injection into the planetary atmosphere is dominated by the density distribution rather than velocity profile of the stellar wind. In this context, we predict a maximum exoplanetary radio emission of 12 mJy at 40 MHz in this system, assuming the crossing of a high-density streamer during periastron passage. Furthermore, in combination with the analysis performed in Alvarado-Gómez et al. (2016), we obtain for the first time a fully simulated mass loss-activity relation, which is compared and discussed in the context of the relation based on astrospheric detections proposed by Wood et al. (2005a). Finally, we provide a characterisation of the 3D properties of the stellar wind of these systems, at the inner edges of their habitable zones. △ Less

Submitted 28 July, 2016; originally announced July 2016.

Comments: 19 pages, 12 Figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 594, A95 (2016)

Solar-type dynamo behaviour in fully convective stars without a tachocline

Authors: Nicholas J. Wright, Jeremy J. Drake

Abstract: In solar-type stars (with radiative cores and convective envelopes), the magnetic field powers star spots, flares and other solar phenomena, as well as chromospheric and coronal emission at ultraviolet to X-ray wavelengths. The dynamo responsible for generating the field depends on the shearing of internal magnetic fields by differential rotation. The shearing has long been thought to take place i… ▽ More In solar-type stars (with radiative cores and convective envelopes), the magnetic field powers star spots, flares and other solar phenomena, as well as chromospheric and coronal emission at ultraviolet to X-ray wavelengths. The dynamo responsible for generating the field depends on the shearing of internal magnetic fields by differential rotation. The shearing has long been thought to take place in a boundary layer known as the tachocline between the radiative core and the convective envelope. Fully convective stars do not have a tachocline and their dynamo mechanism is expected to be very different, although its exact form and physical dependencies are not known. Here we report observations of four fully convective stars whose X-ray emission correlates with their rotation periods in the same way as in Sun-like stars. As the X-ray activity - rotation relationship is a well-established proxy for the behaviour of the magnetic dynamo, these results imply that fully convective stars also operate a solar-type dynamo. The lack of a tachocline in fully convective stars therefore suggests that this is not a critical ingredient in the solar dynamo and supports models in which the dynamo originates throughout the convection zone. △ Less

Submitted 26 July, 2016; originally announced July 2016.

Comments: 6 pages, 1 figure. Accepted for publication in Nature (28 July 2016). Author's version, including Methods

The Missing Magnetic Morphology Term in Stellar Rotation Evolution

Authors: Cecilia Garraffo, Jeremy J. Drake, Ofer Cohen

Abstract: This study examines the relationship between magnetic field complexity and mass and angular momentum losses. Observations of open clusters have revealed a bimodal distribution of the rotation periods of solar-like stars that has proven difficult to explain under the existing rubric of magnetic braking. Recent studies suggest that magnetic complexity can play an important role in controlling stella… ▽ More This study examines the relationship between magnetic field complexity and mass and angular momentum losses. Observations of open clusters have revealed a bimodal distribution of the rotation periods of solar-like stars that has proven difficult to explain under the existing rubric of magnetic braking. Recent studies suggest that magnetic complexity can play an important role in controlling stellar spin-down rates. However, magnetic morphology is still neglected in most rotation evolution models due to the difficulty of properly accounting for its effects on wind driving and angular momentum loss. Using state-of-the-art magnetohydrodynamical magnetized wind simulations we study the effect that different distributions of the magnetic flux at different levels of geometrical complexity have on mass and angular momentum loss rates. Angular momentum loss rates depend strongly on the level of complexity of the field but are independent of the way this complexity is distributed. We deduce the analytical terms representing the magnetic field morphology dependence of mass and angular momentum loss rates. We also define a parameter that best represents complexity for real stars. As a test, we use these analytical methods to estimate mass and angular momentum loss rates for $8$ stars with observed magnetograms and compare them to the simulated results. Magnetic field complexity provides a natural physical basis for stellar rotation evolution models requiring a rapid transition between weak and strong spin-down modes. △ Less

Submitted 20 July, 2016; originally announced July 2016.

Comments: 13 pages, 6 figures, accepted for publication in A&A

Parallel electric fields are inefficient drivers of energetic electrons in magnetic reconnection

Authors: J. T. Dahlin, J. F. Drake, M. Swisdak

Abstract: We present two-dimensional kinetic simulations, with a broad range of initial guide fields, that isolate the role of parallel electric fields ($E_\parallel$) in energetic electron production during collisionless magnetic reconnection. In the strong guide field regime, $E_\parallel$ drives essentially all of the electron energy gain, yet fails to generate an energetic component. We suggest that thi… ▽ More We present two-dimensional kinetic simulations, with a broad range of initial guide fields, that isolate the role of parallel electric fields ($E_\parallel$) in energetic electron production during collisionless magnetic reconnection. In the strong guide field regime, $E_\parallel$ drives essentially all of the electron energy gain, yet fails to generate an energetic component. We suggest that this is due to the weak energy scaling of particle acceleration from $E_\parallel$ compared to that of a Fermi-type mechanism responsible for energetic electron production in the weak guide-field regime. This result has important implications for energetic electron production in astrophysical systems and reconnection-driven dissipation in turbulence. △ Less

Submitted 20 December, 2016; v1 submitted 13 July, 2016; originally announced July 2016.