A prevailing theory for the interstellar production of complex organic molecules (COMs) involves formation on warm dust-grain surfaces, via the diffusion and reaction of radicals produced through grain-surface photodissociation of stable molecules. However, some gas-phase O-bearing COMs, notably acetaldehyde(CH$_3$CHO), methyl formate(CH$_3$OCHO), and dimethyl ether(CH$_3$OCH$_3$), are now observed at very low temperatures, challenging the warm scenario. Here, we introduce a selection of new non-diffusive mechanisms into an astrochemical model, to account for the failure of the standard diffusive picture and to provide a more generalized scenario of COM formation on interstellar grains. New generic rate formulations are provided for cases where: (i) radicals are formed by reactions occurring close to another reactant, producing an immediate follow-on reaction; (ii) radicals are formed in an excited state, allowing them to overcome activation barriers to react with nearby stable molecules; (iii) radicals are formed through photo-dissociation close to a reaction partner, followed by immediate reaction. Each process occurs without the diffusion of large radicals. The new mechanisms significantly enhance cold COM abundances, successfully reproducing key observational results for prestellar core L1544. H-abstraction from grain-surface COMs, followed by recombination, plays a crucial role in amplifying chemical desorption into the gas phase. The UV-induced chemistry produces significant COM abundances in the bulk ices, which are retained on the grains and may persist to later stages. O$_2$ is also formed strongly in the mantle though photolysis, suggesting cometary O$_2$ could indeed be interstellar.

The recent discovery of an Earth-sized planet (TOI-700 d) in the habitable zone of an early-type M-dwarf by the Transiting Exoplanet Survey Satellite constitutes an important advance. In this Letter, we assess the feasibility of this planet to retain an atmosphere -- one of the chief ingredients for surface habitability -- over long timescales by employing state-of-the-art magnetohydrodynamic models to simulate the stellar wind and the associated rates of atmospheric escape. We take two major factors into consideration, namely, the planetary atmospheric composition and magnetic field. In all cases, we determine that the atmospheric ion escape rates are potentially a few orders of magnitude higher than the inner Solar system planets, but TOI-700 d is nevertheless capable of retaining a $1$ bar atmosphere over gigayear timescales for certain regions of the parameter space. The simulations show that the unmagnetized TOI-700 d with a 1 bar Earth-like atmosphere could be stripped away rather quickly ($<$ 1 gigayear), while the unmagnetized TOI-700 d with a 1 bar CO$_2$-dominated atmosphere could persist for many billions of years; we find that the magnetized Earth-like case falls in between these two scenarios. We also discuss the prospects for detecting radio emission of the planet (thereby constraining its magnetic field) and discerning the presence of an atmosphere.

Adopting the Standard Halo Model (SHM) of an isotropic Maxwellian velocity distribution for dark matter (DM) particles in the Galaxy, the most stringent current constraints on their spin-dependent scattering cross-section with nucleons come from the IceCube neutrino observatory and the PICO-60 C$_3$F$_8$ superheated bubble chamber experiments. The former is sensitive to high energy neutrinos from the self-annihilation of DM particles captured in the Sun, while the latter looks for nuclear recoil events from DM scattering off nucleons. Although slower DM particles are more likely to be captured by the Sun, the faster ones are more likely to be detected by PICO. Recent N-body simulations suggest significant deviations from the SHM for the smooth halo component of the DM, while observations hint at a dominant fraction of the local DM being in substructures. We use the method of Ferrer et al. (2015) to exploit the complementarity between the two approaches and derive conservative constraints on DM-nucleon scattering. Our results constrain $\sigma_{\mathrm{SD}} \lesssim 3 \times 10^{-39} \mathrm{cm}^2$ (6 $ \times 10^{-38} \mathrm{cm}^2$) at $\gtrsim 90\%$ C.L. for a DM particle of mass 1~TeV annihilating into $\tau^+ \tau^-$ ($b\bar{b}$) with a local density of $\rho_{\mathrm{DM}} = 0.3~\mathrm{ GeV/cm}^3$. The constraints scale inversely with $\rho_{\mathrm{DM}}$ and are independent of the DM velocity distribution.

