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

Studies of exoplanet demographics require large samples and precise constraints on exoplanet host stars. Using the homogeneous Kepler stellar properties derived using Gaia Data Release 2 by Berger et al. (2020), we re-compute Kepler planet radii and incident fluxes and investigate their distributions with stellar mass and age. We measure the stellar mass dependence of the planet radius valley to be $d \log R_{\mathrm{p}}$/$d \log M_\star = 0.26^{+0.21}_{-0.16}$, consistent with the slope predicted by a planet mass dependence on stellar mass ($0.24-0.35$) and core-powered mass-loss (0.33). We also find first evidence of a stellar age dependence of the planet populations straddling the radius valley. Specifically, we determine that the fraction of super-Earths ($1-1.8 \mathrm{R_\oplus}$) to sub-Neptunes ($1.8-3.5 \mathrm{R_\oplus}$) increases from $0.61 \pm 0.09$ at young ages (< 1 Gyr) to $1.00 \pm 0.10$ at old ages (> 1 Gyr), consistent with the prediction by core-powered mass-loss that the mechanism shaping the radius valley operates over Gyr timescales. Additionally, we find a tentative decrease in the radii of relatively cool ($F_{\mathrm{p}} < 150 \mathrm{F_\oplus}$) sub-Neptunes over Gyr timescales, which suggests that these planets may possess H/He envelopes instead of higher mean molecular weight atmospheres. We confirm the existence of planets within the hot sub-Neptunian "desert" ($2.2 < R_{\mathrm{p}} < 3.8 \mathrm{R_\oplus}$, $F_{\mathrm{p}} > 650 \mathrm{F_\oplus}$) and show that these planets are preferentially orbiting more evolved stars compared to other planets at similar incident fluxes. In addition, we identify candidates for cool ($F_{\mathrm{p}} < 20 \mathrm{F_\oplus}$) inflated Jupiters, present a revised list of habitable zone candidates, and find that the ages of single- and multiple-transiting planet systems are statistically indistinguishable.

An accurate and precise Kepler Stellar Properties Catalog is essential for the interpretation of the Kepler exoplanet survey results. Previous Kepler Stellar Properties Catalogs have focused on reporting the best-available parameters for each star, but this has required combining data from a variety of heterogeneous sources. We present the Gaia-Kepler Stellar Properties Catalog, a set of stellar properties of 186,301 Kepler stars, homogeneously derived from isochrones and broadband photometry, Gaia Data Release 2 parallaxes, and spectroscopic metallicities, where available. Our photometric effective temperatures, derived from $g-K_s$ colors, are calibrated on stars with interferometric angular diameters. Median catalog uncertainties are 112 K for $T_{\mathrm{eff}}$, 0.05 dex for $\log g$, 4% for $R_\star$, 7% for $M_\star$, 13% for $\rho_\star$, 10% for $L_\star$, and 56% for stellar age. These precise constraints on stellar properties for this sample of stars will allow unprecedented investigations into trends in stellar and exoplanet properties as a function of stellar mass and age. In addition, our homogeneous parameter determinations will permit more accurate calculations of planet occurrence and trends with stellar properties.

The rotational evolution of cool dwarfs is poorly constrained after around 1-2 Gyr due to a lack of precise ages and rotation periods for old main-sequence stars. In this work we use velocity dispersion as an age proxy to reveal the temperature-dependent rotational evolution of low-mass Kepler dwarfs, and demonstrate that kinematic ages could be a useful tool for calibrating gyrochronology in the future. We find that a linear gyrochronology model, calibrated to fit the period-Teff relationship of the Praesepe cluster, does not apply to stars older than around 1 Gyr. Although late-K dwarfs spin more slowly than early-K dwarfs when they are young, at old ages we find that late-K dwarfs rotate at the same rate or faster than early-K dwarfs of the same age. This result agrees qualitatively with semi-empirical models that vary the rate of surface-to-core angular momentum transport as a function of time and mass. It also aligns with recent observations of stars in the NGC 6811 cluster, which indicate that the surface rotation rates of K dwarfs go through an epoch of inhibited evolution. We find that the oldest Kepler stars with measured rotation periods are late-K and early-M dwarfs, indicating that these stars maintain spotted surfaces and stay magnetically active longer than more massive stars. Finally, based on their kinematics, we confirm that many rapidly rotating GKM dwarfs are likely to be synchronized binaries.

We obtained high-resolution infrared spectroscopy and short-cadence photometry of the 600-800 Myr Praesepe star K2-100 during transits of its 1.67-day planet. This Neptune-size object, discovered by the NASA K2 mission, is an interloper in the "desert" of planets with similar radii on short period orbits. Our observations can be used to understand its origin and evolution by constraining the orbital eccentricity by transit fitting, measuring the spin-orbit obliquity by the Rossiter-McLaughlin effect, and detecting any extended, escaping hydrogen-helium envelope with the 10830A line of neutral helium in the 2s3S triplet state. Transit photometry with 1-min cadence was obtained by the K2 satellite during Campaign 18 and transit spectra were obtained with the IRD spectrograph on the Subaru telescope. While the elevated activity of K2-100 prevented us from detecting the Rossiter-McLaughlin effect, the new photometry combined with revised stellar parameters allowed us to constrain the eccentricity to e < 0.15/0.28 with 90%/99% confidence. We modeled atmospheric escape as an isothermal, spherically symmetric Parker wind, with photochemistry driven by UV radiation that we estimate by combining the observed spectrum of the active Sun with calibrations from observations of K2-100 and similar young stars in the nearby Hyades cluster. Our non-detection (<5.7mA) of a transit-associated He I line limits mass loss of a solar-composition atmosphere through a T<10000K wind to <0.3Me/Gyr. Either K2-100b is an exceptional desert-dwelling planet, or its mass loss is occurring at a lower rate over a longer interval, consistent with a core accretion-powered scenario for escape.

Current operational forecasts of solar eruptions are made by human experts using a combination of qualitative shape-based classification systems and historical data about flaring frequencies. In the past decade, there has been a great deal of interest in crafting machine-learning (ML) flare-prediction methods to extract underlying patterns from a training set---e.g., a set of solar magnetogram images, each characterized by features derived from the magnetic field and labeled as to whether it was an eruption precursor. These patterns, captured by various methods (neural nets, support vector machines, etc.), can then be used to classify new images. A major challenge with any ML method is the \textit{featurization} of the data: pre-processing the raw images to extract higher-level properties, such as characteristics of the magnetic field, that can streamline the training and use of these methods. It is key to choose features that are informative, from the standpoint of the task at hand. To date, the majority of ML-based solar eruption methods have used physics-based magnetic and electric field features such as the total unsigned magnetic flux, the gradients of the fields, the vertical current density, etc. In this paper, we extend the relevant feature set to include characteristics of the magnetic field that are based purely on the geometry and topology of 2D magnetogram images and show that this improves the prediction accuracy of a neural-net based flare-prediction method.

Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{\beta\beta}\simeq$ 2.46 MeV. For the XENON1T dual-phase time projection chamber, we demonstrate that the relative energy resolution at 1 $\sigma/\mu$ is as low as (0.79$\pm$0.02) % in its one-ton fiducial mass, and for single-site interactions at $Q_{\beta\beta}$. We achieve this by using a signal correction method to rectify the saturation effects of the signal readout system. The very good energy resolution from keV to MeV energies demonstrated in XENON1T opens up new windows for the xenon dual-phase dark matter detectors to simultaneously search for other rare events.

This is the second of three papers describing an `absolute' calibration of the GONG magnetograph using an end-to-end simulation of its measurement process. In the first paper, we described the GONG instrument and our `end-to-end' simulation of its measurement process. In this paper, we consider the theory of calibration, and magnetograph comparison in general, identifying some of the significant issues and pitfalls. The calibration of a magnetograph is a function of whether or not it preserves flux, independent of its spatial resolution. However, we find that the one-dimensional comparison methods most often used for magnetograph calibration and comparison will show dramatic differences between two magnetograms with differing spatial resolution, even if they both preserve flux. Some of the apparent disagreement between magnetograms found in the literature are likely a result of these instrumental resolution differences rather than any intrinsic calibration differences. To avoid them, spatial resolution must be carefully matched prior to comparing magnetograms or making calibration curves. In the third paper, we apply the lessons learned here to absolute calibration of GONG using our `end-to-end' measurement simulation.

Exoplanet catalogs produced by surveys suffer from a lack of completeness (not every planet is detected) and less than perfect reliability (not every planet in the catalog is a true planet), particularly near the survey's detection limit. Exoplanet occurrence rate studies based on such a catalog must be corrected for completeness and reliability. The final Kepler data release, DR25, features a uniformly vetted planet candidate catalog and data products that facilitate corrections. We present a new probabilistic approach to the characterization of Kepler completeness and reliability, making full use of the Kepler DR25 products. We illustrate the impact of completeness and reliability corrections with a Poisson-likelihood occurrence rate method, using a recent stellar properties catalog that incorporates Gaia stellar radii and essentially uniform treatment of the stellar population. Correcting for reliability has a significant impact: the exoplanet occurrence rate for orbital period and radius within 20% of Earth's around GK dwarf stars, corrected for reliability, is 0.015+0.011-0.007, whereas not correcting results in 0.034+0.018-0.012 - correcting for reliability reduces this occurrence rate by more than a factor of two. We further show that using Gaia-based vs. DR25 stellar properties impacts the same occurrence rate by a factor of two. We critically examine the the DR25 catalog and the assumptions behind our occurrence rate method. We propose several ways in which confidence in both the Kepler catalog and occurrence rate calculations can be improved. This work provides an example of how the community can use the DR25 completeness and reliability products.

