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

The Cosmic Linear Anisotropy Solving System (CLASS) IV: Efficient implementation of non-cold relics

We present a new flexible, fast and accurate way to implement massive neutrinos, warm dark matter and any other non-cold dark matter relics in Boltzmann codes. For whatever analytical or numerical form of the phase-space distribution function, the optimal sampling in momentum space compatible with a given level of accuracy is automatically found by comparing quadrature methods. The perturbation integration is made even faster by switching to an approximate viscous fluid description inside the Hubble radius, which differs from previous approximations discussed in the literature. When adding one massive neutrino to the minimal cosmological model, CLASS becomes just 1.5 times slower, instead of about 5 times in other codes (for fixed accuracy requirements). We illustrate the flexibility of our approach by considering a few examples of standard or non-standard neutrinos, as well as warm dark matter models.Comment: 23 pages, 8 figures, 3 tables. Matches published version. Code available at http://class-code.ne

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

Real structured singular value synthesis using the scaled Popov criterion

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/77238/1/AIAA-21537-608.pd

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

Critical Casimir forces and adsorption profiles in the presence of a chemically structured substrate

Motivated by recent experiments with confined binary liquid mixtures near demixing, we study the universal critical properties of a system, which belongs to the Ising universality class, in the film geometry. We employ periodic boundary conditions in the two lateral directions and fixed boundary conditions on the two confining surfaces, such that one of them has a spatially homogeneous adsorption preference while the other one exhibits a laterally alternating adsorption preference, resembling locally a single chemical step. By means of Monte Carlo simulations of an improved Hamiltonian, so that the leading scaling corrections are suppressed, numerical integration, and finite-size scaling analysis we determine the critical Casimir force and its universal scaling function for various values of the aspect ratio of the film. In the limit of a vanishing aspect ratio the critical Casimir force of this system reduces to the mean value of the critical Casimir force for laterally homogeneous ++ and +- boundary conditions, corresponding to the surface spins on the two surfaces being fixed to equal and opposite values, respectively. We show that the universal scaling function of the critical Casimir force for small but finite aspect ratios displays a linear dependence on the aspect ratio which is solely due to the presence of the lateral inhomogeneity. We also analyze the order-parameter profiles at criticality and their universal scaling function which allows us to probe theoretical predictions and to compare with experimental data.Comment: revised version, section 5.2 expanded; 53 pages, 12 figures, iopart clas

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

Computational assessment of insulin secretion and insulin sensitivity from 2-h oral glucose tolerance tests for clinical use for type 2 diabetes

In type 2 diabetes mellitus, glucose homeostasis is tightly maintained through insulin secretion and insulin sensitivity. Therefore, finding an accurate method to assess insulin secretion and sensitivity using clinically available data would enhance the quality of diabetic medical care. In an effort to find such a method, we developed a computational approach to derive indices of these factors using a 2-h oral glucose tolerance test (OGTT). To evaluate our method, clinical data from subjects who received an OGTT and a glucose clamp test were examined. Our insulin secretion index was significantly correlated with an analogous index obtained from a hyperglycemic clamp test (r = 0.90, n = 46, p < 0.001). Our insulin sensitivity index sensitivity was also significantly correlated with an analogous index obtained from a hyperinsulinemic-euglycemic clamp test (r = 0.56, n = 79, p < 0.001). These results suggest that our method can potentially provide an accurate and convenient tool toward improving the management of diabetes in clinical practice by assessing insulin secretion and insulin sensitivity