Alfv\'en-wave-driven 3D magnetohydrodynamic (MHD) models, which are increasingly used to predict stellar wind properties, contain unconstrained parameters and rely on low-resolution stellar magnetograms. We explore the effects of the input Alfv\'en wave energy flux and the surface magnetogram on the wind properties predicted by the Alfv\'en Wave Solar Model (AWSoM). We lowered the resolution of two solar magnetograms during solar cycle maximum and minimum using spherical harmonic decomposition. The Alfv\'en wave energy was altered based on non-thermal velocities determined from a far ultraviolet (FUV) spectrum of the solar twin 18 Sco. Additionally, low-resolution magnetograms of three solar analogues were obtained using Zeeman Doppler imaging (ZDI). Finally, the simulated wind properties were compared to Advanced Composition Explorer (ACE) observations. AWSoM simulations using well constrained input parameters taken from solar observations can reproduce the observed solar wind mass and angular momentum loss rates. The resolution of the magnetogram has a small impact on the wind properties and only during cycle maximum. However, variation in Alfv\'en wave energy influences the wind properties irrespective of the solar cycle activity level. Furthermore, solar wind simulations carried out using the low-resolution magnetogram of the three stars instead of the solar magnetogram could lead to an order of a magnitude difference in the simulated wind properties. The choice in Alfv\'en energy has a stronger influence on the wind output compared to the magnetogram resolution. The influence could be even stronger for stars whose input boundary conditions are not as well constrained as those of the Sun. Unsurprisingly, replacing the solar magnetogram with a stellar magnetogram could lead to completely inaccurate solar wind properties, and should be avoided in solar and stellar wind simulations.

Solar coronal dimmings have been observed extensively in the past two decades and are believed to have close association with coronal mass ejections (CMEs). Recent study found that coronal dimming is the only signature that could differentiate powerful ares that have CMEs from those that do not. Therefore, dimming might be one of the best candidates to observe the stellar CMEs on distant Sun-like stars. In this study, we investigate the possibility of using coronal dimming as a proxy to diagnose stellar CMEs. By simulating a realistic solar CME event and corresponding coronal dimming using a global magnetohydrodynamics model (AWSoM: Alfven-wave Solar Model), we first demonstrate the capability of the model to reproduce solar observations. We then extend the model for simulating stellar CMEs by modifying the input magnetic flux density as well as the initial magnetic energy of the CME flux rope. Our result suggests that with improved instrument sensitivity, it is possible to detect the coronal dimming signals induced by the stellar CMEs.

A deep learning network, Long-Short Term Memory (LSTM) network, is used in this work to predict whether the maximum flare class an active region (AR) will produce in the next 24 hours is class $\Gamma$. We considered $\Gamma$ are $\ge M$, $\ge C$ and any flare class. The essence of using LSTM, which is a recurrent neural network, is its capability to capture temporal information of the data samples. The input features are time sequences of 20 magnetic parameters from SHARPs - Space-weather HMI Active Region Patches. We analyzed active regions from June 2010 to Dec 2018, using the Geostationary Operational Environmental Satellite (GOES) X-ray flare catalogs and label the data samples with identified ARs in the GOES X-ray flare catalogs. Our results (i) shows consistent skill scores with recently published results using LSTMs and better than the previous work using single time input (eg. DeFN) (ii) The skill scores from the model show essential differences when different years of data was chosen for training and testing.

As a part of NASA's Heliophysics System Observatory (HSO) fleet of satellites,the Solar Dynamics Observatory (SDO) has continuously monitored the Sun since2010. Ultraviolet (UV) and Extreme UV (EUV) instruments in orbit, such asSDO's Atmospheric Imaging Assembly (AIA) instrument, suffer time-dependent degradation which reduces instrument sensitivity. Accurate calibration for (E)UV instruments currently depends on periodic sounding rockets, which are infrequent and not practical for heliophysics missions in deep space. In the present work, we develop a Convolutional Neural Network (CNN) that auto-calibrates SDO/AIA channels and corrects sensitivity degradation by exploiting spatial patterns in multi-wavelength observations to arrive at a self-calibration of (E)UV imaging instruments. Our results remove a major impediment to developing future HSOmissions of the same scientific caliber as SDO but in deep space, able to observe the Sun from more vantage points than just SDO's current geosynchronous orbit.This approach can be adopted to perform autocalibration of other imaging systems exhibiting similar forms of degradation

Understanding and monitoring the complex and dynamic processes of the Sun is important for a number of human activities on Earth and in space. For this reason, NASA's Solar Dynamics Observatory (SDO) has been continuously monitoring the multi-layered Sun's atmosphere in high-resolution since its launch in 2010, generating terabytes of observational data every day. The synergy between machine learning and this enormous amount of data has the potential, still largely unexploited, to advance our understanding of the Sun and extend the capabilities of heliophysics missions. In the present work, we show that deep learning applied to SDO data can be successfully used to create a high-fidelity virtual telescope that generates synthetic observations of the solar corona by image translation. Towards this end we developed a deep neural network, structured as an encoder-decoder with skip connections (U-Net), that reconstructs the Sun's image of one instrument channel given temporally aligned images in three other channels. The approach we present has the potential to reduce the telemetry needs of SDO, enhance the capabilities of missions that have less observing channels, and transform the concept development of future missions.

