Intermediate energy single-neutron removal from $^{31}$Mg has been employed to investigate the transition into the N=20 island of inversion. Levels up to 5~MeV excitation energy in $^{30}$Mg were populated and spin-parity assignments were inferred from the corresponding longitudinal momentum distributions and $\gamma$-ray decay scheme. Comparison with eikonal-model calculations also permitted spectroscopic factors to be deduced. Surprisingly, the 0$^{+}_{2}$ level in $^{30}$Mg was found to have a strength much weaker than expected in the conventional picture of a predominantly $2p - 2h$ intruder configuration having a large overlap with the deformed $^{31}$Mg ground state. In addition, negative parity levels were identified for the first time in $^{30}$Mg, one of which is located at low excitation energy. The results are discussed in the light of shell-model calculations employing two newly developed approaches with markedly different descriptions of the structure of $^{30}$Mg. It is concluded that the cross-shell effects in the region of the island of inversion at Z=12 are considerably more complex than previously thought and that $np - nh$ configurations play a major role in the structure of $^{30}$Mg.

We report results from the first search for sterile neutrinos mixing with active neutrinos through a reduction in the rate of neutral-current interactions over a baseline of 810\,km between the NOvA detectors. Analyzing a 14-kton detector equivalent exposure of 6.05$\times$10$^{20}$ protons-on-target in the NuMI beam at Fermilab, we observe 95 neutral-current candidates at the Far Detector compared with $83.5 \pm 9.7 \mbox{(stat.)} \pm 9.4 \mbox{(syst.)}$ events predicted assuming mixing only occurs between active neutrino species. No evidence for $\nu_{\mu} \rightarrow \nu_{s}$ transitions is found. Interpreting these results within a 3+1 model, we place constraints on the mixing angles $\theta_{24}<20.8^{\circ}$ and $\theta_{34}<31.2^{\circ}$ at the 90% C.L. for $0.05~eV^2\leq \Delta m^2_{41}\leq 0.5~eV^2$, the range of mass splittings that produce no significant oscillations over the Near Detector baseline.

Results are reported from an improved measurement of $\nu_\mu \rightarrow \nu_e$ transitions by the NOvA experiment. Using an exposure equivalent to $6.05\times10^{20}$ protons-on-target 33 $\nu_e$ candidates were observed with a background of $8.2\pm0.8$ (syst.). Combined with the latest NOvA $\nu_\mu$ disappearance data and external constraints from reactor experiments on $\sin^22\theta_{13}$, the hypothesis of inverted mass hierarchy with $\theta_{23}$ in the lower octant is disfavored at greater than $93\%$ C.L. for all values of $\delta_{CP}$.

This Letter reports new results on muon neutrino disappearance from NOvA, using a 14 kton detector equivalent exposure of $6.05\times10^{20}$ protons-on-target from the NuMI beam at the Fermi National Accelerator Laboratory. The measurement probes the muon-tau symmetry hypothesis that requires maximal mixing ($\theta_{23} = \pi/4$). Assuming the normal mass hierarchy, we find $\Delta m^2 = (2.67 \pm 0.11)\times 10^{-3}$ eV$^2$ and $\sin^2 \theta_{23}$ at the two statistically degenerate values $0.404^{+0.030}_{-0.022}$ and $0.624^{+0.022}_{-0.030}$, both at the 68% confidence level. Our data disfavor the maximal mixing scenario with 2.6 $\sigma$ significance.

The Landau collision integral is an accurate model for the small-angle dominated Coulomb collisions in fusion plasmas. We investigate a high order accurate, fully conservative, finite element discretization of the nonlinear multi-species Landau integral with adaptive mesh refinement using the PETSc library (www.mcs.anl.gov/petsc). We develop algorithms and techniques to efficiently utilize emerging architectures with an approach that minimizes memory usage and movement and is suitable for vector processing. The Landau collision integral is vectorized with Intel AVX-512 intrinsics and the solver sustains as much as 22% of the theoretical peak flop rate of the Second Generation Intel Xeon Phi, Knights Landing, processor.

