Source author record

J. Kim

J. Kim appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

160works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

160 published item(s)

preprint2026arXiv

CHEX-MATE: Relationship between X-ray and millimetre inferences of galaxy cluster temperature profiles

Thermodynamic profiles from X-ray and millimetre observations of galaxy clusters are often compared under the simplifying assumptions of smooth, spherically symmetric intracluster medium. These approximations lead to expected discrepancies in the inferred profiles, which can provide insights about the cluster structure or cosmology. Motivated by this, we present a joint XMM-\textit{Newton} and \textit{Planck} analysis of 116 CHEX-MATE clusters to measure $η_T = T_X/T_{SZ,X}$, the ratio between spectroscopic X-ray temperatures and a temperature proxy derived from Sunyaev-Zel'dovich (SZ) pressures and X-ray densities. We considered relativistic corrections to the thermal SZ signal and implemented X-ray absorption by Galactic molecular hydrogen. The $η_T$ distribution has a mean of $1.01 \pm 0.03$, with average changes of $8.1\%$ and $2.7\%$ when relativistic corrections and molecular hydrogen absorption are not included, respectively. The $η_T$ distribution is positively skewed, with the scatter mostly affected by cluster morphology: relaxed clusters are closer to unity and less scattered than mixed and disturbed systems. We find little or no correlation with redshift, mass, or temperature.

preprint2026arXiv

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1σ$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments.

preprint2026arXiv

Euclid preparation. Galaxy 2-point correlation function modelling in redshift space

The Euclid satellite will measure spectroscopic redshifts for tens of millions of emission-line galaxies. In the context of Stage-IV surveys, the 3-dimensional clustering of galaxies plays a key role in providing cosmological constraints. In this paper, we conduct a model comparison for the multipole moments of the galaxy 2-point correlation function (2PCF) in redshift space. We test state-of-the-art models, in particular the effective field theory of large-scale structure (EFT), one based on the velocity difference generating function (VDG$_{\infty}$), and different variants of Lagrangian perturbation theory (LPT) models, such as convolutional LPT (CLPT) and its effective-field-theory extension (CLEFT). We analyse the first three even multipoles of the 2PCF in the Flagship 1 simulation, which consists of four snapshots at $z\in\{0.9,1.2,1.5,1.8\}$. We study both template-fitting and full-shape approaches and find that with the template-fitting approach, only the VDG$_{\infty}$ model is able to reach a minimum fitting scale of $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ at $z=0.9$ without biasing the recovered parameters. Indeed, the EFT model becomes inaccurate already at $s_{\rm min}=30\,h^{-1}\,{\rm Mpc}$. Conversely, in the full-shape analysis, the CLEFT and VDG$_{\infty}$ models perform similarly well, but only the CLEFT model can reach $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ while the VDG$_{\infty}$ model is unbiased down to $s_{\rm min}=25\,h^{-1}\,{\rm Mpc}$ at the lowest redshift. Overall, in order to achieve the accuracy required by Euclid, non-perturbative modelling such as in the VDG$_{\infty}$ or CLEFT models should be considered. At $z=1.8$, the CLPT model is sufficient to describe the data with high figure of merit. This comparison selects baseline models that perform best in ideal conditions and sets the stage for an optimal analysis of Euclid data in configuration space.

preprint2026arXiv

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

We model intrinsic alignments (IA) in Euclid's Flagship simulation to investigate its impact on Euclid's weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples.

preprint2026arXiv

Euclid: Galaxy SED reconstruction in the PHZ processing function: impact on the PSF and the role of medium-band filters

Weak lensing surveys require accurate correction for the point spread function (PSF) when measuring galaxy shapes. For a diffraction-limited PSF, as arises in space-based missions, this correction depends on each galaxy SED. In the Euclid mission, galaxy SED reconstruction, a tasks of the photometric-redshift processing function (PHZ PF), relies on broad- and medium-band ancillary photometry. The limited wavelength sampling of the Euclid VIS passband and signal-to-noise ratio may affect the reconstruction accuracy and translate into biases in the weak lensing measurements. In this study, we present the methodology, which is employed in the Euclid PHZ PF, for reconstructing galaxy SEDs at 55 wavelengths, sampling the VIS passband every 10 nm, and we assess whether it fulfils the accuracy requirements imposed on the Euclid PSF model. We employ both physics- and data-driven methods, focusing on a new approach of template-based flux correction and Gaussian processes, and we introduce an SED metric whose bias propagates into PSF quadrupole moment errors. Our findings demonstrate that Gaussian processes and template fitting meet the requirements only in specific, but complementary, redshift intervals. We therefore propose a hybrid approach, which leverages both methods. This solution proves to be effective in meeting the Euclid accuracy requirements for most of the redshift range of the survey. Finally, we investigate the impact on the SED reconstruction of a new set of 16 evenly-spaced medium-band filters for the Subaru telescope, providing quasi-spectroscopic coverage of the VIS passband. This study shows promising results, ensuring accurate SED reconstruction and meeting the mission PSF requirements. This work thus provides not only the methodological foundation of galaxy SED reconstruction in the Euclid PHZ PF, but also a roadmap for future improvements using a new medium-band survey.

preprint2026arXiv

Reconstruction of atmospheric neutrinos in DUNE's horizontal-drift far-detector module

This paper reports on the capabilities in reconstructing and identifying atmospheric neutrino interactions in one of the Deep Underground Neutrino Experiment's (DUNE) far detector modules, a liquid argon time projection chamber (LArTPC) with horizontal drift (FD-HD) of ionization electrons. The reconstruction is based upon the workflow developed for DUNE's long-baseline oscillation analysis, with some necessary machine-learning models' retraining and the addition of features relevant only to atmospheric neutrinos such as the neutrino direction reconstruction. Where relevant, the impact of the detection of the charged particles of the hadronic system is emphasized, and comparisons are carried out between the case when lepton-only information is considered in the reconstruction (as is the case for many neutrino oscillation experiments), versus when all particles identified in the LArTPC were included. Three neutrino direction reconstruction methods have been developed and studied for the atmospheric analyses: using lepton-only information, using all reconstructed particles, and using only correlations from reconstructed hits. The results indicate that incorporating more than just lepton information significantly improves the resolution of both neutrino direction and energy reconstruction. The angle reconstruction algorithms developed in this work result in no strong dependence on particle direction for reconstruction efficiencies or neutrino flavor identification. This comprehensive review of the reconstruction of atmospheric neutrinos in DUNE's FD-HD LArTPC is the first step towards developing a first neutrino oscillation sensitivity analysis, which will ready DUNE for its first measurements.

preprint2025arXiv

CHEX-MATE: towards a consistent universal pressure profile and cluster mass reconstruction

In a self-similar paradigm of structure formation, the thermal pressure of the hot intra-cluster gas follows a universal distribution once the profile of each cluster is normalised based on the proper mass and redshift dependencies. The reconstruction of such a universal pressure profile requires an individual estimate of the mass of each cluster. In this context, we present a method to jointly fit, for the first time, the universal pressure profile and individual cluster $M_{500}$ masses over a sample of galaxy clusters, properly accounting for correlations between the profile shape and amplitude, and masses scaling the individual profiles. We demonstrate the power of the method and show that a consistent exploitation of the universal pressure profile and cluster mass estimates when modelling the thermal pressure in clusters is necessary to avoid biases. In particular, the method, informed by a cluster mass scale, outputs individual cluster masses with same accuracy and better precision than input masses. Using data from the «Cluster HEritage project with XMM-Newton: Mass Assembly and Thermodynamics at the Endpoint of structure formation», we investigate a sample of $\sim 25$ galaxy clusters spanning mass and redshift ranges of $2 \lesssim M_{500}/10^{14} \; \mathrm{M}_{\odot} \lesssim 14$ and $0.07 < z < 0.6$.

preprint2023arXiv

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more.

preprint2022arXiv

All-sky search for gravitational wave emission from scalar boson clouds around spinning black holes in LIGO O3 data

This paper describes the first all-sky search for long-duration, quasi-monochromatic gravitational-wave signals emitted by ultralight scalar boson clouds around spinning black holes using data from the third observing run of Advanced LIGO. We analyze the frequency range from 20~Hz to 610~Hz, over a small frequency derivative range around zero, and use multiple frequency resolutions to be robust towards possible signal frequency wanderings. Outliers from this search are followed up using two different methods, one more suitable for nearly monochromatic signals, and the other more robust towards frequency fluctuations. We do not find any evidence for such signals and set upper limits on the signal strain amplitude, the most stringent being $\approx10^{-25}$ at around 130~Hz. We interpret these upper limits as both an "exclusion region" in the boson mass/black hole mass plane and the maximum detectable distance for a given boson mass, based on an assumption of the age of the black hole/boson cloud system.

preprint2022arXiv

Constraints from LIGO O3 data on gravitational-wave emission due to r-modes in the glitching pulsar PSR J0537-6910

We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537-6910 using data from the LIGO-Virgo Collaboration observing run O3. PSR J0537-6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role in the spin evolution of this pulsar. Theoretical models confirm this possibility and predict emission at a level that can be probed by ground-based detectors. In order to explore this scenario, we search for r-mode emission in the epochs between glitches by using a contemporaneous timing ephemeris obtained from NICER data. We do not detect any signals in the theoretically expected band of 86-97 Hz, and report upper limits on the amplitude of the gravitational waves. Our results improve on previous amplitude upper limits from r-modes in J0537-6910 by a factor of up to 3 and place stringent constraints on theoretical models for r-mode driven spin-down in PSR J0537-6910, especially for higher frequencies at which our results reach below the spin-down limit defined by energy conservation.

preprint2022arXiv

First all-flavor search for transient neutrino emission using 3-years of IceCube DeepCore data

Since the discovery of a flux of high-energy astrophysical neutrinos, searches for their origins have focused primarily at TeV-PeV energies. Compared to sub-TeV searches, high-energy searches benefit from an increase in the neutrino cross section, improved angular resolution on the neutrino direction, and a reduced background from atmospheric neutrinos and muons. However, the focus on high energy does not preclude the existence of sub-TeV neutrino emission where IceCube retains sensitivity. Here we present the first all-flavor search from IceCube for transient emission of low-energy neutrinos, focusing on the energy region of 5.6-100 GeV using three years of data obtained with the IceCube-DeepCore detector. We find no evidence of transient neutrino emission in the data, thus leading to a constraint on the volumetric rate of astrophysical transient sources in the range of $\sim 705-2301\, \text{Gpc}^{-3}\, \text{yr}^{-1}$ for sources following a subphotospheric energy spectrum with a mean energy of 100 GeV and a bolometric energy of $10^{52}$ erg.

preprint2022arXiv

First joint observation by the underground gravitational-wave detector, KAGRA, with GEO600

We report the results of the first joint observation of the KAGRA detector with GEO600. KAGRA is a cryogenic and underground gravitational-wave detector consisting of a laser interferometer with three-kilometer arms, and located in Kamioka, Gifu, Japan. GEO600 is a British--German laser interferometer with 600 m arms, and located near Hannover, Germany. GEO600 and KAGRA performed a joint observing run from April 7 to 20, 2020. We present the results of the joint analysis of the GEO--KAGRA data for transient gravitational-wave signals, including the coalescence of neutron-star binaries and generic unmodeled transients. We also perform dedicated searches for binary coalescence signals and generic transients associated with gamma-ray burst events observed during the joint run. No gravitational-wave events were identified. We evaluate the minimum detectable amplitude for various types of transient signals and the spacetime volume for which the network is sensitive to binary neutron-star coalescences. We also place lower limits on the distances to the gamma-ray bursts analysed based on the non-detection of an associated gravitational-wave signal for several signal models, including binary coalescences. These analyses demonstrate the feasibility and utility of KAGRA as a member of the global gravitational-wave detector network.

preprint2022arXiv

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets.

preprint2022arXiv

Quality Control of Mass-Produced GEM Detectors for the CMS GE1/1 Muon Upgrade

The series of upgrades to the Large Hadron Collider, culminating in the High Luminosity Large Hadron Collider, will enable a significant expansion of the physics program of the CMS experiment. However, the accelerator upgrades will also make the experimental conditions more challenging, with implications for detector operations, triggering, and data analysis. The luminosity of the proton-proton collisions is expected to exceed $2-3\times10^{34}$~cm$^{-2}$s$^{-1}$ for Run 3 (starting in 2022), and it will be at least $5\times10^{34}$~cm$^{-2}$s$^{-1}$ when the High Luminosity Large Hadron Collider is completed for Run 4. These conditions will affect muon triggering, identification, and measurement, which are critical capabilities of the experiment. To address these challenges, additional muon detectors are being installed in the CMS endcaps, based on Gas Electron Multiplier technology. For this purpose, 161 large triple-Gas Electron Multiplier detectors have been constructed and tested. Installation of these devices began in 2019 with the GE1/1 station and will be followed by two additional stations, GE2/1 and ME0, to be installed in 2023 and 2026, respectively. The assembly and quality control of the GE1/1 detectors were distributed across several production sites around the world. We motivate and discuss the quality control procedures that were developed to standardize the performance of the detectors, and we present the final results of the production. Out of 161 detectors produced, 156 detectors passed all tests, and 144 detectors are now installed in the CMS experiment. The various visual inspections, gas tightness tests, intrinsic noise rate characterizations, and effective gas gain and response uniformity tests allowed the project to achieve this high success rate.

preprint2022arXiv

Search for anisotropic gravitational-wave backgrounds using data from Advanced LIGO and Advanced Virgo's first three observing runs

We report results from searches for anisotropic stochastic gravitational-wave backgrounds using data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. For the first time, we include Virgo data in our analysis and run our search with a new efficient pipeline called {\tt PyStoch} on data folded over one sidereal day. We use gravitational-wave radiometry (broadband and narrow band) to produce sky maps of stochastic gravitational-wave backgrounds and to search for gravitational waves from point sources. A spherical harmonic decomposition method is employed to look for gravitational-wave emission from spatially-extended sources. Neither technique found evidence of gravitational-wave signals. Hence we derive 95\% confidence-level upper limit sky maps on the gravitational-wave energy flux from broadband point sources, ranging from $F_{α, Θ} < {\rm (0.013 - 7.6)} \times 10^{-8} {\rm erg \, cm^{-2} \, s^{-1} \, Hz^{-1}},$ and on the (normalized) gravitational-wave energy density spectrum from extended sources, ranging from $Ω_{α, Θ} < {\rm (0.57 - 9.3)} \times 10^{-9} \, {\rm sr^{-1}}$, depending on direction ($Θ$) and spectral index ($α$). These limits improve upon previous limits by factors of $2.9 - 3.5$. We also set 95\% confidence level upper limits on the frequency-dependent strain amplitudes of quasimonochromatic gravitational waves coming from three interesting targets, Scorpius X-1, SN 1987A and the Galactic Center, with best upper limits range from $h_0 < {\rm (1.7-2.1)} \times 10^{-25},$ a factor of $\geq 2.0$ improvement compared to previous stochastic radiometer searches.

preprint2022arXiv

Search for continuous gravitational wave emission from the Milky Way center in O3 LIGO--Virgo data

We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO--Virgo run in the detector frequency band $[10,2000]\rm~Hz$ have been used. No significant detection was found and 95$\%$ confidence level upper limits on the signal strain amplitude were computed, over the full search band, with the deepest limit of about $7.6\times 10^{-26}$ at $\simeq 142\rm~Hz$. These results are significantly more constraining than those reported in previous searches. We use these limits to put constraints on the fiducial neutron star ellipticity and r-mode amplitude. These limits can be also translated into constraints in the black hole mass -- boson mass plane for a hypothetical population of boson clouds around spinning black holes located in the GC.

preprint2022arXiv

Search for continuous gravitational waves from 20 accreting millisecond X-ray pulsars in O3 LIGO data

Results are presented of searches for continuous gravitational waves from 20 accreting millisecond X-ray pulsars with accurately measured spin frequencies and orbital parameters, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. The search algorithm uses a hidden Markov model, where the transition probabilities allow the frequency to wander according to an unbiased random walk, while the $\mathcal{J}$-statistic maximum-likelihood matched filter tracks the binary orbital phase. Three narrow sub-bands are searched for each target, centered on harmonics of the measured spin frequency. The search yields 16 candidates, consistent with a false alarm probability of 30% per sub-band and target searched. These candidates, along with one candidate from an additional target-of-opportunity search done for SAX J1808.4$-$3658, which was in outburst during one month of the observing run, cannot be confidently associated with a known noise source. Additional follow-up does not provide convincing evidence that any are a true astrophysical signal. When all candidates are assumed non-astrophysical, upper limits are set on the maximum wave strain detectable at 95% confidence, $h_0^{95\%}$. The strictest constraint is $h_0^{95\%} = 4.7\times 10^{-26}$ from IGR J17062$-$6143. Constraints on the detectable wave strain from each target lead to constraints on neutron star ellipticity and $r$-mode amplitude, the strictest of which are $ε^{95\%} = 3.1\times 10^{-7}$ and $α^{95\%} = 1.8\times 10^{-5}$ respectively. This analysis is the most comprehensive and sensitive search of continuous gravitational waves from accreting millisecond X-ray pulsars to date.

preprint2022arXiv

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.

preprint2022arXiv

Searching for low-mass dark matter via Migdal effect in COSINE-100

We report on the search for weakly interacting massive particle (WIMP) dark matter candidates in the galactic halo that interact with sodium and iodine nuclei in the COSINE-100 experiment and produce energetic electrons that accompany recoil nuclei via the the Migdal effect. The WIMP mass sensitivity of previous COSINE-100 searches that relied on the detection of ionization signals produced by target nuclei recoiling from elastic WIMP-nucleus scattering was restricted to WIMP masses above $\sim$5 GeV/$c^2$ by the detectors' 1 keVee energy-electron-equivalent threshold. The search reported here looks for recoil signals enhanced by the Migdal electrons that are ejected during the scattering process. This is particularly effective for the detection of low-mass WIMP scattering from the crystals' sodium nuclei in which a relatively larger fraction of the WIMP's energy is transferred to the nucleus recoil energy and the excitation of its orbital electrons. In this analysis, the low-mass WIMP search window of the COSINE-100 experiment is extended to WIMP mass down to 200 MeV/$c^2$. The low-mass WIMP sensitivity will be further improved by lowering the analysis threshold based on a multivariable analysis technique. We consider the influence of these improvements and recent developments in detector performance to re-evaluate sensitivities for the future COSINE-200 experiment. With a 0.2 keVee analysis threshold and high light-yield NaI(Tl) detectors (22 photoelectrons/keVee), the COSINE-200 experiment can explore low-mass WIMPs down to 20 MeV/$c^2$ and probe previously unexplored regions of parameter space.

preprint2022arXiv

Vectorial probing of electric and magnetic transitions in variable optical environments and vice-versa

We use europium doped single crystalline NaYF$_4$ nanorods for probing the electric and magnetic contributions to the local density of optical states (LDOS). Reciprocically, we determine intrinsic properties of the emitters (oscillator strength, quantum yield) by comparing their measured and simulated optical responses in front of a mirror. We first experimentally determine the specifications of the nanoprobe (orientation and oscillator strength of the electric and magnetic dipoles moments) and show significant orientation sensitivity of the branching ratios associated with electric and magnetic transitions. In a second part, we measure the modification of the LDOS in front of a gold mirror in a Drexhage's experiment. We discuss the role of the electric and magnetic LDOS on the basis of numerical simulations, taking into account the orientation of the dipolar emitters. We demonstrate that they behave like degenerated dipoles sensitive to polarized partial LDOS.

preprint2021arXiv

All-sky search for long-duration gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

After the detection of gravitational waves from compact binary coalescences, the search for transient gravitational-wave signals with less well-defined waveforms for which matched filtering is not well-suited is one of the frontiers for gravitational-wave astronomy. Broadly classified into "short" $ \lesssim 1~$\,s and "long" $ \gtrsim 1~$\,s duration signals, these signals are expected from a variety of astrophysical processes, including non-axisymmetric deformations in magnetars or eccentric binary black hole coalescences. In this work, we present a search for long-duration gravitational-wave transients from Advanced LIGO and Advanced Virgo's third observing run from April 2019 to March 2020. For this search, we use minimal assumptions for the sky location, event time, waveform morphology, and duration of the source. The search covers the range of $2~\text{--}~ 500$~s in duration and a frequency band of $24 - 2048$ Hz. We find no significant triggers within this parameter space; we report sensitivity limits on the signal strength of gravitational waves characterized by the root-sum-square amplitude $h_{\mathrm{rss}}$ as a function of waveform morphology. These $h_{\mathrm{rss}}$ limits improve upon the results from the second observing run by an average factor of 1.8.

preprint2021arXiv

All-sky search for short gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