The primary advantage of moderately superheated bubble chamber detectors is their simultaneous sensitivity to nuclear recoils from WIMP dark matter and insensitivity to electron recoil backgrounds. A comprehensive analysis of PICO gamma calibration data demonstrates for the first time that electron recoils in C$_3$F$_8$ scale in accordance with a new nucleation mechanism, rather than one driven by a hot-spike as previously supposed. Using this semi-empirical model, bubble chamber nucleation thresholds may be tuned to be sensitive to lower energy nuclear recoils while maintaining excellent electron recoil rejection. The PICO-40L detector will exploit this model to achieve thermodynamic thresholds as low as 2.8 keV while being dominated by single-scatter events from coherent elastic neutrino-nucleus scattering of solar neutrinos. In one year of operation, PICO-40L can improve existing leading limits from PICO on spin-dependent WIMP-proton coupling by nearly an order of magnitude for WIMP masses greater than 3 GeV c$^{-2}$ and will have the ability to surpass all existing non-xenon bounds on spin-independent WIMP-nucleon coupling for WIMP masses from 3 to 40 GeV c$^{-2}$.

The current progress in the detection of terrestrial type exoplanets has opened a new avenue in the characterization of exoplanetary atmospheres and in the search for biosignatures of life with the upcoming ground-based and space missions. To specify the conditions favorable for the origin, development and sustainment of life as we know it in other worlds, we need to understand the nature of astrospheric, atmospheric and surface environments of exoplanets in habitable zones around G-K-M dwarfs including our young Sun. Global environment is formed by propagated disturbances from the planet-hosting stars in the form of stellar flares, coronal mass ejections, energetic particles, and winds collectively known as astrospheric space weather. Its characterization will help in understanding how an exoplanetary ecosystem interacts with its host star, as well as in the specification of the physical, chemical and biochemical conditions that can create favorable and/or detrimental conditions for planetary climate and habitability along with evolution of planetary internal dynamics over geological timescales. A key linkage of (astro) physical, chemical, and geological processes can only be understood in the framework of interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary and Earth sciences. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets will significantly expand the current definition of the habitable zone to the biogenic zone and provide new observational strategies for searching for signatures of life. The major goal of this paper is to describe and discuss the current status and recent progress in this interdisciplinary field and to provide a new roadmap for the future development of the emerging field of exoplanetary science and astrobiology.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

The EUV (100-912 {\AA}) is a spectral region notoriously difficult to observe due to attenuation by neutral hydrogen gas in the interstellar medium. Despite this, hundreds to thousands of nearby stars of different spectral types and magnetic activity levels are accessible in the EUV range. The EUV probes interesting and complicated regions in the stellar atmosphere like the lower corona and transition region that are inaccessible from other spectral regions. In this white paper we describe how direct EUV observations, which require a dedicated grazing-incidence observatory, cannot yet be accurately substituted with models and theory. Exploring EUV emission from cool dwarf stars in the time domain can make a major contribution to understanding stellar outer atmospheres and magnetism, and offers the clearest path toward detecting coronal mass ejections on stars other than the Sun.

The planetary effective surface temperature alone is insufficient to characterize exoplanet atmospheres and their stability or evolution. Considering the star-planet system as a whole is necessary, and a critical component of the system is the photoionizing stellar extreme ultraviolet emission (EUV; 100-912 {\AA}). EUV photons drive atmospheric mass loss through thermal and nonthermal processes, and an accurate accounting of the EUV energy deposition in a planet's energy budget is essential, especially for terrestrial habitable zone planets and close-in gaseous planets. Direct EUV observations of exoplanet host stars would require a new, dedicated observatory. Archival observations from the $\textit{EUVE}$ satellite, models, and theory alone are insufficient to accurately characterize EUV spectra of the majority of exoplanet host stars, especially for low-mass stars.

In this paper we present a curated dataset from the NASA Solar Dynamics Observatory (SDO) mission in a format suitable for machine learning research. Beginning from level 1 scientific products we have processed various instrumental corrections, downsampled to manageable spatial and temporal resolutions, and synchronized observations spatially and temporally. We illustrate the use of this dataset with two example applications: forecasting future EVE irradiance from present EVE irradiance and translating HMI observations into AIA observations. For each application we provide metrics and baselines for future model comparison. We anticipate this curated dataset will facilitate machine learning research in heliophysics and the physical sciences generally, increasing the scientific return of the SDO mission. This work is a direct result of the 2018 NASA Frontier Development Laboratory Program. Please see the appendix for access to the dataset.

Final results are reported from operation of the PICO-60 C$_3$F$_8$ dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. The chamber was operated at thermodynamic thresholds as low as 1.2 keV without loss of stability. A new blind 1404-kg-day exposure at 2.45 keV threshold was acquired with approximately the same expected total background rate as the previous 1167-kg-day exposure at 3.3 keV. This increased exposure is enabled in part by a new optical tracking analysis to better identify events near detector walls, permitting a larger fiducial volume. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 2.5 $\times$ 10$^{-41}$ cm$^2$ for a 25 GeV WIMP, and improve on previous PICO results for 3-5 GeV WIMPs by an order of magnitude.