We have imaged GM Aur with HST, detected its disk in scattered light at 1400A and 1650A, and compared these with observations at 3300A, 5550A, 1.1 microns, and 1.6 microns. The scattered light increases at shorter wavelengths. The radial surface brightness profile at 3300A shows no evidence of the 24AU radius cavity that has been previously observed in sub-mm observations. Comparison with dust grain opacity models indicates the surface of the entire disk is populated with sub-micron grains. We have compiled an SED from 0.1 microns to 1 mm, and used it to constrain a model of the star+disk system that includes the sub-mm cavity using the Monte Carlo Radiative Transfer code by Barbara Whitney. The best-fit model image indicates that the cavity should be detectable in the F330W bandpass if the cavity has been cleared of both large and small dust grains, but we do not detect it. The lack of an observed cavity can be explained by the presence of sub-microns grains interior to the sub-mm cavity wall. We suggest one explanation for this which could be due to a planet of mass <9 Jupiter masses interior to 24 AU. A unique cylindrical structure is detected in the FUV data from the Advanced Camera for Surveys/Solar Blind Channel. It is aligned along the system semi-minor axis, but does not resemble an accretion-driven jet. The structure is limb-brightened and extends 190 +/- 35 AU above the disk midplane. The inner radius of the limb-brightening is 40 +/- 10 AU, just beyond the sub-millimeter cavity wall.

The neutron-unbound nucleus 13Be was populated with a nucleon-exchange reaction from a 71 MeV/u secondary 13B beam. The decay energy spectrum was reconstructed using invariant mass spectroscopy based on 12Be fragments in coincidence with neutrons. The data could be described with an s-wave resonance at E = 0.73(9) MeV with a width of Gamma = 1.98(34) MeV and a d-wave resonance at E = 2.56(13) MeV with a width of Gamma = 2.29(73) MeV. The observed spectral shape is consistent with previous one-proton removal reaction measurements from 14B.

We investigate ten years of user support emails in the large-scale solver library PETSc in order to identify changes in user requests. For this purpose we assign each email thread to one or several categories describing the type of support request. We find that despite several changes in hardware architecture as well programming models, the relative share of emails for the individual categories does not show a notable change over time. This is particularly remarkable as the total communication volume has increased four-fold in the considered time frame, indicating a considerable growth of the user base. Our data also demonstrates that user support cannot be substituted with what is often referred to as 'better documentation' and that the involvement of core developers in user support is essential.

Fingering convection (or thermohaline convection) is a weak yet important kind of mixing that occurs in stably-stratified stellar radiation zones in the presence of an inverse mean-molecular-weight gradient. Brown et al. (2013) recently proposed a new model for mixing by fingering convection, which contains no free parameter, and was found to fit the results of direct numerical simulations in almost all cases. Notably, however, they found that mixing was substantially enhanced above their predicted values in the few cases where large-scale gravity waves, followed by thermo-compositional layering, grew spontaneously from the fingering convection. This effect is well-known in the oceanographic context, and is attributed to the excitation of the so-called "collective instability". In this work, we build on the results of Brown et al. (2013) and of Traxler et al. (2011b) to determine the conditions under which the collective instability may be expected. We find that it is only relevant in stellar regions which have a relatively large Prandtl number (the ratio of the kinematic viscosity to the thermal diffusivity), $O(10^{-3})$ or larger. This implies that the collective instability cannot occur in main sequence stars, where the Prandtl number is always much smaller than this (except in the outer layers of surface convection zones where fingering is irrelevant anyway). It could in principle be excited in regions of high electron degeneracy, during He core flash, or in the interiors of white dwarfs. We discuss the implications of our findings for these objects, both from a theoretical and from an observational point of view.

Decision trees are a popular technique in statistical data classification. They recursively partition the feature space into disjoint sub-regions until each sub-region becomes homogeneous with respect to a particular class. The basic Classification and Regression Tree (CART) algorithm partitions the feature space using axis parallel splits. When the true decision boundaries are not aligned with the feature axes, this approach can produce a complicated boundary structure. Oblique decision trees use oblique decision boundaries to potentially simplify the boundary structure. The major limitation of this approach is that the tree induction algorithm is computationally expensive. In this article we present a new decision tree algorithm, called HHCART. The method utilizes a series of Householder matrices to reflect the training data at each node during the tree construction. Each reflection is based on the directions of the eigenvectors from each classes' covariance matrix. Considering axis parallel splits in the reflected training data provides an efficient way of finding oblique splits in the unreflected training data. Experimental results show that the accuracy and size of the HHCART trees are comparable with some benchmark methods in the literature. The appealing feature of HHCART is that it can handle both qualitative and quantitative features in the same oblique split.