This is the first of three papers describing an `absolute' calibration of the GONG magnetograph using an end-to-end simulation of its measurement process. The input to this simulation is a MURaM 3D MHD photospheric simulation and the output is the corresponding set of simulated data numbers which would be recorded by the GONG detectors. These simulated data numbers are then used to produce `synthetic magnetograms' which can be compared with the simulation inputs. This paper describes the GONG instrument, the MURaM datacube, our instrument simulator, and calculation of synthetic magnetograms, setting the stage for the subsequent two papers. These will first lay groundwork for calibration (and magnetogram comparison in general), then apply them to calibration of GONG using the simulation results.

This is the last of three papers describing an `absolute' calibration of the GONG magnetograph using and end-to-end simulation of its measurement process. The simulation begins with a MURaM 3D MHD datacube and ends with a `synthetic magnetogram' of the corresponding magnetic field values as they would be observed by GONG. We determine a calibration by comparing the synthetic magnetic field measurements with the MURaM magnetic field values that produced them. The previous two papers have described the GONG measurement process (both instrument and data processing), our simulation of it, and the theory of magnetogram comparison and calibration. In this paper, we address some final points on calibration, combine all of this work into a set of calibration curves, and consider the results. We also review the results of the previous two papers for locality of reference. Our calibration indicates that GONG magnetograms underestimate weak flux by a factor of over 2 near disk center, but that factor decreases to about 1 as the line-of-sight approaches the limb. A preliminary investigation of the generalizability of these results suggests other instruments will be affected in a similar way. We also find that some differences in previous magnetograph comparisons are artifacts of instrumental resolution which do not reflect an intrinsic calibration difference, and the measurements are more similar than sometimes thought. These results are directly applicable to question of solar wind prediction model accuracies, particularly in the search for the cause of the common discrepancy between predicted solar wind magnetic flux at 1 AU and values measured in situ by current satellite missions.

Chang'E 4 is the first mission to the far side of the Moon and consists of a lander, a rover, and a relay spacecraft. Lander and rover were launched at 18:23 UTC on December 7, 2018 and landed in the von K\'arm\'an crater at 02:26 UTC on January 3, 2019. Here we describe the Lunar Lander Neutron \& Dosimetry experiment (LND) which is part of the Chang'E 4 Lander scientific payload. Its chief scientific goal is to obtain first active dosimetric measurements on the surface of the Moon. LND also provides observations of fast neutrons which are a result of the interaction of high-energy particle radiation with the lunar regolith and of their thermalized counterpart, thermal neutrons, which are a sensitive indicator of subsurface water content.

Direct dark matter detection experiments based on a liquid xenon target are leading the search for dark matter particles with masses above $\sim$ 5 GeV/c$^2$, but have limited sensitivity to lighter masses because of the small momentum transfer in dark matter-nucleus elastic scattering. However, there is an irreducible contribution from inelastic processes accompanying the elastic scattering, which leads to the excitation and ionization of the recoiling atom (the Migdal effect) or the emission of a Bremsstrahlung photon. In this letter, we report on a probe of low-mass dark matter with masses down to about 85 MeV/c$^2$ by looking for electronic recoils induced by the Migdal effect and Bremsstrahlung, using data from the XENON1T experiment. Besides the approach of detecting both scintillation and ionization signals, we exploit an approach that uses ionization signals only, which allows for a lower detection threshold. This analysis significantly enhances the sensitivity of XENON1T to light dark matter previously beyond its reach.

We report constraints on light dark matter (DM) models using ionization signals in the XENON1T experiment. We mitigate backgrounds with strong event selections, rather than requiring a scintillation signal, leaving an effective exposure of $(22 \pm 3)$ tonne-days. Above $\sim\!0.4\,\mathrm{keV}_\mathrm{ee}$, we observe $<1 \, \text{event}/(\text{tonne} \times \text{day} \times \text{keV}_\text{ee})$, which is more than one thousand times lower than in similar searches with other detectors. Despite observing a higher rate at lower energies, no DM or CEvNS detection may be claimed because we cannot model all of our backgrounds. We thus exclude new regions in the parameter spaces for DM-nucleus scattering for DM masses $m_\chi$ within $3-6\,\mathrm{GeV}/\mathrm{c}^2$, DM-electron scattering for $m_\chi > 30\,\mathrm{MeV}/\mathrm{c}^2$, and absorption of dark photons and axion-like particles for $m_\chi$ within $0.186 - 1 \, \mathrm{keV}/\mathrm{c}^2$.

We analyze the full $\mathit{Kepler}$ phase curve of KOI-964, a binary system consisting of a hot white dwarf on an eclipsing orbit around an A-type host star. Using all 18 quarters of long-cadence photometry, we carry out a joint light curve fit and obtain improved phase curve amplitudes, occultation depths, orbital parameters, and transit ephemeris over the previous results of Carter et al. 2011. A periodogram of the residuals from the phase curve fit reveals an additional stellar variability signal from the host star with a characteristic period of $0.620276\pm0.000011$ days and a full amplitude of $24\pm2$ ppm. We also present new Keck/HIRES radial velocity observations which we use to measure the orbit and obtain a mass ratio of $q=0.106\pm0.012$. Combining this measurement with the results of a stellar isochrone analysis, we find that the masses of the host star and white dwarf companion are $2.23\pm0.12\,M_{\odot}$ and $0.236^{+0.028}_{-0.027}\,M_{\odot}$, respectively. The effective temperatures of the two components are $9940^{+260}_{-230}$ K and $15080\pm400$ K, respectively, and we determine the age of the system to be $0.21^{+0.11}_{-0.08}$ Gyr. We use the measured system properties to compute predicted phase curve amplitudes and find that while the measured Doppler boosting and mutual illumination components agree well with theory, the ellipsoidal distortion amplitude is significantly underestimated. We detail possible explanations for this discrepancy, including interactions between the dynamical tide of the host star and the tidal bulge and possible non-synchronous rotation of the host star.

We infer the number of planets-per-star as a function of orbital period and planet size using $Kepler$ archival data products with updated stellar properties from the $Gaia$ Data Release 2. Using hierarchical Bayesian modeling and Hamiltonian Monte Carlo, we incorporate planet radius uncertainties into an inhomogeneous Poisson point process model. We demonstrate that this model captures the general features of the outcome of the planet formation process around GK stars, and provides an infrastructure to use the $Kepler$ results to constrain analytic planet distribution models. We report an increased mean and variance in the marginal posterior distributions for the number of planets per $GK$ star when including planet radius measurement uncertainties. We estimate the number of planets-per-$GK$ star between 0.75 and 2.5 $R_{\oplus}$ and 50 to 300 day orbital periods to have a $68\%$ credible interval of $0.49$ to $0.77$ and a posterior mean of $0.63$. This posterior has a smaller mean and a larger variance than the occurrence rate calculated in this work and in Burke et al. (2015) for the same parameter space using the $Q1-Q16$ (previous $Kepler$ planet candidate and stellar catalog), and a larger mean and variance than when using the $DR25$ (latest $Kepler$ planet candidate and stellar catalog). We find that the accuracy and precision of our hierarchical Bayesian model posterior distributions are less sensitive to the total number of planets in the sample, and more so on the characteristics of the catalog completeness and reliability and the span of the planet parameter space.

We report detection of quasi-periodic (1.5 day) dimming of HD 240779, the solar-mass primary in a 5" visual binary (also TIC 284730577), by the Transiting Exoplanet Survey Satellite. This dimming, as has been shown for other "dipper" stars, is likely due to occultation by circumstellar dust. The barycentric space motion, lithium abundance, rotation, and chromospheric emission of the stars in this system point to an age of ~125 Myr, and possible membership in the AB Doradus moving group. As such it occupies an important but poorly explored intermediate regime of stars with transient dimming between young stellar objects in star forming regions and main sequence stars, and between UX Orionis-type Ae/Be stars and M-type "dippers". HD 240779, but not its companion BD+10714B, has WISE-detected excess infrared emission at 12 and 22 microns indicative of circumstellar dust. We propose that infrared emission is produced by collisions of planetesimals during clearing of a residual disk at the end of rocky planet formation, and that quasi-periodic dimming is produced by the rapid disintegration of a 100 km planetesimal near the silicate evaporation radius. Further studies of this and similar systems will illuminate a poorly understood final phase of rocky planet formation like that which produced the inner Solar System.