The LAMOST-Kepler (LK-) project was initiated to use the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) to make spectroscopic follow-up observations for the targets in the field of the Kepler mission. The Kepler field is divided into 14 subfields that are adapted to the LAMOST circular field with a diameter of 5 degrees. During the regular survey phase of LAMOST, the LK-project took data from 2012 June to 2017 June and covered all the 14 subfields at least twice. In particular, we describe in this paper the second Data Release of the LK-project, including all spectra acquired through 2015 May to 2017 June together with the first round observations of the LK-project from 2012 June to 2014 September. The LK-project now counts 227 870 spectra of 156 390 stars, among which we have derived atmospheric parameters (log g, T eff and [Fe/H]) and heliocentric radial velocity (RV) for 173 971 spectra of 126 172 stars. These parameters were obtained with the most recent version of the LAMOST Stellar Parameter Pipeline v 2.9.7. Nearly one half, namely 76 283 targets, are observed both by LAMOST and Kepler telescopes. These spectra, establishing a large spectroscopy library, will be useful for the entire astronomical community, particularly for planetary science and stellar variability on Kepler targets.

Recent detections of high-energy $\gamma$-rays from behind-the-limb (BTL) solar flares by the \emph{Fermi $\gamma$-ray Space Telescope} pose a puzzle and challenge on the particle acceleration and transport mechanisms. In such events, the $\gamma$-ray emission region is located away from the BTL flare site by up to tens of degrees in heliogrpahic longitude. It is thus hypothesized that particles are accelerated at the shock driven by the coronal mass ejection (CME) and then travel from the shock downstream back to the front side of the Sun to produce the observed $\gamma$-rays. To test this scenario, we performed data-driven, global magnetohydrodynamics simulations of the CME associated with a well-observed BTL flare on 2014 September 1. We found that part of the CME-driven shock develops magnetic connectivity with the $\gamma$-ray emission region, facilitating transport of particles back to the Sun. Moreover, the observed increase in $\gamma$-ray flux is temporally correlated with (1) the increase of the shock compression ratio and (2) the presence of a quasi-perpendicular shock over the area that is magnetically connected to the $\gamma$-ray emitting region, both conditions favoring the diffusive shock acceleration (DSA) of particles. These results support the above hypothesis and can help resolve another puzzle, i.e., long-duration (up to 20 hours) $\gamma$-rays flares. We suggest that, in addition to DSA, stochastic acceleration by plasma turbulence may also play a role, especially in the shock downstream region and during the early stage when the shock Alfv\'{e}n Mach number is small.

We report SDO/AIA observations of an extraordinary global extreme ultraviolet (EUV) wave triggered by the X8.2+ flare-CME eruption on 2017 September 10. This was one of the best EUV waves ever observed with modern instruments, yet it was likely the last one of such magnitudes of Solar Cycle 24 as the Sun heads toward the minimum. Its remarkable characteristics include the following. (1) The wave was observed, for the first time, to traverse the full-Sun corona over the entire visible solar disk and off-limb circumference, manifesting a truly global nature, owing to its exceptionally large amplitude, e.g., with EUV enhancements by up to 300% at 1.1 Rsun from the eruption. (2) This leads to strong transmissions (in addition to commonly observed reflections) in and out of both polar coronal holes, which are usually devoid of EUV waves. It has elevated wave speeds >2000 km/s within them, consistent with the expected higher fast-mode magnetosonic wave speeds. The coronal holes essentially serve as new "radiation centers" for the waves being refracted out of them, which then travel toward the equator and collide head-on, causing additional EUV enhancements. (3) The wave produces significant compressional heating to local plasma upon its impact, indicated by long-lasting EUV intensity changes and differential emission measure increases at higher temperatures (e.g., log T=6.2) accompanied by decreases at lower temperatures (e.g., log T=6.0). These characteristics signify the potential of such EUV waves for novel magnetic and thermal diagnostics of the solar corona {\it on global scales}.

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.

The search of life in the Universe is a fundamental problem of astrobiology and a major priority for NASA. A key area of major progress since the NASA Astrobiology Strategy 2015 (NAS15) has been a shift from the exoplanet discovery phase to a phase of characterization and modeling of the physics and chemistry of exoplanetary atmospheres, and the development of observational strategies for the search for life in the Universe by combining expertise from four NASA science disciplines including heliophysics, astrophysics, planetary science and Earth science. The NASA Nexus for Exoplanetary System Science (NExSS) has provided an efficient environment for such interdisciplinary studies. Solar flares, coronal mass ejections and solar energetic particles produce disturbances in interplanetary space collectively referred to as space weather, which interacts with the Earth upper atmosphere and causes dramatic impact on space and ground-based technological systems. Exoplanets within close in habitable zones around M dwarfs and other active stars are exposed to extreme ionizing radiation fluxes, thus making exoplanetary space weather (ESW) effects a crucial factor of habitability. In this paper, we describe the recent developments and provide recommendations in this interdisciplinary effort with the focus on the impacts of ESW on habitability, and the prospects for future progress in searching for signs of life in the Universe as the outcome of the NExSS workshop held in Nov 29 - Dec 2, 2016, New Orleans, LA. This is one of five Life Beyond the Solar System white papers submitted by NExSS to the National Academy of Sciences in support of the Astrobiology Science Strategy for the Search for Life in the Universe.