Context: Solar flare hard X-rays (HXRs) are thought to be produced by nonthermal coronal electrons stopping in the chromosphere, or remaining trapped in the corona. The collisional thick target model (CTTM) predicts that sources produced by harder power-law injection spectra should appear further down the legs or footpoints of a flare loop. Therefore, hardening of the injected power-law electron spectrum during flare onset should be concurrent with a descending hard X-ray source. Aims: To test this implication of the CTTM by comparing its predicted HXR source locations with those derived from observations of a solar flare which exhibits a nonthermally-dominated spectrum before the peak in HXRs, known as an early impulsive event. Methods: HXR images and spectra of an early impulsive C-class flare were obtained using the Ramaty High-Energy Solar Spectroscopic Imager (RHESSI). Images were reconstructed to produce HXR source height evolutions for three energy bands. Spatially-integrated spectral analysis was performed to isolate nonthermal emission, and to determine the power-law index of the electron injection spectrum. The observed height-time evolutions were then fit with CTTM-based simulated heights for each energy. Results: A good match between model and observed source heights was reached, requiring a density model that agreed well with previous studies of flare loop densities. Conclusions: The CTTM has been used to produce a descent of model HXR source heights that compares well with observations of this event. Based on this interpretation, downward motion of nonthermal sources should indeed occur in any flare where there is spectral hardening in the electron distribution during a flare. However, this would often be masked by thermal emission associated with flare plasma pre-heating.

Proton removal reactions from a secondary 22N beam were utilized to populate unbound states in neutron-rich carbon isotopes. Neutrons were measured with the Modular Neutron Array (MoNA) in coincidence with carbon fragments. A resonance with a decay energy of 76(14) keV was observed in the system 18C+n corresponding to a state in 19C at an excitation energy of 653(95)keV. This resonance could correspond to the first 5/2+ state which was recently speculated to be unbound in order to describe 1n and 2n removal cross section measurements from 20C.

A search for the neutron-unbound nucleus $^{21}$C was performed via the single proton removal reaction from a beam of 22 N at 68 MeV/u. Neutrons were detected with the Modular Neutron Array (MoNA) in coincidence with $^{20}$C fragments. No evidence for a low-lying state was found, and the reconstructed $^{20}$C+n decay energy spectrum could be described with an s-wave line shape with a scattering length limit of |as| < 2.8 fm, consistent with shell model predictions. A comparison with a renormalized zero-range three-body model suggests that $^{22}$C is bound by less than 70 keV.

The top quark is the heaviest known elementary particle, with a mass about 40 times larger than the mass of its isospin partner, the bottom quark. It decays almost 100% of the time to a $W$ boson and a bottom quark. Using top-antitop pairs at the Tevatron proton-antiproton collider, the CDF and {\dzero} collaborations have measured the top quark's mass in different final states for integrated luminosities of up to 5.8 fb$^{-1}$. This paper reports on a combination of these measurements that results in a more precise value of the mass than any individual decay channel can provide. It describes the treatment of the systematic uncertainties and their correlations. The mass value determined is $173.18 \pm 0.56 \thinspace ({\rm stat}) \pm 0.75 \thinspace ({\rm syst})$ GeV or $173.18 \pm 0.94$ GeV, which has a precision of $\pm 0.54%$, making this the most precise determination of the top quark mass.

We study perfect state transfer of quantum walks on signed graphs. Our aim is to show that negative edges are useful for perfect state transfer. Specific results we prove include: (1) The signed join of a negative 2-clique with any positive (n,3)-regular graph has perfect state transfer even if the unsigned join does not. Curiously, the perfect state transfer time improves as n increases. (2) A signed complete graph has perfect state transfer if its positive subgraph is a regular graph with perfect state transfer and its negative subgraph is periodic. This shows that signing is useful for creating perfect state transfer since no complete graph (except for the 2-clique) has perfect state transfer. (3) The double-cover of a signed graph has perfect state transfer if the positive subgraph has perfect state transfer and the negative subgraph is periodic. Here, signing is useful for constructing unsigned graphs with perfect state transfer. Furthermore, we study perfect state transfer on a family of signed graphs called the exterior powers which is derived from a many-fermion quantum walk on graphs.

