Increasing the Reliability of Adaptive Quadrature Using Explicit Interpolants

We present two new adaptive quadrature routines. Both routines differ from previously published algorithms in many aspects, most significantly in how they represent the integrand, how they treat non-numerical values of the integrand, how they deal with improper divergent integrals and how they estimate the integration error. The main focus of these improvements is to increase the reliability of the algorithms without significantly impacting their efficiency. Both algorithms are implemented in Matlab and tested using both the "families" suggested by Lyness and Kaganove and the battery test used by Gander and Gautschi and Kahaner. They are shown to be more reliable, albeit in some cases less efficient, than other commonly-used adaptive integrators.Comment: 32 pages, submitted to ACM Transactions on Mathematical Softwar

Point Contact Spectroscopy of Superconducting Gap Anisotropy in Nickel Borocarbide Compound LuNi2B2C

Point contacts are used to investigate the anisotropy of the superconducting energy gap in LuNi2B2C in the ab plane and along the c axis. It is shown that the experimental curves should be described assuming that the superconducting gap is non-uniformly distributed over the Fermi surface. The largest and the smallest gaps have been estimated by two-gap fitting models. It is found that the largest contribution to the point-contact conductivity in the c direction is made by a smaller gap and, in the ab plane by a larger gap. The deviation from the one-gap BCS model is pronounced in the temperature dependence of the gap in both directions. The temperature range, where the deviation occurs, is for the c direction approximately 1.5 times more than in the ab plane. The \Gamma parameter, allowing quantitatively estimate the gap anisotropy by one-gap fitting, in c direction is also about 1.5 times greater than in the ab plane. Since it is impossible to describe satisfactorily such gap distribution either by the one- or two-gap models, a continuous, dual-maxima model of gap distribution over the Fermi surface should be used to describe superconductivity in this material.Comment: 10 pages, 14 Figs, accepted in PR

Relaxation rates and collision integrals for Bose-Einstein condensates

Near equilibrium, the rate of relaxation to equilibrium and the transport properties of excitations (bogolons) in a dilute Bose-Einstein condensate (BEC) are determined by three collision integrals, $\mathcal{G}^{12}$, $\mathcal{G}^{22}$, and $\mathcal{G}^{31}$. All three collision integrals conserve momentum and energy during bogolon collisions, but only $\mathcal{G}^{22}$ conserves bogolon number. Previous works have considered the contribution of only two collision integrals, $\mathcal{G}^{22}$ and $\mathcal{G}^{12}$. In this work, we show that the third collision integral $\mathcal{G}^{31}$ makes a significant contribution to the bogolon number relaxation rate and needs to be retained when computing relaxation properties of the BEC. We provide values of relaxation rates in a form that can be applied to a variety of dilute Bose-Einstein condensates.Comment: 18 pages, 4 figures, accepted by Journal of Low Temperature Physics 7/201

Thermoluminescence of zircon: a kinetic model

The mineral zircon, ZrSiO4, belongs to a class of promising materials for geochronometry by means of thermoluminescence (TL) dating. The development of a reliable and reproducible method for TL dating with zircon requires detailed knowledge of the processes taking place during exposure to ionizing radiation, long-term storage, annealing at moderate temperatures and heating at a constant rate (TL measurements). To understand these processes one needs a kinetic model of TL. This paper is devoted to the construction of such amodel. The goal is to study the qualitative behaviour of the system and to determine the parameters and processes controlling TL phenomena of zircon. The model considers the following processes: (i) Filling of electron and hole traps at the excitation stage as a function of the dose rate and the dose for both (low dose rate) natural and (high dose rate) laboratory irradiation. (ii) Time dependence of TL fading in samples irradiated under laboratory conditions. (iii) Short time annealing at a given temperature. (iv) Heating of the irradiated sample to simulate TL experiments both after laboratory and natural irradiation. The input parameters of the model, such as the types and concentrations of the TL centres and the energy distributions of the hole and electron traps, were obtained by analysing the experimental data on fading of the TL-emission spectra of samples from different geological locations. Electron paramagnetic resonance (EPR) data were used to establish the nature of the TL centres. Glow curves and 3D TL emission spectra are simulated and compared with the experimental data on time-dependent TL fading. The saturation and annealing behaviour of filled trap concentrations has been considered in the framework of the proposed kinetic model and comparedwith the EPR data associated with the rare-earth ions Tb3+ and Dy3+, which play a crucial role as hole traps and recombination centres. Inaddition, the behaviour of some of the SiOmn− centres has been compared with simulation results.