The presence of an atmosphere over sufficiently long timescales is widely perceived as one of the most prominent criteria associated with planetary surface habitability. We address the crucial question as to whether the seven Earth-sized planets transiting the recently discovered ultracool dwarf star TRAPPIST-1 are capable of retaining their atmospheres. To this effect, we carry out numerical simulations to characterize the stellar wind of TRAPPIST-1 and the atmospheric ion escape rates for all the seven planets. We also estimate the escape rates analytically and demonstrate that they are in good agreement with the numerical results. We conclude that the outer planets of the TRAPPIST-1 system are capable of retaining their atmospheres over billion-year timescales. The consequences arising from our results are also explored in the context of abiogenesis, biodiversity, and searches for future exoplanets. In light of the many unknowns and assumptions involved, we recommend that these conclusions must be interpreted with due caution.

New results are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. As in previous PICO bubble chambers, PICO-60 C$_3$F$_8$ exhibits excellent electron recoil and alpha decay rejection, and the observed multiple-scattering neutron rate indicates a single-scatter neutron background of less than 1 event per month. A blind analysis of an efficiency-corrected 1167-kg-day exposure at a 3.3-keV thermodynamic threshold reveals no single-scattering nuclear recoil candidates, consistent with the predicted background. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 3.4 $\times$ 10$^{-41}$ cm$^2$ for a 30-GeV$\thinspace$c$^{-2}$ WIMP, more than one order of magnitude improvement from previous PICO results.

A 30-g xenon bubble chamber, operated at Northwestern University in June and November 2016, has for the first time observed simultaneous bubble nucleation and scintillation by nuclear recoils in a superheated liquid. This chamber is instrumented with a CCD camera for near-IR bubble imaging, a solar-blind photomultiplier tube to detect 175-nm xenon scintillation light, and a piezoelectric acoustic transducer to detect the ultrasonic emission from a growing bubble. The time of nucleation determined from the acoustic signal is used to correlate specific scintillation pulses with bubble-nucleating events. We report on data from this chamber for thermodynamic "Seitz" thresholds from 4.2 to 15.0 keV. The observed single- and multiple-bubble rates when exposed to a $^{252}$Cf neutron source indicate that, for an 8.3-keV thermodynamic threshold, the minimum nuclear recoil energy required to nucleate a bubble is $19\pm6$ keV (1$\sigma$ uncertainty). This is consistent with the observed scintillation spectrum for bubble-nucleating events. We see no evidence for bubble nucleation by gamma rays at any of the thresholds studied, setting a 90% C.L. upper limit of $6.3\times10^{-7}$ bubbles per gamma interaction at a 4.2-keV thermodynamic threshold. This indicates stronger gamma discrimination than in CF$_3$I bubble chambers, supporting the hypothesis that scintillation production suppresses bubble nucleation by electron recoils while nuclear recoils nucleate bubbles as usual. These measurements establish the noble-liquid bubble chamber as a promising new technology for the detection of weakly interacting massive particle dark matter and coherent elastic neutrino-nucleus scattering.

We present an updated global model of the solar corona, including the transition region. We simulate the realistic tree-dimensional (3D) magnetic field using the data from the photospheric magnetic field measurements and assume the magnetohydrodynamic (MHD) Alfv\'en wave turbulence and its non-linear dissipation to be the only source for heating the coronal plasma and driving the solar wind. In closed field regions the dissipation efficiency in a balanced turbulence is enhanced. In the coronal holes we account for a reflection of the outward propagating waves, which is accompanied by generation of weaker counter-propagating waves. The non-linear cascade rate degrades in strongly imbalanced turbulence, thus resulting in colder coronal holes. The distinctive feature of the presented model is the description of the low corona as almost-steady-state low-beta plasma motion and heat flux transfer along the magnetic field lines. We trace the magnetic field lines through each grid point of the lower boundary of the global corona model, chosen at some heliocentric distance, $R=R_{b}\sim1.1\ R_\odot$ well above the transition region. One can readily solve the plasma parameters along the magnetic field line from 1D equations for the plasma motion and heat transport together with the Alfv\'en wave propagation, which adequately describe physics within the heliocentric distances range, $R_{\odot}<R<R_{b}$, in the low solar corona. By interfacing this threaded-field-lines model with the full MHD global corona model at $r=R_{b}$, we find the global solution and achieve a faster-than-real-time performance of the model on $\sim200$ cores.