We present a measurement of the helicity of the W boson produced in top quark decays using ttbar decays in the l+jets and dilepton final states selected from a sample of 5.4 fb^-1 of collisions recorded using the D0 detector at the Fermilab Tevatron ppbar collider. We measure the fractions of longitudinal and right-handed W bosons to be f_0 = 0.669 +- 0.102 [ +- 0.078 (stat.) +- 0.065 (syst.)] and f_+ = 0.023 +- 0.053 [+- 0.041 (stat.) +- 0.034 (syst.)], respectively. This result is consistent at the 98% level with the standard model. A measurement with f_0 fixed to the value from the standard model yields f_+ = 0.010 +- 0.037 [+- 0.022 (stat.) +- 0.030 (syst.) ].

We present a measurement of the elastic differential cross section $d\sigma(p\bar{p}\rightarrow p\bar{p})/dt$ as a function of the four-momentum-transfer squared t. The data sample corresponds to an integrated luminosity of $\approx 31 nb^{-1}$ collected with the D0 detector using dedicated Tevatron $p\bar{p} $ Collider operating conditions at sqrt(s) = 1.96 TeV and covers the range $0.26 <|t|< 1.2 GeV^2$. For $|t|<0.6 GeV^2$, d\sigma/dt is described by an exponential function of the form $Ae^{-b|t|}$ with a slope parameter $ b = 16.86 \pm 0.10(stat) \pm 0.20(syst) GeV^{-2}$. A change in slope is observed at $|t| \approx 0.6 GeV^2$, followed by a more gradual |t| dependence with increasing values of |t|.

We report the combination of recent measurements of the helicity of the W boson from top quark decay by the CDF and D0 collaborations, based on data samples corresponding to integrated luminosities of 2.7 - 5.4 fb^-1 of ppbar collisions collected during Run II of the Fermilab Tevatron Collider. Combining measurements that simultaneously determine the fractions of W bosons with longitudinal (f0) and right-handed (f+) helicities, we find f0 = 0.722 \pm 0.081 [\pm 0.062 (stat.) \pm 0.052 (syst.)] and f+ = -0.033 \pm 0.046 [\pm 0.034 (stat.) \pm 0.031 (syst.)]. Combining measurements where one of the helicity fractions is fixed to the value expected in the standard model, we find f0 = 0.682 \pm 0.057 [\pm 0.035 (stat.) \pm 0.046 (syst.)] and f+ = -0.015\pm0.035 [\pm 0.018 (stat.) \pm 0.030 (syst.)]. The results are consistent with standard model expectations.

We present a search for the pair production of scalar top quarks ($\tilde{t}_{1}$), the lightest supersymmetric partners of the top quarks, in $p\bar{p}$ collisions at a center-of-mass energy of 1.96 TeV, using data corresponding to an integrated luminosity of {7.3 $fb^{-1}$} collected with the \dzero experiment at the Fermilab Tevatron Collider. Each scalar top quark is assumed to decay into a $b$ quark, a charged lepton, and a scalar neutrino ($\tilde{\nu}$). We investigate final states arising from $\tilde{t}_{1} \bar{\tilde{t}_{1}} \rightarrow b\bar{b}\mu\tau \tilde{\nu} \tilde{\nu}$ and $\tilde{t}_{1} \bar{\tilde{t}_{1}} \rightarrow b\bar{b}\tau\tau \tilde{\nu} \tilde{\nu}$. With no significant excess of events observed above the background expected from the standard model, we set exclusion limits on this production process in the ($m_{\tilde{t}_{1}}$,$m_{\tilde{\nu}}$) plane.

