Search

This paper describes a new $\nu_e$ identification method specifically designed to improve the low-energy ($< 30\,\mathrm{GeV}$) $\nu_e$ identification efficiency attained by enlarging the emulsion film scanning volume with the next generation emulsion readout system. A relative increase of 25-70% in the $\nu_e$ low-energy region is expected, leading to improvements in the OPERA sensitivity to neutrino oscillations in the framework of the 3 + 1 model. The method is applied to a subset of data where the detection efficiency increase is expected to be more relevant, and one additional $\nu_e$ candidate is found. The analysis combined with the $\nu_\tau$ appearance results improves the upper limit on $\sin^2 2\theta_{\mu e}$ to 0.016 at 90% C.L. in the MiniBooNE allowed region $\Delta m^2_{41} \sim 0.3\,\mathrm{eV}^2$., Comment: 13 pages, 6 figures

This study conducted an analysis of muons, protons, and charged pions emitted from $\nu_{\mu}$ charged-current interactions on iron using a nuclear emulsion detector. The emulsion detector with a 65$\,$kg iron target was exposed to a neutrino beam corresponding to 4.0$\times$10$^{19}$ protons on target with a mean neutrino energy of 1.49$\,$GeV. The measurements were performed at a momentum threshold of 200 (50)$\,$MeV/$c$ for protons (pions), which are the lowest momentum thresholds attempted up to now. The measured quantities are the multiplicities, emission angles, and momenta of the muons, protons, and charged pions. In addition to these inclusive measurements, exclusive measurements such as the muon-proton emission-angle correlations of specific channels and the opening angle between the protons of CC0$\pi$2$p$ events were performed. The data were compared to Monte Carlo (MC) predictions and some significant differences were observed. The results of the study demonstrate the capability of detailed measurements of neutrino-nucleus interactions using a nuclear emulsion detector to improve neutrino interaction models., Comment: 21 pages, 25 figures

Context. The discovery of an extrasolar planet with an ocean has crucial importance in the search for life beyond Earth. The polarimetric detection of specularly reflected light from a smooth liquid surface is anticipated theoretically, though the polarimetric signature of Earth's oceans has not yet been conclusively detected in disk-integrated planetary light. Aims. We aim to detect and measure the polarimetric signature of the Earth's oceans. Methods. We conducted near-infrared polarimetry for lunar Earthshine and collected data on 32 nights with a variety of ocean fractions in the Earthshine-contributing region. Results. A clear positive correlation was revealed between the polarization degree and ocean fraction. We found hourly variations in polarization in accordance with rotational transition of the ocean fraction. The ratios of the variation to the typical polarization degree were as large as ~0.2-1.4. Conclusions. Our observations provide plausible evidence of the polarimetric signature attributed to Earth's oceans. Near-infrared polarimetry may be considered a prospective technique in the search for exoplanetary oceans., Comment: Published in A&A (20 pages, 14 figures, 3 tables). Note added in proof

We have carried out $\nu_{\mu}$ charged-current interaction measurement on iron using an emulsion detector exposed to the T2K neutrino beam in the J-PARC neutrino facility. The data samples correspond to 4.0$\times$10$^{19}$ protons on target, and the neutrino mean energy is 1.49$\,$GeV. The emulsion detector is suitable for precision measurements of charged particles produced in neutrino-iron interactions with a low momentum threshold thanks to thin-layered structure and sub-$\mu$m spatial resolution. The charged particles are successfully detected, and their multiplicities are measured using the emulsion detector. The cross section was measured to be $\sigma^{\mathrm{Fe}}_{\mathrm{CC}} = (1.28 \pm 0.11({\mathrm{stat.}})^{+0.12}_{-0.11}({\mathrm{syst.}})) \times 10^{-38} \, {\mathrm{cm}}^{2}/{\mathrm{nucleon}}$. The cross section in a limited kinematic phase space of induced muons, $\theta_{\mu} < 45^{\circ}$ and $p_{\mu} > 400 \, {\rm MeV}/c$, on iron was $\sigma^{\mathrm{Fe}}_{\mathrm{CC \hspace{1mm} phase \hspace{0.5mm} space}} = (0.84 \pm 0.07({\mathrm{stat.}})^{+0.07}_{-0.06}({\mathrm{syst.}})) \times 10^{-38} \, {\mathrm{cm}}^{2}/{\mathrm{nucleon}}$. The cross-section results are consistent with previous values obtained via different techniques using the same beamline, and they are well reproduced by current neutrino interaction models. These results demonstrate the capability of the detector towards the detailed measurements of the neutrino-nucleus interactions around the 1$\,$GeV energy region., Comment: 27 pages, 17 figures

This paper reports the track multiplicity and kinematics of muons, charged pions, and protons from charged-current inclusive $\bar{\nu}_{\mu}$ and $\nu_{\mu}$ interactions on a water target, measured using a nuclear emulsion detector in the NINJA experiment. A 3-kg water target was exposed to the T2K antineutrino-enhanced beam with a mean energy of 1.3 GeV. Owing to the high-granularity of the nuclear emulsion, protons with momenta down to 200 MeV/$c$ from the neutrino-water interactions were detected. We find good agreement between the observed data and model predictions for all kinematic distributions other than the number of charged pions. These results demonstrate the capability of measurements with nuclear emulsion to improve neutrino interaction models., Comment: 22 pages, 20 figures

An event topology with two secondary vertices compatible with the decay of short-lived particles was found in the analysis of neutrino interactions in the Opera target. The observed topology is compatible with tau neutrino charged current (CC) interactions with charm production and neutrino neutral current (NC) interactions with $c\overline{c}$ pair production. However, other processes can mimic this topology. A dedicated analysis was implemented to identify the underlying process. A Monte Carlo simulation was developed and complementary procedures were introduced in the kinematic reconstruction. A multivariate analysis technique was used to achieve an optimal separation of signal from background. Most likely, this event is a $\nu_{\tau}$ CC interaction with charm production, the tau and charm particle decaying into 1 prong and 2 prongs, respectively. The significance of this observation is evaluated., Comment: 10 pages, 7 figures; Editors: Marco Roda, Chiara Sirignano and Gabriele Sirri

