Moon Shadow by Cosmic Rays under the Influence of Geomagnetic Field and Search for Antiprotons at Multi-TeV Energies

We have observed the shadowing of galactic cosmic ray flux in the direction of the moon, the so-called moon shadow, using the Tibet-III air shower array operating at Yangbajing (4300 m a.s.l.) in Tibet since 1999. Almost all cosmic rays are positively charged; for that reason, they are bent by the geomagnetic field, thereby shifting the moon shadow westward. The cosmic rays will also produce an additional shadow in the eastward direction of the moon if cosmic rays contain negatively charged particles, such as antiprotons, with some fraction. We selected 1.5 x10^{10} air shower events with energy beyond about 3 TeV from the dataset observed by the Tibet-III air shower array and detected the moon shadow at $\sim 40 \sigma$ level. The center of the moon was detected in the direction away from the apparent center of the moon by 0.23$^\circ$ to the west. Based on these data and a full Monte Carlo simulation, we searched for the existence of the shadow produced by antiprotons at the multi-TeV energy region. No evidence of the existence of antiprotons was found in this energy region. We obtained the 90% confidence level upper limit of the flux ratio of antiprotons to protons as 7% at multi-TeV energies.Comment: 13pages,4figures; Accepted for publication in Astroparticle Physic

Young and Intermediate-age Distance Indicators