The technique of invariant mass spectroscopy has been used to measure, for the first time, the ground state energy of neutron-unbound $^{28}\textrm{F},$ determined to be a resonance in the $^{27}\textrm{F} + n$ continuum at $2\underline{2}0 (\underline{5}0)$ keV. States in $^{28}\textrm{F}$ were populated by the reactions of a 62 MeV/u $^{29}\textrm{Ne}$ beam impinging on a 288 $\textrm{mg/cm}^2$ beryllium target. The measured $^{28}\textrm{F}$ ground state energy is in good agreement with USDA/USDB shell model predictions, indicating that $pf$ shell intruder configurations play only a small role in the ground state structure of $^{28}\textrm{F}$ and establishing a low-$Z$ boundary of the island of inversion for N=19 isotones.

We measure the top quark mass (mt) in ppbar collisions at a center of mass energy of 1.96 TeV using dilepton ttbar->W+bW-bbar->l+nubl-nubarbbar events, where l denotes an electron, a muon, or a tau that decays leptonically. The data correspond to an integrated luminosity of 5.4 fb-1 collected with the D0 detector at the Fermilab Tevatron Collider. We obtain mt = 174.0 +- 1.8(stat) +- 2.4(syst) GeV, which is in agreement with the current world average mt = 173.3 +- 1.1 GeV. This is currently the most precise measurement of mt in the dilepton channel.

Neutron-unbound resonant states of 11Be were populated in neutron knock-out reactions from 12Be and identified by 10Be-n coincidence measurements. A resonance in the decay-energy spectrum at 80(2) keV was attributed to a highly excited unbound state in 11Be at 3.949(2) MeV decaying to the 2+ excited state in 10Be. A knockout cross section of 15(3) mb was inferred for this 3.949(2) MeV state suggesting a spectroscopic factor near unity for this 0p3/2- level, consistent with the detailed shell model calculations.

We report the result of a search for the pair production of the lightest supersymmetric partner of the top quark ($\tilde{t}_1$) in $p\bar{p}$ collisions at a center-of-mass energy of 1.96 TeV at the Fermilab Tevatron collider corresponding to an integrated luminosity of 5.4 fb$^{-1}$. The scalar top quarks are assumed to decay into a $b$ quark, a charged lepton, and a scalar neutrino ($\tilde{\nu}$), and the search is performed in the electron plus muon final state. No significant excess of events above the standard model prediction is detected, and improved exclusion limits at the 95% C.L. are set in the the ($M_{\tilde{t}_1}$,$M_{\tilde{\nu}}$) mass plane.

We extract the total width of the top quark, Gamma_t, from the partial decay width Gamma(t -> W b) measured using the t-channel cross section for single top quark production and from the branching fraction B(t -> W b) measured in ttbar events using up to 2.3 fb^-1 of integrated luminosity collected by the D0 Collaboration at the Tevatron ppbar Collider. The result is Gamma_t = 1.99 +0.69 -0.55 GeV, which translates to a top-quark lifetime of tau_t = (3.3 +1.3 -0.9) x 10^-25 s. Assuming a high mass fourth generation b' quark and unitarity of the four-generation quark-mixing matrix, we set the first upper limit on |Vtb'| < 0.63 at 95% C.L.

We have detected [Fe II] 17.94 um and 24.52 um emission from a sample of M supergiants using TEXES on the IRTF. These low opacity emission lines are resolved at R = 50, 000 and provide new diagnostics of the dynamics and thermodynamics of the stellar wind acceleration zone. The [Fe II] lines, from the first excited term, are sensitive to the warm plasma where energy is deposited into the extended atmosphere to form the chromosphere and wind outflow. These diagnostics complement previous KAO and ISO observations which were sensitive to the cooler and more extended circumstellar envelopes. The turbulent velocities, Vturb is about 12 to 13 km/s, observed in the [Fe II] forbidden lines are found to be a common property of our sample, and are less than that derived from the hotter chromospheric C II] 2325 Angstrom lines observed in alpha Ori, where Vturb is about 17 to 19 km/s. For the first time, we have dynamically resolved the motions of the dominant cool atmospheric component discovered in alpha Ori from multi-wavelength radio interferometry by Lim et al. (1998). Surprisingly, the emission centroids are quite Gaussian and at rest with respect to the M supergiants. These constraints combined with model calculations of the infrared emission line fluxes for alpha Ori imply that the warm material has a low outflow velocity and is located close to the star. We have also detected narrow [Fe I] 24.04 um emission that confirms that Fe II is the dominant ionization state in alpha Ori's extended atmosphere.