The radii and orbital periods of 4000+ confirmed/candidate exoplanets have been precisely measured by the Kepler mission. The radii show a bimodal distribution, with two peaks corresponding to smaller planets (likely rocky) and larger intermediate-size planets, respectively. While only the masses of the planets orbiting the brightest stars can be determined by ground-based spectroscopic observations, these observations allow calculation of their average densities placing constraints on the bulk compositions and internal structures. Yet an important question about the composition of planets ranging from 2 to 4 Earth radii still remains. They may either have a rocky core enveloped in a H2-He gaseous envelope (gas dwarfs) or contain a significant amount of multi-component, H2O-dominated ices/fluids (water worlds). Planets in the mass range of 10-15 Earth masses, if half-ice and half-rock by mass, have radii of 2.5 Earth radii, which exactly match the second peak of the exoplanet radius bimodal distribution. Any planet in the 2-4 Earth radii range requires a gas envelope of at most a few mass percentage points, regardless of the core composition. To resolve the ambiguity of internal compositions, we use a growth model and conduct Monte Carlo simulations to demonstrate that many intermediate-size planets are water worlds.

We report the first experimental results on spin-dependent elastic weakly interacting massive particle (WIMP) nucleon scattering from the XENON1T dark matter search experiment. The analysis uses the full ton year exposure of XENON1T to constrain the spin-dependent proton-only and neutron-only cases. No significant signal excess is observed, and a profile likelihood ratio analysis is used to set exclusion limits on the WIMP-nucleon interactions. This includes the most stringent constraint to date on the WIMP-neutron cross section, with a minimum of $6.3\times10^{-42}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level. The results are compared with those from collider searches and used to exclude new parameter space in an isoscalar theory with an axial-vector mediator.

The discovery of thousands of planetary systems by Kepler has demonstrated that planets are ubiquitous. However, a major challenge has been the confirmation of Kepler planet candidates, many of which still await confirmation. One of the most enigmatic examples is KOI 4.01, Kepler's first discovered planet candidate detection (as KOI 1.01, 2.01, and 3.01 were known prior to launch). Here we present the confirmation and characterization of KOI 4.01 (now Kepler-1658), using a combination of asteroseismology and radial velocities. Kepler-1658 is a massive, evolved subgiant (Mstar = 1.45 +/- 0.06 Msun, Rstar = 2.89 +/- 0.12 Rsun) hosting a massive (Mp = 5.88 +/- 0.47 MJ, Rp = 1.07 +/- 0.05 RJ) hot Jupiter that orbits every 3.85 days. Kepler-1658 joins a small population of evolved hosts with short-period (<=100 days) planets and is now the closest known planet in terms of orbital period to an evolved star. Because of its uniqueness and short orbital period, Kepler-1658 is a new benchmark system for testing tidal dissipation and hot Jupiter formation theories. Using all 4 years of Kepler data, we constrain the orbital decay rate to be Pdot <= -0.42 s/yr, corresponding to a strong observational limit of Qstar >= 4.826 x 10^3 for the tidal quality factor in evolved stars. With an effective temperature Teff ~6200 K, Kepler-1658 sits close to the spin-orbit misalignment boundary at ~6250 K, making it a prime target for follow-up observations to better constrain its obliquity and to provide insight into theories for hot Jupiter formation and migration.

We present first results on the scalar WIMP-pion coupling from 1 t$\times$yr of exposure with the XENON1T experiment. This interaction is generated when the WIMP couples to a virtual pion exchanged between the nucleons in a nucleus. In contrast to most non-relativistic operators, these pion-exchange currents can be coherently enhanced by the total number of nucleons, and therefore may dominate in scenarios where spin-independent WIMP-nucleon interactions are suppressed. Moreover, for natural values of the couplings, they dominate over the spin-dependent channel due to their coherence in the nucleus. Using the signal model of this new WIMP-pion channel, no significant excess is found, leading to an upper limit cross section of $6.4\times10^{-46}$ cm$^2$ (90 % confidence level) at 30 GeV/c$^2$ WIMP mass.

The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed.

One bottleneck for the exploitation of data from the $Kepler$ mission for stellar astrophysics and exoplanet research has been the lack of precise radii and evolutionary states for most of the observed stars. We report revised radii of 177,911 $Kepler$ stars derived by combining parallaxes from $Gaia$ Data Release 2 with the DR25 $Kepler$ Stellar Properties Catalog. The median radius precision is $\approx$ 8%, a typical improvement by a factor of 4-5 over previous estimates for typical $Kepler$ stars. We find that $\approx$ 67% ($\approx$ 120,000) of all $Kepler$ targets are main-sequence stars, $\approx$ 21% ($\approx$ 37,000) are subgiants, and $\approx$ 12% ($\approx$ 21,000) are red giants, demonstrating that subgiant contamination is less severe than some previous estimates and that Kepler targets are mostly main-sequence stars. Using the revised stellar radii, we recalculate the radii for 2123 confirmed and 1922 candidate exoplanets. We confirm the presence of a gap in the radius distribution of small, close-in planets, but find that the gap is mostly limited to incident fluxes $>$ 200$F_\oplus$ and its location may be at a slightly larger radius (closer to $\approx$ 2$R_\oplus$) when compared to previous results. Further, we find several confirmed exoplanets occupying a previously-described "hot super-Earth desert" at high irradiance, show the relation between gas-giant planet radius and incident flux, and establish a bona-fide sample of eight confirmed planets and 30 planet candidates with $R_{\mathrm{p}}$ $<$ 2$R_\oplus$ in circumstellar "habitable zones" (incident fluxes between 0.25--1.50 $F_\oplus$). The results presented here demonstrate the potential for transformative characterization of stellar and exoplanet populations using $Gaia$ data.

We report on a search for Weakly Interacting Massive Particles (WIMPs) using 278.8 days of data collected with the XENON1T experiment at LNGS. XENON1T utilizes a liquid xenon time projection chamber with a fiducial mass of $(1.30 \pm 0.01)$ t, resulting in a 1.0 t$\times$yr exposure. The energy region of interest, [1.4, 10.6] $\mathrm{keV_{ee}}$ ([4.9, 40.9] $\mathrm{keV_{nr}}$), exhibits an ultra-low electron recoil background rate of $(82\substack{+5 \\ -3}\textrm{ (sys)}\pm3\textrm{ (stat)})$ events/$(\mathrm{t}\times\mathrm{yr}\times\mathrm{keV_{ee}})$. No significant excess over background is found and a profile likelihood analysis parameterized in spatial and energy dimensions excludes new parameter space for the WIMP-nucleon spin-independent elastic scatter cross-section for WIMP masses above 6 GeV/c${}^2$, with a minimum of $4.1\times10^{-47}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level.

We present a single-dish mapping algorithm with a number of advantages over traditional techniques. (1) Our algorithm makes use of weighted modeling, instead of weighted averaging, to interpolate between signal measurements. This smooths the data, but without blurring the data beyond instrumental resolution. Techniques that rely on weighted averaging blur point sources sometimes as much as 40%. (2) Our algorithm makes use of local, instead of global, modeling to separate astronomical signal from instrumental and/or environmental signal drift along the telescope's scans. Other techniques, such as basket weaving, model this drift with simple functional forms (linear, quadratic, etc.) across the entirety of scans, limiting their ability to remove such contaminants. (3) Our algorithm makes use of a similar, local modeling technique to separate astronomical signal from radio-frequency interference (RFI), even if only continuum data are available. (4) Unlike other techniques, our algorithm does not require data to be collected on a rectangular grid or regridded before processing. (5) Data from any number of observations, overlapping or not, may be appended and processed together. (6) Any pixel density may be selected for the final image. We present our algorithm, and evaluate it using both simulated and real data. We are integrating it into the image-processing library of the Skynet Robotic Telescope Network, which includes optical telescopes spanning four continents, and now also Green Bank Observatory's 20-meter diameter radio telescope in West Virginia. Skynet serves hundreds of professional users, and additionally tens of thousands of students, of all ages. Default data products are generated on the fly, but will soon be customizable after the fact.

Sigma clipping is commonly used in astronomy for outlier rejection, but the number of standard deviations beyond which one should clip data from a sample ultimately depends on the size of the sample. Chauvenet rejection is one of the oldest, and simplest, ways to account for this, but, like sigma clipping, depends on the sample's mean and standard deviation, neither of which are robust quantities: Both are easily contaminated by the very outliers they are being used to reject. Many, more robust measures of central tendency, and of sample deviation, exist, but each has a tradeoff with precision. Here, we demonstrate that outlier rejection can be both very robust and very precise if decreasingly robust but increasingly precise techniques are applied in sequence. To this end, we present a variation on Chauvenet rejection that we call "robust" Chauvenet rejection (RCR), which uses three decreasingly robust/increasingly precise measures of central tendency, and four decreasingly robust/increasingly precise measures of sample deviation. We show this sequential approach to be very effective for a wide variety of contaminant types, even when a significant -- even dominant -- fraction of the sample is contaminated, and especially when the contaminants are strong. Furthermore, we have developed a bulk-rejection variant, to significantly decrease computing times, and RCR can be applied both to weighted data, and when fitting parameterized models to data. We present aperture photometry in a contaminated, crowded field as an example. RCR may be used by anyone at https://skynet.unc.edu/rcr, and source code is available there as well.