We present the analysis of the binary-lens microlensing event OGLE-2013-BLG-0911. The best-fit solutions indicate the binary mass ratio of q~0.03 which differs from that reported in Shvartzvald+2016. The event suffers from the well-known close/wide degeneracy, resulting in two groups of solutions for the projected separation normalized by the Einstein radius of s~0.15 or s~7. The finite source and the parallax observations allow us to measure the lens physical parameters. The lens system is an M-dwarf orbited by a massive Jupiter companion at very close (M_{host}=0.30^{+0.08}_{-0.06} M_{Sun}, M_{comp}=10.1^{+2.9}_{-2.2} M_{Jup}, a_{exp}=0.40^{+0.05}_{-0.04} au) or wide (M_{host}=0.28^{+0.10}_{-0.08} M_{Sun}, M_{comp}=9.9^{+3.8}_{-3.5}M_{Jup}, a_{exp}=18.0^{+3.2}_{-3.2} au) separation. Although the mass ratio is slightly above the planet-brown dwarf (BD) mass-ratio boundary of q=0.03 which is generally used, the median physical mass of the companion is slightly below the planet-BD mass boundary of 13M_{Jup}. It is likely that the formation mechanisms for BDs and planets are different and the objects near the boundaries could have been formed by either mechanism. It is important to probe the distribution of such companions with masses of ~13M_{Jup} in order to statistically constrain the formation theories for both BDs and massive planets. In particular, the microlensing method is able to probe the distribution around low-mass M-dwarfs and even BDs which is challenging for other exoplanet detection methods., Comment: 22 pages, 10 figures, Accepted for publication in The Astronomical Journal

The OPERA experiment has conclusively observed the appearance of tau neutrinos in the muon neutrino CNGS beam. Exploiting the OPERA detector capabilities, it was possible to isolate high purity samples of $\nu_{e}$, $\nu_{\mu}$ and $\nu_{\tau}$ charged current weak neutrino interactions, as well as neutral current weak interactions. In this Letter, the full dataset is used for the first time to test the three-flavor neutrino oscillation model and to derive constraints on the existence of a light sterile neutrino within the framework of the $3+1$ neutrino model. For the first time, tau and electron neutrino appearance channels are jointly used to test the sterile neutrino hypothesis. A significant fraction of the sterile neutrino parameter space allowed by LSND and MiniBooNE experiments is excluded at 90% C.L. In particular, the best-fit values obtained by MiniBooNE combining neutrino and antineutrino data are excluded at 3.3 $\sigma$ significance., Comment: 7 pages, 4 figures; Editors: Budimir Kli\v{c}ek and Matteo Tenti

OPERA is a long-baseline experiment designed to search for $\nu_{\mu}\to\nu_{\tau}$ oscillations in appearance mode. It was based at the INFN Gran Sasso laboratory (LNGS) and took data from 2008 to 2012 with the CNGS neutrino beam from CERN. After the discovery of $\nu_\tau$ appearance in 2015, with $5.1\sigma$ significance, the criteria to select $\nu_\tau$ candidates have been extended and a multivariate approach has been used for events identification. In this way the statistical uncertainty in the measurement of the oscillation parameters and of $\nu_\tau$ properties has been improved. Results are reported., Comment: 7 pages, 1 figure, conference: The 15th International Workshop on Tau Lepton Physics (tau2018)AB

The OPERA experiment discovered muon neutrino into tau neutrino oscillations in appearance mode, detecting tau leptons by means of nuclear emulsion films. The apparatus was also endowed with electronic detectors with tracking capability, such as scintillator strips and resistive plate chambers. Because of its location, in the underground Gran Sasso laboratory, under 3800 m.w.e., the OPERA detector has also been used as an observatory for TeV muons produced by cosmic rays in the atmosphere. In this paper the measurement of the single muon flux modulation and of its correlation with the seasonal variation of the atmospheric temperature are reported.

We report an orbital characterization of GJ1108Aab that is a low-mass binary system in pre-main-sequence phase. Via the combination of astrometry using adaptive optics and radial velocity measurements, an eccentric orbital solution of $e$=0.63 is obtained, which might be induced by the Kozai-Lidov mechanism with a widely separated GJ1108B system. Combined with several observed properties, we confirm the system is indeed young. Columba is the most probable moving group, to which the GJ1108A system belongs, although its membership to the group has not been established. If the age of Columba is assumed for GJ1108A, the dynamical masses of both GJ1108Aa and GJ1108Ab ($M_{\rm dynamical,GJ1108Aa}=0.72\pm0.04 M_{\odot}$ and $M_{\rm dynamical,GJ1108Ab}=0.30\pm0.03 M_{\odot}$) are more massive than what an evolutionary model predicts based on the age and luminosities. We consider the discrepancy in mass comparison can attribute to an age uncertainty; the system is likely older than stars in Columba, and effects that are not implemented in classical models such as accretion history and magnetic activity are not preferred to explain the mass discrepancy. We also discuss the performance of the evolutionary model by compiling similar low-mass objects in evolutionary state based on the literature. Consequently, it is suggested that the current model on average reproduces the mass of resolved low-mass binaries without any significant offsets., Comment: Accepted in ApJ

