Fast smoothing parameter separation in multidimensional generalized P-splines: the SAP algorithm

A new computational algorithm for estimating the smoothing parameters of a multidimensional penalized spline generalized linear model with anisotropic penalty is presented. This new proposal is based on the mixed model representation of a multidimensional P-spline, in which the smoothing parameter for each covariate is expressed in terms of variance components. On the basis of penalized quasi-likelihood methods, closed-form expressions for the estimates of the variance components are obtained. This formulation leads to an efficient implementation that considerably reduces the computational burden. The proposed algorithm can be seen as a generalization of the algorithm by Schall (1991)-for variance components estimation-to deal with non-standard structures of the covariance matrix of the random effects. The practical performance of the proposed algorithm is evaluated by means of simulations, and comparisons with alternative methods are made on the basis of the mean square error criterion and the computing time. Finally, we illustrate our proposal with the analysis of two real datasets: a two dimensional example of historical records of monthly precipitation data in USA and a three dimensional one of mortality data from respiratory disease according to the age at death, the year of death and the month of death.The authors would like to express their gratitude for the support received in the form of the Spanish Ministry of Economy and Competitiveness grants MTM2011-28285-C02-01 and MTM2011-28285-C02-02. The research of Dae-Jin Lee was funded by an NIH grant for the Superfund Metal Mixtures, Biomarkers and Neurodevelopment project 1PA2ES016454-01A2

Ti6Al4V metal cutting chip formation experiments and modelling over a wide range of cutting speeds

Measured forces, chip geometry and tool temperatures from machining a mill annealed Ti6Al4V at cutting speeds mainly from 1 to 100 m/min, but in some cases down to 0.1 m/min, are reported, as well as mechanical testing of the material. Finite element simulations with inputs the measured flow stress, and subsequently a small high temperature strain hardening recovery correction, and a failure model calibrated from the cutting tests at speeds from 1 to 10 m/min, give satisfactory agreement with the higher speed tests once surface strain hardening and damage from the previous pass of the tool are taken into account. This paper’s originality is firstly to show that more complicated flow stress models involving large strain softening are not needed provided shear failure is included; and secondly its failure model: this proposes a non-zero failed shear stress depending on local pressure and temperature. The simulations provide relations between tool mechanical and thermal loading and cutting conditions to aid process improvement

Holographic dark energy with time varying c2c^2 parameter

We consider the holographic dark energy model in which the model parameter $c^2$ evolves slowly with time. First we calculate the evolution of EoS parameter as well as the deceleration parameter in this generalized version of holographic dark energy (GHDE). Depending on the parameter $c^2$, the phantom regime can be achieved earlier or later compare with original version of holographic dark energy. The evolution of energy density of GHDE model is investigated in terms of parameter $c^2$. We also show that the time-dependency of $c^2$ can effect on the transition epoch from decelerated phase to accelerated expansion. Finally, we perform the statefinder diagnostic for GHDE model and show that the evolutionary trajectories of the model in $s-r$ plane are strongly depend on the parameter $c^2$.Comment: 16 pages, 4 figures, accepted by Astrophys Space Sc

Au+Au Reactions at the AGS: Experiments E866 and E917

Particle production and correlation functions from Au+Au reactions have been measured as a function of both beam energy (2-10.7AGeV) and impact parameter. These results are used to probe the dynamics of heavy-ion reactions, confront hadronic models over a wide range of conditions and to search for the onset of new phenomena.Comment: 12 pages, 14 figures, Talk presented at Quark Matter '9

Horizontal Branch Stars: The Interplay between Observations and Theory, and Insights into the Formation of the Galaxy