We perform and analyze results of a global magnetohydrodyanmic (MHD) simulation of the fast coronal mass ejection (CME) that occurred on 2011 March 7. The simulation is made using the newly developed Alfv\'en Wave Solar Model (AWSoM), which describes the background solar wind starting from the upper chromosphere and extends to 24 R$_{\odot}$. Coupling AWSoM to an inner heliosphere (IH) model with the Space Weather Modeling Framework (SWMF) extends the total domain beyond the orbit of Earth. Physical processes included in the model are multi-species thermodynamics, electron heat conduction (both collisional and collisionless formulations), optically thin radiative cooling, and Alfv\'en-wave turbulence that accelerates and heats the solar wind. The Alfv\'en-wave description is physically self-consistent, including non-Wentzel-Kramers-Brillouin (WKB) reflection and physics-based apportioning of turbulent dissipative heating to both electrons and protons. Within this model, we initiate the CME by using the Gibson-Low (GL) analytical flux rope model and follow its evolution for days, in which time it propagates beyond STEREO A. A detailed comparison study is performed using remote as well as \textit{in situ} observations. Although the flux rope structure is not compared directly due to lack of relevant ejecta observation at 1 AU in this event, our results show that the new model can reproduce many of the observed features near the Sun (e.g., CME-driven extreme ultraviolet (EUV) waves, deflection of the flux rope from the coronal hole, "double-front" in the white light images) and in the heliosphere (e.g., shock propagation direction, shock properties at STEREO A).

We present a first-principles-based coronal mass ejection (CME) model suitable for both scientific and operational purposes by combining a global magnetohydrodynamics (MHD) solar wind model with a flux rope-driven CME model. Realistic CME events are simulated self-consistently with high fidelity and forecasting capability by constraining initial flux rope parameters with observational data from GONG, SOHO/LASCO, and STEREO/COR. We automate this process so that minimum manual intervention is required in specifying the CME initial state. With the newly developed data-driven Eruptive Event Generator Gibson-Low (EEGGL), we present a method to derive Gibson-Low (GL) flux rope parameters through a handful of observational quantities so that the modeled CMEs can propagate with the desired CME speeds near the Sun. A test result with CMEs launched with different Carrington rotation magnetograms are shown. Our study shows a promising result for using the first-principles-based MHD global model as a forecasting tool, which is capable of predicting the CME direction of propagation, arrival time, and ICME magnetic field at 1 AU (see companion paper by Jin et al. 2016b).

We analyzed both HCN J=1-0 and HNC J=1-0 line profiles to study the inflow motions in different evolutionary stages of massive star formation: 54 infrared dark clouds (IRDCs), 69 high-mass protostellar object (HMPOs), and 54 ultra-compact HII regions (UCHIIs). The inflow asymmetry in HCN spectra seems to be prevalent throughout all the three evolutionary phases, with IRDCs showing the largest excess in blue profile. In the case of HNC spectra, the prevalence of blue sources does not appear, excepting for IRDCs. We suggest that this line is not appropriate to trace inflow motion in evolved stages of massive star formation because the abundance of HNC decreases at high temperatures. This result spotlights the importance of considering chemistry in the dynamics study of massive star-forming regions. The fact that the IRDCs show the highest blue excess in both transitions indicates that the most active inflow occurs in the early phase of star formation, i.e., the IRDC phase rather than in the later phases. However, mass is still inflowing onto some UCHIIs. We also found that the absorption dips of the HNC spectra in 6 out of 7 blue sources are red-shifted relative to their systemic velocities. These red-shifted absorption dips may indicate global collapse candidates, although mapping observations with better resolution are needed to examine this feature in more detail.

New data are reported from a second run of the 2-liter PICO-2L C$_3$F$_8$ bubble chamber with a total exposure of 129$\,$kg-days at a thermodynamic threshold energy of 3.3$\,$keV. These data show that measures taken to control particulate contamination in the superheated fluid resulted in the absence of the anomalous background events observed in the first run of this bubble chamber. One single nuclear-recoil event was observed in the data, consistent both with the predicted background rate from neutrons and with the observed rate of unambiguous multiple-bubble neutron scattering events. The chamber exhibits the same excellent electron-recoil and alpha decay rejection as was previously reported. These data provide the most stringent direct detection constraints on weakly interacting massive particle (WIMP)-proton spin-dependent scattering to date for WIMP masses $<$ 50$\,$GeV/c$^2$.

