The choice of a suitable sampler for benthic macroinvertebrates in deep rivers

Both chemical and biological methods are used to assess the water quality of rivers. Many standard physical and chemical methods are now established, but biological procedures of comparable accuracy and versatility are still lacking. This is unfortunate because the biological assessment of water quality has several advantages over physical and chemical analyses. Several groups of organisms have been used to assess water quality in rivers and these include Bacteria, Protozoa, Algae, macrophytes, macroinvertebrates and fish. Hellawell (1978) provides an excellent review of the advantages and disadvantages of these groups, and concludes that macroinvertebrates are the most useful for monitoring water quality. Although macroinvertebrates are relatively easy to sample in shallow water (depth 1m). The present paper first considers different types of samplers with emphasis on immediate samplers, and then discusses some problems in choosing a suitable sampler for benthic macroinvertebrates in deep rivers

Weak convergence of finite element approximations of linear stochastic evolution equations with additive noise II. Fully discrete schemes

We present an abstract framework for analyzing the weak error of fully discrete approximation schemes for linear evolution equations driven by additive Gaussian noise. First, an abstract representation formula is derived for sufficiently smooth test functions. The formula is then applied to the wave equation, where the spatial approximation is done via the standard continuous finite element method and the time discretization via an I-stable rational approximation to the exponential function. It is found that the rate of weak convergence is twice that of strong convergence. Furthermore, in contrast to the parabolic case, higher order schemes in time, such as the Crank-Nicolson scheme, are worthwhile to use if the solution is not very regular. Finally we apply the theory to parabolic equations and detail a weak error estimate for the linearized Cahn-Hilliard-Cook equation as well as comment on the stochastic heat equation

Two-magnon Raman scattering in insulating cuprates: Modifications of the effective Raman operator

Calculations of Raman scattering intensities in spin 1/2 square-lattice Heisenberg model, using the Fleury-Loudon-Elliott theory, have so far been unable to describe the broad line shape and asymmetry of the two magnon peak found experimentally in the cuprate materials. Even more notably, the polarization selection rules are violated with respect to the Fleury-Loudon-Elliott theory. There is comparable scattering in $B_{1g}$ and $A_{1g}$ geometries, whereas the theory would predict scattering in only $B_{1g}$ geometry. We review various suggestions for this discrepency and suggest that at least part of the problem can be addressed by modifying the effective Raman Hamiltonian, allowing for two-magnon states with arbitrary total momentum. Such an approach based on the Sawatzsky-Lorenzana theory of optical absorption assumes an important role of phonons as momentum sinks. It leaves the low energy physics of the Heisenberg model unchanged but substantially alters the Raman line-shape and selection rules, bringing the results closer to experiments.Comment: 7 pages, 6 figures, revtex. Contains some minor revisions from previous versio

Special fast diffusion with slow asymptotics. Entropy method and flow on a Riemannian manifold

We consider the asymptotic behaviour of positive solutions $u(t,x)$ of the fast diffusion equation $u_t=\Delta (u^{m}/m)={\rm div} (u^{m-1}\nabla u)$ posed for x\in\RR^d, $t>0$, with a precise value for the exponent $m=(d-4)/(d-2)$. The space dimension is $d\ge 3$ so that $m<1$, and even $m=-1$ for $d=3$. This case had been left open in the general study \cite{BBDGV} since it requires quite different functional analytic methods, due in particular to the absence of a spectral gap for the operator generating the linearized evolution. The linearization of this flow is interpreted here as the heat flow of the Laplace-Beltrami operator of a suitable Riemannian Manifold (\RR^d,{\bf g}), with a metric ${\bf g}$ which is conformal to the standard \RR^d metric. Studying the pointwise heat kernel behaviour allows to prove {suitable Gagliardo-Nirenberg} inequalities associated to the generator. Such inequalities in turn allow to study the nonlinear evolution as well, and to determine its asymptotics, which is identical to the one satisfied by the linearization. In terms of the rescaled representation, which is a nonlinear Fokker--Planck equation, the convergence rate turns out to be polynomial in time. This result is in contrast with the known exponential decay of such representation for all other values of $m$.Comment: 37 page