We review HB stars in a broad astrophysical context, including both variable and non-variable stars. A reassessment of the Oosterhoff dichotomy is presented, which provides unprecedented detail regarding its origin and systematics. We show that the Oosterhoff dichotomy and the distribution of globular clusters (GCs) in the HB morphology-metallicity plane both exclude, with high statistical significance, the possibility that the Galactic halo may have formed from the accretion of dwarf galaxies resembling present-day Milky Way satellites such as Fornax, Sagittarius, and the LMC. A rediscussion of the second-parameter problem is presented. A technique is proposed to estimate the HB types of extragalactic GCs on the basis of integrated far-UV photometry. The relationship between the absolute V magnitude of the HB at the RR Lyrae level and metallicity, as obtained on the basis of trigonometric parallax measurements for the star RR Lyrae, is also revisited, giving a distance modulus to the LMC of (m-M)_0 = 18.44+/-0.11. RR Lyrae period change rates are studied. Finally, the conductive opacities used in evolutionary calculations of low-mass stars are investigated. [ABRIDGED]Comment: 56 pages, 22 figures. Invited review, to appear in Astrophysics and Space Scienc

Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter

Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz'min energy threshold, \nobreak{$6\times 10^{19}$eV}. The anisotropy was measured by the fraction of arrival directions that are less than $3.1^\circ$ from the position of an active galactic nucleus within 75 Mpc (using the V\'eron-Cetty and V\'eron $12^{\rm th}$ catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is $(38^{+7}_{-6})$, compared with $21$ expected for isotropic cosmic rays. This is down from the early estimate of $(69^{+11}_{-13})$. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation.Comment: Accepted for publication in Astroparticle Physics on 31 August 201

The Fluorescence Detector of the Pierre Auger Observatory

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions.Comment: 53 pages. Submitted to Nuclear Instruments and Methods in Physics Research Section

Advanced functionality for radio analysis in the Offline software framework of the Pierre Auger Observatory

The advent of the Auger Engineering Radio Array (AERA) necessitates the development of a powerful framework for the analysis of radio measurements of cosmic ray air showers. As AERA performs "radio-hybrid" measurements of air shower radio emission in coincidence with the surface particle detectors and fluorescence telescopes of the Pierre Auger Observatory, the radio analysis functionality had to be incorporated in the existing hybrid analysis solutions for fluoresence and surface detector data. This goal has been achieved in a natural way by extending the existing Auger Offline software framework with radio functionality. In this article, we lay out the design, highlights and features of the radio extension implemented in the Auger Offline framework. Its functionality has achieved a high degree of sophistication and offers advanced features such as vectorial reconstruction of the electric field, advanced signal processing algorithms, a transparent and efficient handling of FFTs, a very detailed simulation of detector effects, and the read-in of multiple data formats including data from various radio simulation codes. The source code of this radio functionality can be made available to interested parties on request.Comment: accepted for publication in NIM A, 13 pages, minor corrections to author list and references in v

Search for First Harmonic Modulation in the Right Ascension Distribution of Cosmic Rays Detected at the Pierre Auger Observatory

We present the results of searches for dipolar-type anisotropies in different energy ranges above $2.5\times 10^{17}$ eV with the surface detector array of the Pierre Auger Observatory, reporting on both the phase and the amplitude measurements of the first harmonic modulation in the right-ascension distribution. Upper limits on the amplitudes are obtained, which provide the most stringent bounds at present, being below 2% at 99% $C.L.$ for EeV energies. We also compare our results to those of previous experiments as well as with some theoretical expectations.Comment: 28 pages, 11 figure

Search for a W' boson decaying to a bottom quark and a top quark in pp collisions at sqrt(s) = 7 TeV

Results are presented from a search for a W' boson using a dataset corresponding to 5.0 inverse femtobarns of integrated luminosity collected during 2011 by the CMS experiment at the LHC in pp collisions at sqrt(s)=7 TeV. The W' boson is modeled as a heavy W boson, but different scenarios for the couplings to fermions are considered, involving both left-handed and right-handed chiral projections of the fermions, as well as an arbitrary mixture of the two. The search is performed in the decay channel W' to t b, leading to a final state signature with a single lepton (e, mu), missing transverse energy, and jets, at least one of which is tagged as a b-jet. A W' boson that couples to fermions with the same coupling constant as the W, but to the right-handed rather than left-handed chiral projections, is excluded for masses below 1.85 TeV at the 95% confidence level. For the first time using LHC data, constraints on the W' gauge coupling for a set of left- and right-handed coupling combinations have been placed. These results represent a significant improvement over previously published limits.Comment: Submitted to Physics Letters B. Replaced with version publishe