The OPERA experiment was designed to study $\nu_\mu\to\nu_\tau$ oscillations in appearance mode in the CNGS neutrino beam. In this letter we report the final analysis of the full data sample collected between 2008 and 2012, corresponding to $17.97\cdot 10^{19}$ protons on target. Selection criteria looser than in previous analyses have produced ten $\nu_\tau$ candidate events, thus reducing the statistical uncertainty in the measurement of the oscillation parameters and of $\nu_\tau$ properties. A multivariate approach for event identification has been applied to the candidate events and the discovery of $\nu_\tau$ appearance is confirmed with an improved significance level of 6.1 $\sigma$. $\Delta m^2_{23}$ has been measured, in appearance mode, with an accuracy of 20%. The measurement of $\nu_\tau$ CC cross-section, for the first time with a negligible contamination from $\bar{\nu}_\tau$, and the first direct observation of the $\nu_\tau$ lepton number are also reported., Comment: 7 pages, 3 figures

We report the discovery and the analysis of the short timescale binary-lens microlensing event, MOA-2015-BLG-337. The lens system could be a planetary system with a very low mass host, around the brown dwarf/planetary mass boundary, or a brown dwarf binary. We found two competing models that explain the observed light curves with companion/host mass ratios of q~0.01 and ~0.17, respectively. From the measurement of finite source effects in the best-fit planetary model, we find a relatively small angular Einstein radius of theta_E ~ 0.03 mas which favors a low mass lens. We conduct a Bayesian analysis to obtain the probability distribution of the lens properties. The results for the planetary models strongly depend on the minimum mass, M_min, in the assumed mass function. In summary, there are two solutions of the lens system: (1) a brown dwarf/planetary mass boundary object orbited by a super-Neptune (the planetary model with M_min=0.001 M_sun) and (2) a brown dwarf binary (the binary model). If the planetary models is correct, this system can be one of a new class of planetary system, having a low host mass and also a planetary mass ratio (q <0.03) between the host and its companion. The discovery of the event is important for the study of planetary formation in very low mass objects. In addition, it is important to consider all viable solutions in these kinds of ambiguous events in order for the future comprehensive statistical analyses of planetary/binary microlensing events., Comment: 17 pages, 7 figures, Accepted for publication in AJ

The OPERA experiment has discovered the tau neutrino appearance in the CNGS muon neutrino beam, in agreement with the 3 neutrino flavour oscillation hypothesis. The OPERA neutrino interaction target, made of Emulsion Cloud Chamber, was particularly efficient in the reconstruction of electromagnetic showers. Moreover, thanks to the very high granularity of the emulsion films, showers induced by electrons can be distinguished from those induced by $\pi^0$s, thus allowing the detection of charged current interactions of electron neutrinos. In this paper the results of the search for electron neutrino events using the full dataset are reported. An improved method for the electron neutrino energy estimation is exploited. Data are compatible with the 3 neutrino flavour mixing model expectations and are used to set limits on the oscillation parameters of the 3+1 neutrino mixing model, in which an additional mass eigenstate $m_{4}$ is introduced. At high $\Delta m^{2}_{41}$ $( \gtrsim 0.1~\textrm{eV}^{2})$, an upper limit on $\sin^2 2\theta_{\mu e}$ is set to 0.021 at 90% C.L. and $\Delta m^2_{41} \gtrsim 4 \times 10^{-3}~\textrm{eV}^{2}$ is excluded for maximal mixing in appearance mode., Comment: Editors: M. Tenti and S. Vasina

We report the discovery and the analysis of the planetary microlensing event, OGLE-2013-BLG-1761. There are some degenerate solutions in this event because the planetary anomaly is only sparsely sampled. But the detailed light curve analysis ruled out all stellar binary models and shows that the lens to be a planetary system. There is the so-called close/wide degeneracy in the solutions with the planet/host mass ratio of $q \sim (7.5 \pm 1.5) \times 10^{-3}$ and $q \sim (9.3 \pm 2.9) \times 10^{-3}$ with the projected separation in Einstein radius units of $s = 0.95$ (close) and $s = 1.19$ (wide), respectively. The microlens parallax effect is not detected but the finite source effect is detected. Our Bayesian analysis indicates that the lens system is located at $D_{\rm L}=6.9_{-1.2}^{+1.0} \ {\rm kpc}$ away from us and the host star is an M/K-dwarf with the mass of $M_{\rm L}=0.33_{-0.18}^{+0.32} \ M_{\odot}$ orbited by a super-Jupiter mass planet with the mass of $m_{\rm P}=2.8_{-1.5}^{+2.5} \ M_{\rm Jup}$ at the projected separation of $a_{\perp}=1.8_{-0.5}^{+0.5} \ {\rm AU}$. The preference of the large lens distance in the Bayesian analysis is due to the relatively large observed source star radius. The distance and other physical parameters can be constrained by the future high resolution imaging by ground large telescopes or HST. If the estimated lens distance is correct, this planet provides another sample for testing the claimed deficit of planets in the Galactic bulge., Comment: 22pages, 7figures, 2tables, submitted to AAS Journals

We describe the first ever implementation of an emulsion multi-stage shifter in an accelerator neutrino experiment. The system was installed in the neutrino monitor building in J-PARC as a part of a test experiment T60 and stable operation was maintained for a total of 126.6 days. By applying time information to emulsion films, various results were obtained. Time resolutions of 5.3 to 14.7 s were evaluated in an operation spanning 46.9 days (time resolved numbers of 3.8--1.4$\times10^{5}$). By using timing and spatial information, a reconstruction of coincident events that consisted of high multiplicity events and vertex events, including neutrino events was performed. Emulsion events were matched to events observed by INGRID, one of near detectors of the T2K experiment, with high reliability (98.5\%) and hybrid analysis was established via use of the multi-stage shifter. The results demonstrate that the multi-stage shifter is feasible for use in neutrino experiments.