Finite Element Convergence for the Joule Heating Problem with Mixed Boundary Conditions

We prove strong convergence of conforming finite element approximations to the stationary Joule heating problem with mixed boundary conditions on Lipschitz domains in three spatial dimensions. We show optimal global regularity estimates on creased domains and prove a priori and a posteriori bounds for shape regular meshes.Comment: Keywords: Joule heating problem, thermistors, a posteriori error analysis, a priori error analysis, finite element metho

Signatures of Thermal Dilepton Radiation at RHIC

The properties of thermal dilepton production from heavy-ion collisions in the RHIC energy regime are evaluated for invariant masses ranging from 0.5 to 3 GeV. Using an expanding thermal fireball to model the evolution through both quark-gluon and hadronic phases various features of the spectra are addressed. In the low-mass region, due to an expected large background, the focus is on possible medium modifications of the narrow resonance structures from $\omega$ and $\phi$ mesons, whereas in the intermediate-mass region the old idea of identifying QGP radiation is reiterated including effects of chemical under-saturation in the early stages of central Au+Au collisions.Comment: 17 pages ReVTeX including 16 figure

Optimal designs for rational function regression

We consider optimal non-sequential designs for a large class of (linear and nonlinear) regression models involving polynomials and rational functions with heteroscedastic noise also given by a polynomial or rational weight function. The proposed method treats D-, E-, A-, and $\Phi_p$-optimal designs in a unified manner, and generates a polynomial whose zeros are the support points of the optimal approximate design, generalizing a number of previously known results of the same flavor. The method is based on a mathematical optimization model that can incorporate various criteria of optimality and can be solved efficiently by well established numerical optimization methods. In contrast to previous optimization-based methods proposed for similar design problems, it also has theoretical guarantee of its algorithmic efficiency; in fact, the running times of all numerical examples considered in the paper are negligible. The stability of the method is demonstrated in an example involving high degree polynomials. After discussing linear models, applications for finding locally optimal designs for nonlinear regression models involving rational functions are presented, then extensions to robust regression designs, and trigonometric regression are shown. As a corollary, an upper bound on the size of the support set of the minimally-supported optimal designs is also found. The method is of considerable practical importance, with the potential for instance to impact design software development. Further study of the optimality conditions of the main optimization model might also yield new theoretical insights.Comment: 25 pages. Previous version updated with more details in the theory and additional example

Geometric partial differential equations: Theory, numerics and applications

This workshop concentrated on partial differential equations involving stationary and evolving surfaces in which geometric quantities play a major role. Mutual interest in this emerging field stimulated the interaction between analysis, numerical solution, and applications

On a diffuse interface model for tumour growth with non-local interactions and degenerate mobilities

We study a non-local variant of a diffuse interface model proposed by Hawkins--Darrud et al. (2012) for tumour growth in the presence of a chemical species acting as nutrient. The system consists of a Cahn--Hilliard equation coupled to a reaction-diffusion equation. For non-degenerate mobilities and smooth potentials, we derive well-posedness results, which are the non-local analogue of those obtained in Frigeri et al. (European J. Appl. Math. 2015). Furthermore, we establish existence of weak solutions for the case of degenerate mobilities and singular potentials, which serves to confine the order parameter to its physically relevant interval. Due to the non-local nature of the equations, under additional assumptions continuous dependence on initial data can also be shown.Comment: 28 page

Fertility, Living Arrangements, Care and Mobility

There are four main interconnecting themes around which the contributions in this book are based. This introductory chapter aims to establish the broad context for the chapters that follow by discussing each of the themes. It does so by setting these themes within the overarching demographic challenge of the twenty-first century – demographic ageing. Each chapter is introduced in the context of the specific theme to which it primarily relates and there is a summary of the data sets used by the contributors to illustrate the wide range of cross-sectional and longitudinal data analysed