With the global view and high-cadence observations from SDO/AIA and STEREO, many spatially separated solar eruptive events appear to be coupled. However, the mechanisms for "sympathetic" events are still largely unknown. In this study, we investigate the impact of an erupting flux rope on surrounding solar structures through large-scale magnetic coupling. We build a realistic environment of the solar corona on 2011 February 15 using a global magnetohydrodynamics (MHD) model and initiate coronal mass ejections (CMEs) in active region (AR) 11158 by inserting Gibson-Low analytical flux ropes. We show that a CME's impact on the surrounding structures depends not only on the magnetic strength of these structures and their distance to the source region, but also on the interaction between the CME with the large-scale magnetic field. Within the CME expansion domain where the flux rope field directly interacts with the solar structures, expansion-induced reconnection often modifies the overlying field, thereby increasing the decay index. This effect may provide a primary coupling mechanism underlying the sympathetic eruptions. The magnitude of the impact is found to depend on the orientation of the erupting flux rope, with the largest impacts occurring when the flux rope is favorably oriented for reconnecting with the surrounding regions. Outside the CME expansion domain, the influence of the CME is mainly through field line compression or post-eruption relaxation. Based on our numerical experiments, we discuss a way to quantify the eruption impact, which could be useful for forecasting purposes.

New data are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 36.8 kg of CF$_3$I and located in the SNOLAB underground laboratory. PICO-60 is the largest bubble chamber to search for dark matter to date. With an analyzed exposure of 92.8 livedays, PICO-60 exhibits the same excellent background rejection observed in smaller bubble chambers. Alpha decays in PICO-60 exhibit frequency-dependent acoustic calorimetry, similar but not identical to that reported recently in a C$_3$F$_8$ bubble chamber. PICO-60 also observes a large population of unknown background events, exhibiting acoustic, spatial, and timing behaviors inconsistent with those expected from a dark matter signal. These behaviors allow for analysis cuts to remove all background events while retaining $48.2\%$ of the exposure. Stringent limits on weakly interacting massive particles interacting via spin-dependent proton and spin-independent processes are set, and most interpretations of the DAMA/LIBRA modulation signal as dark matter interacting with iodine nuclei are ruled out.

New data are reported from the operation of a 2-liter C$_3$F$_8$ bubble chamber in the 2100 meter deep SNOLAB underground laboratory, with a total exposure of 211.5 kg-days at four different recoil energy thresholds ranging from 3.2 keV to 8.1 keV. These data show that C3F8 provides excellent electron recoil and alpha rejection capabilities at very low thresholds, including the first observation of a dependence of acoustic signal on alpha energy. Twelve single nuclear recoil event candidates were observed during the run. The candidate events exhibit timing characteristics that are not consistent with the hypothesis of a uniform time distribution, and no evidence for a dark matter signal is claimed. These data provide the most sensitive direct detection constraints on WIMP-proton spin-dependent scattering to date, with significant sensitivity at low WIMP masses for spin-independent WIMP-nucleon scattering.

Using the HCN and HNC J=1--0 line observations, the abundance ratio of HCN/HNC has been estimated for different evolutionary stages of massive star formation: Infrared dark clouds (IRDCs), High-mass protostellar object (HMPOs), and Ultra-compact HII regions (UCHIIs). IRDCs were divided into `quiescent IRDC cores' and `active IRDC cores', depending on star formation activity. The HCN/HNC ratio is known to be higher at active and high temperature regions related to ongoing star formation, compared to cold and quiescent regions. Our observations toward 8 quiescent IRDC cores, 16 active IRDC cores, 23 HMPOs, and 31 UCHIIs show consistent results; the ratio is 0.97~($\pm~0.10$), 2.65~($\pm~0.88$), 4.17~($\pm~1.03$) and 8.96~($\pm~3.32$) in these respective evolutionary stages, increasing from quiescent IRDC cores to UCHIIs. The change of the HCN/HNC abundance ratio, therefore, seems directly associated with the evolutionary stages of star formation, which have different temperatures. One suggested explanation for this trend is the conversion of HNC to HCN, which occurs effectively at higher temperatures. To test the explanation, we performed a simple chemical model calculation. In order to fit the observed results, the energy barrier of the conversion must be much lower than the value provided by theoretical calculations.

We present a new version of the Alfven Wave Solar Model (AWSoM), a global model from the upper chromosphere to the corona and the heliosphere. The coronal heating and solar wind acceleration are addressed with low-frequency Alfven wave turbulence. The injection of Alfven wave energy at the inner boundary is such that the Poynting flux is proportional to the magnetic field strength. The three-dimensional magnetic field topology is simulated using data from photospheric magnetic field measurements. This model does not impose open-closed magnetic field boundaries; those develop self-consistently. The physics includes: (1) The model employs three different temperatures, namely the isotropic electron temperature and the parallel and perpendicular ion temperatures. The firehose, mirror, and ion-cyclotron instabilities due to the developing ion temperature anisotropy are accounted for. (2) The Alfven waves are partially reflected by the Alfven speed gradient and the vorticity along the field lines. The resulting counter-propagating waves are responsible for the nonlinear turbulent cascade. The balanced turbulence due to uncorrelated waves near the apex of the closed field lines and the resulting elevated temperatures are addressed. (3) To apportion the wave dissipation to the three temperatures, we employ the results of the theories of linear wave damping and nonlinear stochastic heating. (4) We have incorporated the collisional and collisionless electron heat conduction. We compare the simulated multi-wavelength EUV images of CR2107 with the observations from STEREO/EUVI and SDO/AIA instruments. We demonstrate that the reflection due to strong magnetic fields in proximity of active regions intensifies the dissipation and observable emission sufficiently.