This paper presents the results of a search for generic short-duration gravitational-wave transients in data from the third observing run of Advanced LIGO and Advanced Virgo. Transients with durations of milliseconds to a few seconds in the 24--4096 Hz frequency band are targeted by the search, with no assumptions made regarding the incoming signal direction, polarization or morphology. Gravitational waves from compact binary coalescences that have been identified by other targeted analyses are detected, but no statistically significant evidence for other gravitational wave bursts is found. Sensitivities to a variety of signals are presented. These include updated upper limits on the source rate-density as a function of the characteristic frequency of the signal, which are roughly an order of magnitude better than previous upper limits. This search is sensitive to sources radiating as little as $\sim$10$^{-10} M_{\odot} c^2$ in gravitational waves at $\sim$70 Hz from a distance of 10~kpc, with 50\% detection efficiency at a false alarm rate of one per century. The sensitivity of this search to two plausible astrophysical sources is estimated: neutron star f-modes, which may be excited by pulsar glitches, as well as selected core-collapse supernova models.

preprint2021arXiv

All-sky, all-frequency directional search for persistent gravitational-waves from Advanced LIGO's and Advanced Virgo's first three observing runs

We present the first results from an all-sky all-frequency (ASAF) search for an anisotropic stochastic gravitational-wave background using the data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. Upper limit maps on broadband anisotropies of a persistent stochastic background were published for all observing runs of the LIGO-Virgo detectors. However, a broadband analysis is likely to miss narrowband signals as the signal-to-noise ratio of a narrowband signal can be significantly reduced when combined with detector output from other frequencies. Data folding and the computationally efficient analysis pipeline, {\tt PyStoch}, enable us to perform the radiometer map-making at every frequency bin. We perform the search at 3072 {\tt{HEALPix}} equal area pixels uniformly tiling the sky and in every frequency bin of width $1/32$~Hz in the range $20-1726$~Hz, except for bins that are likely to contain instrumental artefacts and hence are notched. We do not find any statistically significant evidence for the existence of narrowband gravitational-wave signals in the analyzed frequency bins. Therefore, we place $95\%$ confidence upper limits on the gravitational-wave strain for each pixel-frequency pair, the limits are in the range $(0.030 - 9.6) \times10^{-24}$. In addition, we outline a method to identify candidate pixel-frequency pairs that could be followed up by a more sensitive (and potentially computationally expensive) search, e.g., a matched-filtering-based analysis, to look for fainter nearly monochromatic coherent signals. The ASAF analysis is inherently independent of models describing any spectral or spatial distribution of power. We demonstrate that the ASAF results can be appropriately combined over frequencies and sky directions to successfully recover the broadband directional and isotropic results.

preprint2021arXiv

Constraint on the Matter-Antimatter Symmetry-Violating Phase in Neutrino Oscillations

The current laws of physics do not explain the observed imbalance of matter and antimatter in the universe. Sakharov proposed that an explanation would require the violation of CP symmetry between matter and antimatter. The only CP violation observed so far is in the weak interactions of quarks, and it is too small to explain the matter-antimatter imbalance of the universe. It has been shown that CP violation in the lepton sector could generate the matter-antimatter disparity through the process called leptogenesis. The quantum mixing of neutrinos, the neutral leptons in the Standard Model, provides a potential source of CP violation through a complex phase dCP, which may have consequences for theoretical models of leptogenesis. This CP violation can be measured in muon neutrino to electron neutrino oscillations and the corresponding antineutrino oscillations, which are experimentally accessible with accelerator-produced beams as established by the T2K experiment. Until now, the value of dCP has not been significantly constrained by neutrino oscillation experiments. Here the T2K collaboration reports a measurement that favors large enhancement of the neutrino oscillation probability, excluding values of dCP which result in a large enhancement of the observed anti-neutrino oscillation probability at three standard deviations (3 sigma). The 3 sigma confidence level interval for dCP, which is cyclic and repeats every 2pi, is [-3.41,-0.03] for the so-called normal mass ordering, and [-2.54,-0.32] for the inverted mass ordering. Our results show an indication of CP violation in the lepton sector. Herein we establish methods for sensitive searches for matter-antimatter asymmetry in neutrino oscillations using accelerator-produced neutrino beams. Future measurements with larger data samples will determine whether the leptonic CP violation is larger than the quark sector CP violation.

preprint2021arXiv

Detection of a particle shower at the Glashow resonance with IceCube

The Glashow resonance describes the resonant formation of a $W^-$ boson during the interaction of a high-energy electron antineutrino with an electron, peaking at an antineutrino energy of 6.3 petaelectronvolts (PeV) in the rest frame of the electron. Whereas this energy scale is out of reach for currently operating and future planned particle accelerators, natural astrophysical phenomena are expected to produce antineutrinos with energies beyond the PeV scale. Here we report the detection by the IceCube neutrino observatory of a cascade of high-energy particles (a particle shower) consistent with being created at the Glashow resonance. A shower with an energy of $6.05 \pm 0.72$ PeV (determined from Cherenkov radiation in the Antarctic Ice Sheet) was measured. Features consistent with the production of secondary muons in the particle shower indicate the hadronic decay of a resonant $W^-$ boson, confirm that the source is astrophysical and provide improved directional localization. The evidence of the Glashow resonance suggests the presence of electron antineutrinos in the astrophysical flux, while also providing further validation of the standard model of particle physics. Its unique signature indicates a method of distinguishing neutrinos from antineutrinos, thus providing a way to identify astronomical accelerators that produce neutrinos via hadronuclear or photohadronic interactions, with or without strong magnetic fields. As such, knowledge of both the flavour (that is, electron, muon or tau neutrinos) and charge (neutrino or antineutrino) will facilitate the advancement of neutrino astronomy.

preprint2021arXiv

Emergent interlayer magnetic order via strain-induced orthorhombic distortion in the 5d Mott insulator Sr2IrO4

We report a La2CuO4-like interlayer antiferromagnetic order in Sr2IrO4 films with large orthorhombic distortion (> 1.5%). The biaxial lattice strain in epitaxial heterostructures of Sr2IrO4/Ca3Ru2O7 lowers the crystal symmetry of Sr2IrO4 from tetragonal (C4) to orthorhombic (C2), guiding the Ir 5d Jeff = 1/2 pseudospin moment parallel to the elongated b-axis via magnetic anisotropy. From resonant X-ray scattering experiments, we observed an antiferromagnetic order in the orthorhombic Sr2IrO4 film whose interlayer stacking pattern is inverted from that of the tetragonal Sr2IrO4 crystal. This interlayer stacking is similar to that of the orthorhombic La2CuO4, implying that the asymmetric interlayer exchange interaction along a and b-directions exceeds the anisotropic interlayer pseudo-dipolar interaction. Our result suggests that strain-induced distortion can provide a delicate knob for tuning the long-range magnetic order in quasi-two-dimensional systems by evoking the competition between the interlayer exchange coupling and the pseudo-dipolar interaction.

preprint2021arXiv

IceCube Data for Neutrino Point-Source Searches Years 2008-2018

IceCube has performed several all-sky searches for point-like neutrino sources using track-like events, including a recent time-integrated analysis using 10 years of IceCube data. This paper accompanies the public data release of these neutrino candidates detected by IceCube between April 6, 2008 and July 8, 2018. The selection includes through-going tracks, primarily due to muon neutrino candidates, that reach the detector from all directions, as well as neutrino track events that start within the instrumented volume. An updated selection and reconstruction for data taken after April 2012 slightly improves the sensitivity of the sample. While more than 80% of the sample overlaps between the old and new versions, differing events can lead to changes relative to the previous 7 year event selection. An a posteriori estimate of the significance of the 2014-2015 TXS flare is reported with an explanation of observed discrepancies with previous results. This public data release, which includes 10 years of data and binned detector response functions for muon neutrino signal events, shows improved sensitivity in generic time-integrated point source analyses and should be preferred over previous releases.

preprint2021arXiv

Modeling the triple-GEM detector response to background particles for the CMS Experiment

An estimate of environmental background hit rate on triple-GEM chambers is performed using Monte Carlo (MC) simulation and compared to data taken by test chambers installed in the CMS experiment (GE1/1) during Run-2 at the Large Hadron Collider (LHC). The hit rate is measured using data collected with proton-proton collisions at 13 TeV and a luminosity of 1.5$\times10^{34}$ cm$^{-2}$ s$^{-1}$. The simulation framework uses a combination of the FLUKA and Geant4 packages to obtain the hit rate. FLUKA provides the radiation environment around the GE1/1 chambers, which is comprised of the particle flux with momentum direction and energy spectra ranging from $10^{-11}$ to $10^{4}$ MeV for neutrons, $10^{-3}$ to $10^{4}$ MeV for $γ$'s, $10^{-2}$ to $10^{4}$ MeV for $e^{\pm}$, and $10^{-1}$ to $10^{4}$ MeV for charged hadrons. Geant4 provides an estimate of detector response (sensitivity) based on an accurate description of detector geometry, material composition and interaction of particles with the various detector layers. The MC simulated hit rate is estimated as a function of the perpendicular distance from the beam line and agrees with data within the assigned uncertainties of 10-14.5%. This simulation framework can be used to obtain a reliable estimate of background rates expected at the High Luminosity LHC.

preprint2021arXiv

Overview of KAGRA: Calibration, detector characterization, physical environmental monitors, and the geophysics interferometer

KAGRA is a newly built gravitational wave observatory, a laser interferometer with a 3 km arm length, located in Kamioka, Gifu, Japan. In this series of articles, we present an overview of the baseline KAGRA, for which we finished installing the designed configuration in 2019. This article describes the method of calibration (CAL) used for reconstructing gravitational wave signals from the detector outputs, as well as the characterization of the detector (DET). We also review the physical environmental monitors (PEM) system and the geophysics interferometer (GIF). Both are used for characterizing and evaluating the data quality of the gravitational wave channel. They play important roles in utilizing the detector output for gravitational wave searches. These characterization investigations will be even more important in the near future, once gravitational wave detection has been achieved, and in using KAGRA in the gravitational wave astronomy era.

preprint2021arXiv

Search for Gravitational Waves Associated with Gamma-Ray Bursts Detected by Fermi and Swift During the LIGO-Virgo Run O3b

We search for gravitational-wave signals associated with gamma-ray bursts detected by the Fermi and Swift satellites during the second half of the third observing run of Advanced LIGO and Advanced Virgo (1 November 2019 15:00 UTC-27 March 2020 17:00 UTC).We conduct two independent searches: a generic gravitational-wave transients search to analyze 86 gamma-ray bursts and an analysis to target binary mergers with at least one neutron star as short gamma-ray burst progenitors for 17 events. We find no significant evidence for gravitational-wave signals associated with any of these gamma-ray bursts. A weighted binomial test of the combined results finds no evidence for sub-threshold gravitational wave signals associated with this GRB ensemble either. We use several source types and signal morphologies during the searches, resulting in lower bounds on the estimated distance to each gamma-ray burst. Finally, we constrain the population of low luminosity short gamma-ray bursts using results from the first to the third observing runs of Advanced LIGO and Advanced Virgo. The resulting population is in accordance with the local binary neutron star merger rate.

preprint2021arXiv

Search for subsolar-mass binaries in the first half of Advanced LIGO and Virgo's third observing run