Ancient numerical daemons of conceptual hydrological modeling 2. Impact of time stepping schemes on model analysis and prediction

Despite the widespread use of conceptual hydrological models in environmental research and operations, they remain frequently implemented using numerically unreliable methods. This paper considers the impact of the time stepping scheme on model analysis (sensitivity analysis, parameter optimization, and Markov chain Monte Carlo-based uncertainty estimation) and prediction. It builds on the companion paper (Clark and Kavetski, 2010), which focused on numerical accuracy, fidelity, and computational efficiency. Empirical and theoretical analysis of eight distinct time stepping schemes for six different hydrological models in 13 diverse basins demonstrates several critical conclusions. (1) Unreliable time stepping schemes, in particular, fixed-step explicit methods, suffer from troublesome numerical artifacts that severely deform the objective function of the model. These deformations are not rare isolated instances but can arise in any model structure, in any catchment, and under common hydroclimatic conditions. (2) Sensitivity analysis can be severely contaminated by numerical errors, often to the extent that it becomes dominated by the sensitivity of truncation errors rather than the model equations. (3) Robust time stepping schemes generally produce "better behaved" objective functions, free of spurious local optima, and with sufficient numerical continuity to permit parameter optimization using efficient quasi Newton methods. When implemented within a multistart framework, modern Newton-type optimizers are robust even when started far from the optima and provide valuable diagnostic insights not directly available from evolutionary global optimizers. (4) Unreliable time stepping schemes lead to inconsistent and biased inferences of the model parameters and internal states. (5) Even when interactions between hydrological parameters and numerical errors provide "the right result for the wrong reason" and the calibrated model performance appears adequate, unreliable time stepping schemes make the model unnecessarily fragile in predictive mode, undermining validation assessments and operational use. Erroneous or misleading conclusions of model analysis and prediction arising from numerical artifacts in hydrological models are intolerable, especially given that robust numerics are accepted as mainstream in other areas of science and engineering. We hope that the vivid empirical findings will encourage the conceptual hydrological community to close its Pandora's box of numerical problems, paving the way for more meaningful model application and interpretation. Copyright 2010 by the American Geophysical Union.Dmitri Kavetski and Martyn P. Clar

Multidecadal accumulation of anthropogenic and remineralized dissolved inorganic carbon along the Extended Ellett Line in the northeast Atlantic Ocean

Marine carbonate chemistry measurements have been carried out annually since 2009 during UK research cruises along the Extended Ellett Line (EEL), a hydrographic transect in the northeast Atlantic Ocean. The EEL intersects several water masses that are key to the global thermohaline circulation, and therefore the cruises sample a region in which it is critical to monitor secular physical and biogeochemical changes. We have combined results from these EEL cruises with existing quality-controlled observational data syntheses to produce a hydrographic time series for the EEL from 1981 to 2013. This reveals multidecadal increases in dissolved inorganic carbon (DIC) throughout the water column, with a near-surface maximum rate of 1.80 ± 0.45 µmol kg−1 yr−1. Anthropogenic CO2 accumulation was assessed, using simultaneous changes in apparent oxygen utilization (AOU) and total alkalinity (TA) as proxies for the biogeochemical processes that influence DIC. The stable carbon isotope composition of DIC (δ13CDIC) was also determined and used as an independent test of our method. We calculated a volume-integrated anthropogenic CO2 accumulation rate of 2.8 ± 0.4 mg C m−3 yr−1 along the EEL, which is about double the global mean. The anthropogenic CO2 component accounts for only 31 ± 6% of the total DIC increase. The remainder is derived from increased organic matter remineralization, which we attribute to the lateral redistribution of water masses that accompanies subpolar gyre contraction. Output from a general circulation ecosystem model demonstrates that spatiotemporal heterogeneity in the observations has not significantly biased our multidecadal rate of change calculations and indicates that the EEL observations have been tracking distal changes in the surrounding North Atlantic and Nordic Seas