Precise neutrino--nucleus interaction measurements in the sub-multi GeV region are important to reduce the systematic uncertainty in future neutrino oscillation experiments. Furthermore, the excess of ${\nu_e}$ interactions, as a possible interpretation of the existence of a sterile neutrino has been observed in such an energy region. The nuclear emulsion technique can measure all the final state particles with low energy threshold for a variety of targets (Fe, C, H${_2}$O, and so on). Its sub-$\mu$m position resolution allows measurements of the ${\nu_e}$ cross-section with good electron/gamma separation capability. We started a new experiment at J-PARC to study sub-multi GeV neutrino interactions by introducing the nuclear emulsion technique. The J-PARC T60 experiment has been implemented as a first step of such a project. Systematic neutrino event analysis with full scanning data in the nuclear emulsion detector was performed for the first time. The first neutrino event detection and its analysis is described in this paper., Comment: 26 pages, 29 figures, 4 table, prepared for submission to PTEP

We present high signal-to-noise ratio, precise $YJH$ photometry and $Y$ band (\gpiwave~$\mu$m) spectroscopy of HD 1160 B, a young substellar companion discovered from the Gemini NICI Planet Finding Campaign, using the Subaru Coronagraphic Extreme Adaptive Optics instrument and the Gemini Planet Imager. HD 1160 B has typical mid-M dwarf-like infrared colors and a spectral type of M5.5$^{+1.0}_{-0.5}$, where the blue edge of our $Y$ band spectrum rules out earlier spectral types. Atmospheric modeling suggests HD 1160 B having an effective temperature of 3000--3100 $K$, a surface gravity of log $g$ = 4--4.5, a radius of~\bestfitradius~$R_{\rm J}$, and a luminosity of log $L$/$L_{\odot} = -2.76 \pm 0.05$. Neither the primary's Hertzspring-Russell diagram position nor atmospheric modeling of HD 1160 B show evidence for a sub-solar metallicity. The interpretation of the HD 1160 B depends on which stellar system components are used to estimate an age. Considering HD 1160 A, B and C jointly, we derive an age of 80--125 Myr, implying that HD 1160 B straddles the hydrogen-burning limit (70--90 $M_{\rm J}$). If we consider HD 1160 A alone, younger ages (20--125 Myr) and a brown dwarf-like mass (35--90 $M_{\rm J}$) are possible. Interferometric measurements of the primary, a precise GAIA parallax, and moderate resolution spectroscopy can better constrain the system's age and how HD 1160 B fits within the context of (sub)stellar evolution., Comment: 50 pages, 6 Tables, 13 Figures; Accepted for publication in The Astrophysical Journal

We report the direct imaging detection of a low-mass companion to a young, moderately active star V450 And, that was previously identified with the radial velocity method. The companion was found in high-contrast images obtained with the Subaru Telescope equipped with the HiCIAO camera and AO188 adaptive optics system. From the public ELODIE and SOPHIE archives we extracted available high-resolution spectra and radial velocity (RV) measurements, along with RVs from the Lick planet search program. We combined our multi-epoch astrometry with these archival, partially unpublished RVs, and found that the companion is a low-mass star, not a brown dwarf, as previously suggested. We found the best-fitting dynamical masses to be $m_1=1.141_{-0.091}^{+0.037}$ and $m_2=0.279^{+0.023}_{-0.020}$ M$_\odot$. We also performed spectral analysis of the SOPHIE spectra with the iSpec code. The Hipparcos time-series photometry shows a periodicity of $P=5.743$ d, which is also seen in SOPHIE spectra as an RV modulation of the star A. We interpret it as being caused by spots on the stellar surface, and the star to be rotating with the given period. From the rotation and level of activity, we found that the system is $380^{+220}_{-100}$ Myr old, consistent with an isochrone analysis ($220^{+2120}_{-90}$ Myr). This work may serve as a test case for future studies of low-mass stars, brown dwarfs and exoplanets by combination of RV and direct imaging data., Comment: 15 pages, 9 figures, 7 tables, to appear in ApJ

The short pulse equation was introduced by Schaefer--Wayne (2004) for modeling the propagation of ultrashort optical pulses. While it can describe a wide range of solutions, its ultrashort pulse solutions with a few cycles, which the conventional nonlinear Schroedinger equation does not possess, have drawn much attention. In such a region, existing numerical methods turn out to require very fine numerical mesh, and accordingly are computationally expensive. In this paper, we establish a new efficient numerical method by combining the idea of the hodograph transformation and the structure-preserving numerical methods. The resulting scheme is a self-adaptive moving mesh scheme that can successfully capture not only the ultrashort pulses but also exotic solutions such as loop soliton solutions.

Dust trapping accelerates the coagulation of dust particles, and thus it represents an initial step toward the formation of planetesimals. We report $H$-band (1.6 um) linear polarimetric observations and 0.87 mm interferometric continuum observations toward a transitional disk around LkHa 330. As results, a pair of spiral arms were detected in the $H$-band emission and an asymmetric (potentially arm-like) structure was detected in the 0.87 mm continuum emission. We discuss the origin of the spiral arm and the asymmetric structure, and suggest that a massive unseen planet is the most plausible explanation. The possibility of dust trapping and grain growth causing the asymmetric structure was also investigated through the opacity index (beta) by plotting the observed SED slope between 0.87 mm from our SMA observation and 1.3 mm from literature. The results imply that grains are indistinguishable from ISM-like dust in the east side ($beta = 2.0 pm 0.5$), but much smaller in the west side $beta = 0.7^{+0.5}_{-0.4}$, indicating differential dust size distribution between the two sides of the disk. Combining the results of near-infrared and submillimeter observations, we conjecture that the spiral arms exist at the upper surface and an asymmetric structure resides in the disk interior. Future observations at centimeter wavelengths and differential polarization imaging in other bands (Y to K) with extreme AO imagers are required to understand how large dust grains form and to further explore the dust distribution in the disk., Comment: 10 pages, 4 figures, AJ accepted for publication