We describe, analyze and validate the recently developed Alfv\'en Wave Solar Model (AWSoM), a 3D global model starting from the top of the chromosphere and extending into interplanetary space (up to 1-2 AU). This model solves the extended two temperature magnetohydrodynamics equations coupled to a wave kinetic equation for low frequency Alfv\'en waves. In this picture, heating and acceleration of the plasma are due to wave dissipation and wave pressure gradients, respectively. The dissipation process is described by a fully developed turbulent cascade of counter-propagating waves. We adopt a unified approach for calculating the wave dissipation in both open and closed magnetic field lines, allowing for a self-consistent treatment of any magnetic topology. Wave dissipation is the only heating mechanism assumed in the model, and no geometric heating functions are invoked. Electron heat conduction and radiative cooling are also included. We demonstrate that the large-scale, steady-state (in the co-rotating frame) properties of the solar environment are reproduced, using three adjustable parameters: the Poynting flux of chromospheric Alfv\'en waves, the perpendicular correlation length of the turbulence, and a pseudo-reflection coefficient. We compare model results for Carrington Rotation 2063 (November-December 2007) to remote observations in the EUV and X-ray ranges from STEREO, SOHO and Hinode spacecraft, as well as to in-situ measurements performed by Ulysses. The model results are in good agreement with observations. This is the first global model capable of simultaneously reproducing the multi-wavelength observations of the lower corona and the wind structure beyond Earth's orbit.

We present a new global model of the solar corona, including the low corona, the transition region and the top of chromosphere. The realistic 3D magnetic field is simulated using the data from the photospheric magnetic field measurements. The distinctive feature of the new model is incorporating the MHD Alfven wave turbulence. We assume this turbulence and its non-linear dissipation to be the only momentum and energy source for heating the coronal plasma and driving the solar wind. The difference between the turbulence dissipation efficiency in coronal holes and that in closed field regions is because the non-linear cascade rate degrades in strongly anisotropic turbulence in coronal holes (no inward propagating wave), thus resulting in colder coronal holes with the bi-modal solar wind originating from them. The detailed presentation of the theoretical model is illustrated with the synthetic images for multi-wavelength EUV emission compared with the observations from SDO AIA and Stereo EUVI instruments for the Carrington rotation 2107.

We study the correlation heights, which indicate the formation height of Extreme Ultraviolet (EUV) lines in an active region using observations from the EUV Imaging Spectrometer (EIS) and Solar Optical Telescope (SOT) on board \emph{Hinode}. The nonlinear force-free field (NLFFF) optimization method is adopted to extrapolate the 3D magnetic fields to higher layers. Three subregions with different characteristics are selected in the active region for this study. The results show that the formation heights in different subregions vary with their different magnetic fields or velocity patterns. After solving the line blending problem between the He {\sc \romannumeral 2} 256.32 \AA and Si {\sc \romannumeral 10} 256.37 \AA lines by the double Gaussian curve fitting, we find that the transition region lies higher in a strong magnetic area. In a pre-flare heating area there possibly exist multithermal loops as implied by comparing the Doppler velocity and the magnetic field on the solar disk.

A general analysis of Meissner effect and spin susceptibility of a uniform superconductor in an asymmetric two-component fermion system is presented in nonrelativistic field theory approach. We found that, the pairing mechanism dominates the magnetization property of superconductivity, and the asymmetry enhances the paramagnetism of the system. At the turning point from BCS to breached pairing superconductivity, the Meissner mass squared and spin susceptibility are divergent at zero temperature. In the breached pairing state induced by chemical potential difference and mass difference between the two kinds of fermions, the system goes from paramagnetism to diamagnetism, when the mass ratio of the two species increases.

The electromagnetic(EM) interactions between charged protons on the correlations of nucleons are discussed by introducing the Anderson-Higgs mechanism of broken U(1) EM symmetry into the relativistic nuclear theory with a parametric photon mass. The non-saturating Coulomb force contribution is emphasized on the equation of state of nuclear matter with charge symmetry breaking(CSB) at finite temperature and the breached $^1S_0$ pairing correlations of proton-proton and neutron-neutron. The universal properties given by an order parameter field with a non-zero vacuum expectation value (VEV) nearby phase transition are explored within the mean field theory(MFT) level. This mechanism can be extended to the charged or charge neutralized strongly coupling multi-components system for the discussion of binding or pairing issues.