We present the spectral analysis of 21 blue supergiant stars of spectral type late B to early A within the Local Group dwarf galaxy IC 1613 based on VLT-FORS2 low resolution spectra. Combining our results with studies of early B type blue supergiants we report a wide bi-modal distribution of metallicities with two peaks around [Z] $\sim -0.50$ dex and [Z] $\sim -0.85$ dex. The bi-modal distribution correlates with spatial location, when compared with column densities of neutral hydrogen in IC 1613. While the low [Z] objects appear in regions of relatively high ISM HI column densities or close to them, the high [Z] supergiants are found in the central HI hole almost devoid of hydrogen. This suggests varied chemical evolution histories of the young stellar populations in IC 1613. Utilizing the Flux-Weighted Gravity-Luminosity Relation (FGLR), we determine IC 1613's distance modulus as $m - M$ = 24.39 $\pm$ 0.11 mag. This value is in agreement within previous distance measurements using the near-infrared period-luminosity relationship of Cepheids and the tip of the red giant branch.

In this paper, we describe the XENON100 data analyses used to assess the target-intrinsic background sources radon ($^{222}$Rn), thoron ($^{220}$Rn) and krypton ($^{85}$Kr). We detail the event selections of high-energy alpha particles and decay-specific delayed coincidences. We derive distributions of the individual radionuclides inside the detector and quantify their abundances during the main three science runs of the experiment over a period of $\sim$ 4 years, from January 2010 to January 2014. We compare our results to external measurements of radon emanation and krypton concentrations where we find good agreement. We report an observed reduction in concentrations of radon daughters that we attribute to the plating-out of charged ions on the negatively biased cathode.

The lithium doublet at 6708 \AA\ provides an age diagnostic for main sequence FGK dwarfs. We measured the abundance of lithium in 1305 stars with detected transiting planets from the Kepler Mission using high-resolution spectroscopy. Our catalog of lithium measurements from this sample have a range of abundance from A(Li) = 3.11 $\pm$ 0.07 to an upper limit of $-$0.84 dex. For a magnitude-limited sample that comprises 960 of the 1305 stars, our Keck-HIRES spectra have a median S/N = 45 per pixel at $\sim$6700 \AA\ with spectral resolution $\frac{\lambda}{\Delta \lambda}$ = $R$ = 55,000. We identify 80 young stars that have A(Li) values greater than the Hyades at their respective effective temperatures; these stars are younger than $\sim$650 Myr old, the approximate age of the Hyades. We then compare the distribution of A(Li) with planet size, multiplicity, orbital period, and insolation flux. We find larger planets preferentially in younger systems, with an A-D two-sided test p-value = 0.002, a $>3\sigma$ confidence that the older and younger planet samples do not come from the same parent distribution. This is consistent with planet inflation/photoevaporation at early ages. The other planet parameters ($Kepler$ planet multiplicity, orbital period, and insolation flux) are uncorrelated with age.

We report on the response of liquid xenon to low energy electronic recoils below 15 keV from beta decays of tritium at drift fields of 92 V/cm, 154 V/cm and 366 V/cm using the XENON100 detector. A data-to-simulation fitting method based on Markov Chain Monte Carlo is used to extract the photon yields and recombination fluctuations from the experimental data. The photon yields measured at the two lower fields are in agreement with those from literature; additional measurements at a higher field of 366 V/cm are presented. The electronic and nuclear recoil discrimination as well as its dependence on the drift field and photon detection efficiency are investigated at these low energies. The results provide new measurements in the energy region of interest for dark matter searches using liquid xenon.

We report the first dark matter search results from XENON1T, a $\sim$2000-kg-target-mass dual-phase (liquid-gas) xenon time projection chamber in operation at the Laboratori Nazionali del Gran Sasso in Italy and the first ton-scale detector of this kind. The blinded search used 34.2 live days of data acquired between November 2016 and January 2017. Inside the (1042$\pm$12) kg fiducial mass and in the [5, 40] $\mathrm{keV}_{\mathrm{nr}}$ energy range of interest for WIMP dark matter searches, the electronic recoil background was $(1.93 \pm 0.25) \times 10^{-4}$ events/(kg $\times$ day $\times \mathrm{keV}_{\mathrm{ee}}$), the lowest ever achieved in a dark matter detector. A profile likelihood analysis shows that the data is consistent with the background-only hypothesis. We derive the most stringent exclusion limits on the spin-independent WIMP-nucleon interaction cross section for WIMP masses above 10 GeV/c${}^2$, with a minimum of 7.7 $\times 10^{-47}$ cm${}^2$ for 35-GeV/c${}^2$ WIMPs at 90% confidence level.

We analyze solar quiescent prominence bubble characteristics and instability dynamics using Hinode/Solar Optical Telescope (SOT) data. We measure bubble expansion rate, prominence downflows, and the profile of the boundary layer brightness and thickness as a function of time. The largest bubble analyzed rises into the prominence with a speed of about 1.3 km/s until it is destabilized by a localized shear flow on the boundary. Boundary layer thickness grows gradually as prominence downflows deposit plasma onto the bubble with characteristic speeds of 20 - 35 km/s. Lateral downflows initiate from the thickened boundary layer with characteristic speeds of 25 - 50 km/s, "draining" the layer of plasma. Strong shear flow across one bubble boundary leads to an apparent coupled Kelvin-Helmholtz Rayleigh-Taylor (KH-RT) instability. We measure shear flow speeds above the bubble of 10 km/s and infer interior bubble flow speeds on the order of 100 km/s. Comparing the measured growth rate of the instability to analytic expressions, we infer a magnetic flux density across the bubble boundary of ~10^{-3} T (10 gauss) at an angle of ~70 degrees to the prominence plane. The results are consistent with the hypothesis that prominence bubbles are caused by magnetic flux that emerges below a prominence, setting up the conditions for RT, or combined KH-RT, instability flows that transport flux, helicity, and hot plasma upward into the overlying coronal magnetic flux rope.

We present the first search for dark matter-induced delayed coincidence signals in a dual-phase xenon time projection chamber, using the 224.6 live days of the XENON100 science run II. This very distinct signature is predicted in the framework of magnetic inelastic dark matter which has been proposed to reconcile the modulation signal reported by the DAMA/LIBRA collaboration with the null results from other direct detection experiments. No candidate event has been found in the region of interest and upper limits on the WIMP's magnetic dipole moment are derived. The scenarios proposed to explain the DAMA/LIBRA modulation signal by magnetic inelastic dark matter interactions of WIMPs with masses of 58.0 GeV/c$^2$ and 122.7 GeV/c$^2$ are excluded at 3.3 $\sigma$ and 9.3 $\sigma$, respectively.

We present results of searches for vector and pseudo-scalar bosonic super-WIMPs, which are dark matter candidates with masses at the keV-scale, with the XENON100 experiment. XENON100 is a dual-phase xenon time projection chamber operated at the Laboratori Nazionali del Gran Sasso. A profile likelihood analysis of data with an exposure of 224.6 live days $\times$ 34\,kg showed no evidence for a signal above the expected background. We thus obtain new and stringent upper limits in the $(8-125)$\,keV/c$^2$ mass range, excluding couplings to electrons with coupling constants of $g_{ae} > 3\times10^{-13}$ for pseudo-scalar and $\alpha'/\alpha > 2\times10^{-28}$ for vector super-WIMPs, respectively. These limits are derived under the assumption that super-WIMPs constitute all of the dark matter in our galaxy.

The XENON1T experiment at the Laboratori Nazionali del Gran Sasso (LNGS) is the first WIMP dark matter detector operating with a liquid xenon target mass above the ton-scale. Out of its 3.2t liquid xenon inventory, 2.0t constitute the active target of the dual-phase time projection chamber. The scintillation and ionization signals from particle interactions are detected with low-background photomultipliers. This article describes the XENON1T instrument and its subsystems as well as strategies to achieve an unprecedented low background level. First results on the detector response and the performance of the subsystems are also presented.

We report on WIMP search results in the XENON100 detector using a non-relativistic effective field theory approach. The data from science run II (34 kg $\times$ 224.6 live days) was re-analyzed, with an increased recoil energy interval compared to previous analyses, ranging from $(6.6 - 240)~\mathrm{keV_\mathrm{nr}}$. The data is found to be compatible with the background-only hypothesis. We present 90% confidence level exclusion limits on the coupling constants of WIMP-nucleon effective operators using a binned profile likelihood method. We also consider the case of inelastic WIMP scattering, where incident WIMPs may up-scatter to a higher mass state, and set exclusion limits on this model as well.

We present the first constraints on the spin-dependent, inelastic scattering cross section of Weakly Interacting Massive Particles (WIMPs) on nucleons from XENON100 data with an exposure of 7.64$\times$10$^3$\,kg\,day. XENON100 is a dual-phase xenon time projection chamber with 62\,kg of active mass, operated at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy and designed to search for nuclear recoils from WIMP-nucleus interactions. Here we explore inelastic scattering, where a transition to a low-lying excited nuclear state of $^{129}$Xe is induced. The experimental signature is a nuclear recoil observed together with the prompt de-excitation photon. We see no evidence for such inelastic WIMP-$^{129}$Xe interactions. A profile likelihood analysis allows us to set a 90\% C.L. upper limit on the inelastic, spin-dependent WIMP-nucleon cross section of $3.3 \times 10^{-38}$\,cm$^{2}$ at 100\,GeV/c$^2$. This is the most constraining result to date, and sets the pathway for an analysis of this interaction channel in upcoming, larger dual-phase xenon detectors.