We report the discovery of a microlensing planet OGLE-2012-BLG-0950Lb with the planet/host mass ratio of $q \sim 2 \times 10^{-4}$. A long term distortion detected in both MOA and OGLE light curve can be explained by the microlens parallax due to the Earth's orbital motion around the Sun. Although the finite source effect is not detected, we obtain the lens flux by the high resolution Keck AO observation. Combining the microlens parallax and the lens flux reveal the nature of the lens: a planet with mass of $M_{\rm p} = 35^{+17}_{-9} M_{\oplus}$ is orbiting around a M-dwarf with mass of $M_{\rm host} = 0.56^{+0.12}_{-0.16} M_{\odot}$ with a planet-host projected separation of $r_{\perp} =2.7^{+0.6}_{-0.7}$ AU located at $D_{\rm L} = 3.0^{+0.8}_{-1.1}$ kpc from us. This is the first mass measurement from only microlens parallax and the lens flux without the finite source effect. In the coming space observation-era with $Spitzer$, $K2$, $Euclid$, and $WFIRST$, we expect many such events for which we will not be able to measure any finite source effect. This work demonstrates an ability of mass measurements in such events., Comment: 40 pages, 10 figures, 4 Tables. Accepted to AJ

We present high-contrast angular differential imaging (ADI) observations of the debris disk around HD 32297 in H-band, as well as the first polarimetric images for this system in polarized differential imaging (PDI) mode with Subaru/HICIAO. In ADI, we detect the nearly edge-on disk at >5sigma levels from ~0.45 arcsec to ~1.7 arcsec (50-192 AU) from the star and recover the spine deviation from the midplane already found in previous works. We also find for the first time imaging and surface brightness (SB) indications for the presence of a gapped structure on both sides of the disk at distances of ~0.75 arcsec (NE side) and ~0.65 arcsec (SW side). Global forward-modelling work delivers a best-fit model disk and well-fitting parameter intervals that essentially match previous results, with high-forward scattering grains and a ring located at 110 AU. However, this single ring model cannot account for the gapped structure seen in our SB profiles. We create simple double ring models and achieve a satisfactory fit with two rings located at 60 and 95 AU, respectively, low-forward scattering grains and very sharp inner slopes. In polarized light we retrieve the disk extending from ~0.25-1.6 arcsec, although the central region is quite noisy and high S/N are only found in the range ~0.75-1.2 arcsec. The disk is polarized in the azimuthal direction, as expected, and the departure from the midplane is also clearly observed. Evidence for a gapped scenario is not found in the PDI data. We obtain a linear polarization degree of the grains that increases from ~10% at 0.55 arcsec to ~25% at 1.6 arcsec. The maximum is found at scattering angles of ~90degrees, either from the main components of the disk or from dust grains blown out to larger radii., Comment: Accepted for publication in Astronomy and Astrophysics

We report the discovery of a planet by the microlensing method, OGLE-2012-BLG-0724Lb. Although the duration of the planetary signal for this event was one of the shortest seen for a planetary event, the anomaly was well covered thanks to high cadence observations taken by the survey groups OGLE and MOA. By analyzing the light curve, this planetary system is found to have a mass ratio $q=(1.58\pm0.15)\times10^{-3}$. By conducting a Bayesian analysis, we estimate that the host star is an M-dwarf star with a mass of $M_{\rm L}=0.29_{-0.16}^{+0.33} \ M_{\odot}$ located at $D_{\rm L}=6.7_{-1.2}^{+1.1} \ {\rm kpc}$ away from the Earth and the companion's mass is $m_{\rm P}=0.47_{-0.26}^{+0.54} \ M_{\rm Jup}$. The projected planet-host separation is $a_{\perp}=1.6_{-0.3}^{+0.4} \ {\rm AU}$. Because the lens-source relative proper motion is relatively high, future high resolution images would detect the lens host star and determine the lens properties uniquely. This system is likely a Saturn-mass exoplanet around an M-dwarf and such systems are commonly detected by gravitational microlensing. This adds an another example of a possible pileup of sub-Jupiters $(0.2 < m_{\rm P}/M_{\rm Jup} < 1)$ in contrast to a lack of Jupiters ($\sim 1 - 2 \ M_{\rm Jup}$) around M-dwarfs, supporting the prediction by core accretion models that Jupiter-mass or more massive planets are unlikely to form around M-dwarfs., Comment: 22 pages, 7 figures, 2 tables, accepted for publication in ApJ

The OPERA experiment was designed to search for $\nu_{\mu} \rightarrow \nu_{\tau}$ oscillations in appearance mode, i.e. by detecting the $\tau$-leptons produced in charged current $\nu_{\tau}$ interactions. The experiment took data from 2008 to 2012 in the CERN Neutrinos to Gran Sasso beam. The observation of $\nu_{\mu} \rightarrow \nu_{\tau}$ appearance, achieved with four candidate events in a sub-sample of the data, was previously reported. In this paper, a fifth $\nu_{\tau}$ candidate event, found in an enlarged data sample, is described. Together with a further reduction of the expected background, the candidate events detected so far allow assessing the discovery of $\nu_{\mu}\rightarrow \nu_{\tau}$ oscillations in appearance mode with a significance larger than 5 $\sigma$., Comment: 7 pages, 3 figures

The Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) instrument is a multipurpose high-contrast imaging platform designed for the discovery and detailed characterization of exoplanetary systems and serves as a testbed for high-contrast imaging technologies for ELTs. It is a multi-band instrument which makes use of light from 600 to 2500nm allowing for coronagraphic direct exoplanet imaging of the inner 3 lambda/D from the stellar host. Wavefront sensing and control are key to the operation of SCExAO. A partial correction of low-order modes is provided by Subaru's facility adaptive optics system with the final correction, including high-order modes, implemented downstream by a combination of a visible pyramid wavefront sensor and a 2000-element deformable mirror. The well corrected NIR (y-K bands) wavefronts can then be injected into any of the available coronagraphs, including but not limited to the phase induced amplitude apodization and the vector vortex coronagraphs, both of which offer an inner working angle as low as 1 lambda/D. Non-common path, low-order aberrations are sensed with a coronagraphic low-order wavefront sensor in the infrared (IR). Low noise, high frame rate, NIR detectors allow for active speckle nulling and coherent differential imaging, while the HAWAII 2RG detector in the HiCIAO imager and/or the CHARIS integral field spectrograph (from mid 2016) can take deeper exposures and/or perform angular, spectral and polarimetric differential imaging. Science in the visible is provided by two interferometric modules: VAMPIRES and FIRST, which enable sub-diffraction limited imaging in the visible region with polarimetric and spectroscopic capabilities respectively. We describe the instrument in detail and present preliminary results both on-sky and in the laboratory., Comment: Accepted for publication, 20 pages, 10 figures