Distance measurements beyond geometrical and semi-geometrical methods, rely mainly on standard candles. As the name suggests, these objects have known luminosities by virtue of their intrinsic proprieties and play a major role in our understanding of modern cosmology. The main caveats associated with standard candles are their absolute calibration, contamination of the sample from other sources and systematic uncertainties. The absolute calibration mainly depends on their chemical composition and age. To understand the impact of these effects on the distance scale, it is essential to develop methods based on different sample of standard candles. Here we review the fundamental properties of young and intermediate-age distance indicators such as Cepheids, Mira variables and Red Clump stars and the recent developments in their application as distance indicators.Comment: Review article, 63 pages (28 figures), Accepted for publication in Space Science Reviews (Chapter 3 of a special collection resulting from the May 2016 ISSI-BJ workshop on Astronomical Distance Determination in the Space Age

Exploration of the equilibrium operating space for NSTX-Upgrade

This paper explores a range of high-performance equilibrium scenarios available in the NSTX-Upgrade device [J.E. Menard, submitted for publication to Nuclear Fusion]. NSTX-Upgrade is a substantial upgrade to the existing NSTX device [M. Ono, et al., Nuclear Fusion 40, 557 (2000)], with significantly higher toroidal field and solenoid capabilities, and three additional neutral beam sources with significantly larger current drive efficiency. Equilibria are computed with freeboundary TRANSP, allowing a self consistent calculation of the non-inductive current drive sources, the plasma equilibrium, and poloidal field coil current, using the realistic device geometry. The thermal profiles are taken from a variety of existing NSTX discharges, and different assumptions for the thermal confinement scalings are utilized. The no-wall and idealwall n=1 stability limits are computed with the DCON code. The central and minimum safety factors are quite sensitive to many parameters: they generally increases with large outer plasmawall gaps and higher density, but can have either trend with the confinement enhancement factor. In scenarios with strong central beam current drive, the inclusion of non-classical fast ion diffusion raises qmin, decreases the pressure peaking, and generally improves the global stability, at the expense of a reduction in the non-inductive current drive fraction; cases with less beam current drive are largely insensitive to additional fast ion diffusion. The non-inductive current level is quite sensitive to the underlying confinement and profile assumptions. For instance, for BT=1.0 T and Pinj=12.6 MW, the non-inductive current level varies from 875 kA with ITER-98y,2 thermal confinement scaling and narrow thermal profiles to 1325 kA for an ST specific scaling expression and broad profiles. This sensitivity should facilitate the determination of the correct scaling of transport with current and field to use for future fully non-inductive ST devices. Scenarios are presented which can be sustained for 8-10 seconds, or (20-30)τCR, at βN=3.8-4.5, facilitating, for instance, the study of disruption avoidance for very long pulse. Scenarios have been documented which can operate with βT~25% and equilibrated qmin>1. The value of qmin can be controlled at either fixed non-inductive fraction of 100% or fixed plasma current, by varying which beam sources are used, opening the possibility for feedback qmin control. In terms of quantities like collisionality, neutron emission, non-inductive fraction, or stored energy, these scenarios represent a significant performance extension compared to NSTX and other present spherical torii

Peeling-ballooning stability of tokamak plasmas with applied 3D magnetic fields

The poloidal harmonics of the toroidal normal modes of an unstable axisymmetric tokamak plasma are employed as basis functions for the minimisation of the 3D energy functional. This approach presents a natural extension of the perturbative method considered in Anastopoulos Tzanis et al (2019 Nucl. Fusion 59 126028). This variational formulation is applied to the stability of tokamak plasmas subject to external non-axisymmetric magnetic fields. A comparison of the variational and perturbative methods shows that for D-shaped, high β N plasmas, the coupling of normal modes becomes strong at experimentally relevant applied 3D fields, leading to violation of the assumptions that justify a perturbative analysis. The variational analysis employed here addresses strong coupling, minimising energy with respect to both toroidal and poloidal Fourier coefficients. In general, it is observed that ballooning unstable modes are further destabilised by the applied 3D fields and field-aligned localisation of the perturbation takes place, as local ballooning theory suggests. For D-shaped high β N plasmas, relevant to experimental cases, it is observed that the existence of intermediate n unstable peeling-ballooning modes, where a maximum in the growth rate spectrum typically occurs, leads to a destabilising synergistic coupling that strongly degrades the stability of the 3D system

Framework for assessing and improving the performance of recursive digital filters for baseflow estimation with application to the Lyne and Hollick filter

Baseflow is often regarded as the streamflow component derived predominantly from groundwater discharge. The estimation of baseflow is important for water supply, water allocation, investigation of contamination impacts, low flow hydrology and flood hydrology. Baseflow is commonly estimated using graphical methods, recursive digital filters (RDFs), tracer based methods, and conceptual models. Of all of these methods, RDFs are the most commonly used, due to their relatively easy and efficient implementation. This paper presents a generic framework for assessing and improving the performance of RDFs for baseflow estimation for catchments with different characteristics and subject to different hydrological conditions. As part of the framework, a fully integrated surface water/groundwater (SW/GW) model is used to obtain estimates of streamflow and baseflow for catchments with different properties, such as soil types and rainfall patterns. An RDF is then applied to the simulated streamflow to assess how well the baseflow obtained using the filter matches the baseflow obtained using the fully integrated SW/GW model. In order to improve the performance of the filter, the user-defined parameter(s) controlling filter operation can be adjusted in order to obtain the best match between the baseflow obtained using the filter and that obtained using the fully integrated SW/GW model (i.e. through calibration). The proposed framework is tested by applying it to a common SW/GW benchmarking problem, the tilted V-catchment, for a range of soil properties. HydroGeoSphere (HGS) is used to develop the fully integrated SW/GW model and the Lyne and Hollick (LH) filter is used as the RDF. The performance of the LH filter is assessed using the commonly used value of the filter parameter of 0.925, as well as calibrated filter parameter values. The results obtained show that the performance of the LH filter is affected significantly by the saturated hydraulic conductivity (Ks) of the soil and that calibrated LH filter parameter can result in significant improvements in filter performance. © 2012 Elsevier Ltd.L. Li, H.R. Maier, M.F. Lambert, C.T. Simmons, D. Partingto

On the mechanisms governing gas penetration into a tokamak plasma during a massive gas injection

A new 1D radial fluid code, IMAGINE, is used to simulate the penetration of gas into a tokamak plasma during a massive gas injection (MGI). The main result is that the gas is in general strongly braked as it reaches the plasma, due to mechanisms related to charge exchange and (to a smaller extent) recombination. As a result, only a fraction of the gas penetrates into the plasma. Also, a shock wave is created in the gas which propagates away from the plasma, braking and compressing the incoming gas. Simulation results are quantitatively consistent, at least in terms of orders of magnitude, with experimental data for a D 2 MGI into a JET Ohmic plasma. Simulations of MGI into the background plasma surrounding a runaway electron beam show that if the background electron density is too high, the gas may not penetrate, suggesting a possible explanation for the recent results of Reux et al in JET (2015 Nucl. Fusion 55 093013)