Based on five years of observations with the 40-foot telescope at Green Bank Observatory (GBO), Reichart & Stephens (2000) found that the radio source Cassiopeia A had either faded more slowly between the mid-1970s and late 1990s than Baars et al. (1977) had found it to be fading between the late 1940s and mid-1970s, or that it had rebrightened and then resumed fading sometime between the mid-1970s and mid-1990s, in L band (1.4 GHz). Here, we present 15 additional years of observations of Cas A and Cyg A with the 40-foot in L band, and three and a half additional years of observations of Cas A, Cyg A, Tau A, and Vir A with GBO's recently refurbished 20-meter telescope in L and X (9 GHz) bands. We also present a more sophisticated analysis of the 40-foot data, and a reanalysis of the Baars et al. (1977) data, which reveals small, but non-negligible differences. We find that overall, between the late 1950s and late 2010s, Cas A faded at an average rate of $0.670 \pm 0.019$ %/yr in L band, consistent with Reichart & Stephens (2000). However, we also find, at the 6.3$\sigma$ credible level, that it did not fade at a constant rate. Rather, Cas A faded at a faster rate through at least the late 1960s, rebrightened (or at least faded at a much slower rate), and then resumed fading at a similarly fast rate by, at most, the late 1990s. Given these differences from the original Baars et al. (1977) analysis, and given the importance of their fitted spectral and temporal models for flux-density calibration in radio astronomy, we update and improve on these models for all four of these radio sources. In doing so, we additionally find that Tau A is fading at a rate of $0.102^{+0.042}_{-0.043}$ %/yr in L band.

We report on WIMP search results of the XENON100 experiment, combining three runs summing up to 477 live days from January 2010 to January 2014. Data from the first two runs were already published. A blind analysis was applied to the last run recorded between April 2013 and January 2014 prior to combining the results. The ultra-low electromagnetic background of the experiment, ~$5 \times 10^{-3}$ events/(keV$_{\mathrm{ee}}\times$kg$\times$day) before electronic recoil rejection, together with the increased exposure of 48 kg $\times$ yr improves the sensitivity. A profile likelihood analysis using an energy range of (6.6 - 43.3) keV$_{\mathrm{nr}}$ sets a limit on the elastic, spin-independent WIMP-nucleon scattering cross section for WIMP masses above 8 GeV/$c^2$, with a minimum of 1.1 $\times 10^{-45}$ cm$^2$ at 50 GeV/$c^2$ and 90% confidence level. We also report updated constraints on the elastic, spin-dependent WIMP-nucleon cross sections obtained with the same data. We set upper limits on the WIMP-neutron (proton) cross section with a minimum of 2.0 $\times 10^{-40}$ cm$^2$ (52$\times 10^{-40}$ cm$^2$) at a WIMP mass of 50 GeV/$c^2$, at 90% confidence level.

We report on a search for electronic recoil event rate modulation signatures in the XENON100 data accumulated over a period of 4 years, from January 2010 to January 2014. A profile likelihood method, which incorporates the stability of the XENON100 detector and the known electronic recoil background model, is used to quantify the significance of periodicity in the time distribution of events. There is a weak modulation signature at a period of $431^{+16}_{-14}$ days in the low energy region of $(2.0-5.8)$ keV in the single scatter event sample, with a global significance of $1.9\,\sigma$, however no other more significant modulation is observed. The expected annual modulation of a dark matter signal is not compatible with this result. Single scatter events in the low energy region are thus used to exclude the DAMA/LIBRA annual modulation as being due to dark matter electron interactions via axial vector coupling at $5.7\,\sigma$.

We perform a low-mass dark matter search using an exposure of 30\,kg$\times$yr with the XENON100 detector. By dropping the requirement of a scintillation signal and using only the ionization signal to determine the interaction energy, we lowered the energy threshold for detection to 0.7\,keV for nuclear recoils. No dark matter detection can be claimed because a complete background model cannot be constructed without a primary scintillation signal. Instead, we compute an upper limit on the WIMP-nucleon scattering cross section under the assumption that every event passing our selection criteria could be a signal event. Using an energy interval from 0.7\,keV to 9.1\,keV, we derive a limit on the spin-independent WIMP-nucleon cross section that excludes WIMPs with a mass of 6\,GeV/$c^2$ above $1.4 \times 10^{-41}$\,cm$^2$ at 90\% confidence level.

DARk matter WImp search with liquid xenoN (DARWIN) will be an experiment for the direct detection of dark matter using a multi-ton liquid xenon time projection chamber at its core. Its primary goal will be to explore the experimentally accessible parameter space for Weakly Interacting Massive Particles (WIMPs) in a wide mass-range, until neutrino interactions with the target become an irreducible background. The prompt scintillation light and the charge signals induced by particle interactions in the xenon will be observed by VUV sensitive, ultra-low background photosensors. Besides its excellent sensitivity to WIMPs above a mass of 5 GeV/c2, such a detector with its large mass, low-energy threshold and ultra-low background level will also be sensitive to other rare interactions. It will search for solar axions, galactic axion-like particles and the neutrinoless double-beta decay of 136-Xe, as well as measure the low-energy solar neutrino flux with <1% precision, observe coherent neutrino-nucleus interactions, and detect galactic supernovae. We present the concept of the DARWIN detector and discuss its physics reach, the main sources of backgrounds and the ongoing detector design and R&D efforts.

The XENON1T experiment is currently in the commissioning phase at the Laboratori Nazionali del Gran Sasso, Italy. In this article we study the experiment's expected sensitivity to the spin-independent WIMP-nucleon interaction cross section, based on Monte Carlo predictions of the electronic and nuclear recoil backgrounds. The total electronic recoil background in $1$ tonne fiducial volume and ($1$, $12$) keV electronic recoil equivalent energy region, before applying any selection to discriminate between electronic and nuclear recoils, is $(1.80 \pm 0.15) \cdot 10^{-4}$ ($\rm{kg} \cdot day \cdot keV)^{-1}$, mainly due to the decay of $^{222}\rm{Rn}$ daughters inside the xenon target. The nuclear recoil background in the corresponding nuclear recoil equivalent energy region ($4$, $50$) keV, is composed of $(0.6 \pm 0.1)$ ($\rm{t} \cdot y)^{-1}$ from radiogenic neutrons, $(1.8 \pm 0.3) \cdot 10^{-2}$ ($\rm{t} \cdot y)^{-1}$ from coherent scattering of neutrinos, and less than $0.01$ ($\rm{t} \cdot y)^{-1}$ from muon-induced neutrons. The sensitivity of XENON1T is calculated with the Profile Likelihood Ratio method, after converting the deposited energy of electronic and nuclear recoils into the scintillation and ionization signals seen in the detector. We take into account the systematic uncertainties on the photon and electron emission model, and on the estimation of the backgrounds, treated as nuisance parameters. The main contribution comes from the relative scintillation efficiency $\mathcal{L}_\mathrm{eff}$, which affects both the signal from WIMPs and the nuclear recoil backgrounds. After a $2$ y measurement in $1$ t fiducial volume, the sensitivity reaches a minimum cross section of $1.6 \cdot 10^{-47}$ cm$^2$ at m$_\chi$=$50$ GeV/$c^2$.

Spectroscopic observations of prominence eruptions associated with coronal mass ejections (CMEs), although relatively rare, can provide valuable plasma and 3D geometry diagnostics. We report the first observations by the Interface Region Imaging Spectrograph (IRIS) mission of a spectacular fast CME/prominence eruption associated with an equivalent X1.6 flare on 2014 May 9. The maximum plane-of-sky and Doppler velocities of the eruption are 1200 and 460 km/s, respectively. There are two eruption components separated by ~200 km/s in Doppler velocity: a primary, bright component and a secondary, faint component, suggesting a hollow, rather than solid, cone-shaped distribution of material. The eruption involves a left-handed helical structure undergoing counter-clockwise (viewed top-down) unwinding motion. There is a temporal evolution from upward eruption to downward fallback with less-than-free-fall speeds and decreasing nonthermal line widths. We find a wide range of Mg II k/h line intensity ratios (less than ~2 expected for optically-thin thermal emission): the lowest ever-reported median value of 1.17 found in the fallback material and a comparably high value of 1.63 in nearby coronal rain and intermediate values of 1.53 and 1.41 in the two eruption components. The fallback material exhibits a strong ($> 5 \sigma$) linear correlation between the k/h ratio and the Doppler velocity as well as the line intensity. We demonstrate that Doppler dimming of scattered chromospheric emission by the erupted material can potentially explain such characteristics.