We present a new Subaru/HiCIAO high-contrast H-band polarized intensity (PI) image of a nearby transitional disk associated with TW Hydrae. The scattered light from the disk was detected from 0.2" to 1.5" (11 - 81 AU) and the PI image shows a clear axisymmetric depression in polarized intensity at ~ 0.4" (~ 20 AU) from the central star, similar to the ~ 80 AU gap previously reported from HST images. Azimuthal polarized intensity profile also shows the disk beyond 0.2" is almost axisymmetric. We discuss two possible scenarios explaining the origin of the polarized intensity depression: 1) a gap structure may exist at ~ 20 AU from the central star because of shallow slope seen in the polarized intensity profile, and 2) grain growth may be occurring in the inner region of the disk. Multi-band observations at NIR and millimeter/sub-millimeter wavelengths play a complementary role in investigating dust opacity and may help reveal the origin of the gap more precisely., Comment: 8 pages, 2 figures, ApJL accepted for publication

The OPERA experiment, exposed to the CERN to Gran Sasso $\nu_\mu$ beam, collected data from 2008 to 2012. Four oscillated $\nu_\tau$ Charged Current interaction candidates have been detected in appearance mode, which are consistent with $\nu_\mu \to \nu_\tau$ oscillations at the atmospheric $\Delta m^2$ within the "standard" three-neutrino framework. In this paper, the OPERA $\nu_\tau$ appearance results are used to derive limits on the mixing parameters of a massive sterile neutrino., Comment: 11 pages, 4 figures; reference to Planck result updated in the Introduction. Submitted to JHEP

We previously developed an automatic track scanning system which enables the detection of large-angle nuclear fragments in the nuclear emulsion films of the OPERA experiment. As a next step, we have investigated this system's track recognition capability for large-angle minimum ionizing particles $(1.0 \leq |tan \theta| \leq 3.5)$. This paper shows that, for such tracks, the system has a detection efficiency of 95$\%$ or higher and reports the achieved angular accuracy of the automatically recognized tracks. This technology is of general purpose and will likely contribute not only to various analyses in the OPERA experiment, but also to future experiments, e.g. on low-energy neutrino and hadron interactions, or to future research on cosmic rays using nuclear emulsions carried by balloons., Comment: 11 pages, 10 figures, accepted by JINST

The formation scenario of a gapped disk, i.e., transitional disk, and its asymmetry is still under debate. Proposed scenarios such as disk-planet interaction, photoevaporation, grain growth, anticyclonic vortex, eccentricity, and their combinations would result in different radial distributions of the gas and the small (sub-$\mu$m size) and large (millimeter size) dust grains as well as asymmetric structures in a disk. Optical/near-infrared (NIR) imaging observations and (sub-)millimeter interferometry can trace small and large dust grains, respectively; therefore multi-wavelength observations could help elucidate the origin of complicated structures of a disk. Here we report SMA observations of the dust continuum at 1.3~mm and $^{12}$CO~$J=2\rightarrow1$ line emission of the pre-transitional protoplanetary disk around the solar-mass star PDS~70. PDS~70, a weak-lined T Tauri star, exhibits a gap in the scattered light from its disk with a radius of $\sim$65~AU at NIR wavelengths. However, we found a larger gap in the disk with a radius of $\sim$80~AU at 1.3~mm. Emission from all three disk components (the gas and the small and large dust grains) in images exhibits a deficit in brightness in the central region of the disk, in particular, the dust-disk in small and large dust grains has asymmetric brightness. The contrast ratio of the flux density in the dust continuum between the peak position to the opposite side of the disk reaches 1.4. We suggest the asymmetries and different gap-radii of the disk around PDS~70 are potentially formed by several (unseen) accreting planets inducing dust filtration., Comment: Accepted by ApJ

The OPERA experiment is searching for nu_mu -> nu_tau oscillations in appearance mode i.e. via the direct detection of tau leptons in nu_tau charged current interactions. The evidence of nu_mu -> nu_tau appearance has been previously reported with three nu_tau candidate events using a sub-sample of data from the 2008-2012 runs. We report here a fourth nu_tau candidate event, with the tau decaying into a hadron, found after adding the 2012 run events without any muon in the final state to the data sample. Given the number of analysed events and the low background, nu_mu -> nu_tau oscillations are established with a significance of 4.2sigma., Comment: Submitted to Progress of Theoretical and Experimental Physics (PTEP)

The OPERA long-baseline neutrino-oscillation experiment has observed the direct appearance of $\nu_\tau$ in the CNGS $\nu_\mu$ beam. Two large muon magnetic spectrometers are used to identify muons produced in the $\tau$ leptonic decay and in $\nu_\mu^{CC}$ interactions by measuring their charge and momentum. Besides the kinematic analysis of the $\tau$ decays, background resulting from the decay of charmed particles produced in $\nu_\mu^{CC}$ interactions is reduced by efficiently identifying the muon track. A new method for the charge sign determination has been applied, via a weighted angular matching of the straight track-segments reconstructed in the different parts of the dipole magnets. Results obtained for Monte Carlo and real data are presented. Comparison with a method where no matching is used shows a significant reduction of up to 40\% of the fraction of wrongly determined charges., Comment: 10 pages. Improvements in the text