We report on a search for compact binary coalescences where at least one binary component has a mass between 0.2 $M_\odot$ and 1.0 $M_\odot$ in Advanced LIGO and Advanced Virgo data collected between 1 April 2019 1500 UTC and 1 October 2019 1500 UTC. We extend previous analyses in two main ways: we include data from the Virgo detector and we allow for more unequal mass systems, with mass ratio $q \geq 0.1$. We do not report any gravitational-wave candidates. The most significant trigger has a false alarm rate of 0.14 $\mathrm{yr}^{-1}$. This implies an upper limit on the merger rate of subsolar binaries in the range $[220-24200] \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$, depending on the chirp mass of the binary. We use this upper limit to derive astrophysical constraints on two phenomenological models that could produce subsolar-mass compact objects. One is an isotropic distribution of equal-mass primordial black holes. Using this model, we find that the fraction of dark matter in primordial black holes is $f_\mathrm{PBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 6\%$. The other is a dissipative dark matter model, in which fermionic dark matter can collapse and form black holes. The upper limit on the fraction of dark matter black holes depends on the minimum mass of the black holes that can be formed: the most constraining result is obtained at $M_\mathrm{min}=1 M_\odot$, where $f_\mathrm{DBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 0.003\%$. These are the tightest limits on spinning subsolar-mass binaries to date.

preprint2021arXiv

Searches for neutrinos from cosmic-ray interactions in the Sun using seven years of IceCube data

Cosmic-ray interactions with the solar atmosphere are expected to produce particle showers which in turn produce neutrinos from weak decays of mesons. These solar atmospheric neutrinos (SA$ν$s) have never been observed experimentally. A detection would be an important step in understanding cosmic-ray propagation in the inner solar system and the dynamics of solar magnetic fields. SA$ν$s also represent an irreducible background to solar dark matter searches and a detection would allow precise characterization of this background. Here, we present the first experimental search based on seven years of data collected from May 2010 to May 2017 in the austral winter with the IceCube Neutrino Observatory. An unbinned likelihood analysis is performed for events reconstructed within 5 degrees of the center of the Sun. No evidence for a SA$ν$ flux is observed. After inclusion of systematic uncertainties, we set a 90\% upper limit of $1.02^{+0.20}_{-0.18}\cdot10^{-13}$~$\mathrm{GeV^{-1}cm^{-2}s^{-1}}$ at 1 TeV.

preprint2021arXiv

The Nearby Evolved Stars Survey II: Constructing a volume-limited sample and first results from the James Clerk Maxwell Telescope

The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $μ$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust.

preprint2020arXiv

A search for IceCube events in the direction of ANITA neutrino candidates

During the first three flights of the Antarctic Impulsive Transient Antenna (ANITA) experiment, the collaboration detected several neutrino candidates. Two of these candidate events were consistent with an ultra-high-energy up-going air shower and compatible with a tau neutrino interpretation. A third neutrino candidate event was detected in a search for Askaryan radiation in the Antarctic ice, although it is also consistent with the background expectation. The inferred emergence angle of the first two events is in tension with IceCube and ANITA limits on isotropic cosmogenic neutrino fluxes. Here, we test the hypothesis that these events are astrophysical in origin, possibly caused by a point source in the reconstructed direction. Given that any ultra-high-energy tau neutrino flux traversing the Earth should be accompanied by a secondary flux in the TeV-PeV range, we search for these secondary counterparts in seven years of IceCube data using three complementary approaches. In the absence of any significant detection, we set upper limits on the neutrino flux from potential point sources. We compare these limits to ANITA's sensitivity in the same direction and show that an astrophysical explanation of these anomalous events under standard model assumptions is severely constrained regardless of source spectrum.

preprint2020arXiv

A Search for MeV to TeV Neutrinos from Fast Radio Bursts with IceCube

We present two searches for IceCube neutrino events coincident with 28 fast radio bursts (FRBs) and one repeating FRB. The first improves upon a previous IceCube analysis -- searching for spatial and temporal correlation of events with FRBs at energies greater than roughly 50 GeV -- by increasing the effective area by an order of magnitude. The second is a search for temporal correlation of MeV neutrino events with FRBs. No significant correlation is found in either search, therefore, we set upper limits on the time-integrated neutrino flux emitted by FRBs for a range of emission timescales less than one day. These are the first limits on FRB neutrino emission at the MeV scale, and the limits set at higher energies are an order-of-magnitude improvement over those set by any neutrino telescope.

preprint2020arXiv

An experimental proof that resistance-switching memories are not memristors

It has been suggested that all resistive-switching memory cells are memristors. The latter are hypothetical, ideal devices whose resistance, as originally formulated, depends only on the net charge that traverses them. Recently, an unambiguous test has been proposed [J. Phys. D: Appl. Phys. {\bf 52}, 01LT01 (2019)] to determine whether a given physical system is indeed a memristor or not. Here, we experimentally apply such a test to both in-house fabricated Cu-SiO2 and commercially available electrochemical metallization cells. Our results unambiguously show that electrochemical metallization memory cells are not memristors. Since the particular resistance-switching memories employed in our study share similar features with many other memory cells, our findings refute the claim that all resistance-switching memories are memristors. They also cast doubts on the existence of ideal memristors as actual physical devices that can be fabricated experimentally. Our results then lead us to formulate two memristor impossibility conjectures regarding the impossibility of building a model of physical resistance-switching memories based on the memristor model.

preprint2020arXiv

ANTARES and IceCube Combined Search for Neutrino Point-like and Extended Sources in the Southern Sky

A search for point-like and extended sources of cosmic neutrinos using data collected by the ANTARES and IceCube neutrino telescopes is presented. The data set consists of all the track-like and shower-like events pointing in the direction of the Southern Sky included in the nine-year ANTARES point-source analysis, combined with the through-going track-like events used in the seven-year IceCube point-source search. The advantageous field of view of ANTARES and the large size of IceCube are exploited to improve the sensitivity in the Southern Sky by a factor $\sim$2 compared to both individual analyses. In this work, the Southern Sky is scanned for possible excesses of spatial clustering, and the positions of preselected candidate sources are investigated. In addition, special focus is given to the region around the Galactic Centre, whereby a dedicated search at the location of SgrA* is performed, and to the location of the supernova remnant RXJ 1713.7-3946. No significant evidence for cosmic neutrino sources is found and upper limits on the flux from the various searches are presented.

preprint2020arXiv

Characteristics of the diffuse astrophysical electron and tau neutrino flux with six years of IceCube high energy cascade data

We report on the first measurement of the astrophysical neutrino flux using particle showers (cascades) in IceCube data from 2010 -- 2015. Assuming standard oscillations, the astrophysical neutrinos in this dedicated cascade sample are dominated ($\sim 90 \%$) by electron and tau flavors. The flux, observed in the sensitive energy range from $16\,\mathrm{TeV}$ to $2.6\,\mathrm{PeV}$, is consistent with a single power-law model as expected from Fermi-type acceleration of high energy particles at astrophysical sources. We find the flux spectral index to be $γ=2.53\pm0.07$ and a flux normalization for each neutrino flavor of $ϕ_{astro} = 1.66^{+0.25}_{-0.27}$ at $E_{0} = 100\, \mathrm{TeV}$, in agreement with IceCube's complementary muon neutrino results and with all-neutrino flavor fit results. In the measured energy range we reject spectral indices $γ\leq2.28$ at $\ge3σ$ significance level. Due to high neutrino energy resolution and low atmospheric neutrino backgrounds, this analysis provides the most detailed characterization of the neutrino flux at energies below $\sim100\,{\rm{TeV}}$ compared to previous IceCube results. Results from fits assuming more complex neutrino flux models suggest a flux softening at high energies and a flux hardening at low energies (p-value $\ge 0.06$). The sizable and smooth flux measured below $\sim 100\,{\rm{TeV}}$ remains a puzzle. In order to not violate the isotropic diffuse gamma-ray background as measured by the Fermi-LAT, it suggests the existence of astrophysical neutrino sources characterized by dense environments which are opaque to gamma-rays.

preprint2020arXiv

Constraints on Neutrino Emission from Nearby Galaxies Using the 2MASS Redshift Survey and IceCube

The distribution of galaxies within the local universe is characterized by anisotropic features. Observatories searching for the production sites of astrophysical neutrinos can take advantage of these features to establish directional correlations between a neutrino dataset and overdensities in the galaxy distribution in the sky. The results of two correlation searches between a seven-year time-integrated neutrino dataset from the IceCube Neutrino Observatory, and the 2MASS Redshift Survey (2MRS) catalog are presented here. The first analysis searches for neutrinos produced via interactions between diffuse intergalactic Ultra-High Energy Cosmic Rays (UHECRs) and the matter contained within galaxies. The second analysis searches for low-luminosity sources within the local universe, which would produce subthreshold multiplets in the IceCube dataset that directionally correlate with galaxy distribution. No significant correlations were observed in either analyses. Constraints are presented on the flux of neutrinos originating within the local universe through diffuse intergalactic UHECR interactions, as well as on the density of standard candle sources of neutrinos at low luminosities.

preprint2020arXiv

Development of an analysis to probe the neutrino mass ordering with atmospheric neutrinos using three years of IceCube DeepCore data

The Neutrino Mass Ordering (NMO) remains one of the outstanding questions in the field of neutrino physics. One strategy to measure the NMO is to observe matter effects in the oscillation pattern of atmospheric neutrinos above $\sim 1\,\mathrm{GeV}$, as proposed for several next-generation neutrino experiments. Moreover, the existing IceCube DeepCore detector can already explore this type of measurement. We present rthe development and application of two independent analyses to search for the signature of the NMO with three years of DeepCore data. These analyses include a full treatment of systematic uncertainties and a statistically-rigorous method to determine the significance for the NMO from a fit to the data. Both analyses show that the dataset is fully compatible with both mass orderings. For the more sensitive analysis, we observe a preference for Normal Ordering with a $p$-value of $p_\mathrm{IO} = 15.3\%$ and $\mathrm{CL}_\mathrm{s}=53.3\%$ for the Inverted Ordering hypothesis, while the experimental results from both analyses are consistent within their uncertainties. Since the result is independent of the value of $δ_\mathrm{CP}$ and obtained from energies $E_ν\gtrsim 5\,\mathrm{GeV}$, it is complementary to recent results from long-baseline experiments. These analyses set the groundwork for the future of this measurement with more capable detectors, such as the IceCube Upgrade and the proposed PINGU detector.

preprint2020arXiv

Effective Elastic Wave Characteristics of Composite Media

We derive exact expressions for effective elastodynamic properties of two-phase composites in the long-wavelength (quasistatic) regime via homogenized constitutive relations that are local in space. This is accomplished by extending the "strong-contrast" expansion formalism that was previously applied to the static problem. These strong-contrast expansions explicitly incorporate complete microstructural information of the composite via an infinite set of $n$-point correlation functions. Utilizing the rapid-convergence properties of these series expansions (even for extreme contrast ratios), we extract accurate approximations that depend on the microstructure via the spectral density, which is easy to compute or measure for any composite. We also investigate the predictive power of modifications of such approximation formulas postulated elsewhere [J. Kim and S. Torquato, Proc. Nat. Acad. Sci. ${\bf 117}$, 8764 (2020)] to extend their applicability ${\it beyond ~the ~quasistatic ~regime.}$ The accuracy of these nonlocal microstructure-dependent approximations is validated by comparison to full-waveform simulation results for certain models of dispersions. We apply our formulas to a variety of models of nonhyperuniform and hyperuniform disordered composites. We demonstrate that hyperuniform systems are less lossy than their nonhyperuniform counterparts in the quasistatic regime, and stealthy hyperuniform media can be perfectly transparent for a wide range of wavenumber. Finally, we discuss how to utilize our approximations for engineering composites with prescribed elastic wave characteristics.

preprint2020arXiv

IceCube Search for High-Energy Neutrino Emission from TeV Pulsar Wind Nebulae

Pulsar wind nebulae (PWNe) are the main gamma-ray emitters in the Galactic plane. They are diffuse nebulae that emit nonthermal radiation. Pulsar winds, relativistic magnetized outflows from the central star, shocked in the ambient medium produce a multiwavelength emission from the radio through gamma rays. Although the leptonic scenario is able to explain most PWNe emission, a hadronic contribution cannot be excluded. A possible hadronic contribution to the high-energy gamma-ray emission inevitably leads to the production of neutrinos. Using 9.5 yr of all-sky IceCube data, we report results from a stacking analysis to search for neutrino emission from 35 PWNe that are high-energy gamma-ray emitters. In the absence of any significant correlation, we set upper limits on the total neutrino emission from those PWNe and constraints on hadronic spectral components.

preprint2020arXiv

IceCube Search for Neutrinos Coincident with Compact Binary Mergers from LIGO-Virgo's First Gravitational-Wave Transient Catalog

Using the IceCube Neutrino Observatory, we search for high-energy neutrino emission coincident with compact binary mergers observed by the LIGO and Virgo gravitational wave (GW) detectors during their first and second observing runs. We present results from two searches targeting emission coincident with the sky localization of each gravitational wave event within a 1000 second time window centered around the reported merger time. One search uses a model-independent unbinned maximum likelihood analysis, which uses neutrino data from IceCube to search for point-like neutrino sources consistent with the sky localization of GW events. The other uses the Low-Latency Algorithm for Multi-messenger Astrophysics, which incorporates astrophysical priors through a Bayesian framework and includes LIGO-Virgo detector characteristics to determine the association between the GW source and the neutrinos. No significant neutrino coincidence is seen by either search during the first two observing runs of the LIGO-Virgo detectors. We set upper limits on the time-integrated neutrino emission within the 1000 second window for each of the 11 GW events. These limits range from 0.02-0.7 $\mathrm{GeV~cm^{-2}}$. We also set limits on the total isotropic equivalent energy, $E_{\mathrm{iso}}$, emitted in high-energy neutrinos by each GW event. These limits range from 1.7 $\times$ 10$^{51}$ - 1.8 $\times$ 10$^{55}$ erg. We conclude with an outlook for LIGO-Virgo observing run O3, during which both analyses are running in real time.

preprint2020arXiv

Interstrip Capacitances of the Readout Board used in Large Triple-GEM Detectors for the CMS Muon Upgrade

We present analytical calculations, Finite Element Analysis modeling, and physical measurements of the interstrip capacitances for different potential strip geometries and dimensions of the readout boards for the GE2/1 triple-Gas Electron Multiplier detector in the CMS muon system upgrade. The main goal of the study is to find configurations that minimize the interstrip capacitances and consequently maximize the signal-to-noise ratio for the detector. We find agreement at the 1.5--4.8% level between the two methods of calculations and on the average at the 17% level between calculations and measurements. A configuration with halved strip lengths and doubled strip widths results in a measured 27--29% reduction over the original configuration while leaving the total number of strips unchanged. We have now adopted this design modification for all eight module types of the GE2/1 detector and will produce the final detector with this new strip design.

preprint2020arXiv

Measurements of the Absolute Branching Fractions of $B^\pm \to K^\pm X_{c\bar c}$

A study of the two body decays $B^\pm\rightarrow X_{c\bar c}K^\pm$, where X$_{c\bar c}$ refers to one charmonium state, is reported by BaBar collaboration using a data sample of 424 fb$^{-1}$. The absolute determination of branching fractions for these decays are significantly improved compared to previous BaBaR measurements. Evidence is found for the decay $B^+\rightarrow X(3872)K^+$ at the $3σ$ level. The absolute branching fraction ${\cal B}(B^+\rightarrow X(3872)K^+) = (2.1\pm0.6({\rm stat})\pm0.3({\rm syst}))\times 10^{-4}$ is measured for the first time. It follows that ${\cal B}(X(3872)\rightarrow J/ψπ^+π^-)=(4.1\pm1.3)\%$, supporting the hypothesis of a molecular component for this resonance.

preprint2020arXiv

Measurements of the Time-Dependent Cosmic-Ray Sun Shadow with Seven Years of IceCube Data -- Comparison with the Solar Cycle and Magnetic Field Models

Observations of the time-dependent cosmic-ray Sun shadow have been proven as a valuable diagnostic for the assessment of solar magnetic field models. In this paper, seven years of IceCube data are compared to solar activity and solar magnetic field models. A quantitative comparison of solar magnetic field models with IceCube data on the event rate level is performed for the first time. Additionally, a first energy-dependent analysis is presented and compared to recent predictions. We use seven years of IceCube data for the Moon and the Sun and compare them to simulations on data rate level. The simulations are performed for the geometrical shadow hypothesis for the Moon and the Sun and for a cosmic-ray propagation model governed by the solar magnetic field for the case of the Sun. We find that a linearly decreasing relationship between Sun shadow strength and solar activity is preferred over a constant relationship at the 6.4sigma level. We test two commonly used models of the coronal magnetic field, both combined with a Parker spiral, by modeling cosmic-ray propagation in the solar magnetic field. Both models predict a weakening of the shadow in times of high solar activity as it is also visible in the data. We find tensions with the data on the order of $3σ$ for both models, assuming only statistical uncertainties. The magnetic field model CSSS fits the data slightly better than the PFSS model. This is generally consistent with what is found previously by the Tibet AS-gamma Experiment, a deviation of the data from the two models is, however, not significant at this point. Regarding the energy dependence of the Sun shadow, we find indications that the shadowing effect increases with energy during times of high solar activity, in agreement with theoretical predictions.

preprint2020arXiv

Nickel Foam as a Substrate for III-nitride Nanowire Growth

This article presents the use of flexible metal foam substrates for the growth of III-nitride nanowire light emitters to tackle the inherent limitations of thin-film light emitting diodes as well as fabrication and application issues of traditional substrates. A dense packing of gallium nitride nanowires were grown on a nickel foam substrate. The nanowires grew predominantly along the a-plane direction, normal to the local surface of the nickel foam. Strong luminescence was observed from undoped GaN and InGaN quantum well light emitting diode nanowires.

preprint2020arXiv

Opportunities and Challenges in MOCVD of β-Ga2O3 for Power Electronic Devices

Recent breakthroughs in bulk crystal growth of beta-Ga2O3 by the edge-defined film-fed technique has led to the commercialization of large-area beta-Ga2O3 substrates. Standard epitaxy approaches are being utilized to develop various thin-film beta-Ga2O3 based devices including lateral transistors. This article will discuss the challenges for metal organic chemical vapor deposition (MOCVD) of beta-Ga2O3 and the design criteria for use of this material system in power electronic device structures.

preprint2020arXiv

Overview of KAGRA: Detector design and construction history

KAGRA is a newly built gravitational-wave telescope, a laser interferometer comprising arms with a length of 3\,km, located in Kamioka, Gifu, Japan. KAGRA was constructed under the ground and it is operated using cryogenic mirrors that help in reducing the seismic and thermal noise. Both technologies are expected to provide directions for the future of gravitational-wave telescopes. In 2019, KAGRA finished all installations with the designed configuration, which we call the baseline KAGRA. In this occasion, we present an overview of the baseline KAGRA from various viewpoints in a series of of articles. In this article, we introduce the design configurations of KAGRA with its historical background.

preprint2020arXiv

Planck 2018 results. V. CMB power spectra and likelihoods

This paper describes the 2018 Planck CMB likelihoods, following a hybrid approach similar to the 2015 one, with different approximations at low and high multipoles, and implementing several methodological and analysis refinements. With more realistic simulations, and better correction and modelling of systematics, we can now make full use of the High Frequency Instrument polarization data. The low-multipole 100x143 GHz EE cross-spectrum constrains the reionization optical-depth parameter $τ$ to better than 15% (in combination with with the other low- and high-$\ell$ likelihoods). We also update the 2015 baseline low-$\ell$ joint TEB likelihood based on the Low Frequency Instrument data, which provides a weaker $τ$ constraint. At high multipoles, a better model of the temperature-to-polarization leakage and corrections for the effective calibrations of the polarization channels (polarization efficiency or PE) allow us to fully use the polarization spectra, improving the constraints on the $Λ$CDM parameters by 20 to 30% compared to TT-only constraints. Tests on the modelling of the polarization demonstrate good consistency, with some residual modelling uncertainties, the accuracy of the PE modelling being the main limitation. Using our various tests, simulations, and comparison between different high-$\ell$ implementations, we estimate the consistency of the results to be better than the 0.5$σ$ level. Minor curiosities already present before (differences between $\ell$<800 and $\ell$>800 parameters or the preference for more smoothing of the $C_\ell$ peaks) are shown to be driven by the TT power spectrum and are not significantly modified by the inclusion of polarization. Overall, the legacy Planck CMB likelihoods provide a robust tool for constraining the cosmological model and represent a reference for future CMB observations. (Abridged)

preprint2020arXiv

Planck 2018 results. VII. Isotropy and Statistics of the CMB

Analysis of the Planck 2018 data set indicates that the statistical properties of the cosmic microwave background (CMB) temperature anisotropies are in excellent agreement with previous studies using the 2013 and 2015 data releases. In particular, they are consistent with the Gaussian predictions of the $Λ$CDM cosmological model, yet also confirm the presence of several so-called "anomalies" on large angular scales. The novelty of the current study, however, lies in being a first attempt at a comprehensive analysis of the statistics of the polarization signal over all angular scales, using either maps of the Stokes parameters, $Q$ and $U$, or the $E$-mode signal derived from these using a new methodology (which we describe in an appendix). Although remarkable progress has been made in reducing the systematic effects that contaminated the 2015 polarization maps on large angular scales, it is still the case that residual systematics (and our ability to simulate them) can limit some tests of non-Gaussianity and isotropy. However, a detailed set of null tests applied to the maps indicates that these issues do not dominate the analysis on intermediate and large angular scales (i.e., $\ell \lesssim 400$). In this regime, no unambiguous detections of cosmological non-Gaussianity, or of anomalies corresponding to those seen in temperature, are claimed. Notably, the stacking of CMB polarization signals centred on the positions of temperature hot and cold spots exhibits excellent agreement with the $Λ$CDM cosmological model, and also gives a clear indication of how Planck provides state-of-the-art measurements of CMB temperature and polarization on degree scales.

preprint2020arXiv

Planck intermediate results. LV. Reliability and thermal properties of high-frequency sources in the Second Planck Catalogue of Compact Sources

We describe an extension of the most recent version of the Planck Catalogue of Compact Sources (PCCS2), produced using a new multi-band Bayesian Extraction and Estimation Package (BeeP). BeeP assumes that the compact sources present in PCCS2 at 857 GHz have a dust-like spectral energy distribution, which leads to emission at both lower and higher frequencies, and adjusts the parameters of the source and its SED to fit the emission observed in Planck's three highest frequency channels at 353, 545, and 857 GHz, as well as the IRIS map at 3000 GHz. In order to reduce confusion regarding diffuse cirrus emission, BeeP's data model includes a description of the background emission surrounding each source, and it adjusts the confidence in the source parameter extraction based on the statistical properties of the spatial distribution of the background emission. BeeP produces the following three new sets of parameters for each source: (a) fits to a modified blackbody (MBB) thermal emission model of the source; (b) SED-independent source flux densities at each frequency considered; and (c) fits to an MBB model of the background in which the source is embedded. BeeP also calculates, for each source, a reliability parameter, which takes into account confusion due to the surrounding cirrus. We define a high-reliability subset (BeeP/base), containing 26 083 sources (54.1 per cent of the total PCCS2 catalogue), the majority of which have no information on reliability in the PCCS2. The results of the BeeP extension of PCCS2, which are made publicly available via the PLA, will enable the study of the thermal properties of well-defined samples of compact Galactic and extra-galactic dusty sources.

preprint2020arXiv

Planck intermediate results. LVI. Detection of the CMB dipole through modulation of the thermal Sunyaev-Zeldovich effect: Eppur si muove II

The largest temperature anisotropy in the cosmic microwave background (CMB) is the dipole, which has been measured with increasing accuracy for more than three decades, particularly with the Planck satellite. The simplest interpretation of the dipole is that it is due to our motion with respect to the rest frame of the CMB. Since current CMB experiments infer temperature anisotropies from angular intensity variations, the dipole modulates the temperature anisotropies with the same frequency dependence as the thermal Sunyaev-Zeldovich (tSZ) effect. We present the first, and significant, detection of this signal in the tSZ maps and find that it is consistent with direct measurements of the CMB dipole, as expected. The signal contributes power in the tSZ maps, which is modulated in a quadrupolar pattern, and we estimate its contribution to the tSZ bispectrum, noting that it contributes negligible noise to the bispectrum at relevant scales.

preprint2020arXiv

Precision measurement of the ${\cal B}(Υ(3S)\toτ^+τ^-)/{\cal B}(Υ(3S)\toμ^+μ^-)$ ratio

We report on a precision measurement of the ratio ${\cal R}_{τμ}^{Υ(3S)} = {\cal B}(Υ(3S)\toτ^+τ^-)/{\cal B}(Υ(3S)\toμ^+μ^-)$ using data collected with the BaBar detector at the SLAC PEP-II $e^+e^-$ collider. The measurement is based on a 28 fb$^{-1}$ data sample collected at a center-of-mass energy of 10.355 GeV corresponding to a sample of 122 million $Υ(3S)$ mesons. The ratio is measured to be ${\cal R}_{τμ}^{Υ(3S)} = 0.966 \pm 0.008_\mathrm{stat} \pm 0.014_\mathrm{syst}$ and is in agreement with the Standard Model prediction of 0.9948 within 2 standard deviations. The uncertainty in ${\cal R}_{τμ}^{Υ(3S)}$ is almost an order of magnitude smaller than the only previous measurement.

preprint2020arXiv

Search for heavy neutrinos with the T2K near detector ND280

This paper reports on the search for heavy neutrinos with masses in the range $140 < M_N < 493$ MeV/c$^2$ using the off-axis near detector ND280 of the T2K experiment. These particles can be produced from kaon decays in the standard neutrino beam and then subsequently decay in ND280. The decay modes under consideration are $N \to \ell^{\pm}_α π^{\mp}$ and $N \to \ell^+_α \ell^-_β ν(\barν)$ ($α,β=e,μ$). A search for such events has been made using the Time Projection Chambers of ND280, where the background has been reduced to less than two events in the current dataset in all channels. No excess has been observed in the signal region. A combined Bayesian statistical approach has been applied to extract upper limits on the mixing elements of heavy neutrinos to electron-, muon- and tau- flavoured currents ($U_e^2$, $U_μ^2$, $U_τ^2$) as a function of the heavy neutrino mass, e.g. $U_e^2 < 10^{-9}$ at $90\%$ C.L. for a mass of $390$ MeV/c$^2$. These constraints are competitive with previous experiments.

preprint2020arXiv

Search for lepton-flavor violating decays $D^{0}\rightarrow X^{0}e^{\pm}μ^{\mp}$

We present a search for seven lepton-flavor-violating neutral charm decays of the type $D^{0}\rightarrow X^{0} e^{\pm} μ^{\mp}$, where $X^{0}$ represents a $π^{0}$, $K^{0}_{\rm S}$, $\bar{K^{*0}}$, $ρ^{0}$, $ϕ$, $ω$, or $η$ meson. The analysis is based on $468$ fb$^{-1}$ of $e^+e^-$ annihilation data collected at or close to the $Υ(4S)$ resonance with the BaBar detector at the SLAC National Accelerator Laboratory. No significant signals are observed, and we establish 90\% confidence level upper limits on the branching fractions in the range $(5.0 - 22.5)\times 10^{-7}$. The limits are between one and two orders of magnitude more stringent than previous measurements.

preprint2020arXiv

Search for PeV Gamma-Ray Emission from the Southern Hemisphere with 5 Years of Data from the IceCube Observatory

The measurement of diffuse PeV gamma-ray emission from the Galactic plane would provide information about the energy spectrum and propagation of Galactic cosmic rays, and the detection of a point-like source of PeV gamma rays would be strong evidence for a Galactic source capable of accelerating cosmic rays up to at least a few PeV. This paper presents several un-binned maximum likelihood searches for PeV gamma rays in the Southern Hemisphere using 5 years of data from the IceTop air shower surface detector and the in-ice array of the IceCube Observatory. The combination of both detectors takes advantage of the low muon content and deep shower maximum of gamma-ray air showers, and provides excellent sensitivity to gamma rays between $\sim$0.6 PeV and 100 PeV. Our measurements of point-like and diffuse Galactic emission of PeV gamma rays are consistent with background, so we constrain the angle-integrated diffuse gamma-ray flux from the Galactic Plane at 2 PeV to $2.61 \times 10^{-19}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$ at 90% confidence, assuming an E$^{-3}$ spectrum, and we estimate 90% upper limits on point-like emission at 2 PeV between 10$^{-21}$ - 10$^{-20}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$ for an E$^{-2}$ spectrum, depending on declination. Furthermore, we exclude unbroken power-law emission up to 2 PeV for several TeV gamma-ray sources observed by H.E.S.S., and calculate upper limits on the energy cutoffs of these sources at 90% confidence. We also find no PeV gamma rays correlated with neutrinos from IceCube's high-energy starting event sample. These are currently the strongest constraints on PeV gamma-ray emission.

preprint2020arXiv

Search for rare or forbidden decays of the $D^{0}$ meson

We present a search for nine lepton-number-violating and three lepton-flavor-violating neutral charm decays of the type $D^0\rightarrow h^{\prime -} h^{-}\ell^{\prime +} \ell^{+}$ and $D^0\rightarrow h^{\prime -} h^{+}\ell^{\prime\pm} \ell^{\mp}$, where $h$ and $h^{\prime}$ represent a $K$ or $π$ meson and $\ell$ and $\ell^{\prime}$ an electron or muon. The analysis is based on $468$ fb$^{-1}$ of $e^+e^-$ annihilation data collected at or close to the $Y(4S)$ resonance with the BaBar detector at the SLAC National Accelerator Laboratory. No significant signal is observed for any of the twelve modes and we establish 90% confidence level upper limits on the branching fractions in the range $(1.0 - 30.6)\times 10^{-7}$. The limits are between one and three orders of magnitude times more stringent than previous measurements.

preprint2020arXiv

Semimetal to semiconductor transition in $\text{Bi}/\text{TiO}_{2}$ core/shell nanowires

We demonstrate the full thermoelectric and structural characterization of individual bismuth-based (Bi-based) core/shell nanowires. The influence of strain on the temperature dependence of the electrical conductivity, the absolute Seebeck coefficient and the thermal conductivity of bismuth/titanium dioxide ($\text{Bi}/\text{TiO}_{2}$) nanowires with different diameters is investigated and compared to bismuth (Bi) and bismuth/tellurium (Bi/Te) nanowires and bismuth bulk. Scattering at surfaces, crystal defects and interfaces between the core and the shell reduces the electrical conductivity to less than $5\,\%$ and the thermal conductivity to less than $25\,\%$ to $50\,\%$ of the bulk value at room temperature. On behalf of a compressive strain, $\text{Bi}/\text{TiO}_{2}$ core/shell nanowires show a decreasing electrical conductivity with decreasing temperature opposed to that of Bi and Bi/Te nanowires. We find that the compressive strain induced by the $\text{TiO}_{2}$ shell can lead to a band opening of bismuth increasing the absolute Seebeck coefficient by $10\,\%$ to $30\,\%$ compared to bulk at room temperature. In the semiconducting state, the activation energy is determined to $\left|41.3\pm0.2\right|\,\text{meV}$. We show that if the strain exceeds the elastic limit the semimetallic state is recovered due to the lattice relaxation.

preprint2020arXiv

Visualizing Orbital Content of Electronic Bands in Anisotropic 2D Semiconducting ReSe$_{2}$

Many properties of layered materials change as they are thinned from their bulk forms down to single layers, with examples including indirect-to-direct band gap transition in 2H semiconducting transition metal dichalcogenides as well as thickness-dependent changes in the valence band structure in post-transition metal monochalcogenides and black phosphorus. Here, we use angle-resolved photoemission spectroscopy to study the electronic band structure of monolayer ReSe$_{2}$, a semiconductor with a distorted 1T structure and in-plane anisotropy. By changing the polarization of incoming photons, we demonstrate that for ReSe$_{2}$, in contrast to the 2H materials, the out-of-plane transition metal $d_{z^{2}}$ and chalcogen $p_{z}$ orbitals do not contribute significantly to the top of the valence band which explains the reported weak changes in the electronic structure of this compound as a function of layer number. We estimate a band gap of 1.7 eV in pristine ReSe$_{2}$ using scanning tunneling spectroscopy and explore the implications on the gap following surface-doping with potassium. A lower bound of 1.4 eV is estimated for the gap in the fully doped case, suggesting that doping-dependent many-body effects significantly affect the electronic properties of ReSe$_{2}$. Our results, supported by density functional theory calculations, provide insight into the mechanisms behind polarization-dependent optical properties of rhenium dichalcogenides and highlight their place amongst two-dimensional crystals.

preprint2019arXiv

A Search for Neutrino Point-Source Populations in 7 Years of IceCube Data with Neutrino-count Statistics

The presence of a population of point sources in a dataset modifies the underlying neutrino-count statistics from the Poisson distribution. This deviation can be exactly quantified using the non-Poissonian template fitting technique, and in this work we present the first application this approach to the IceCube high-energy neutrino dataset. Using this method, we search in 7 years of IceCube data for point-source populations correlated with the disk of the Milky Way, the Fermi bubbles, the Schlegel, Finkbeiner, and Davis dust map, or with the isotropic extragalactic sky. No evidence for such a population is found in the data using this technique, and in the absence of a signal we establish constraints on population models with source count distribution functions that can be described by a power-law with a single break. The derived limits can be interpreted in the context of many possible source classes. In order to enhance the flexibility of the results, we publish the full posterior from our analysis, which can be used to establish limits on specific population models that would contribute to the observed IceCube neutrino flux.

preprint2019arXiv

Benchmarking an 11-qubit quantum computer

The field of quantum computing has grown from concept to demonstration devices over the past 20 years. Universal quantum computing offers efficiency in approaching problems of scientific and commercial interest, such as factoring large numbers, searching databases, simulating intractable models from quantum physics, and optimizing complex cost functions. Here, we present an 11-qubit fully-connected, programmable quantum computer in a trapped ion system composed of 13 $^{171}$Yb$^{+}$ ions. We demonstrate average single-qubit gate fidelities of 99.5$\%$, average two-qubit-gate fidelities of 97.5$\%$, and state preparation and measurement errors of 0.7$\%$. To illustrate the capabilities of this universal platform and provide a basis for comparison with similarly-sized devices, we compile the Bernstein-Vazirani (BV) and Hidden Shift (HS) algorithms into our native gates and execute them on the hardware with average success rates of 78$\%$ and 35$\%$, respectively. These algorithms serve as excellent benchmarks for any type of quantum hardware, and show that our system outperforms all other currently available hardware.

preprint2019arXiv

Combined sensitivity to the neutrino mass ordering with JUNO, the IceCube Upgrade, and PINGU

The ordering of the neutrino mass eigenstates is one of the fundamental open questions in neutrino physics. While current-generation neutrino oscillation experiments are able to produce moderate indications on this ordering, upcoming experiments of the next generation aim to provide conclusive evidence. In this paper we study the combined performance of the two future multi-purpose neutrino oscillation experiments JUNO and the IceCube Upgrade, which employ two very distinct and complementary routes towards the neutrino mass ordering. The approach pursued by the $20\,\mathrm{kt}$ medium-baseline reactor neutrino experiment JUNO consists of a careful investigation of the energy spectrum of oscillated $\barν_e$ produced by ten nuclear reactor cores. The IceCube Upgrade, on the other hand, which consists of seven additional densely instrumented strings deployed in the center of IceCube DeepCore, will observe large numbers of atmospheric neutrinos that have undergone oscillations affected by Earth matter. In a joint fit with both approaches, tension occurs between their preferred mass-squared differences $ Δm_{31}^{2}=m_{3}^{2}-m_{1}^{2} $ within the wrong mass ordering. In the case of JUNO and the IceCube Upgrade, this allows to exclude the wrong ordering at $>5σ$ on a timescale of 3--7 years --- even under circumstances that are unfavorable to the experiments' individual sensitivities. For PINGU, a 26-string detector array designed as a potential low-energy extension to IceCube, the inverted ordering could be excluded within 1.5 years (3 years for the normal ordering) in a joint analysis.

preprint2019arXiv

Design and Performance of the first IceAct Demonstrator at the South Pole

In this paper we describe the first results of a compact imaging air-Cherenkov telescope, IceAct, operating in coincidence with the IceCube Neutrino Observatory (IceCube) at the geographic South Pole. An array of IceAct telescopes (referred to as the IceAct project) is under consideration as part of the IceCube-Gen2 extension to IceCube. Surface detectors in general will be a powerful tool in IceCube-Gen2 for distinguishing astrophysical neutrinos from the dominant backgrounds of cosmic-ray induced atmospheric muons and neutrinos: the IceTop array is already in place as part of IceCube, but has a high energy threshold. Although the duty cycle will be lower for the IceAct telescopes than the present IceTop tanks, the IceAct telescopes may prove to be more effective at lowering the detection threshold for air showers. Additionally, small imaging air-Cherenkov telescopes in combination with IceTop, the deep IceCube detector or other future detector systems might improve measurements of the composition of the cosmic ray energy spectrum. In this paper we present measurements of a first 7-pixel imaging air Cherenkov telescope demonstrator, proving the capability of this technology to measure air showers at the South Pole in coincidence with IceTop and the deep IceCube detector.

preprint2019arXiv

Efficient propagation of systematic uncertainties from calibration to analysis with the SnowStorm method in IceCube

Efficient treatment of systematic uncertainties that depend on a large number of nuisance parameters is a persistent difficulty in particle physics experiments. Where low-level effects are not amenable to simple parameterization or re-weighting, analyses often rely on discrete simulation sets to quantify the effects of nuisance parameters on key analysis observables. Such methods may become computationally untenable for analyses requiring high statistics Monte Carlo with a large number of nuisance degrees of freedom, especially in cases where these degrees of freedom parameterize the shape of a continuous distribution. In this paper we present a method for treating systematic uncertainties in a computationally efficient and comprehensive manner using a single simulation set with multiple and continuously varied nuisance parameters. This method is demonstrated for the case of the depth-dependent effective dust distribution within the IceCube Neutrino Telescope.

preprint2019arXiv

First Measurement of the Charged Current $\overlineν_μ$ Double Differential Cross Section on a Water Target without Pions in the final state

This paper reports the first differential measurement of the charged-current $\overlineν_μ$ interaction cross section on water with no pions in the final state. The unfolded flux-averaged measurement using the T2K off-axis near detector is given in double differential bins of $μ^+$ momentum and angle. The integrated cross section in a restricted phase space is $σ=\left(1.11\pm0.18\right)\times10^{-38}$ cm$^{2}$ per water molecule. Comparisons with several nuclear models are also presented.

preprint2019arXiv

Measurement of neutrino and antineutrino neutral-current quasielastic-like interactions on oxygen by detecting nuclear de-excitation $γ$-rays

Neutrino- and antineutrino-oxygen neutral-current quasielastic-like interactions are measured at Super-Kamiokande using nuclear de-excitation $γ$-rays to identify signal-like interactions in data from a $14.94 \ (16.35)\times 10^{20}$ protons-on-target exposure of the T2K neutrino (antineutrino) beam. The measured flux-averaged cross sections on oxygen nuclei are $\langle σ_{ν{\rm -NCQE}} \rangle = 1.70 \pm 0.17 ({\rm stat.}) ^{+ {\rm 0.51}}_{- {\rm 0.38}} ({\rm syst.}) \times 10^{-38} \ {\rm cm^2/oxygen}$ with a flux-averaged energy of 0.82 GeV and $\langle σ_{\barν {\rm -NCQE}} \rangle = 0.98 \pm 0.16 ({\rm stat.}) ^{+ {\rm 0.26}}_{- {\rm 0.19}} ({\rm syst.}) \times 10^{-38} \ {\rm cm^2/oxygen}$ with a flux-averaged energy of 0.68 GeV, for neutrinos and antineutrinos, respectively. These results are the most precise to date, and the antineutrino result is the first cross section measurement of this channel. They are compared with various theoretical predictions. The impact on evaluation of backgrounds to searches for supernova relic neutrinos at present and future water Cherenkov detectors is also discussed.

preprint2019arXiv

Measurement of the muon neutrino charged-current single $π^+$ production on hydrocarbon using the T2K off-axis near detector ND280

We report the measurements of single and double differential cross section of muon neutrino charged-current interactions on carbon with a single positively charged pion in the final state at the T2K off-axis near detector using $5.56\times10^{20}$ protons on target. The analysis uses data control samples for the background subtraction and the cross section signal, defined as a single negatively charged muon and a single positively charged pion exiting from the target nucleus, is extracted using an unfolding method. The model dependent cross section, integrated over the T2K off-axis neutrino beam spectrum peaking at $0.6$~GeV, is measured to be $σ= (11.76 \pm 0.44 \text{(stat)} \pm 2.39 \text{(syst)}) \times 10^{-40} \text{cm}^2$~$\text{nucleon}^{-1}$. Various differential cross sections are measured, including the first measurement of the Adler angles for single charged pion production in neutrino interactions with heavy nuclei target.

preprint2019arXiv

Planck 2018 results. I. Overview and the cosmological legacy of Planck

The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter LCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (theta_*) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the LCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances.

preprint2019arXiv

Planck 2018 results. VIII. Gravitational lensing

We present measurements of the cosmic microwave background (CMB) lensing potential using the final $\textit{Planck}$ 2018 temperature and polarization data. We increase the significance of the detection of lensing in the polarization maps from $5\,σ$ to $9\,σ$. Combined with temperature, lensing is detected at $40\,σ$. We present an extensive set of tests of the robustness of the lensing-potential power spectrum, and construct a minimum-variance estimator likelihood over lensing multipoles $8 \le L \le 400$. We find good consistency between lensing constraints and the results from the $\textit{Planck}$ CMB power spectra within the $\rm{ΛCDM}$ model. Combined with baryon density and other weak priors, the lensing analysis alone constrains $σ_8 Ω_{\rm m}^{0.25}=0.589\pm 0.020$ ($1\,σ$ errors). Also combining with baryon acoustic oscillation (BAO) data, we find tight individual parameter constraints, $σ_8=0.811\pm0.019$, $H_0=67.9_{-1.3}^{+1.2}\,\text{km}\,\text{s}^{-1}\,\rm{Mpc}^{-1}$, and $Ω_{\rm m}=0.303^{+0.016}_{-0.018}$. Combining with $\textit{Planck}$ CMB power spectrum data, we measure $σ_8$ to better than $1\,\%$ precision, finding $σ_8=0.811\pm 0.006$. We find consistency with the lensing results from the Dark Energy Survey, and give combined lensing-only parameter constraints that are tighter than joint results using galaxy clustering. Using $\textit{Planck}$ cosmic infrared background (CIB) maps we make a combined estimate of the lensing potential over $60\,\%$ of the sky with considerably more small-scale signal. We demonstrate delensing of the $\textit{Planck}$ power spectra, detecting a maximum removal of $40\,\%$ of the lensing-induced power in all spectra. The improvement in the sharpening of the acoustic peaks by including both CIB and the quadratic lensing reconstruction is detected at high significance (abridged).

preprint2019arXiv

Planck 2018 results. X. Constraints on inflation

We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be $n_\mathrm{s}=0.9649\pm 0.0042$ at 68% CL and show no evidence for a scale dependence of $n_\mathrm{s}.$ Spatial flatness is confirmed at a precision of 0.4% at 95% CL with the combination with BAO data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, $r_{0.002}<0.10$, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain $r_{0.002}<0.056$. In the framework of single-field inflationary models with Einstein gravity, these results imply that: (a) slow-roll models with a concave potential, $V" (ϕ) < 0,$ are increasingly favoured by the data; and (b) two different methods for reconstructing the inflaton potential find no evidence for dynamics beyond slow roll. Non-parametric reconstructions of the primordial power spectrum consistently confirm a pure power law. A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectrum, a result further strengthened for certain oscillatory models by a new combined analysis that includes Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadrupolar modulation of the primordial fluctuations. However, the polarization data do not confirm physical models for a scale-dependent dipolar modulation.

preprint2019arXiv

Planck 2018 results. XII. Galactic astrophysics using polarized dust emission

We present 353 GHz full-sky maps of the polarization fraction $p$, angle $ψ$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on the zero level in total intensity. The observed inverse behaviour between $p$ and $S$ is interpreted with models of the polarized sky that include effects from only the topology of the turbulent Galactic magnetic field. Thus, the statistical properties of $p$, $ψ$, and $S$ mostly reflect the structure of the magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map $S \times p$, looking for residual trends. While $p$ decreases by a factor of 3--4 between $N_H=10^{20}$ cm$^{-2}$ and $N_H=2\times 10^{22}$ cm$^{-2}$, $S \times p$ decreases by only about 25%, a systematic trend observed in both the diffuse ISM and molecular clouds. Second, we find no systematic trend of $S \times p$ with the dust temperature, even though in the diffuse ISM lines of sight with high $p$ and low $S$ tend to have colder dust. We also compare Planck data with starlight polarization in the visible at high latitudes. The agreement in polarization angles is remarkable. Two polarization emission-to-extinction ratios that characterize dust optical properties depend only weakly on $N_H$ and converge towards the values previously determined for translucent lines of sight. We determine an upper limit for the polarization fraction in extinction of 13%, compatible with the $p_\mathrm{max}$ observed in emission. These results provide strong constraints for models of Galactic dust in diffuse gas.

preprint2019arXiv

Search for Electron Antineutrino Appearance in a Long-baseline Muon Antineutrino Beam

Electron antineutrino appearance is measured by the T2K experiment in an accelerator-produced antineutrino beam, using additional neutrino beam operation to constrain parameters of the PMNS mixing matrix. T2K observes 15 candidate electron antineutrino events with a background expectation of 9.3 events. Including information from the kinematic distribution of observed events, the hypothesis of no electron antineutrino appearance is disfavored with a significance of 2.40σ and no discrepancy between data and PMNS predictions is found. A complementary analysis that introduces an additional free parameter which allows non-PMNS values of electron neutrino and antineutrino appearance also finds no discrepancy between data and PMNS predictions.

preprint2019arXiv

Search for Sources of Astrophysical Neutrinos Using Seven Years of IceCube Cascade Events

Low background searches for astrophysical neutrino sources anywhere in the sky can be performed using cascade events induced by neutrinos of all flavors interacting in IceCube with energies as low as ~1 TeV. Previously, we showed that even with just two years of data, the resulting sensitivity to sources in the southern sky is competitive with IceCube and ANTARES analyses using muon tracks induced by charge current muon neutrino interactions - especially if the neutrino emission follows a soft energy spectrum or originates from an extended angular region. Here, we extend that work by adding five more years of data, significantly improving the cascade angular resolution, and including tests for point-like or diffuse Galactic emission to which this dataset is particularly well-suited. For many of the signal candidates considered, this analysis is the most sensitive of any experiment. No significant clustering was observed, and thus many of the resulting constraints are the most stringent to date. In this paper we will describe the improvements introduced in this analysis and discuss our results in the context of other recent work in neutrino astronomy.

preprint2019arXiv

Search for transient optical counterparts to high-energy IceCube neutrinos with Pan-STARRS1

In order to identify the sources of the observed diffuse high-energy neutrino flux, it is crucial to discover their electromagnetic counterparts. IceCube began releasing alerts for single high-energy ($E > 60$ TeV) neutrino detections with sky localisation regions of order 1 deg radius in 2016. We used Pan-STARRS1 to follow-up five of these alerts during 2016-2017 to search for any optical transients that may be related to the neutrinos. Typically 10-20 faint ($m < 22.5$ mag) extragalactic transients are found within the Pan-STARRS1 footprints and are generally consistent with being unrelated field supernovae (SNe) and AGN. We looked for unusual properties of the detected transients, such as temporal coincidence of explosion epoch with the IceCube timestamp. We found only one transient that had properties worthy of a specific follow-up. In the Pan-STARRS1 imaging for IceCube-160427A (probability to be of astrophysical origin of $\sim$50 %), we found a SN PS16cgx, located at 10.0' from the nominal IceCube direction. Spectroscopic observations of PS16cgx showed that it was an H-poor SN at z = 0.2895. The spectra and light curve resemble some high-energy Type Ic SNe, raising the possibility of a jet driven SN with an explosion epoch temporally coincident with the neutrino detection. However, distinguishing Type Ia and Type Ic SNe at this redshift is notoriously difficult. Based on all available data we conclude that the transient is more likely to be a Type Ia with relatively weak SiII absorption and a fairly normal rest-frame r-band light curve. If, as predicted, there is no high-energy neutrino emission from Type Ia SNe, then PS16cgx must be a random coincidence, and unrelated to the IceCube-160427A. We find no other plausible optical transient for any of the five IceCube events observed down to a 5$σ$ limiting magnitude of $m \sim 22$ mag, between 1 day and 25 days after detection.

preprint2019arXiv

Time-integrated Neutrino Source Searches with 10 years of IceCube Data

This paper presents the results from point-like neutrino source searches using ten years of IceCube data collected between Apr.~6, 2008 and Jul.~10, 2018. We evaluate the significance of an astrophysical signal from a point-like source looking for an excess of clustered neutrino events with energies typically above $\sim1\,$TeV among the background of atmospheric muons and neutrinos. We perform a full-sky scan, a search within a selected source catalog, a catalog population study, and three stacked Galactic catalog searches. The most significant point in the Northern hemisphere from scanning the sky is coincident with the Seyfert II galaxy NGC 1068, which was included in the source catalog search. The excess at the coordinates of NGC 1068 is inconsistent with background expectations at the level of $2.9\,σ$ after accounting for statistical trials. The combination of this result along with excesses observed at the coordinates of three other sources, including TXS 0506+056, suggests that, collectively, correlations with sources in the Northern catalog are inconsistent with background at 3.3$\,σ$ significance. These results, all based on searches for a cumulative neutrino signal integrated over the ten years of available data, motivate further study of these and similar sources, including time-dependent analyses, multimessenger correlations, and the possibility of stronger evidence with coming upgrades to the detector.

preprint2018arXiv

Planck 2018 results. II. Low Frequency Instrument data processing

We present a final description of the data-processing pipeline for the Planck, Low Frequency Instrument (LFI), implemented for the 2018 data release. Several improvements have been made with respect to the previous release, especially in the calibration process and in the correction of instrumental features such as the effects of nonlinearity in the response of the analogue-to-digital converters. We provide a brief pedagogical introduction to the complete pipeline, as well as a detailed description of the important changes implemented. Self-consistency of the pipeline is demonstrated using dedicated simulations and null tests. We present the final version of the LFI full sky maps at 30, 44, and 70 GHz, both in temperature and polarization, together with a refined estimate of the Solar dipole and a final assessment of the main LFI instrumental parameters.

preprint2018arXiv

Planck 2018 results. XI. Polarized dust foregrounds

The study of polarized dust emission has become entwined with the analysis of the cosmic microwave background (CMB) polarization. We use new Planck maps to characterize Galactic dust emission as a foreground to the CMB polarization. We present Planck EE, BB, and TE power spectra of dust polarization at 353 GHz for six nested sky regions covering from 24 to 71 % of the sky. We present power-law fits to the angular power spectra, yielding evidence for statistically significant variations of the exponents over sky regions and a difference between the values for the EE and BB spectra. The TE correlation and E/B power asymmetry extend to low multipoles that were not included in earlier Planck polarization papers. We also report evidence for a positive TB dust signal. Combining data from Planck and WMAP, we determine the amplitudes and spectral energy distributions (SEDs) of polarized foregrounds, including the correlation between dust and synchrotron polarized emission, for the six sky regions as a function of multipole. This quantifies the challenge of the component separation procedure required for detecting the reionization and recombination peaks of primordial CMB B modes. The SED of polarized dust emission is fit well by a single-temperature modified blackbody emission law from 353 GHz to below 70 GHz. For a dust temperature of 19.6 K, the mean spectral index for dust polarization is $β_{\rm d}^{P} = 1.53\pm0.02 $. By fitting multi-frequency cross-spectra, we examine the correlation of the dust polarization maps across frequency. We find no evidence for decorrelation. If the Planck limit for the largest sky region applies to the smaller sky regions observed by sub-orbital experiments, then decorrelation might not be a problem for CMB experiments aiming at a primordial B-mode detection limit on the tensor-to-scalar ratio $r\simeq0.01$ at the recombination peak.

preprint2018arXiv

Search for Multi-messenger Sources of Gravitational Waves and High-energy Neutrinos with Advanced LIGO during its first Observing Run, ANTARES and IceCube

Astrophysical sources of gravitational waves, such as binary neutron star and black hole mergers or core-collapse supernovae, can drive relativistic outflows, giving rise to non-thermal high-energy emission. High-energy neutrinos are signatures of such outflows. The detection of gravitational waves and high-energy neutrinos from common sources could help establish the connection between the dynamics of the progenitor and the properties of the outflow. We searched for associated emission of gravitational waves and high-energy neutrinos from astrophysical transients with minimal assumptions using data from Advanced LIGO from its first observing run O1, and data from the ANTARES and IceCube neutrino observatories from the same time period. We focused on candidate events whose astrophysical origin could not be determined from a single messenger. We found no significant coincident candidate, which we used to constrain the rate density of astrophysical sources dependent on their gravitational wave and neutrino emission processes.

preprint2017arXiv

Construction of KAGRA: an Underground Gravitational Wave Observatory

Major construction and initial-phase operation of a second-generation gravitational-wave detector KAGRA has been completed. The entire 3-km detector is installed underground in a mine in order to be isolated from background seismic vibrations on the surface. This allows us to achieve a good sensitivity at low frequencies and high stability of the detector. Bare-bones equipment for the interferometer operation has been installed and the first test run was accomplished in March and April of 2016 with a rather simple configuration. The initial configuration of KAGRA is named {\it iKAGRA}. In this paper, we summarize the construction of KAGRA, including the study of the advantages and challenges of building an underground detector and the operation of the iKAGRA interferometer together with the geophysics interferometer that has been constructed in the same tunnel.

preprint2017arXiv

Search for High-energy Neutrinos from Binary Neutron Star Merger GW170817 with ANTARES, IceCube, and the Pierre Auger Observatory

The Advanced LIGO and Advanced Virgo observatories recently discovered gravitational waves from a binary neutron star inspiral. A short gamma-ray burst (GRB) that followed the merger of this binary was also recorded by the Fermi Gamma-ray Burst Monitor (Fermi-GBM), and the Anticoincidence Shield for the Spectrometer for the International Gamma-Ray Astrophysics Laboratory (INTEGRAL), indicating particle acceleration by the source. The precise location of the event was determined by optical detections of emission following the merger. We searched for high-energy neutrinos from the merger in the GeV--EeV energy range using the ANTARES, IceCube, and Pierre Auger Observatories. No neutrinos directionally coincident with the source were detected within $\pm500$ s around the merger time. Additionally, no MeV neutrino burst signal was detected coincident with the merger. We further carried out an extended search in the direction of the source for high-energy neutrinos within the 14-day period following the merger, but found no evidence of emission. We used these results to probe dissipation mechanisms in relativistic outflows driven by the binary neutron star merger. The non-detection is consistent with model predictions of short GRBs observed at a large off-axis angle.

preprint2017arXiv

The status of KAGRA underground cryogenic gravitational wave telescope

KAGRA is a 3-km interferometric gravitational wave telescope located in the Kamioka mine in Japan. It is the first km-class gravitational wave telescope constructed underground to reduce seismic noise, and the first km-class telescope to use cryogenic cooling of test masses to reduce thermal noise. The construction of the infrastructure to house the interferometer in the tunnel, and the initial phase operation of the interferometer with a simple 3-km Michelson configuration have been completed. The first cryogenic operation is expected in 2018, and the observing runs with a full interferometer are expected in 2020s. The basic interferometer configuration and the current status of KAGRA are described.

preprint2016arXiv

A First Targeted Search for Gravitational-Wave Bursts from Core-Collapse Supernovae in Data of First-Generation Laser Interferometer Detectors

We present results from a search for gravitational-wave bursts coincident with a set of two core-collapse supernovae observed between 2007 and 2011. We employ data from the Laser Interferometer Gravitational-wave Observatory (LIGO), the Virgo gravitational-wave observatory, and the GEO 600 gravitational-wave observatory. The targeted core-collapse supernovae were selected on the basis of (1) proximity (within approximately 15 Mpc), (2) tightness of observational constraints on the time of core collapse that defines the gravitational-wave search window, and (3) coincident operation of at least two interferometers at the time of core collapse. We find no plausible gravitational-wave candidates. We present the probability of detecting signals from both astrophysically well-motivated and more speculative gravitational-wave emission mechanisms as a function of distance from Earth, and discuss the implications for the detection of gravitational waves from core-collapse supernovae by the upgraded Advanced LIGO and Virgo detectors.

preprint2016arXiv

All-sky search for long-duration gravitational wave transients with LIGO

We present the results of a search for long-duration gravitational wave transients in two sets of data collected by the LIGO Hanford and LIGO Livingston detectors between November 5, 2005 and September 30, 2007, and July 7, 2009 and October 20, 2010, with a total observational time of 283.0 days and 132.9 days, respectively. The search targets gravitational wave transients of duration 10 - 500 s in a frequency band of 40 - 1000 Hz, with minimal assumptions about the signal waveform, polarization, source direction, or time of occurrence. All candidate triggers were consistent with the expected background; as a result we set 90% confidence upper limits on the rate of long-duration gravitational wave transients for different types of gravitational wave signals. For signals from black hole accretion disk instabilities, we set upper limits on the source rate density between $3.4 \times 10^{-5}$ - $9.4 \times 10^{-4}$ Mpc$^{-3}$ yr$^{-1}$ at 90% confidence. These are the first results from an all-sky search for unmodeled long-duration transient gravitational waves.

preprint2016arXiv

An improved analysis of GW150914 using a fully spin-precessing waveform model

This paper presents updated estimates of source parameters for GW150914, a binary black-hole coalescence event detected by the Laser Interferometer Gravitational-wave Observatory (LIGO) on September 14, 2015 [1]. Reference presented parameter estimation [2] of the source using a 13-dimensional, phenomenological precessing-spin model (precessing IMRPhenom) and a 11-dimensional nonprecessing effective-one-body (EOB) model calibrated to numerical-relativity simulations, which forces spin alignment (nonprecessing EOBNR). Here we present new results that include a 15-dimensional precessing-spin waveform model (precessing EOBNR) developed within the EOB formalism. We find good agreement with the parameters estimated previously [2], and we quote updated component masses of $35^{+5}_{-3}\mathrm{M}_\odot$ and $30^{+3}_{-4}\mathrm{M}_\odot$ (where errors correspond to 90% symmetric credible intervals). We also present slightly tighter constraints on the dimensionless spin magnitudes of the two black holes, with a primary spin estimate $0.65$ and a secondary spin estimate $0.75$ at 90% probability. Reference [2] estimated the systematic parameter-extraction errors due to waveform-model uncertainty by combining the posterior probability densities of precessing IMRPhenom and nonprecessing EOBNR. Here we find that the two precessing-spin models are in closer agreement, suggesting that these systematic errors are smaller than previously quoted.

preprint2016arXiv

An invisible non-volatile solid-state memory

Information technologies require entangling data stability with encryption for a next generation of secure data storage. Current magnetic memories, ranging from low-density stripes up to high-density hard drives, can ultimately be detected using routinely available probes or manipulated by external magnetic perturbations. Antiferromagnetic resistors feature unrivalled robustness but the stable resistive states reported scarcely differ by more than a fraction of a percent at room temperature. Here we show that the metamagnetic (ferromagnetic to antiferromagnetic) transition in intermetallic Fe0.50Rh0.50 can be electrically controlled in a magnetoelectric heterostructure to reveal or cloak a given ferromagnetic state. From an aligned ferromagnetic phase, magnetic states are frozen into the antiferromagnetic phase by the application of an electric field, thus eliminating the stray field and likewise making it insensitive to external magnetic field. Application of a reverse electric field reverts the antiferromagnetic state to the original ferromagnetic state. Our work demonstrates the building blocks of a feasible, extremely stable, non-volatile, electrically addressable, low-energy dissipation, magnetoelectric multiferroic memory.

preprint2016arXiv

Co-Designing a Scalable Quantum Computer with Trapped Atomic Ions

The first generation of quantum computers are on the horizon, fabricated from quantum hardware platforms that may soon be able to tackle certain tasks that cannot be performed or modelled with conventional computers. These quantum devices will not likely be universal or fully programmable, but special-purpose processors whose hardware will be tightly co-designed with particular target applications. Trapped atomic ions are a leading platform for first generation quantum computers, but are also fundamentally scalable to more powerful general purpose devices in future generations. This is because trapped ion qubits are atomic clock standards that can be made identical to a part in 10^15, and their quantum circuit connectivity can be reconfigured through the use of external fields, without modifying the arrangement or architecture of the qubits themselves. In this article we show how a modular quantum computer of any size can be engineered from ion crystals, and how the wiring between ion trap qubits can be tailored to a variety of applications and quantum computing protocols.

preprint2016arXiv

Comprehensive All-sky Search for Periodic Gravitational Waves in the Sixth Science Run LIGO Data

We report on a comprehensive all-sky search for periodic gravitational waves in the frequency band 100-1500 Hz and with a frequency time derivative in the range of $[-1.18, +1.00]\times 10^{-8}$ Hz/s. Such a signal could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the data from the Initial LIGO sixth science run and covers a larger parameter space with respect to any past search. A Loosely Coherent detection pipeline was applied to follow up weak outliers in both Gaussian (95% recovery rate) and non-Gaussian (75% recovery rate) bands. No gravitational wave signals were observed, and upper limits were placed on their strength. Our smallest upper limit on worst-case (linearly polarized) strain amplitude $h_0$ is ${9.7}\times 10^{-25}$ near 169 Hz, while at the high end of our frequency range we achieve a worst-case upper limit of ${5.5}\times 10^{-24}$. Both cases refer to all sky locations and entire range of frequency derivative values.

preprint2016arXiv

Design Study of a Superconducting Gantry for Carbon Beam Therapy

This paper describes the design study of a gantry for a carbon beam. The designed gantry is compact such that its size is comparable to the size of the proton gantry. This is possible by introducing superconducting double helical coils for dipole magnets. The gantry optics is designed in such a way that it provides rotation-invariant optics and variable beam size as well as point-to-parallel scanning of a beam. For large-aperture magnet, three-dimensional magnetic field distribution is obtained by invoking a computer code, and a number of particles are tracked by integrating equations of motion numerically together with three-dimensional interpolation. The beam-shape distortion due to the fringe field is reduced to an acceptable level by optimizing the coil windings with the help of genetic algorithm. Higher-order transfer coefficients are calculated and shown to be reduced greatly with appropriate optimization of the coil windings.

preprint2016arXiv

Directly comparing GW150914 with numerical solutions of Einstein's equations for binary black hole coalescence

We compare GW150914 directly to simulations of coalescing binary black holes in full general relativity, accounting for all the spin-weighted quadrupolar modes, and separately accounting for all the quadrupolar and octopolar modes. Consistent with the posterior distributions reported in LVC_PE[1] (at 90% confidence), we find the data are compatible with a wide range of nonprecessing and precessing simulations. Followup simulations performed using previously-estimated binary parameters most resemble the data. Comparisons including only the quadrupolar modes constrain the total redshifted mass Mz \in [64 - 82M_\odot], mass ratio q = m2/m1 \in [0.6,1], and effective aligned spin χ_eff \in [-0.3, 0.2], where χ_{eff} = (S1/m1 + S2/m2) \cdot\hat{L} /M. Including both quadrupolar and octopolar modes, we find the mass ratio is even more tightly constrained. Simulations with extreme mass ratios and effective spins are highly inconsistent with the data, at any mass. Several nonprecessing and precessing simulations with similar mass ratio and χ_{eff} are consistent with the data. Though correlated, the components' spins (both in magnitude and directions) are not significantly constrained by the data. For nonprecessing binaries, interpolating between simulations, we reconstruct a posterior distribution consistent with previous results. The final black hole's redshifted mass is consistent with Mf,z between 64.0 - 73.5M_\odot and the final black hole's dimensionless spin parameter is consistent with af = 0.62 - 0.73. As our approach invokes no intermediate approximations to general relativity and can strongly reject binaries whose radiation is inconsistent with the data, our analysis provides a valuable complement to LVC_PE[1].

preprint2016arXiv

Enhanced nonlinear refractive index in epsilon-near-zero materials

New propagation regimes for light arise from the ability to tune the dielectric permittivity to extremely low values. Here we demonstrate a universal approach based on the low linear permittivity values attained in the epsilon-near-zero (ENZ) regime for enhancing the nonlinear refractive index, which enables remarkable light-induced changes of the material properties. Experiments performed on Al-doped ZnO (AZO) thin films show a six-fold increase of the Kerr nonlinear refractive index ($n_2$) at the ENZ wavelength, located in the 1300 nm region. This in turn leads to ultrafast light-induced refractive index changes of the order of unity, thus representing a new paradigm for nonlinear optics.

preprint2016arXiv

GW150914: First results from the search for binary black hole coalescence with Advanced LIGO

On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-wave Observatory (LIGO) simultaneously observed the binary black hole merger GW150914. We report the results of a matched-filter search using relativistic models of compact-object binaries that recovered GW150914 as the most significant event during the coincident observations between the two LIGO detectors from September 12 to October 20, 2015. GW150914 was observed with a matched filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203 000 years, equivalent to a significance greater than 5.1 σ.

preprint2016arXiv

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

preprint2016arXiv

Measurement of coherent $π^{+}$ production in low energy neutrino-Carbon scattering

We report the first measurement of the flux-averaged cross section for charged current coherent $π^{+}$ production on carbon for neutrino energies less than 1.5 GeV to a restricted final state phase space region in the T2K near detector, ND280. Comparisons are made with predictions from the Rein-Sehgal coherent production model and the model by Alvarez-Ruso {\it et al.}, the latter representing the first implementation of an instance of the new class of microscopic coherent models in a neutrino interaction Monte Carlo event generator. This results contradicts the null results reported by K2K and SciBooNE in a similar neutrino energy region.

preprint2016arXiv

Measurement of double-differential muon neutrino charged-current interactions on C$_8$H$_8$ without pions in the final state using the T2K off-axis beam

We report the measurement of muon neutrino charged-current interactions on carbon without pions in the final state at the T2K beam energy using 5.734$\times10^{20}$ protons on target. For the first time the measurement is reported as a flux-integrated, double-differential cross-section in muon kinematic variables ($\cosθ_μ$, $p_μ$), without correcting for events where a pion is produced and then absorbed by final state interactions. Two analyses are performed with different selections, background evaluations and cross-section extraction methods to demonstrate the robustness of the results against biases due to model-dependent assumptions. The measurements compare favorably with recent models which include nucleon-nucleon correlations but, given the present precision, the measurement does not solve the degeneracy between different models. The data also agree with Monte Carlo simulations which use effective parameters that are tuned to external data to describe the nuclear effects. The total cross-section in the full phase space is $σ= (0.417 \pm 0.047 \text{(syst)} \pm 0.005 \text{(stat)})\times 10^{-38} \text{cm}^2$ $\text{nucleon}^{-1}$ and the cross-section integrated in the region of phase space with largest efficiency and best signal-over-background ratio ($\cosθ_μ>0.6$ and $p_μ> 200$ MeV) is $σ= (0.202 \pm 0.0359 \text{(syst)} \pm 0.0026 \text{(stat)}) \times 10^{-38} \text{cm}^2$ $\text{nucleon}^{-1}$.

preprint2016arXiv

Measurement of Muon Antineutrino Oscillations with an Accelerator-Produced Off-Axis Beam

T2K reports its first measurements of the parameters governing the disappearance of $\barν_μ$ in an off-axis beam due to flavor change induced by neutrino oscillations. The quasimonochromatic $\barν_μ$ beam, produced with a peak energy of 0.6 GeV at J-PARC, is observed at the far detector Super-Kamiokande, 295 km away, where the $\barν_μ$ survival probability is expected to be minimal. Using a dataset corresponding to $4.01 \times 10^{20}$ protons on target, $34$ fully contained $μ$-like events were observed. The best-fit oscillation parameters are $\sin^2 (\barθ_{23}) = 0.45$ and $|Δ\bar{m}^2_{32}| = 2.51 \times 10^{-3}$ eV$^2$ with 68% confidence intervals of 0.38 - 0.64 and 2.26 - 2.80 $\times 10^{-3}$ eV$^2$ respectively. These results are in agreement with existing antineutrino parameter measurements and also with the $ν_μ$ disappearance parameters measured by T2K.

preprint2016arXiv

Optimization of the Multi-turn Injection Efficiency for Medical Synchrotron

We present a method for optimization of the multi-turn injection effciency for medical synchrotron. We show that for given injection energy the injection efficiency can be greatly enhanced by choosing transverse tunes appropriately as well as optimizing the injection bump and the number of turns required for beam injection. We verify our study by applying the method to the Korea Heavy Ion Medical Accelerator (KHIMA) synchrotron which is currently built at the campus of Dongnam Institute of Radiological and Medical Sciences (DIRAMS) in Busan, Korea. First the frequency map analysis is performed with the help of ELEGANT and ACCSIM codes. The tunes which yield the good injection efficiency are then selected. With these tunes the injection bump and the number of turns required for injection are then optimized by tracking a number of particles up to one thousand turns after injection beyond which there is no further beam loss. Results for optimization of the injection efficiency for proton ion are presented.

preprint2016arXiv

Planck 2015 results. XIX. Constraints on primordial magnetic fields

We compute and investigate four types of imprint of a stochastic background of primordial magnetic fields (PMFs) on the cosmic microwave background (CMB) anisotropies: the impact of PMFs on the CMB spectra; the effect on CMB polarization induced by Faraday rotation; the impact of PMFs on the ionization history; magnetically-induced non-Gaussianities; and the magnetically-induced breaking of statistical isotropy. Overall, Planck data constrain the amplitude of PMFs to less than a few nanogauss. In particular, individual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are $B_{1\,\mathrm{Mpc}}< 4.4$ nG (where $B_{1\,\mathrm{Mpc}}$ is the comoving field amplitude at a scale of 1 Mpc) at 95% confidence level, assuming zero helicity, and $B_{1\,\mathrm{Mpc}}< 5.6$ nG for a maximally helical field.For nearly scale-invariant PMFs we obtain $B_{1\,\mathrm{Mpc}}<2.0$ nG and $B_{1\,\mathrm{Mpc}}<0.9$ nG if the impact of PMFs on the ionization history of the Universe is included. From the analysis of magnetically-induced non-Gaussianity we obtain three different values, corresponding to three applied methods, all below 5 nG. The constraint from the magnetically-induced passive-tensor bispectrum is $B_{1\,\mathrm{Mpc}}< 2.8$ nG. A search for preferred directions in the magnetically-induced passive bispectrum yields $B_{1\,\mathrm{Mpc}}< 4.5$ nG, whereas the the compensated-scalar bispectrum gives $B_{1\,\mathrm{Mpc}}< 3$ nG. The analysis of the Faraday rotation of CMB polarization by PMFs uses the Planck power spectra in $EE$ and $BB$ at 70 GHz and gives $B_{1\,\mathrm{Mpc}}< 1380$ nG. In our final analysis, we consider the harmonic-space correlations produced by Alfvén waves, finding no significant evidence for the presence of these waves. Together, these results comprise a comprehensive set of constraints on possible PMFs with Planck data.

preprint2016arXiv

Planck 2015 results. XVI. Isotropy and statistics of the CMB

We test the statistical isotropy and Gaussianity of the cosmic microwave background (CMB) anisotropies using observations made by the Planck satellite. Our results are based mainly on the full Planck mission for temperature, but also include some polarization measurements. In particular, we consider the CMB anisotropy maps derived from the multi-frequency Planck data by several component-separation methods. For the temperature anisotropies, we find excellent agreement between results based on these sky maps over both a very large fraction of the sky and a broad range of angular scales, establishing that potential foreground residuals do not affect our studies. Tests of skewness, kurtosis, multi-normality, N-point functions, and Minkowski functionals indicate consistency with Gaussianity, while a power deficit at large angular scales is manifested in several ways, for example low map variance. The results of a peak statistics analysis are consistent with the expectations of a Gaussian random field. The "Cold Spot" is detected with several methods, including map kurtosis, peak statistics, and mean temperature profile. We thoroughly probe the large-scale dipolar power asymmetry, detecting it with several independent tests, and address the subject of a posteriori correction. Tests of directionality suggest the presence of angular clustering from large to small scales, but at a significance that is dependent on the details of the approach. We perform the first examination of polarization data, finding the morphology of stacked peaks to be consistent with the expectations of statistically isotropic simulations. Where they overlap, these results are consistent with the Planck 2013 analysis based on the nominal mission data and provide our most thorough view of the statistics of the CMB fluctuations to date.

preprint2016arXiv

Planck 2015 results. XVII. Constraints on primordial non-Gaussianity

The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (KSW), binned, and modal - we obtain consistent values for the local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone fNL^local=2.5+\-5.7, fNL^equil=-16+\-70 and fNL^ortho=-34+\-33(68%CL). Combining temperature and polarization data we obtain fNL^local=0.8+\-5.0, fNL^equil=-4+\-43 and fNL^ortho=-26+\-21 (68%CL). The results are based on cross-validation of these estimators on simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with Minkowski functionals based measurements. The effect of time-domain de-glitching systematics on the bispectrum is negligible. In spite of these test outcomes we conservatively label the results including polarization data as preliminary, owing to a known mismatch of the noise model in simulations and the data. Beyond fNL estimates, we present model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally dependent vector models. We also present a wide survey of scale-dependent oscillatory bispectra, and we look for isocurvature NG. Our constraint on the local primordial trispectrum amplitude is gNL^local=(-9.0+\-7.7)x10^4 (68%CL), and we perform an analysis of additional trispectrum shapes. The global picture is one of consistency with the premises of the LambdaCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.[abridged]

preprint2016arXiv

Properties of the Binary Black Hole Merger GW150914

On September 14, 2015, the Laser Interferometer Gravitational-wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in general relativity. GW150914 was produced by a nearly equal mass binary black hole of $36^{+5}_{-4} M_\odot$ and $29^{+4}_{-4} M_\odot$; for each parameter we report the median value and the range of the 90% credible interval. The dimensionless spin magnitude of the more massive black hole is bound to be $<0.7$ (at 90% probability). The luminosity distance to the source is $410^{+160}_{-180}$ Mpc, corresponding to a redshift $0.09^{+0.03}_{-0.04}$ assuming standard cosmology. The source location is constrained to an annulus section of $610$ deg$^2$, primarily in the southern hemisphere. The binary merges into a black hole of $62^{+4}_{-4} M_\odot$ and spin $0.67^{+0.05}_{-0.07}$. This black hole is significantly more massive than any other inferred from electromagnetic observations in the stellar-mass regime.

preprint2016arXiv

Search for transient gravitational waves in coincidence with short duration radio transients during 2007-2013

We present an archival search for transient gravitational-wave bursts in coincidence with 27 single pulse triggers from Green Bank Telescope pulsar surveys, using the LIGO, Virgo and GEO interferometer network. We also discuss a check for gravitational-wave signals in coincidence with Parkes Fast Radio Bursts using similar methods. Data analyzed in these searches were collected between 2007 and 2013. Possible sources of emission of both short-duration radio signals and transient gravitational-wave emission include starquakes on neutron stars, binary coalescence of neutron stars, and cosmic string cusps. While no evidence for gravitational-wave emission in coincidence with these radio transients was found, the current analysis serves as a prototype for similar future searches using more sensitive second-generation interferometers.

preprint2016arXiv

Sensitivity of the T2K accelerator-based neutrino experiment with an Extended run to $20\times10^{21}$ POT

Recent measurements at the T2K experiment indicate that CP violation in neutrino mixing may be observed in the future by long-baseline neutrino oscillation experiments. We explore the physics program of an extension to the currently approved T2K running of $7.8\times 10^{21}$ protons-on-target to $20\times 10^{21}$ protons-on-target,aiming at initial observation of CP violation with 3$\,σ$ or higher significance for the case of maximum CP violation. With accelerator and beam line upgrades, as well as analysis improvements, this program would occur before the next generation of long-baseline neutrino oscillation experiments that are expected to start operation in 2026.

preprint2016arXiv

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

preprint2016arXiv

Supplement: The Rate of Binary Black Hole Mergers Inferred from Advanced LIGO Observations Surrounding GW150914

Supplemental information for a Letter reporting the rate of binary black hole (BBH) coalescences inferred from 16 days of coincident Advanced LIGO observations surrounding the transient gravitational wave signal GW150914. In that work we reported various rate estimates whose 90\% credible intervals fell in the range $2$--$600 \, \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$. Here we give details of our method and computations, including information about our search pipelines, a derivation of our likelihood function for the analysis, a description of the astrophysical search trigger distribution expected from merging BBHs, details on our computational methods, a description of the effects and our model for calibration uncertainty, and an analytic method of estimating our detector sensitivity that is calibrated to our measurements.

preprint2016arXiv

The Rate of Binary Black Hole Mergers Inferred from Advanced LIGO Observations Surrounding GW150914

A transient gravitational-wave signal, GW150914, was identified in the twin Advanced LIGO detectors on September 14, 2015 at 09:50:45 UTC. To assess the implications of this discovery, the detectors remained in operation with unchanged configurations over a period of 39 d around the time of the signal. At the detection statistic threshold corresponding to that observed for GW150914, our search of the 16 days of simultaneous two-detector observational data is estimated to have a false alarm rate (FAR) of $< 4.9 \times 10^{-6} \, \mathrm{yr}^{-1}$, yielding a $p$-value for GW150914 of $< 2 \times 10^{-7}$. Parameter estimation followup on this trigger identifies its source as a binary black hole (BBH) merger with component masses $(m_1, m_2) = \left(36^{+5}_{-4},29^{+4}_{-4}\right) \, M_\odot$ at redshift $z = 0.09^{+0.03}_{-0.04}$ (median and 90\% credible range). Here we report on the constraints these observations place on the rate of BBH coalescences. Considering only GW150914, assuming that all BBHs in the Universe have the same masses and spins as this event, imposing a search FAR threshold of 1 per 100 years, and assuming that the BBH merger rate is constant in the comoving frame, we infer a 90% credible range of merger rates between $2$--$53 \, \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$ (comoving frame). Incorporating all search triggers that pass a much lower threshold while accounting for the uncertainty in the astrophysical origin of each trigger, we estimate a higher rate, ranging from $13$--$600 \, \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$ depending on assumptions about the BBH mass distribution. All together, our various rate estimates fall in the conservative range $2$--$600 \, \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$.

preprint2016arXiv

Upper limits on the rates of binary neutron star and neutron-star--black-hole mergers from Advanced LIGO's first observing run

We report here the non-detection of gravitational waves from the merger of binary neutron star systems and neutron-star--black-hole systems during the first observing run of Advanced LIGO. In particular we searched for gravitational wave signals from binary neutron star systems with component masses $\in [1,3] M_{\odot}$ and component dimensionless spins $< 0.05$. We also searched for neutron-star--black-hole systems with the same neutron star parameters, black hole mass $\in [2,99] M_{\odot}$ and no restriction on the black hole spin magnitude. We assess the sensitivity of the two LIGO detectors to these systems, and find that they could have detected the merger of binary neutron star systems with component mass distributions of $1.35\pm0.13 M_{\odot}$ at a volume-weighted average distance of $\sim$ 70Mpc, and for neutron-star--black-hole systems with neutron star masses of $1.4M_\odot$ and black hole masses of at least $5M_\odot$, a volume-weighted average distance of at least $\sim$ 110Mpc. From this we constrain with 90% confidence the merger rate to be less than 12,600 Gpc$^{-3}$yr$^{-1}$ for binary-neutron star systems and less than 3,600 Gpc$^{-3}$yr$^{-1}$ for neutron-star--black-hole systems. We find that if no detection of neutron-star binary mergers is made in the next two Advanced LIGO and Advanced Virgo observing runs we would place significant constraints on the merger rates. Finally, assuming a rate of $10^{+20}_{-7}$Gpc$^{-3}$yr$^{-1}$ short gamma ray bursts beamed towards the Earth and assuming that all short gamma-ray bursts have binary-neutron-star (neutron-star--black-hole) progenitors we can use our 90% confidence rate upper limits to constrain the beaming angle of the gamma-ray burst to be greater than ${2.3^{+1.7}_{-1.1}}^{\circ}$ (${4.3^{+3.1}_{-1.9}}^{\circ}$).

preprint2015arXiv

An Overview of the 2014 ALMA Long Baseline Campaign

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

preprint2015arXiv

Design and RF measurements of a 5 GHz 500 kW window for the ITER LHCD system

CEA/IRFM is conducting R\&D efforts in order to validate the critical RF components of the 5 GHz ITER LHCD system, which is expected to transmit 20 MW of RF power to the plasma. Two 5 GHz 500 kW BeO pill-box type window prototypes have been manufactured in 2012 by the PMB Company, in close collaboration with CEA/IRFM. Both windows have been validated at low power, showing good agreement between measured and modeling, with a return loss better than 32 dB and an insertion loss below 0.05 dB. This paper reports on the window RF design and the low power measurements. The high power tests up to 500kW have been carried out in March 2013 in collaboration with NFRI. Results of these tests are also reported. In the current ITER LHCD design, 20 MW Continuous Wave (CW) of Radio-Frequency power at 5 GHz are expected to be generated and transmitted to the plasma. In order to separate the vacuum vessel pressure from the cryostat waveguide pressure, forty eight 5 GHz 500kW CW windows are to be assembled on the waveguides at the equatorial port flange. For nuclear safety reasons, forty eight additional windows could be located in the cryostat section, to separate and monitor the cryostat waveguide pressure from the exterior transmission line pressure. These windows are identified as being one of the main critical components for the ITER LHCD system since first ITER LHCD studies [1] [2] [3] or more recently [4] [5] , and clearly require an important R\&D effort. In this context and even if the LHCD system is not part of the construction baseline, the CEA/IRFM is conducting a R\&D effort in order to validate a design and the performances of these RF windows. In order to begin the assessment of this need, two 5 GHz 500 kW/5 s pill-box type windows prototypes have been manufactured in 2012 by the PMB Company in close collaboration with the CEA/IRFM [6]. The section 2 of this paper reports the RF and mechanical design of a 5 GHz window. Some features of the mechanical design and the experimental RF measurements at low power are reported in section 3. High power results, made in collaboration with NFRI, are detailed in section 4. The development of CW windows is discussed in the conclusion. 2-RF AND MECHANICAL DESIGN The proposed 5 GHz RF window is based on a pill-box design [2] , i.e. a ceramic brazed in portion of a circular waveguide, connected on either side to a rectangular waveguide section. Typical design rules of thumb of such device are circular section diameter about the same size of the diagonal of the rectangular waveguide (cf. FIGURE 1). Without taking into account the ceramic, the circular section length is approximately half a guided wavelength of the circular TE 11 mode, in order for the device to act as a half-wave transformer. Once optimized, taking into account the ceramic, matching is correct only for a narrow band of frequency and is very sensitive to the device dimensions and the ceramic relative permittivity. The heat losses in the ceramic, which have to be extracted by an active water cooling, depends on the inside electric field topology and of ceramic dielectric loss (loss tangent). Undesirable modes due to parasitic resonances can be excited in the ceramic volume, raising the electric field and

preprint2015arXiv

Epsilon-Near-Zero Al-Doped ZnO for Ultrafast Switching at Telecom Wavelengths: Outpacing the Traditional Amplitude-Bandwidth Trade-Off

Transparent conducting oxides have recently gained great attention as CMOS-compatible materials for applications in nanophotonics due to their low optical loss, metal-like behavior, versatile/tailorable optical properties, and established fabrication procedures. In particular, aluminum doped zinc oxide (AZO) is very attractive because its dielectric permittivity can be engineered over a broad range in the near infrared and infrared. However, despite all these beneficial features, the slow (> 100 ps) electron-hole recombination time typical of these compounds still represents a fundamental limitation impeding ultrafast optical modulation. Here we report the first epsilon-near-zero AZO thin films which simultaneously exhibit ultra-fast carrier dynamics (excitation and recombination time below 1 ps) and an outstanding reflectance modulation up to 40% for very low pump fluence levels (< 4 mJ/cm2) at the telecom wavelength of 1.3 μm. The unique properties of the demonstrated AZO thin films are the result of a low temperature fabrication procedure promoting oxygen vacancies and an ultra-high carrier concentration. As a proof-of-concept, an all-optical AZO-based plasmonic modulator achieving 3 dB modulation in 7.5 μm and operating at THz frequencies is numerically demonstrated. Our results overcome the traditional "modulation depth vs. speed" trade-off by at least an order of magnitude, placing AZO among the most promising compounds for tunable/switchable nanophotonics.

preprint2015arXiv

Measurement of absorption and charge exchange of $π^+$ on carbon

The combined cross section for absorption and charge exchange interactions of positively charged pions with carbon nuclei for the momentum range 200 MeV/c to 300 MeV/c have been measured with the DUET experiment at TRIUMF. The uncertainty is reduced by nearly half compared to previous experiments. This result will be a valuable input to existing models to constrain pion interactions with nuclei.

preprint2015arXiv

Measurement of the $ν_μ$ CCQE cross section on carbon with the ND280 detector at T2K

The Charged-Current Quasi-Elastic (CCQE) interaction, $ν_{l} + n \rightarrow l^{-} + p$, is the dominant CC process at $E_ν\sim 1$ GeV and contributes to the signal in accelerator-based long-baseline neutrino oscillation experiments operating at intermediate neutrino energies. This paper reports a measurement by the T2K experiment of the $ν_μ$ CCQE cross section on a carbon target with the off-axis detector based on the observed distribution of muon momentum ($p_μ$) and angle with respect to the incident neutrino beam ($θ_μ$). The flux-integrated CCQE cross section was measured to be $(0.83 \pm 0.12) \times 10^{-38}\textrm{ cm}^{2}$ in good agreement with NEUT MC value of ${0.88 \times 10^{-38}} \textrm{ cm}^{2}$. The energy dependence of the CCQE cross section is also reported. The axial mass, $M_A^{QE}$, of the dipole axial form factor was extracted assuming the Smith-Moniz CCQE model with a relativistic Fermi gas nuclear model. Using the absolute (shape-only) $p_μcosθ_μ$ distribution, the effective $M_A^{QE}$ parameter was measured to be ${1.26^{+0.21}_{-0.18} \textrm{ GeV}/c^{2}}$ (${1.43^{+0.28}_{-0.22} \textrm{ GeV}/c^{2}}$).

preprint2015arXiv

Measurement of the $ν_μ$ charged current quasi-elastic cross-section on carbon with the T2K on-axis neutrino beam

We report a measurement of the $ν_μ$ charged current quasi-elastic cross-sections on carbon in the T2K on-axis neutrino beam. The measured charged current quasi-elastic cross-sections on carbon at mean neutrino energies of 1.94 GeV and 0.93 GeV are $(11.95\pm 0.19(stat.)_{-1.47}^{+1.82} (syst.))\times 10^{-39}\mathrm{cm}^2/\mathrm{neutron}$ and $(10.64\pm 0.37(stat.)_{-1.65}^{+2.03} (syst.))\times 10^{-39}\mathrm{cm}^2/\mathrm{neutron}$, respectively. These results agree well with the predictions of neutrino interaction models. In addition, we investigated the effects of the nuclear model and the multi-nucleon interaction.

preprint2015arXiv

Measurement of the Electron Neutrino Charged-current Interaction Rate on Water with the T2K ND280 pi-zero Detector

This paper presents a measurement of the charged current interaction rate of the electron neutrino beam component of the beam above $1.5$~GeV using the large fiducial mass of the T2K $π^0$ detector. The predominant poriton of the $ν_e$ flux ($\sim$85 %) at these energies comes from kaon decays. The measured ratio of the observed beam interaction rate to the predicted rate in the detector with water targets filled is 0.89 $\pm$ 0.08 (stat.) $\pm$ 0.11 (sys.), and with the water targets emptied is 0.90 $\pm$ 0.09 (stat.) $\pm$ 0.13 (sys.). The ratio obtained for the interactions on water only from an event subtraction method is 0.87 $\pm$ 0.33 (stat.) $\pm$ 0.21 (sys.). This is the first measurement of the interaction rate of electron neutrinos on water, which is particularly of interest to experiments with water Cherenkov detectors.

preprint2015arXiv

Measurement of the muon neutrino inclusive charged-current cross section in the energy range of 1-3 GeV with the T2K INGRID detector

We report a measurement of the $ν_μ$-nucleus inclusive charged current cross section (=$σ^{cc}$) on iron using data from exposed to the J-PARC neutrino beam. The detector consists of 14 modules in total, which are spread over a range of off-axis angles from 0$^\circ$ to 1.1$^\circ$. The variation in the neutrino energy spectrum as a function of the off-axis angle, combined with event topology information, is used to calculate this cross section as a function of neutrino energy. The cross section is measured to be $σ^{cc}(1.1\text{ GeV}) = 1.10 \pm 0.15$ $(10^{-38}\text{cm}^2/\text{nucleon})$, $σ^{cc}(2.0\text{ GeV}) = 2.07 \pm 0.27$ $(10^{-38}\text{cm}^2/\text{nucleon})$, and $σ^{cc}(3.3\text{ GeV}) = 2.29 \pm 0.45$ $(10^{-38}\text{cm}^2/\text{nucleon})$, at energies of 1.1, 2.0, and 3.3 GeV, respectively. These results are consistent with the cross section calculated by the neutrino interaction generators currently used by T2K. More importantly, the method described here opens up a new way to determine the energy dependence of neutrino-nucleus cross sections.

preprint2015arXiv

Measurements of neutrino oscillation in appearance and disappearance channels by the T2K experiment with 6.6E20 protons on target

We report on measurements of neutrino oscillation using data from the T2K long-baseline neutrino experiment collected between 2010 and 2013. In an analysis of muon neutrino disappearance alone, we find the following estimates and 68% confidence intervals for the two possible mass hierarchies: Normal Hierarchy: $\sin^2θ_{23}=0.514^{+0.055}_{-0.056}$ and $Δm^2_{32}=(2.51\pm0.10)\times 10^{-3}$ eV$^2$/c$^4$ Inverted Hierarchy: $\sin^2θ_{23}=0.511\pm0.055$ and $Δm^2_{13}=(2.48\pm0.10)\times 10^{-3}$ eV$^2$/c$^4$ The analysis accounts for multi-nucleon mechanisms in neutrino interactions which were found to introduce negligible bias. We describe our first analyses that combine measurements of muon neutrino disappearance and electron neutrino appearance to estimate four oscillation parameters and the mass hierarchy. Frequentist and Bayesian intervals are presented for combinations of these parameters, with and without including recent reactor measurements. At 90% confidence level and including reactor measurements, we exclude the region: $δ_{CP}=[0.15,0.83]π$ for normal hierarchy and $δ_{CP}=[-0.08,1.09]π$ for inverted hierarchy. The T2K and reactor data weakly favor the normal hierarchy with a Bayes Factor of 2.2. The most probable values and 68% 1D credible intervals for the other oscillation parameters, when reactor data are included, are: $\sin^2θ_{23}=0.528^{+0.055}_{-0.038}$ and $|Δm^2_{32}|=(2.51\pm0.11)\times 10^{-3}$ eV$^2$/c$^4$.

preprint2015arXiv

Neutrino Oscillation Physics Potential of the T2K Experiment

The observation of the recent electron neutrino appearance in a muon neutrino beam and the high-precision measurement of the mixing angle $θ_{13}$ have led to a re-evaluation of the physics potential of the T2K long-baseline neutrino oscillation experiment. Sensitivities are explored for CP violation in neutrinos, non-maximal $\sin^22θ_{23}$, the octant of $θ_{23}$, and the mass hierarchy, in addition to the measurements of $δ_{CP}$, $\sin^2θ_{23}$, and $Δm^2_{32}$, for various combinations of $ν$-mode and \(\barν\)-mode data-taking. With an exposure of $7.8\times10^{21}$~protons-on-target, T2K can achieve 1-$σ$ resolution of 0.050(0.054) on $\sin^2θ_{23}$ and $0.040(0.045)\times10^{-3}~\rm{eV}^2$ on $Δm^2_{32}$ for 100\%(50\%) neutrino beam mode running assuming $\sin^2θ_{23}=0.5$ and $Δm^2_{32} = 2.4\times10^{-3}$ eV$^2$. T2K will have sensitivity to the CP-violating phase $δ_{\rm{CP}}$ at 90\% C.L. or better over a significant range. For example, if $\sin^22θ_{23}$ is maximal (i.e $θ_{23}$=$45^\circ$) the range is $-115^\circ<δ_{\rm{CP}}<-60^\circ$ for normal hierarchy and $+50^\circ<δ_{\rm{CP}}<+130^\circ$ for inverted hierarchy. When T2K data is combined with data from the NO$ν$A experiment, the region of oscillation parameter space where there is sensitivity to observe a non-zero $δ_{CP}$ is substantially increased compared to if each experiment is analyzed alone.

preprint2015arXiv

Planck 2015 results. I. Overview of products and scientific results

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

preprint2015arXiv

Remotely induced magnetism in a normal metal using a superconducting spin-valve

Superconducting spintronics has emerged in the last decade as a promising new field that seeks to open a new dimension for nanoelectronics by utilizing the internal spin structure of the superconducting Cooper pair as a new degree of freedom. Its basic building blocks are spin-triplet Cooper pairs with equally aligned spins, which are promoted by proximity of a conventional superconductor to a ferromagnetic material with inhomogeneous macroscopic magnetization. Using low-energy muon spin rotation experiments, we find an entirely unexpected novel effect: the appearance of a magnetization in a thin layer of a non-magnetic metal (gold), separated from a ferromagnetic double layer by a 50 nm thick superconducting layer of Nb. The effect can be controlled by either temperature or by using a magnetic field to control the state of the remote ferromagnetic elements and may act as a basic building block for a new generation of quantum interference devices based on the spin of a Cooper pair.

preprint2015arXiv

X-ray scattering study of pyrochlore iridates: crystal structure, electronic and magnetic excitations

We have investigated the structural, electronic, and magnetic properties of the pyrochlore iridates Eu2Ir2O7 and Pr2Ir2O7 using a combination of resonant elastic x-ray scattering, x-ray powder diffraction, and resonant inelastic x-ray scattering (RIXS). The structural parameters of Eu2Ir2O7 have been examined as a function of temperature and applied pressure, with a particular emphasis on regions of the phase diagram where electronic and magnetic phase transitions have been reported. We find no evidence of crystal symmetry change over the range of temperatures (~6 to 300 K) and pressures (~0.1 to 17 GPa) studied. We have also investigated the electronic and magnetic excitations in single crystal samples of Eu2Ir2O7 and Pr2Ir2O7 using high resolution Ir L3-edge RIXS. In spite of very different ground state properties, we find these materials exhibit qualitatively similar excitation spectra, with crystal field excitations at ~3-5 eV, spin-orbit excitations at ~0.5-1 eV, and broad low-lying excitations below ~0.15 eV. In Eu2Ir2O7 we observe highly damped magnetic excitations at ~45 meV, which display significant momentum dependence. We compare these results with recent dynamical structure factor calculations.

preprint2014arXiv

Charmed and light pseudoscalar meson decay constants from four-flavor lattice QCD with physical light quarks

We compute the leptonic decay constants $f_{D^+}$, $f_{D_s}$, and $f_{K^+}$, and the quark-mass ratios $m_c/m_s$ and $m_s/m_l$ in unquenched lattice QCD using the experimentally determined value of $f_{π^+}$ for normalization. We use the MILC highly improved staggered quark (HISQ) ensembles with four dynamical quark flavors---up, down, strange, and charm---and with both physical and unphysical values of the light sea-quark masses. The use of physical pions removes the need for a chiral extrapolation, thereby eliminating a significant source of uncertainty in previous calculations. Four different lattice spacings ranging from $a\approx 0.06$ fm to $0.15$ fm are included in the analysis to control the extrapolation to the continuum limit. Our primary results are $f_{D^+} = 212.6(0.4)({}^{+1.0}_{-1.2})\ \mathrm{MeV}$, $f_{D_s} = 249.0(0.3)({}^{+1.1}_{-1.5})\ \mathrm{MeV}$, and $f_{D_s}/f_{D^+} = 1.1712(10)({}^{+29}_{-32})$, where the errors are statistical and total systematic, respectively. The errors on our results for the charm decay constants and their ratio are approximately two to four times smaller than those of the most precise previous lattice calculations. We also obtain $f_{K^+}/f_{π^+} = 1.1956(10)({}^{+26}_{-18})$, updating our previous result, and determine the quark-mass ratios $m_s/m_l = 27.35(5)({}^{+10}_{-7})$ and $m_c/m_s = 11.747(19)({}^{+59}_{-43})$. When combined with experimental measurements of the decay rates, our results lead to precise determinations of the CKM matrix elements $|V_{us}| = 0.22487(51) (29)(20)(5)$, $|V_{cd}|=0.217(1) (5)(1)$ and $|V_{cs}|= 1.010(5)(18)(6)$, where the errors are from this calculation of the decay constants, the uncertainty in the experimental decay rates, structure-dependent electromagnetic corrections, and, in the case of $|V_{us}|$, the uncertainty in $|V_{ud}|$, respectively.

preprint2014arXiv

Charmed and light pseudoscalar meson decay constants from HISQ simulations

We compute the leptonic decay constants $f_{D^+}$, $f_{D_s}$, and $f_{K^+}$, and the quark-mass ratios $m_c/m_s$ and $m_s/m_l$ in unquenched lattice QCD. We use the MILC highly improved staggered quark (HISQ) ensembles with four dynamical quark flavors. Our primary results are $f_{D^+} = 212.6(0.4)({}^{+1.0}_{-1.2})\ \mathrm{MeV}$, $f_{D_s} = 249.0(0.3)({}^{+1.1}_{-1.5})\ \mathrm{MeV}$, and $f_{D_s}/f_{D^+} = 1.1712(10)({}^{+29}_{-32})$, where the errors are statistical and total systematic, respectively. We also obtain $f_{K^+}/f_{π^+} = 1.1956(10)({}^{+26}_{-18})$, updating our previous result, and determine the quark-mass ratios $m_s/m_l = 27.35(5)({}^{+10}_{-7})$ and $m_c/m_s = 11.747(19)({}^{+59}_{-43})$. When combined with experimental measurements of the decay rates, our results lead to precise determinations of the CKM matrix elements $|V_{us}| = 0.22487(51) (29)(20)(5)$, $|V_{cd}|=0.217(1) (5)(1)$ and $|V_{cs}|= 1.010(5)(18)(6)$, where the errors are from this calculation of the decay constants, the uncertainty in the experimental decay rates, structure-dependent electromagnetic corrections, and, in the case of $|V_{us}|$, the uncertainty in $|V_{ud}|$, respectively.

preprint2014arXiv

Controlled suppression of superconductivity by the generation of polarized Cooper pairs in spin valve structures

Transport measurements are presented on thin-film superconducting spin-valve systems, where the controlled non-collinear arrangement of two ferromagnetic Co layers can be used to influence the superconducting state of Nb. We observe a very clear oscillation of the superconducting transition temperature with the relative orientation of the two ferromagnetic layers. Our measurements allow us to distinguish between the competing influences of domain averaging, stray dipolar fields and the formation of superconducting spin triplets. Domain averaging is shown to lead to a weak enhancement of transition temperature for the anti-parallel configuration of exchange fields, while much larger changes are observed for other configurations, which can be attributed to drainage currents due to spin triplet formation.

preprint2014arXiv

Determination of $|V_{us}|$ from a lattice-QCD calculation of the $K\toπ\ellν$ semileptonic form factor with physical quark masses

We calculate the kaon semileptonic form factor $f_+(0)$ from lattice QCD, working, for the first time, at the physical light-quark masses. We use gauge configurations generated by the MILC collaboration with $N_f=2+1+1$ flavors of sea quarks, which incorporate the effects of dynamical charm quarks as well as those of up, down, and strange. We employ data at three lattice spacings to extrapolate to the continuum limit. Our result, $f_+(0) = 0.9704(32)$, where the error is the total statistical plus systematic uncertainty added in quadrature, is the most precise determination to date. Combining our result with the latest experimental measurements of $K$ semileptonic decays, one obtains the Cabibbo-Kobayashi-Maskawa matrix element $|V_{us}|=0.22290(74)(52)$, where the first error is from $f_+(0)$ and the second one is from experiment. In the first-row test of Cabibbo-Kobayashi-Maskawa unitarity, the error stemming from $|V_{us}|$ is now comparable to that from $|V_{ud}|$.

preprint2014arXiv

Double Fano resonances in a composite metamaterial possessing tripod plasmonic resonances

By embedding four-rod resonators inside double-split ring resonators superlattice, a planar composite metamaterial possessing tripod plasmonic resonances is fabricated. Double Fano resonances are observed where a common subradiant driven oscillator is coupled with two superradiant oscillators. As a classical analogue of four-level tripod atomic system, the transmission spectrum of the composite metamaterial exhibits a double Fano-based coherent effect. Transfer of absorbed power between two superradiant oscillators is controlled by manipulating two coupling strengths conjugated through the polarization angle of a normally incident electromagnetic wave.

preprint2014arXiv

Ferroelectric and Magnetic Domains in LuFe2O4 Observed by Scanning Probe Microscopy

LuFe2O4 is a multiferroic system which exhibits ferroelectricity, charge order, and ferrimagnetic order simultaneously below ~230 K. The ferroelectric domains of LuFe2O4 are imaged with both piezoresponse force microscopy (PFM) and electrostatic force microscopy (EFM), while the magnetic domains are characterized by magnetic force microscopy (MFM). Comparison of PFM and EFM results lead to a conclusion that the ferroelectricity is of electronic origin as opposed to the usual displacive one. Simultaneous characterization of ferroelectric and magnetic domains by EFM and MFM, respectively, on the same surface of LuFe2O4 reveals that both domains have irregular patterns of similar shape, but the length scales are quite different. The domain size is approximately 100 nm for the ferroelectric domains while the magnetic domain size is much larger and gets as large as 1μm. We also demonstrate that the origin of the formation of irregular domains in LuFe2O4 is not extrinsic but intrinsic.

preprint2014arXiv

Individual addressing of trapped $^{171}$Yb$^+$ ion qubits using a MEMS-based beam steering system

The ability to individually manipulate the increasing number of qubits is one of the many challenges towards scalable quantum information processing with trapped ions. Using micro-mirrors fabricated with micro-electromechanical systems (MEMS) technology, we focus laser beams on individual ions in a linear chain and steer the focal point in two dimensions. We demonstrate sequential single qubit gates on multiple $^{171}$Yb$^+$ qubits and characterize the gate performance using quantum state tomography. Our system features negligible crosstalk to neighboring ions ($< 3\times 10^{-4}$), and switching speed comparable to typical single qubit gate times ($<$ 2 $μ$s).

preprint2014arXiv

K semileptonic form factor with HISQ fermions at the physical point

We present results for the form factor $f_+^{K π}(0)$, needed to extract the CKM matrix element $|V_{us}|$ from experimental data on semileptonic $K$ decays, on the HISQ $N_f=2+1+1$ MILC configurations. The HISQ action is also used for the valence sector. The data set used for our final result includes three different values of the lattice spacing and data at the physical light quark masses. We discuss the error budget and how this calculation improves on our previous determination of $f_+^{K π}(0)$ on the asqtad $N_f=2+1$ MILC configurations.

preprint2014arXiv

Measurement of the inclusive $ν_μ$ charged current cross section on iron and hydrocarbon in the T2K on-axis neutrino beam

We report a measurement of the $ν_μ$ inclusive charged current cross sections on iron and hydrocarbon in the T2K on-axis neutrino beam. The measured inclusive charged current cross sections on iron and hydrocarbon averaged over the T2K on-axis flux with a mean neutrino energy of 1.51 GeV are $(1.444\pm0.002(stat.)_{-0.157}^{+0.189}(syst.))\times 10^{-38}\mathrm{cm}^2/\mathrm{nucleon}$, and $(1.379\pm0.009(stat.)_{-0.147}^{+0.178}(syst.))\times 10^{-38}\mathrm{cm}^2/\mathrm{nucleon}$, respectively, and their cross section ratio is $1.047\pm0.007(stat.)\pm0.035(syst.)$. These results agree well with the predictions of the neutrino interaction model, and thus we checked the correct treatment of the nuclear effect for iron and hydrocarbon targets in the model within the measurement precisions.

preprint2014arXiv

Measurement of the Inclusive Electron Neutrino Charged Current Cross Section on Carbon with the T2K Near Detector

The T2K off-axis near detector, ND280, is used to make the first differential cross-section measurements of electron neutrino charged current interactions at energies ~1 GeV as a function of electron momentum, electron scattering angle and four-momentum transfer of the interaction. The total flux-averaged $ν_e$ charged current cross-section on carbon is measured to be $1.11\pm0.09~(stat)\pm0.18~(syst)\times10^{-38} cm^2/nucleon$. The differential and total cross-section measurements agree with the predictions of two leading neutrino interaction generators, NEUT and GENIE. The NEUT prediction is $1.23\times10^{-38} cm^2/nucleon$ and the GENIE prediction is $1.08\times10^{-38} cm^2/nucleon$. The total $ν_e$ charged current cross-section result is also in agreement with data from the Gargamelle experiment.

preprint2014arXiv

Measurement of the intrinsic electron neutrino component in the T2K neutrino beam with the ND280 detector

The T2K experiment has reported the first observation of the appearance of electron neutrinos in a muon neutrino beam. The main and irreducible background to the appearance signal comes from the presence in the neutrino beam of a small intrinsic component of electron neutrinos originating from muon and kaon decays. In T2K, this component is expected to represent 1.2% of the total neutrino flux. A measurement of this component using the near detector (ND280), located 280 m from the target, is presented. The charged current interactions of electron neutrinos are selected by combining the particle identification capabilities of both the time projection chambers and electromagnetic calorimeters of ND280. The measured ratio between the observed electron neutrino beam component and the prediction is 1.01+-0.10 providing a direct confirmation of the neutrino fluxes and neutrino cross section modeling used for T2K neutrino oscillation analyses. Electron neutrinos coming from muons and kaons decay are also separately measured, resulting in a ratio with respect to the prediction of 0.68+-0.30 and 1.10+-0.14, respectively.

preprint2014arXiv

Measurement of the neutrino-oxygen neutral-current interaction cross section by observing nuclear deexcitation $γ$ rays

We report the first measurement of the neutrino-oxygen neutral-current quasielastic (NCQE) cross section. It is obtained by observing nuclear deexcitation $γ$-rays which follow neutrino-oxygen interactions at the Super-Kamiokande water Cherenkov detector. We use T2K data corresponding to $3.01 \times 10^{20}$ protons on target. By selecting only events during the T2K beam window and with well-reconstructed vertices in the fiducial volume, the large background rate from natural radioactivity is dramatically reduced. We observe 43 events in the $4-30$ MeV reconstructed energy window, compared with an expectation of 51.0, which includes an estimated 16.2 background events. The background is primarily nonquasielastic neutral-current interactions and has only 1.2 events from natural radioactivity. The flux-averaged NCQE cross section we measure is $1.55 \times 10^{-38}$ cm$^2$ with a 68\% confidence interval of $(1.22, 2.20) \times 10^{-38}$ cm$^2$ at a median neutrino energy of 630 MeV, compared with the theoretical prediction of $2.01 \times 10^{-38}$ cm$^2$.

preprint2014arXiv

Planck 2013 results. I. Overview of products and scientific results

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

preprint2014arXiv

Planck 2013 results. XXIII. Isotropy and statistics of the CMB

The two fundamental assumptions of the standard cosmological model - that the initial fluctuations are statistically isotropic and Gaussian - are rigorously tested using maps of the cosmic microwave background (CMB) anisotropy from the Planck satellite. Deviations from isotropy have been found and demonstrated to be robust against component separation algorithm, mask choice and frequency dependence. Many of these anomalies were previously observed in the WMAP data, and are now confirmed at similar levels of significance (about 3 sigma). However, we find little evidence for non-Gaussianity, with the exception of a few statistical signatures that seem to be associated with specific anomalies. In particular, we find that the quadrupole-octopole alignment is also connected to a low observed variance of the CMB signal. A power asymmetry is now found to persist to scales corresponding to about l=600, and can be described in the low-l regime by a phenomenological dipole modulation model. However, any primordial power asymmetry is strongly scale-dependent and does not extend to arbitrarily small angular scales. Finally, it is plausible that some of these features may be reflected in the angular power spectrum of the data, which shows a deficit of power on similar scales. Indeed, when the power spectra of two hemispheres defined by a preferred direction are considered separately, one shows evidence for a deficit in power, while its opposite contains oscillations between odd and even modes that may be related to the parity violation and phase correlations also detected in the data. Although these analyses represent a step forward in building an understanding of the anomalies, a satisfactory explanation based on physically motivated models is still lacking.

preprint2014arXiv

Precise Measurement of the Neutrino Mixing Parameter θ_{23} from Muon Neutrino Disappearance in an Off-axis Beam

New data from the T2K neutrino oscillation experiment produce the most precise measurement of the neutrino mixing parameter theta_{23}. Using an off-axis neutrino beam with a peak energy of 0.6 GeV and a data set corresponding to 6.57 x 10^{20} protons on target, T2K has fit the energy-dependent nu_mu oscillation probability to determine oscillation parameters. Marginalizing over the values of other oscillation parameters yields sin^2 (theta_{23}) = 0.514 +0.055/-0.056 (0.511 +- 0.055), assuming normal (inverted) mass hierarchy. The best-fit mass-squared splitting for normal hierarchy is Delta m^2_{32} = (2.51 +- 0.10) x 10^{-3} eV^2/c^4 (inverted hierarchy: Delta m^2_{13} = (2.48 +- 0.10) x 10^{-3} eV^2/c^4). Adding a model of multinucleon interactions that affect neutrino energy reconstruction is found to produce only small biases in neutrino oscillation parameter extraction at current levels of statistical uncertainty.

preprint2014arXiv

Search for short baseline $ν_e$ disappearance with the T2K near detector

The T2K experiment has performed a search for $ν_e$ disappearance due to sterile neutrinos using $5.9 \times 10^{20}$ protons on target for a baseline of $280 m$ in a neutrino beam peaked at about $500 MeV$. A sample of ν_e CC interactions in the off-axis near detector has been selected with a purity of 63\% and an efficiency of 26\%. The p-value for the null hypothesis is 0.085 and the excluded region at 95\% CL is approximately $sin^2 2 θ_{ee} > 0.3$ for $Δm^2_{eff} > 7 eV^2 / c^4$.

preprint2014arXiv

The $D_s$, $D^+$, $B_s$ and $B$ decay constants from $2+1$ flavor lattice QCD

We present a study of the $D$ and $B$ leptonic decay constants on the MILC $N_f=2+1$ asqtad gauge ensembles using asqtad-improved staggered light quarks and clover heavy quarks in the Fermilab interpretation. Our previous analysis \cite{Bazavov:2011aa} computed the decay constants at lattice spacings $a \approx 0.14, 0.11$ and $0.083$ fm. We have extended the simulations to finer $a \approx 0.058$ and $0.043$ fm lattice spacings, and have also increased statistics; this allows us to address many important sources of uncertainty. Technical advances include a two-step two-point fit procedure, better tuning of the heavy quark masses and a better determination of the axial-vector current matching. The present analysis remains blinded, so here we focus on the improvements and their predicted impact on the error budget compared to the prior analysis.

preprint2013arXiv

Charmed and strange pseudoscalar meson decay constants from HISQ simulations

We update our determinations of $f_{D^+}$, $f_{D_s}$, $f_K$, and quark mass ratios from simulations with four flavors of HISQ dynamical quarks. The availability of ensembles with light quarks near their physical mass means that we can extract physical results with only small corrections for valence- and sea-quark mass mistunings instead of a chiral extrapolation. The adjusted valence-quark masses and lattice spacings may be determined from an ensemble-by-ensemble analysis, and the results for the quark mass ratios then extrapolated to the continuum limit. Our central values of the charmed meson decay constants, however, come from an alternative analysis, which uses staggered chiral perturbation theory for the heavy-light mesons, and allows us to incorporate data at unphysical quark masses where statistical errors are often smaller. A jackknife analysis propagated through all of these steps takes account of the correlations among all the quantities used in the analysis. Systematic errors from the finite spatial size and EM effects are estimated by varying the parameters in the analysis, and systematic errors from the assumptions in the continuum extrapolation are estimated from the spread of values from different extrapolations.

preprint2013arXiv

Large Scale Modular Quantum Computer Architecture with Atomic Memory and Photonic Interconnects

The practical construction of scalable quantum computer hardware capable of executing non-trivial quantum algorithms will require the juxtaposition of different types of quantum systems. We analyze a modular ion trap quantum computer architecture with a hierarchy of interactions that can scale to very large numbers of qubits. Local entangling quantum gates between qubit memories within a single register are accomplished using natural interactions between the qubits, and entanglement between separate registers is completed via a probabilistic photonic interface between qubits in different registers, even over large distances. We show that this architecture can be made fault-tolerant, and demonstrate its viability for fault-tolerant execution of modest size quantum circuits.

preprint2013arXiv

Leptonic decay-constant ratio f_{K^+}/f_{pi^+} from lattice QCD with physical light quarks

A calculation of the ratio of leptonic decay constants f_{K^+}/f_{π^+} makes possible a precise determination of the ratio of CKM matrix elements |V_{us}|/|V_{ud}| in the Standard Model, and places a stringent constraint on the scale of new physics that would lead to deviations from unitarity in the first row of the CKM matrix. We compute f_{K^+}/f_{π^+} numerically in unquenched lattice QCD using gauge-field ensembles recently generated that include four flavors of dynamical quarks: up, down, strange, and charm. We analyze data at four lattice spacings a ~ 0.06, 0.09, 0.12, and 0.15 fm with simulated pion masses down to the physical value 135 MeV. We obtain f_{K^+}/f_{π^+} = 1.1947(26)(37), where the errors are statistical and total systematic, respectively. This is our first physics result from our N_f = 2+1+1 ensembles, and the first calculation of f_{K^+}/f_{π^+} from lattice-QCD simulations at the physical point. Our result is the most precise lattice-QCD determination of f_{K^+}/f_{π^+}, with an error comparable to the current world average. When combined with experimental measurements of the leptonic branching fractions, it leads to a precise determination of |V_{us}|/|V_{ud}| = 0.2309(9)(4) where the errors are theoretical and experimental, respectively.

preprint2013arXiv

Measurement of Neutrino Oscillation Parameters from Muon Neutrino Disappearance with an Off-axis Beam

The T2K collaboration reports a precision measurement of muon neutrino disappearance with an off-axis neutrino beam with a peak energy of 0.6 GeV. Near detector measurements are used to constrain the neutrino flux and cross section parameters. The Super-Kamiokande far detector, which is 295 km downstream of the neutrino production target, collected data corresponding to $3.01 \times 10^{20}$ protons on target. In the absence of neutrino oscillations, $205 \pm 17$ (syst.) events are expected to be detected and only 58 muon neutrino event candidates are observed. A fit to the neutrino rate and energy spectrum assuming three neutrino flavors, normal mass hierarchy and $θ_{23}\leq π/4$ yields a best-fit mixing angle $\sin^2(2θ_{23})=1.000$ and mass splitting $|Δm^2_{32}| =2.44 \times 10^{-3}$ eV$^2$/c$^4$. If $θ_{23}\geq π/4$ is assumed, the best-fit mixing angle changes to $\sin^2(2θ_{23})=0.999$ and the mass splitting remains unchanged.

preprint2013arXiv

Measurement of the Inclusive NuMu Charged Current Cross Section on Carbon in the Near Detector of the T2K Experiment

T2K has performed the first measurement of νμ inclusive charged current interactions on carbon at neutrino energies of ~1 GeV where the measurement is reported as a flux-averaged double differential cross section in muon momentum and angle. The flux is predicted by the beam Monte Carlo and external data, including the results from the NA61/SHINE experiment. The data used for this measurement were taken in 2010 and 2011, with a total of 10.8 x 10^{19} protons-on-target. The analysis is performed on 4485 inclusive charged current interaction candidates selected in the most upstream fine-grained scintillator detector of the near detector. The flux-averaged total cross section is <σ_CC>_ϕ=(6.91 +/- 0.13 (stat) +/- 0.84 (syst)) x10^{-39} cm^2/nucleon for a mean neutrino energy of 0.85 GeV.

preprint2013arXiv

Multi-Determinant Wave-functions in Quantum Monte Carlo

Quantum Monte Carlo (QMC) methods have received considerable attention over the last decades due to their great promise for providing a direct solution to the many-body Schrodinger equation in electronic systems. Thanks to their low scaling with number of particles, QMC methods present a compelling competitive alternative for the accurate study of large molecular systems and solid state calculations. In spite of such promise, the method has not permeated the quantum chemistry community broadly, mainly because of the fixed-node error, which can be large and whose control is difficult. In this Perspective, we present a systematic application of large scale multi-determinant expansions in QMC, and report on its impressive performance with first row dimers and the 55 molecules of the G1 test set. We demonstrate the potential of this strategy for systematically reducing the fixed-node error in the wave function and for achieving chemical accuracy in energy predictions. When compared to traditional quantum chemistry methods like MP2, CCSD(T), and various DFT approximations, the QMC results show a marked improvement over all of them. In fact, only the explicitly-correlated CCSD(T) method with a large basis set produces more accurate results. Further developments in trial wave functions and algorithmic improvements appear promising for rendering QMC as the benchmark standard in large electronic systems.

preprint2013arXiv

The T2K Neutrino Flux Prediction

The Tokai-to-Kamioka (T2K) experiment studies neutrino oscillations using an off-axis muon neutrino beam with a peak energy of about 0.6 GeV that originates at the J-PARC accelerator facility. Interactions of the neutrinos are observed at near detectors placed at 280 m from the production target and at the far detector -- Super-Kamiokande (SK) -- located 295 km away. The flux prediction is an essential part of the successful prediction of neutrino interaction rates at the T2K detectors and is an important input to T2K neutrino oscillation and cross section measurements. A FLUKA and GEANT3 based simulation models the physical processes involved in the neutrino production, from the interaction of primary beam protons in the T2K target, to the decay of hadrons and muons that produce neutrinos. The simulation uses proton beam monitor measurements as inputs. The modeling of hadronic interactions is re-weighted using thin target hadron production data, including recent charged pion and kaon measurements from the NA61/SHINE experiment. For the first T2K analyses the uncertainties on the flux prediction are evaluated to be below 15% near the flux peak. The uncertainty on the ratio of the flux predictions at the far and near detectors is less than 2% near the flux peak.

preprint2012arXiv

Can residuals of the Solar system foreground explain low multipole anomalies of the CMB ?

The low multipole anomalies of the Cosmic Microwave Background has received much attention during the last few years. It is still not ascertained whether these anomalies are indeed primordial or the result of systematics or foregrounds. An example of a foreground, which could generate some non-Gaussian and statistically anisotropic features at low multipole range, is the very symmetric Kuiper Belt in the outer solar system. In this paper, expanding upon the methods presented by Maris et al. (2011), we investigate the contributions from the Kuiper Belt objects (KBO) to the WMAP ILC 7 map, whereby we can minimize the contrast in power between even and odd multipoles in the CMB, discussed discussed by Kim & Naselsky (2010). We submit our KBO de-correlated CMB signal to several tests, to analyze its validity, and find that incorporation of the KBO emission can decrease the quadrupole-octupole alignment and parity asymmetry problems, provided that the KBO signals has a non-cosmological dipole modulation, associated with the statistical anisotropy of the ILC 7 map. Additionally, we show that the amplitude of the dipole modulation, within a 2 sigma interval, is in agreement with the corresponding amplitudes, discussed by Lew (2008).

preprint2012arXiv

Faraday Rotation as a diagnostic of Galactic foreground contamination of CMB maps

The contribution from the residuals of the foreground can have a significant impact on the temperature maps of the Cosmic Microwave Background (CMB). Mostly, the focus has been on the galactic plane, when foreground cleaning has taken place. However, in this paper, we will investigate the possible foreground contamination, from sources outside the galactic plane in the CMB maps. We will analyze the correlation between the Faraday rotation map and the CMB temperature map. The Faraday rotation map is dependent on the galactic magnetic field, as well as the thermal electron density, and both may contribute to the CMB temperature. We find that the standard deviation for the mean cross correlation deviate from that of simulations at the 99.9% level. Additionally, a comparison between the CMB temperature extrema and the extremum points of the Faraday rotation is also performed, showing a general overlap between the two. Also we find that the CMB Cold Spot is located at an area of strong negative cross correlation, meaning that it may be explained by a galactic origin. Further, we investigate nearby supernova remnants in the galaxy, traced by the galactic radio loops. These super nova remnants are located at high and low galactic latitude, and thus well outside the galactic plane. We find some correlation between the Faraday Rotation and the CMB temperature, at select radio loops. This indicate, that the galactic foregrounds may affect the CMB, at high galactic latitudes

preprint2012arXiv

First Muon-Neutrino Disappearance Study with an Off-Axis Beam

We report a measurement of muon-neutrino disappearance in the T2K experiment. The 295-km muon-neutrino beam from Tokai to Kamioka is the first implementation of the off-axis technique in a long-baseline neutrino oscillation experiment. With data corresponding to 1.43 10**20 protons on target, we observe 31 fully-contained single muon-like ring events in Super-Kamiokande, compared with an expectation of 104 +- 14 (syst) events without neutrino oscillations. The best-fit point for two-flavor nu_mu -> nu_tau oscillations is sin**2(2 theta_23) = 0.98 and |Δm**2_32| = 2.65 10**-3 eV**2. The boundary of the 90 % confidence region includes the points (sin**2(2 theta_23),|Δm**2_32|) = (1.0, 3.1 10**-3 eV**2), (0.84, 2.65 10**-3 eV**2) and (1.0, 2.2 10**-3 eV**2).

preprint2012arXiv

Is the CMB asymmetry due to the kinematic dipole?

Parity violation found in the Cosmic Microwave Background (CMB) radiation is a crucial clue for the non-standard cosmological model or the possible contamination of various foreground residuals and/or calibration of the CMB data sets. In this paper, we study the directional properties of the CMB parity asymmetry by excluding the $m=0$ modes in the definition of parity parameters. We find that the preferred directions of the parity parameters coincide with the CMB kinematic dipole, which implies that the CMB parity asymmetry may be connected with the possible contamination of the residual dipole component. We also find that such tendency is not only localized at $l=2,3$, but in the extended multipole ranges up to $l\sim 22$.

preprint2012arXiv

Kaon semileptonic decay form factors with HISQ valence quarks

We report on the status of our kaon semileptonic form factor calculations using the highly-improved staggered quark (HISQ) formulation to simulate the valence fermions. We present results for the form factor f_+^{K π}(0) on the asqtad N_f=2+1 MILC configurations, discuss the chiral-continuum extrapolation, and give a preliminary estimate of the total error. We also present a more preliminary set of results for the same form factor but with the sea quarks also simulated with the HISQ action; these results include data at the physical light quark masses. The improvements that we expect to achieve with the use of the HISQ configurations and simulations at the physical quark masses are briefly discussed.

preprint2012arXiv

Morphology of High-Multiplicity Events in Heavy Ion Collisions

We discuss opportunities that may arise from subjecting high-multiplicity events in relativistic heavy ion collisions to an analysis similar to the one used in cosmology for the study of fluctuations of the Cosmic Microwave Background (CMB). To this end, we discuss examples of how pertinent features of heavy ion collisions including global characteristics, signatures of collective flow and event-wise fluctuations are visually represented in a Mollweide projection commonly used in CMB analysis, and how they are statistically analyzed in an expansion over spherical harmonic functions. If applied to the characterization of purely azimuthal dependent phenomena such as collective flow, the expansion coefficients of spherical harmonics are seen to contain redundancies compared to the set of harmonic flow coefficients commonly used in heavy ion collisions. Our exploratory study indicates, however, that these redundancies may offer novel opportunities for a detailed characterization of those event-wise fluctuations that remain after subtraction of the dominant collective flow signatures. By construction, the proposed approach allows also for the characterization of more complex collective phenomena like higher-order flow and other sources of fluctuations, and it may be extended to the characterization of phenomena of non-collective origin such as jets.

preprint2011arXiv

Indication of Electron Neutrino Appearance from an Accelerator-produced Off-axis Muon Neutrino Beam

The T2K experiment observes indications of $ν_μ\rightarrow ν_e$ appearance in data accumulated with $1.43\times10^{20}$ protons on target. Six events pass all selection criteria at the far detector. In a three-flavor neutrino oscillation scenario with $|Δm_{23}^2|=2.4\times10^{-3}$ eV$^2$, $\sin^2 2θ_{23}=1$ and $\sin^2 2θ_{13}=0$, the expected number of such events is 1.5$\pm$0.3(syst.). Under this hypothesis, the probability to observe six or more candidate events is 7$\times10^{-3}$, equivalent to 2.5$σ$ significance. At 90% C.L., the data are consistent with 0.03(0.04)$<\sin^2 2θ_{13}<$ 0.28(0.34) for $δ_{\rm CP}=0$ and a normal (inverted) hierarchy.

preprint2011arXiv

Measurements of the T2K neutrino beam properties using the INGRID on-axis near detector

Precise measurement of neutrino beam direction and intensity was achieved based on a new concept with modularized neutrino detectors. INGRID (Interactive Neutrino GRID) is an on-axis near detector for the T2K long baseline neutrino oscillation experiment. INGRID consists of 16 identical modules arranged in horizontal and vertical arrays around the beam center. The module has a sandwich structure of iron target plates and scintillator trackers. INGRID directly monitors the muon neutrino beam profile center and intensity using the number of observed neutrino events in each module. The neutrino beam direction is measured with accuracy better than 0.4 mrad from the measured profile center. The normalized event rate is measured with 4% precision.

preprint2011arXiv

Pearson's random walk in the space of the CMB phases: evidence for parity asymmetry

The temperature fluctuations of the Cosmic Microwave Background (CMB) are supposed to be distributed randomly in both magnitude and phase, following to the simplest model of inflation. In this paper, we look at the odd and even multipoles of the spherical harmonic decomposition of the CMB, and the different characteristics of these, giving rise to a parity asymmetry. We compare the even and odd multipoles in the CMB power spectrum, and also the even and odd mean angles. We find for the multipoles of the power spectrum, that there is power excess in odd multipoles, compared to even ones, meaning that we have a parity asymmetry. Further, for the phases, we present a random walk for the mean angles, and find a significant separation for even/odd mean angles, especially so for galactic coordinates. This is further tested and confirmed with a directional parity test, comparing the parity asymmetry in galactic and ecliptic coordinates.

preprint2011arXiv

Properties of light pseudoscalars from lattice QCD with HISQ ensembles

We fit lattice-QCD data for light-pseudoscalar masses and decay constants, from HISQ configurations generated by MILC, to SU(3) staggered chiral perturbation theory. At present such fits have rather high values of chi^2/d.o.f., possibly due to the lack of ensembles with lighter-than-physical sea strange-quark masses. We propose solutions to this problem for future work. We also perform simple linear interpolations near the physical point on two ensembles with different lattice spacings, and obtain the preliminary result (f_K / f_pi)^phys = 1.1872(41) in the continuum limit.

preprint2011arXiv

Symmetry of the CMB sky as a new test of its statistical isotropy. Non Cosmological Octupole?

In this article we propose a novel test for statistical anisotropy of the CMB. The test is based on the fact, that the Galactic foregrounds have a remarcably strong symmetry with respect to their antipodal points and with respect to the Galactic plane, while the cosmological signal should not be symmetric or asymmetric under these transitions. We have applied the test for the octupole component of the WMAP ILC 7 map, by looking at a_3,1 and a_3,3, and their ratio to a_3,2 both for real and imaginary values. We find abnormal symmetry of the octupole component at the level of 0.58%, compared to Monte Carlo simulations. By using the analysis of the phases of the octupole we found remarkably strong cross-correlations between the phases of kinematic dipole and ILC 7 octupole, in full agreement with previous results. We further test the multipole range 2<l<100, by investigating the ratio between the l+m=even and l+m=odd parts of power spectra. We compare the results to simulations of a Gaussian random sky, and find significant departure from the statistically isotropic and homogeneous case, for a very broad range of multipoles. We found that for the most prominent peaks of our estimator, the phases of the corresponding harmonics are coherent with phases of the octupole. We believe, our test would be very useful for detections of various types of residuals of the foreground and systematic effects at very a broad range of multipoles 2 < l < 1500-3000 for the forthcoming PLANCK CMB map, before any conclusions about primordial non-Gaussianity and statistical anisotropy of the CMB.

preprint2009arXiv

Algorithmic comparisons of decaying, isothermal, supersonic turbulence

Contradicting results have been reported in the literature with respect to the performance of the numerical techniques employed for the study of supersonic turbulence. We aim at characterising the performance of different particle-based and grid-based techniques on the modelling of decaying supersonic turbulence. Four different grid codes (ENZO, FLASH, TVD, ZEUS) and three different SPH codes (GADGET, PHANTOM, VINE) are compared. We additionally analysed two calculations denoted as PHANTOM A and PHANTOM B using two different implementations of artificial viscosity. Our analysis indicates that grid codes tend to be less dissipative than SPH codes, though details of the techniques used can make large differences in both cases. For example, the Morris & Monaghan viscosity implementation for SPH results in less dissipation (PHANTOM B and VINE versus GADGET and PHANTOM A). For grid codes, using a smaller diffusion parameter leads to less dissipation, but results in a larger bottleneck effect (our ENZO versus FLASH runs). As a general result, we find that by using a similar number of resolution elements N for each spatial direction means that all codes (both grid-based and particle-based) show encouraging similarity of all statistical quantities for isotropic supersonic turbulence on spatial scales k<N/32 (all scales resolved by more than 32 grid cells), while scales smaller than that are significantly affected by the specific implementation of the algorithm for solving the equations of hydrodynamics. At comparable numerical resolution, the SPH runs were on average about ten times more computationally intensive than the grid runs, although with variations of up to a factor of ten between the different SPH runs and between the different grid runs. (abridged)

preprint2008arXiv

Localized electronic states induced by defects and possible origin of ferroelectricity in strontium titanate thin films

Several defect configurations including oxygen vacancies have been investigated as possible origins of the reported room-temperature ferroelectricity of strontium titanate (STO) thin films [Appl. Phys. Letts. 91, 042908 (2007)]. First-principles calculations revealed that the Sr-O-O vacancy complexes create deep localized states in the band gap of SrTiO3 without affecting its insulating property. These results are in agreement with electronic structural changes determined from optical transmission and X-ray absorption measurements. This work opens the way to exploiting oxygen vacancies and their complexes as a source of ferroelectricity in perovskite oxide thin films, including STO.

preprint2002arXiv

Charged Higgs boson in the next-to-minimal supersymmetric standard model with explicit CP violation

The phenomenology of the explicit CP violation in the Higgs sector of the next-to-minimal supersymmetric standard model (NMSSM) is investigated, with emphasis on the charged Higgs boson. The radiative corrections due to both quarks and scalar-quarks of the third generation are taken into account, and the negative result of the search for the Higgs bosons at CERN LEP2, with the discovery limit of 0.1 pb, is imposed as a constraint. It is found that there are parameter regions of the NMSSM where the lightest neutral Higgs boson may even be massless, without being detected at LEP2. This implies that the LEP2 data do not contradict the existence of a massless neutral Higgs boson in the NMSSM. For the charged Higgs boson, the radiative corrections to its mass may be negative in some parameter regions of the NMSSM. The phenomenological lower bound on the radiatively corrected mass of the charged Higgs boson is increased as the CP violation becomes maximal, i.e., as the CP violating phase becomes $π/2$. At the maximal CP violation, its lower bound is about 110 GeV for 5 $\leqslant \tan β\leqslant$ 40. The vacuum expectation value (VEV) of the neutral Higgs singlet is shown to be no smaller than 16 GeV for any parameter values of the NMSSM with explicit CP violation. This value of the lower limit is found to increase up to about 45 GeV as the ratio ($\tan β$) of the VEVs of the two Higgs doublets decreases to smaller values ($\sim$ 2). The discovery limit of the Higgs boson search at LEP2 is found to cover about a half of the kinematically allowed part of the whole parameter space of the NMSSM, and the portion is roughly stable against the CP violating phase.

preprint1996arXiv

Carbon superatom thin films

Assembling clusters on surfaces has emerged as a novel way to grow thin films with targeted properties. In particular, it has been proposed from experimental findings that fullerenes deposited on surfaces could give rise to thin films retaining the bonding properties of the incident clusters. However the microscopic structure of such films is still unclear. By performing quantum molecular dynamics simulations, we show that C_28 fullerenes can be deposited on a surface to form a thin film of nearly defect free molecules, which act as carbon superatoms. Our findings help clarify the structure of disordered small fullerene films and also support the recently proposed hyperdiamond model for solid C_28.