The Interface Region Imaging Spectrograph (IRIS) small explorer spacecraft provides simultaneous spectra and images of the photosphere, chromosphere, transition region, and corona with 0.33-0.4 arcsec spatial resolution, 2 s temporal resolution and 1 km/s velocity resolution over a field-of-view of up to 175 arcsec x 175 arcsec. IRIS was launched into a Sun-synchronous orbit on 27 June 2013 using a Pegasus-XL rocket and consists of a 19-cm UV telescope that feeds a slit-based dual-bandpass imaging spectrograph. IRIS obtains spectra in passbands from 1332-1358, 1389-1407 and 2783-2834 Angstrom including bright spectral lines formed in the chromosphere (Mg II h 2803 Angstrom and Mg II k 2796 Angstrom) and transition region (C II 1334/1335 Angstrom and Si IV 1394/1403 Angstrom). Slit-jaw images in four different passbands (C II 1330, Si IV 1400, Mg II k 2796 and Mg II wing 2830 Angstrom) can be taken simultaneously with spectral rasters that sample regions up to 130 arcsec x 175 arcsec at a variety of spatial samplings (from 0.33 arcsec and up). IRIS is sensitive to emission from plasma at temperatures between 5000 K and 10 MK and will advance our understanding of the flow of mass and energy through an interface region, formed by the chromosphere and transition region, between the photosphere and corona. This highly structured and dynamic region not only acts as the conduit of all mass and energy feeding into the corona and solar wind, it also requires an order of magnitude more energy to heat than the corona and solar wind combined. The IRIS investigation includes a strong numerical modeling component based on advanced radiative-MHD codes to facilitate interpretation of observations of this complex region. Approximately eight Gbytes of data (after compression) are acquired by IRIS each day and made available for unrestricted use within a few days of the observation.

The measurement of the Sun's surface motions with a high spatial and temporal resolution is still a challenge. We wish to validate horizontal velocity measurements all over the visible disk of the Sun from Solar Dynamics Observatory/ Helioseismic and Magnetic Imager (SDO/HMI) data. Horizontal velocity fields are measured by following the proper motions of solar granules using a newly developed version of the Coherent Structure Tracking (CST) code. The comparison of the surface flows measured at high spatial resolution (Hinode, 0.1 arcsec) and low resolution (SDO/HMI, 0.5 arcsec) allows us to determine corrections to be applied to the horizontal velocity measured from HMI white light data. We derive horizontal velocity maps with spatial and temporal resolutions of respectively 2.5 Mm and 30 min. From the two components of the horizontal velocity Vx and Vy measured in the sky plane and the simultaneous line of sight component from SDO/HMI dopplergrams v_D, we derive the spherical velocity components (Vr, Vtheta, Vphi). The azimuthal component Vphi gives the solar differential rotation with a high precision (+-0.037km/s) from a temporal sequence of only three hours. By following the proper motions of the solar granules, we can revisit the dynamics of the solar surface at high spatial and temporal resolutions from hours to months and years with the SDO data.

We report the first analyses of SDO/AIA observations of the formation of a quiescent polar crown prominence in a coronal cavity. The He II 304 \AA\ (log T_{max} ~ 4.8 K) data show both the gradual disappearance of the prominence due to vertical drainage and lateral transport of plasma followed by the formation of a new prominence some 12 hours later. The formation of the prominence is preceded by the appearance of a bright emission "cloud" in the central region of the coronal cavity. The peak brightness of the cloud progressively shifts in time from the Fe XIV 211 \AA\ channel, through the Fe XII 193 \AA\ channel, to the Fe IX 171 \AA\ channel (log T_{max} ~ 6.2, 6.1, 5.8 K, respectively) while simultaneously decreasing in altitude. Filter ratio analysis estimates the initial temperature of the cloud in the cavity to be approximately log T \sim 6.25 K with evidence of cooling over time. The subsequent growth of the prominence is accompanied by darkening of the cavity in the 211 \AA\ channel. The observations imply the possibility of prominence formation via in situ condensation of hot plasma from the coronal cavity, in support of the proposed process of magneto-thermal convection in coronal magnetic flux ropes.

We investigate the properties of acoustic events (AEs), defined as spatially concentrated and short duration energy flux, in the quiet sun using observations of a 2D field of view (FOV) with high spatial and temporal resolution provided by the Solar Optical Telescope (SOT) onboard \textit{Hinode}. Line profiles of Fe \textsc{i} 557.6 nm were recorded by the Narrow band Filter Imager (NFI) on a $82" \times 82"$ FOV during 75 min with a time step of 28.75 s and 0.08$"$ pixel size. Vertical velocities were computed at three atmospheric levels (80, 130 and 180 km) using the bisector technique allowing the determination of energy flux in the range 3-10 mHz using two complementary methods (Hilbert transform and Fourier power spectra). Horizontal velocities were computed using local correlation tracking (LCT) of continuum intensities providing divergences. The net energy flux is upward. In the range 3-10 mHz, a full FOV space and time averaged flux of 2700 W m$^{-2}$ (lower layer 80-130 km) and 2000 W m$^{-2}$ (upper layer 130-180 km) is concentrated in less than 1% of the solar surface in the form of narrow (0.3$"$) AE. Their total duration (including rise and decay) is of the order of $10^{3}$ s. Inside each AE, the mean flux is $1.6 10^{5}$ W m$^{-2}$ (lower layer) and $1.2 10^{5}$ W m$^{-2}$ (upper). Each event carries an average energy (flux integrated over space and time) of $2.5 10^{19}$ J (lower layer) to $1.9 10^{19}$ J (upper). More than $10^{6}$ events could exist permanently on the Sun, with a birth and decay rate of 3500 s$^{-1}$. Most events occur in intergranular lanes, downward velocity regions, and areas of converging motions.

This series of papers investigates the dynamic interior of a quiescent prominence revealed by recent {\it Hinode} and {\it SDO/AIA} high-resolution observations. This first paper is a study of the static equilibrium of the Kippenhahn-Schl\"{u}ter diffuse plasma slab, suspended vertically in a bowed magnetic field, under the frozen-in condition and subject to a theoretical thermal balance among an optically-thin radiation, heating, and field-aligned thermal conduction. The everywhere-analytical solutions to this nonlinear problem are an extremely restricted subset of the physically admissible states of the system. For most values of the total mass frozen into a given bowed field, force-balance and steady energy-transport cannot both be met without a finite fraction of the total mass having collapsed into a cold sheet of zero thickness, within which the frozen-in condition must break down. An exact, resistive hydromagnetic extension of the Kippenhahn-Schl\"{u}ter slab is also presented, resolving the mass-sheet singularity into a finite-thickness layer of steadily-falling dense fluid. Our hydromagnetic result suggests that the narrow, vertical prominence $H_{\alpha}$ threads may be falling across magnetic fields, with optically-thick cores much denser and ionized to much lower degrees than conventionally considered. This implication is discussed in relation to (i) the recent {\it SDO/AIA} observations of quiescent prominences that are massive and yet draining mass everywhere in their interiors, (ii) the canonical range of $5-60 G$ determined from spectral-polarimetric observations of prominence magnetic fields over the years and (iii) the need for a more realistic multi-fluid treatment.

For the first time, the motion of granules (solar plasma on the surface on scales larger than 2.5 Mm) has been followed over the entire visible surface of the Sun, using SDO/HMI white-light data. Horizontal velocity fields are derived from image correlation tracking using a new version of the coherent structure tracking algorithm.The spatial and temporal resolutions of the horizontal velocity map are 2.5 Mm and 30 min respectively . From this reconstruction, using the multi-resolution analysis, one can obtain to the velocity field at different scales with its derivatives such as the horizontal divergence or the vertical component of the vorticity. The intrinsic error on the velocity is ~0.25 km/s for a time sequence of 30 minutes and a mesh size of 2.5 Mm.This is acceptable compared to the granule velocities, which range between 0.3 km/s and 1.8 km/s. A high correlation between velocities computed from Hinode and SDO/HMI has been found (85%). From the data we derive the power spectrum of the supergranulation horizontal velocity field, the solar differential rotation, and the meridional velocity.

Observations of solar prominences reveal a complex, dynamic flow field within them. The flow field within quiescent prominences is characterized by long ``threads'' and dark ``bubbles'' that fall and rise (respectively) in a thin sheet. The flow field in active prominences display more helical motions that travel along the axis of the prominence. We explore the possible dynamics of both of these with the aid of 2.5D MHD simulations. Our model, compressible plasma possesses density and temperature gradients and resides in magnetic field configurations that mimc those of a solar prominence. We present results of various configurations and discuss the nonlinear behavior of the resulting dynamics.

Imaging solar coronal condensation forming prominences was difficult in the past, a situation recently changed by Hinode and SDO. We present the first example observed with SDO/AIA, in which material gradually cools through multiple EUV channels in a transequatorial loop system that confines an earlier eruption. Nine hours later, this leads to eventual condensation at the dip of these loops, forming a moderate-size prominence of ~$10^{14}$ gram, to be compared to the characteristic $10^{15}$ gram mass of a CME. The prominence mass is not static but maintained by condensation at a high estimated rate of $10^{10}$ gram/sec against a comparable, sustained drainage through numerous vertical downflow threads, such that 96% of the total condensation (~$10^{15}$ gram) is drained in approximately one day. The mass condensation and drainage rates temporally correlate with the total prominence mass. The downflow velocity has a narrow Gaussian distribution with a mean of 30 km/s, while the downward acceleration distribution has an exponential drop with a mean of ~1/6 $g_{Sun}$, indicating a significant canceling of gravity, possibly by the Lorentz force. Our observations show that a macroscopic quiescent prominence is microscopically dynamic, involving the passage of a significant mass through it, maintained by a continual mass supply against a comparable mass drainage, which bears important implications for CME initiation mechanisms in which mass unloading is important.