We conduct a statistical analysis of a combined sample of direct imaging data, totalling nearly 250 stars. The stars cover a wide range of ages and spectral types, and include five detections ($\kappa$ And b, two $\sim$60 M$_{\rm J}$ brown dwarf companions in the Pleiades, PZ Tel B, and CD$-$35 2722B). For some analyses we add a currently unpublished set of SEEDS observations, including the detections GJ 504b and GJ 758B. We conduct a uniform, Bayesian analysis of all stellar ages using both membership in a kinematic moving group and activity/rotation age indicators. We then present a new statistical method for computing the likelihood of a substellar distribution function. By performing most of the integrals analytically, we achieve an enormous speedup over brute-force Monte Carlo. We use this method to place upper limits on the maximum semimajor axis of the distribution function derived from radial-velocity planets, finding model-dependent values of $\sim$30--100 AU. Finally, we model the entire substellar sample, from massive brown dwarfs to a theoretically motivated cutoff at $\sim$5 M$_{\rm Jup}$, with a single power law distribution. We find that $p(M, a) \propto M^{-0.65\pm0.60} a^{-0.85\pm0.39}$ (1$\sigma$ errors) provides an adequate fit to our data, with 1.0--3.1\% (68\% confidence) of stars hosting 5--70 $M_{\rm Jup}$ companions between 10 and 100 AU. This suggests that many of the directly imaged exoplanets known, including most (if not all) of the low-mass companions in our sample, formed by fragmentation in a cloud or disk, and represent the low-mass tail of the brown dwarfs., Comment: ApJ in press, replaced with accepted version. Methodology clarified, conclusions mostly unchanged

The OPERA experiment, designed to perform the first observation of $\nu_\mu \rightarrow \nu_\tau$ oscillations in appearance mode through the detection of the $\tau$ leptons produced in $\nu_\tau$ charged current interactions, has collected data from 2008 to 2012. In the present paper, the procedure developed to detect $\tau$ particle decays, occurring over distances of the order of 1 mm from the neutrino interaction point, is described in detail. The results of its application to the search for charmed hadrons are then presented as a validation of the methods for $\nu_\tau$ appearance detection.

The OPERA detector, designed to search for $\nu_{\mu} \to \nu_{\tau}$ oscillations in the CNGS beam, is located in the underground Gran Sasso laboratory, a privileged location to study TeV-scale cosmic rays. For the analysis here presented, the detector was used to measure the atmospheric muon charge ratio in the TeV region. OPERA collected charge-separated cosmic ray data between 2008 and 2012. More than 3 million atmospheric muon events were detected and reconstructed, among which about 110000 multiple muon bundles. The charge ratio $R_{\mu} \equiv N_{\mu^+}/N_{\mu^-}$ was measured separately for single and for multiple muon events. The analysis exploited the inversion of the magnet polarity which was performed on purpose during the 2012 Run. The combination of the two data sets with opposite magnet polarities allowed minimizing systematic uncertainties and reaching an accurate determination of the muon charge ratio. Data were fitted to obtain relevant parameters on the composition of primary cosmic rays and the associated kaon production in the forward fragmentation region. In the surface energy range 1-20 TeV investigated by OPERA, $R_{\mu}$ is well described by a parametric model including only pion and kaon contributions to the muon flux, showing no significant contribution of the prompt component. The energy independence supports the validity of Feynman scaling in the fragmentation region up to $200$ TeV/nucleon primary energy., Comment: 8 pages, 2 figures, 4 tables

The OPERA experiment is designed to search for $\nu_{\mu} \rightarrow \nu_{\tau}$ oscillations in appearance mode i.e. through the direct observation of the $\tau$ lepton in $\nu_{\tau}$ charged current interactions. The experiment has taken data for five years, since 2008, with the CERN Neutrino to Gran Sasso beam. Previously, two $\nu_{\tau}$ candidates with a $\tau$ decaying into hadrons were observed in a sub-sample of data of the 2008-2011 runs. Here we report the observation of a third $\nu_\tau$ candidate in the $\tau^-\to\mu^-$ decay channel coming from the analysis of a sub-sample of the 2012 run. Taking into account the estimated background, the absence of $\nu_{\mu} \rightarrow \nu_{\tau}$ oscillations is excluded at the 3.4 $\sigma$ level., Comment: 9 pages, 5 figures, 1 tables

Most exoplanets detected by direct imaging so far have been characterized by relatively hot (> ~1000 K) and cloudy atmospheres. A surprising feature in some of their atmospheres has been a distinct lack of methane, possibly implying non-equilibrium chemistry. Recently, we reported the discovery of a planetary companion to the Sun-like star GJ 504 using Subaru/HiCIAO within the SEEDS survey. The planet is substantially colder (<600 K) than previously imaged planets, and has indications of fewer clouds, which implies that it represents a new class of planetary atmospheres with expected similarities to late T-type brown dwarfs in the same temperature range. If so, one might also expect the presence of significant methane absorption, which is characteristic of such objects. Here, we report the detection of deep methane absorption in the atmosphere of GJ 504 b, using the Spectral Differential Imaging mode of HiCIAO to distinguish the absorption feature around 1.6 um. We also report updated JHK photometry based on new Ks-band data and a re-analysis of the existing data. The results support the notion that GJ 504 b has atmospheric properties distinct from other imaged exoplanets, and will become a useful reference object for future planets in the same temperature range., Comment: 9 pages, 4 figures, accepted for publication in ApJ Letters