The launch of the Hinode satellite has allowed unprecedented high resolution, stable images of solar quiescent prominences to be taken over extended periods of time. These new images led to the discovery of dark upflows that propagated from the base of prominences developing highly turbulent profiles. As yet, how these flows are driven is not fully understood. To study the physics behind this phenomena we use 3-D magnetohydrodynamic (MHD) simulations to investigate the nonlinear stability of the Kippenhahn-Shl\"{u}ter prominence model to the magnetic Rayleigh-Taylor instability. The model simulates the rise of a buoyant tube inside a quiescent prominence, where the upper boundary between the tube and prominence model is perturbed to excite the interchange of magnetic field lines. We found upflows of constant velocity (maximum found 6\,km\,s$^{-1}$) and a maximum plume width $\approx 1500$\,km which propagate through a height of approximately 6\,Mm in the no guide field case. The case with the strong guide field (initially $B_y=2B_x$) results in a large plume that rises though the prominence model at $\sim 5$\,km\,s$^{-1}$ with width $\sim 900$\,km (resulting in width of 2400\,km when viewed along the axis of the prominence) reaching a height of $\sim 3.1$\,Mm. In both cases nonlinear processes were important for determining plume dynamics.

We present observations of a chromospheric jet and growing "loop" system that show new evidence of a fan-spine topology resulting from magnetic flux emergence. This event, occurring in an equatorial coronal hole on 2007 February 9, was observed by the Hinode Solar Optical Telescope in the Ca II H line in unprecedented detail. The predecessor of the jet is a bundle of fine material threads that extend above the chromosphere and appear to rotate about the bundle axis at ~50 km/s (period <200 s). These rotations or transverse oscillations propagate upward at velocities up to 786 km/s. The bundle first slowly and then rapidly swings up, with the transition occurring at the onset of an A4.9 flare. A loop expands simultaneously in these two phases (velocity: 16-135 km/s). Near the peak of the flare, the loop appears to rupture; simultaneous upward ejecta and mass downflows faster than free-fall appear in one of the loop legs. The material bundle then swings back in a whiplike manner and develops into a collimated jet, which is orientated along the inferred open field lines with transverse oscillations continuing at slower rates. Some material falls back along smooth streamlines, showing no more oscillations. At low altitudes, the streamlines bifurcate at presumably a magnetic null point and bypass an inferred dome, depicting an inverted-Y geometry. These streamlines closely match in space the late Ca II H loop and X-ray flare loop. These observations are consistent with the model that flux emergence in an open-field region leads to magnetic reconnection, forming a jet and fan-spine topology. We propose that the material bundle and collimated jet represent the outer spine in quasi-static and eruptive stages, respectively, and the growing loop is a 2D projection of the 3D fan surface.

Continuous observations were performed of a quiescent prominence with the Solar Optical Telescope (SOT) on board the /emph{Hinode} satellite on 2006 December 23--24. A peculiar slowly-rising column of $/sim10^{4}$ K plasma develops from the lower atmosphere during the observations. The apparent ascent speed of the column is 2 km s$^{-1}$, while the fine structures of the column exhibit much faster motion of up to 20 km s$^{-1}$. The column eventually becomes a faint low-lying prominence. Associated with the appearance of the column, an overlying coronal cavity seen in the X-ray and EUV moves upward at $/sim$5 km s$^{-1}$. We discuss the relationship between these episodes, and suggest that they are due to the emergence of a helical flux rope that undergoes reconnection with lower coronal fields, possibly carrying material into the coronal cavity. Under the assumption of the emerging flux scenario, the lower velocity of 2 km s$^{-1}$ and the higher one of 20 km s$^{-1}$ in the column are attributed to the rising motion of the emerging flux and to the outflow driven by magnetic reconnection between the emerging flux and the pre-existing coronal field, respectively. The present paper gives a coherent explanation of the enigmatic phenomenon of the rising column with the emergence of the helical rope, and its effect on the corona. We discuss the implications that the emergence of such a helical rope has on the dynamo process in the convection zone.

The aim of this work is to give new observational constraints on solar surface flows by determining the horizontal scale dependence of the velocity and intensity fields, as represented by their power spectra, and to offer some theoretical guidelines to interpret these spectra. We use long time series of images taken by SOT/Hinode and reconstruct both horizontal (by granule tracking) and vertical (by Doppler effect) velocity fields in a field of view 75x75Mm^2. At small sub-granulation scales, the kinetic energy spectral density associated with vertical motions exhibits a k^{-10/3}-like spectrum, while the intensity fluctuation spectrum follows a k^{-3} or k^{-17/3}-like spectrum at the two continuum levels investigated (525 and 450 nm respectively). We discuss the physical origin of these scalings and argue that they provide a direct observational signature of buoyancy-driven turbulent dynamics in a strongly thermally diffusive regime. In the mesogranulation range and up to a scale of 25Mm, we find that the vertical velocity field amplitude decreases like L^{-3/2} with the horizontal scale L. This behaviour corresponds to a k^2 spectral power law. Still in the mesoscale range, [2.5, 10]Mm, we find that intensity fluctuations in the blue continuum follow a power law in k^2. We show that granule tracking cannot sample scales below 2.5Mm. We locate the supergranulation energy peak around 30Mm and show that the emergence of a pore erases this spectral peak. Thanks to a scale height estimate, we find that supergranular flows are shallow. (abridged abstract).

We report a chromospheric jet lasting for more than 1 hr observed by Hinode Solar Optical Telescope in unprecedented detail. The ejection occurred in three episodes separated by 12-14 min, with the amount and velocity of material decreasing with time. The upward velocities range from 438 to 33 km/s, while the downward velocities of the material falling back have smaller values (mean: -56 km/s) and a narrower distribution (standard deviation: 14 km/s). The average acceleration inferred from parabolic space-time tracks is 141 m/s^2, a fraction of the solar gravitational acceleration. The jet consists of fine threads (0.5-2.0 arcsec wide), which exhibit coherent, oscillatory transverse motions perpendicular to the jet axis and about a common equilibrium position. These motions propagate upward along the jet, with the maximum phase speed of 744 +/- 11 km/s at the leading front of the jet. The transverse oscillation velocities range from 151 to 26 km/s, amplitudes from 6.0 to 1.9 Mm, and periods from 250 to 536 s. The oscillations slow down with time and cease when the material starts to fall back. The falling material travels along almost straight lines in the original direction of ascent, showing no transverse motions. These observations are consistent with the scenario that the jet involves untwisting helical threads, which rotate about the axis of a single large cylinder and shed magnetic helicity into the upper atmosphere.

There are significant advantages in combining Hinode observations with ground-based instruments that can observe additional spectral diagnostics at higher data rates and with greater flexibility. However, ground-based observations, because of the random effects of weather and seeing as well as the complexities data analysis due to changing instrumental configurations, have traditionally been less efficient than satellite observations in producing useful datasets. Future large ground-based telescopes will need to find new ways to optimize both their operational efficiency and scientific output. We have begun experimenting with service-mode or queue-mode observations at the Dunn Solar Telescope using the Interferometric Bidimensional Spectrometer (IBIS) as part of joint Hinode campaigns. We describe our experiences and the advantages of such an observing mode for solar physics.

The formation and evolution process and magnetic configuration of solar prominences remain unclear. In order to study the formation process of prominences, we examine continuous observations of a prominence in NOAA AR 10953 with the Solar Optical Telescope on the \emph{Hinode} satellite. As reported in our previous Letter, we find a signature suggesting that a helical flux rope emerges from below the photosphere under a pre-existing prominence. Here we investigate more detailed properties and photospheric indications of the emerging helical flux rope, and discuss their relationship to the formation of the prominence. Our main conclusions are: (1) A dark region with absence of strong vertical magnetic fields broadens and then narrows in Ca \textsc{ii} H-line filtergrams. This phenomenon is consistent with the emergence of the helical flux rope as photospheric counterparts. The size of the flux rope is roughly 30,000 km long and 10,000 km wide. The width is larger than that of the prominence. (2) No shear motion or converging flows are detected, but we find diverging flows such as mesogranules along the polarity inversion line. The presence of mesogranules may be related to the emergence of the helical flux rope. (3) The emerging helical flux rope reconnects with magnetic fields of the pre-existing prominence to stabilize the prominence for the next several days. We thus conjecture that prominence coronal magnetic fields emerge in the form of helical flux ropes that contribute to the formation and maintenance of the prominence.

Aims: The interactions of velocity scales on the Sun's surface, from granulation to supergranulation are still not understood, nor are their interaction with magnetic fields. We thus aim at giving a better description of dynamics in the mesoscale range which lies between the two scales mentioned above. Method: We analyse a 48h high-resolution time sequence of the quiet Sun photosphere at the disk center obtained with the Solar Optical Telescope onboard Hinode. The observations, which have a field of view of 100 \arcsec$\times$ 100 \arcsec, typically contain four supergranules. We monitor in detail the motion and evolution of granules as well as those of the radial magnetic field. Results: This analysis allows us to better characterize Trees of Fragmenting Granules issued from repeated fragmentation of granules, especially their lifetime statistics. Using floating corks advected by measured velocity fields, we show their crucial role in the advection of the magnetic field and in the build up of the network. Finally, thanks to the long duration of the time series, we estimate that the turbulent diffusion coefficient induced by horizontal motion is approximately $430 \mathrm{km}^2 \mathrm{s}^{-1}$. Conclusions: These results demonstrate that the long living families contribute to the formation of the magnetic network and suggest that supergranulation could be an emergent length scale building up as small magnetic elements are advected and concentrated by TFG flows. Our estimate for the magnetic diffusion associated with this horizontal motion might provide a useful input for mean-field dynamo models.

The heating of the solar chromosphere and corona is a long-standing puzzle in solar physics. Hinode observations show the ubiquitous presence of chromospheric anemone jets outside sunspots in active regions. They are typically 3 to 7 arc seconds = 2000 to 5000 kilometers long and 0.2 to 0.4 arc second = 150 to 300 kilometers wide, and their velocity is 10 to 20 kilometers per second. These small jets have an inverted Y-shape, similar to the shape of x-ray anemone jets in the corona. These features imply that magnetic reconnection similar to that in the corona is occurring at a much smaller spatial scale throughout the chromosphere and suggest that the heating of the solar chromosphere and corona may be related to small-scale ubiquitous reconnection.

We use sequences of images and magnetograms from Hinode to study magnetic elements in internetwork parts of the quiet solar photosphere. Visual inspection shows the existence of many long-lived (several hours) structures that interact frequently, and may migrate over distances ~7 Mm over a period of a few hours. About a fifth of the elements have an associated bright point in G-band or Ca II H intensity. We apply a hysteresis-based algorithm to identify elements. The algorithm is able to track elements for about 10 min on average. Elements intermittently drop below the detection limit, though the associated flux apparently persists and often reappears some time later. We infer proper motions of elements from their successive positions, and find that they obey a Gaussian distribution with an rms of 1.57+-0.08 km/s. The apparent flows indicate a bias of about 0.2 km/s toward the network boundary. Elements of negative polarity show a higher bias than elements of positive polarity, perhaps as a result of to the dominant positive polarity of the network in the field of view, or because of increased mobility due to their smaller size. A preference for motions in X is likely explained by higher supergranular flow in that direction. We search for emerging bipoles by grouping elements of opposite polarity that appear close together in space and time. We find no evidence supporting Joy's law at arcsecond scales.

On-orbit performance of the Solar Optical Telescope (SOT) aboard Hinode is described with some attentions on its unpredicted aspects. In general, SOT revealed an excellent performance and has been providing outstanding data. Some unexpected features exist, however, in behaviors of the focus position, throughput and structural stability. Most of them are recovered by the daily operation i.e., frequent focus adjustment, careful heater setting and corrections in data analysis. The tunable filter contains air bubbles which degrade the data quality significantly. Schemes for tuning the filter without disturbing the bubbles have been developed and tested, and some useful procedures to obtain Dopplergram and magnetogram are now available. October and March when the orbit of satellite becomes nearly perpendicular to the direction towards the sun provide a favorable condition for continuous runs of the narrow-band filter imager.

Solar prominences are cool 10$^4$ Kelvin plasma clouds supported in the surrounding 10$^6$ Kelvin coronal plasma by as-yet undetermined mechanisms. Observations from \emph{Hinode} show fine-scale threadlike structures oscillating in the plane of the sky with periods of several minutes. We suggest these transverse magnetohydrodynamic waves may represent Alfv\'en waves propagating on coronal magnetic field lines and these may play a role in heating the corona.

Continuous observations were obtained of active region 10953 with the Solar Optical Telescope (SOT) on board the \emph{Hinode} satellite during 2007 April 28 to May 9. A prominence was located over the polarity inversion line (PIL) in the south-east of the main sunspot. These observations provided us with a time series of vector magnetic fields on the photosphere under the prominence. We found four features: (1) The abutting opposite-polarity regions on the two sides along the PIL first grew laterally in size and then narrowed. (2) These abutting regions contained vertically-weak, but horizontally-strong magnetic fields. (3) The orientations of the horizontal magnetic fields along the PIL on the photosphere gradually changed with time from a normal-polarity configuration to a inverse-polarity one. (4) The horizontal-magnetic field region was blueshifted. These indicate that helical flux rope was emerging from below the photosphere into the corona along the PIL under the pre-existing prominence. We suggest that this supply of a helical magnetic flux into the corona is associated with evolution and maintenance of active-region prominences.

The Solar Optical Telescope (SOT) on board Hinode satellite observed an X3.4 class flare on 2006 December 13. Typical two-ribbon structure was observed, not only in the chromospheric CaII H line but also in G-band and FeI 6302A line. The high-resolution, seeing-free images achieved by SOT revealed, for the first time, the sub-arcsec fine structures of the "white light" flare. The G-band flare ribbons on sunspot umbrae showed a sharp leading edge followed by a diffuse inside, as well as previously known core-halo structure. The underlying structures such as umbral dots, penumbral filaments and granules were visible in the flare ribbons. Assuming that the sharp leading edge was directly heated by particle beam and the diffuse parts were heated by radiative back-warming, we estimate the depth of the diffuse flare emission using the intensity profile of the flare ribbon. We found that the depth of the diffuse emission is about 100 km or less from the height of the source of radiative back-warming. The flare ribbons were also visible in the Stokes-V images of FeI 6302A, as a transient polarity reversal. This is probably related to "magnetic transient" reported in the literature. The intensity increase in Stokes-I images indicates that the FeI 6302A line was significantly deformed by the flare, which may cause such a magnetic transient.

The Solar Optical Telescope (SOT) aboard HINODE successfully and continuously observed a formation process of a light bridge in a matured sunspot of the NOAA active region 10923 for several days with high spatial resolution. During its formation, many umbral dots were observed emerging from the leading edges of penumbral filaments, and intruding into the umbra rapidly. The precursor of the light bridge formation was also identified as the relatively slow inward motion of the umbral dots which emerged not near the penumbra, but inside the umbra. The spectro-polarimeter on SOT provided physical conditions in the photosphere around the umbral dots and the light bridges. We found the light bridges and the umbral dots had significantly weaker magnetic fields associated with upflows relative to the core of the umbra, which implies that there was hot gas with weak field strength penetrating from subphotosphere to near the visible surface inside those structures. There needs to be a mechanism to drive the inward motion of the hot gas along the light bridges. We suggest that the emergence and the inward motion are triggered by a buoyant penumbral flux tube as well as the subphotospheric flow crossing the sunspot.

Results from initial helioseismic observations by Solar Optical Telescope onboard Hinode are reported. It has been demonstrated that intensity oscillation data from Broadband Filter Imager can be used for various helioseismic analyses. The k-omega power spectra, as well as corresponding time-distance cross-correlation function that promises high-resolution time-distance analysis below 6-Mm travelling distance, were obtained for G-band and CaII-H data. Subsurface supergranular patterns have been observed from our first time-distance analysis. The results show that the solar oscillation spectrum is extended to much higher frequencies and wavenumbers, and the time-distance diagram is extended to much shorter travel distances and times than they were observed before, thus revealing great potential for high-resolution helioseismic observations from Hinode.

The Altcriss project aims to perform a long term survey of the radiation environment on board the International Space Station. Measurements are being performed with active and passive devices in different locations and orientations of the Russian segment of the station. The goal is to perform a detailed evaluation of the differences in particle fluence and nuclear composition due to different shielding material and attitude of the station. The Sileye-3/Alteino detector is used to identify nuclei up to Iron in the energy range above 60 MeV/n. Several passive dosimeters (TLDs, CR39) are also placed in the same location of Sileye-3 detector. Polyethylene shielding is periodically interposed in front of the detectors to evaluate the effectiveness of shielding on the nuclear component of the cosmic radiation. The project was submitted to ESA in reply to the AO in the Life and Physical Science of 2004 and data taking began in December 2005. Dosimeters and data cards are rotated every six months: up to now three launches of dosimeters and data cards have been performed and have been returned with the end of expedition 12 and 13.

The Sun provides an excellent natural laboratory for nonlinear phenomena. We use motions of magnetic bright points on the solar surface, at the smallest scales yet observed, to study the small scale dynamics of the photospheric plasma. The paths of the bright points are analyzed within a continuous time random walk framework. Their spatial and temporal scaling suggest that the observed motions are the walks of imperfectly correlated tracers on a turbulent fluid flow in the lanes between granular convection cells.

We investigate the possibility of searching for ultra high energy neutrinos in cosmic rays using acoustic techniques in ocean water. The type of information provided by the acoustic detection is complementary to that of other techniques, and the filtering effect of the atmosphere, imposed by the fact that detection only happens if a shower fully develops in water, would provide a clear neutrino identification. We find that it may be possible to implement this technique with very limited resources using existing high frequency underwater hydrophone arrays. We review the expected acoustic signals produced by neutrino-induced showers in water and develop an optimal filtering algorithm able to suppress statistical noise. The algorithm found is computationally appropriate to be used as a trigger for the signal processors available on existing arrays. We estimate the noise rates for a trigger system on a very large size hydrophone array of the US Navy and find that, while a higher density of hydrophones would be desirable, the existing system may already provide useful data.