We previously reported the direct detection of a low mass companion at a projected separation of 55+-2 AU around the B9 type star {\kappa} Andromedae. The properties of the system (mass ratio, separation) make it a benchmark for the understanding of the formation and evolution of gas giant planets and brown dwarfs on wide-orbits. We present new angular differential imaging (ADI) images of the Kappa Andromedae system at 2.146 (Ks), 3.776 (L'), 4.052 (NB 4.05) and 4.78 {\mu}m (M') obtained with Keck/NIRC2 and LBTI/LMIRCam, as well as more accurate near-infrared photometry of the star with the MIMIR instrument. We derive a more accurate J = 15.86 +- 0.21, H = 14.95 +- 0.13, Ks = 14.32 +- 0.09 mag for {\kappa} And b. We redetect the companion in all our high contrast observations. We confirm previous contrasts obtained at Ks and L' band. We derive NB 4.05 = 13.0 +- 0.2 and M' = 13.3 +- 0.3 mag and estimate Log10(L/Lsun) = -3.76 +- 0.06. We build the 1-5 microns spectral energy distribution of the companion and compare it to seven PHOENIX-based atmospheric models in order to derive Teff = 1900+100-200 K. Models do not set constrains on the surface gravity. ``Hot-start" evolutionary models predict masses of 14+25-2 MJup based on the luminosity and temperature estimates, and considering a conservative age range for the system (30+120-10 Myr). ``warm-start" evolutionary tracks constrain the mass to M >= 11 MJup. Therefore, the mass of {\kappa} Andromedae b mostly falls in the brown-dwarf regime, due to remaining uncertainties in age and mass-luminosity models. According to the formation models, disk instability in a primordial disk could account for the position and a wide range of plausible masses of {\kappa} And b., Comment: 20 pages, 16 figures, accepted for publication in Astronomy and Astrophysics on August 6, 2013

The OPERA neutrino experiment is designed to perform the first observation of neutrino oscillations in direct appearance mode in the $\nu_\mu \to \nu_\tau$ channel, via the detection of the $\tau$-leptons created in charged current $\nu_\tau$ interactions. The detector, located in the underground Gran Sasso Laboratory, consists of an emulsion/lead target with an average mass of about 1.2 kt, complemented by electronic detectors. It is exposed to the CERN Neutrinos to Gran Sasso beam, with a baseline of 730 km and a mean energy of 17 GeV. The observation of the first $\nu_\tau$ candidate event and the analysis of the 2008-2009 neutrino sample have been reported in previous publications. This work describes substantial improvements in the analysis and in the evaluation of the detection efficiencies and backgrounds using new simulation tools. The analysis is extended to a sub-sample of 2010 and 2011 data, resulting from an electronic detector-based pre-selection, in which an additional $\nu_\tau$ candidate has been observed. The significance of the two events in terms of a $\nu_\mu \to \nu_\tau$ oscillation signal is of 2.40 $\sigma$., Comment: 33 pages, 12 figures, 5 tables. Prepared for submission to JHEP

Several exoplanets have recently been imaged at wide separations of >10 AU from their parent stars. These span a limited range of ages (<50 Myr) and atmospheric properties, with temperatures of 800--1800 K and very red colors (J - H > 0.5 mag), implying thick cloud covers. Furthermore, substantial model uncertainties exist at these young ages due to the unknown initial conditions at formation, which can lead to an order of magnitude of uncertainty in the modeled planet mass. Here, we report the direct imaging discovery of a Jovian exoplanet around the Sun-like star GJ 504, detected as part of the SEEDS survey. The system is older than all other known directly-imaged planets; as a result, its estimated mass remains in the planetary regime independent of uncertainties related to choices of initial conditions in the exoplanet modeling. Using the most common exoplanet cooling model, and given the system age of 160 [+350, -60] Myr, GJ 504 b has an estimated mass of 4 [+4.5, -1.0] Jupiter masses, among the lowest of directly imaged planets. Its projected separation of 43.5 AU exceeds the typical outer boundary of ~30 AU predicted for the core accretion mechanism. GJ 504 b is also significantly cooler (510 [+30, -20] K) and has a bluer color (J-H = -0.23 mag) than previously imaged exoplanets, suggesting a largely cloud-free atmosphere accessible to spectroscopic characterization. Thus, it has the potential of providing novel insights into the origins of giant planets, as well as their atmospheric properties., Comment: 20 pages, 12 figures, Accepted for publication in ApJ. Minor updates from the version 1

We present results from the first three years of observations of moving group targets in the SEEDS high-contrast imaging survey of exoplanets and disks using the Subaru telescope. We achieve typical contrasts of ~10^5 at 1" and ~10^6 beyond 2" around 63 proposed members of nearby kinematic moving groups. We review each of the kinematic associations to which our targets belong, concluding that five, \beta Pictoris (~20 Myr), AB Doradus (~100 Myr), Columba (~30 Myr), Tucana-Horogium (~30 Myr), TW Hydrae (~10 Myr), are sufficiently well-defined to constrain the ages of individual targets. Somewhat less than half of our targets are high-probability members of one of these moving groups. For all of our targets, we combine proposed moving group membership with other age indicators where available, including Ca II HK emission, X-ray activity, and rotation period, to produce a posterior probability distribution of age. SEEDS observations discovered a substellar companion to one of our targets, \kappa And, a late B star. We do not detect any other substellar companions, but do find seven new close binary systems, of which one still needs to be confirmed. A detailed analysis of the statistics of this sample, and of the companion mass constraints given our age probability distributions and exoplanet cooling models, will be presented in a forthcoming paper., Comment: 26 pages, 4 figures, 6 tables. Replaced with published ApJ version

A first result of the search for \numu $\rightarrow$ \nue oscillations in the OPERA experiment, located at the Gran Sasso Underground Laboratory, is presented. The experiment looked for the appearance of \nue in the CNGS neutrino beam using the data collected in 2008 and 2009. Data are compatible with the non-oscillation hypothesis in the three-flavour mixing model. A further analysis of the same data constrains the non-standard oscillation parameters $\theta_{new}$ and $\Delta m^2_{new}$ suggested by the LSND and MiniBooNE experiments. For large $\Delta m^{2}_{new}$ values ($>$0.1 eV$^{2}$), the OPERA 90% C.L. upper limit on sin$^{2}(2\theta_{new})$ based on a Bayesian statistical method reaches the value $7.2 \times 10^{-3}$.