The statistics of refractive error maps : managing wavefront aberration analysis without Zernike polynomials

The refractive error of a human eye varies across the pupil and therefore may be treated as a random variable. The probability distribution of this random variable provides a means for assessing the main refractive properties of the eye without the necessity of traditional functional representation of wavefront aberrations. To demonstrate this approach, the statistical properties of refractive error maps are investigated. Closed-form expressions are derived for the probability density function (PDF) and its statistical moments for the general case of rotationally-symmetric aberrations. A closed-form expression for a PDF for a general non-rotationally symmetric wavefront aberration is difficult to derive. However, for specific cases, such as astigmatism, a closed-form expression of the PDF can be obtained. Further, interpretation of the distribution of the refractive error map as well as its moments is provided for a range of wavefront aberrations measured in real eyes. These are evaluated using a kernel density and sample moments estimators. It is concluded that the refractive error domain allows non-functional analysis of wavefront aberrations based on simple statistics in the form of its sample moments. Clinicians may find this approach to wavefront analysis easier to interpret due to the clinical familiarity and intuitive appeal of refractive error maps.

Zernike radial slope polynomials for wavefront reconstruction and refraction

Ophthalmic wavefront sensors typically measure wavefront slope, from which wavefront phase is reconstructed. We show that ophthalmic prescriptions (in power-vector format) can be obtained directly from slope measurements without wavefront reconstruction. This is achieved by fitting the measurement data with a new set of orthonormal basis functions called Zernike radial slope polynomials. Coefficients of this expansion can be used to specify the ophthalmic power vector using explicit formulas derived by a variety of methods. Zernike coefficients for wavefront error can be recovered from the coefficients of radial slope polynomials, thereby offering an alternative way to perform wavefront reconstruction.

Numerical simulation of a new two-dimensional variable-order fractional percolation equation in non-homogeneous porous media

Percolation flow problems are discussed in many research fields, such as seepage hydraulics, groundwater hydraulics, groundwater dynamics and fluid dynamics in porous media. Many physical processes appear to exhibit fractional-order behavior that may vary with time, or space, or space and time. The theory of pseudodifferential operators and equations has been used to deal with this situation. In this paper we use a fractional Darcys law with variable order Riemann-Liouville fractional derivatives, this leads to a new variable-order fractional percolation equation. In this paper, a new two-dimensional variable-order fractional percolation equation is considered. A new implicit numerical method and an alternating direct method for the two-dimensional variable-order fractional model is proposed. Consistency, stability and convergence of the implicit finite difference method are established. Finally, some numerical examples are given. The numerical results demonstrate the effectiveness of the methods. This technique can be used to simulate a three-dimensional variable-order fractional percolation equation.

Cell migration and proliferation on homogeneous and non-homogeneous domains : modelling on the scale of individuals and populations

Cell migration is a behaviour critical to many key biological effects, including wound healing, cancerous cell invasion and morphogenesis, the development of an organism from an embryo. However, given that each of these situations is distinctly different and cells are extremely complicated biological objects, interest lies in more basic experiments which seek to remove conflating factors and present a less complex environment within which cell migration can be experimentally examined. These include in vitro studies like the scratch assay or circle migration assay, and ex vivo studies like the colonisation of the hindgut by neural crest cells. The reduced complexity of these experiments also makes them much more enticing as problems to mathematically model, like done here. The primary goal of the mathematical models used in this thesis is to shed light on which cellular behaviours work to generate the travelling waves of invasion observed in these experiments, and to explore how variations in these behaviours can potentially predict differences in this invasive pattern which are experimentally observed when cell types or chemical environment are changed. Relevant literature has already identified the difficulty of distinguishing between these behaviours when using traditional mathematical biology techniques operating on a macroscopic scale, and so here a sophisticated individual-cell-level model, an extension of the Cellular Potts Model (CPM), is been constructed and used to model a scratch assay experiment. This model includes a novel mechanism for dealing with cell proliferations that allowed for the differing properties of quiescent and proliferative cells to be implemented into their behaviour. This model is considered both for its predictive power and used to make comparisons with the travelling waves which result in more traditional macroscopic simulations. These comparisons demonstrate a surprising amount of agreement between the two modelling frameworks, and suggest further novel modifications to the CPM that would allow it to better model cell migration. Considerations of the model’s behaviour are used to argue that the dominant effect governing cell migration (random motility or signal-driven taxis) likely depends on the sort of invasion demonstrated by cells, as easily seen by microscopic photography. Additionally, a scratch assay simulated on a non-homogeneous domain consisting of a ’fast’ and ’slow’ region is also used to further differentiate between these different potential cell motility behaviours. A heterogeneous domain is a novel situation which has not been considered mathematically in this context, nor has it been constructed experimentally to the best of the candidate’s knowledge. Thus this problem serves as a thought experiment used to test the conclusions arising from the simulations on homogeneous domains, and to suggest what might be observed should this non-homogeneous assay situation be experimentally realised. Non-intuitive cell invasion patterns are predicted for diffusely-invading cells which respond to a cell-consumed signal or nutrient, contrasted with rather expected behaviour in the case of random-motility-driven invasion. The potential experimental observation of these behaviours is demonstrated by the individual-cell-level model used in this thesis, which does agree with the PDE model in predicting these unexpected invasion patterns. In the interest of examining such a case of a non-homogeneous domain experimentally, some brief suggestion is made as to how this could be achieved.

Mid- and near-infrared spectroscopic investigation of homogeneous cation distribution in MgxZnyAl(x+y)/2-layered double hydroxide (LDH)

In this report, a detailed FTIR fitting analysis was used to recognize Mg, Zn and Al homogeneous distribution in MgxZnyAl(x+y)/2-Layered double hydroxide (LDH) hydroxyl layer. In detail, OH-Mg2Al:OH-Mg3 ratios decreased from 95.2:4.8 (MIR) and 94.2:5.8 (NIR) to 58.9:41.1 (MIR) and 61.8:38.2 (NIR), when Mg:Al increased from 2.2:1.0 to 4.1:1.0 in MgAl-LDHs. These fitting results were similar with theoretical calculations of 94.3:5.7 and 59.0:41.0. In a further analysis of MgxZnyAl(x+y)/2-LDHs, OH bonded Zn2Mg, Zn2Al, MgZnAl, Mg2Al and Mg2Zn peaks were identified at 3420, 3430, 3445–3450, 3454 and 3545 cm-1, respectively. With the decrease of Mg:Zn from 3:1 to 1:3, metal-hydroxyl bands changed from OH-Mg2Al and MgZnAl (with a ratio of 49.4:50.6) to OH-MgZnAl and Zn2Al (with a ratio of 55.0:45.0). They were also similar with theoretical calculations of 47.6:52.4 and 54.6:45.4. As a result, these results show that there is an ordered cation distribution in MgxZnyAl(x+y)/2-LDH, and FTIR is feasible in recognizing this structure.

Permutation polynomials of the form (xp −x +δ)s +L(x)

Recently, several classes of permutation polynomials of the form (x2 + x + δ)s + x over F2m have been discovered. They are related to Kloosterman sums. In this paper, the permutation behavior of polynomials of the form (xp − x + δ)s + L(x) over Fpm is investigated, where L(x) is a linearized polynomial with coefficients in Fp. Six classes of permutation polynomials on F2m are derived. Three classes of permutation polynomials over F3m are also presented.

Revisiting Polya's summation techniques using a spreadsheet : from addition tables to Bernoulli polynomials

Recurrence relations in mathematics form a very powerful and compact way of looking at a wide range of relationships. Traditionally, the concept of recurrence has often been a difficult one for the secondary teacher to convey to students. Closely related to the powerful proof technique of mathematical induction, recurrences are able to capture many relationships in formulas much simpler than so-called direct or closed formulas. In computer science, recursive coding often has a similar compactness property, and, perhaps not surprisingly, suffers from similar problems in the classroom as recurrences: the students often find both the basic concepts and practicalities elusive. Using models designed to illuminate the relevant principles for the students, we offer a range of examples which use the modern spreadsheet environment to powerfully illustrate the great expressive and computational power of recurrences.

Homogeneous bent functions of degree n in 2n variables do not exist for n>3

We prove that homogeneous bent functions f:GF(2)^2n --> GF(2) of degree n do not exist for n>3. Consequently homogeneous bent functions must have degree 3.

The eight variable homogeneous degree three bent functions

We determine the affine equivalence classes of the eight variable degree three homogeneous bent functions using a new algorithm. Our algorithm applies to general bent functions and can systematically determine the automorphism groups. We provide a partial verification of the enumeration of eight variable degree three homogeneous bent functions obtained by Meng et al. We determine the affine equivalence classes of these functions.

One-step homogeneous C-reactive protein assay for saliva

Background: Cardiovascular disease is the leading cause of death in the world. Human C-reactive protein (CRP) has been used in the risk assessment of coronary events. Human saliva mirrors the body's health and well-being and is non-invasive, easy to collect and ideal for third world countries as well as for large patient screening. The aim was to establish a saliva CRP reference range and to demonstrate the clinical utility of salivary CRP levels in assessing the coronary events in a primary health care setting. Methods: We have used a homogeneous bead based assay to detect CRP levels in human saliva. We have developed a rapid 15 min (vs 90 min), sequential, one-step assay to detect CRP in saliva. Saliva was collected from healthy volunteers (n = 55, ages 20-70 years) as well as from cardiac patients (n = 28, ages 43-86 years). Results: The assay incubation time was optimised from 90 min to 15 mm and generated a positive correlation (n = 29, range 10-2189 pg/mL, r2 = 0.94; Passing Bablok slope 0.885. Intercept 0, p>0.10), meaning we could decrease the incubation time and produce equivalent results with confidence. The mean CRP level in the saliva of healthy human volunteers was 285 pg/mL and in cardiac patients was 1680 pg/mL (p<0.01). Analysis of CRP concentrations in paired serum and saliva samples from cardiac patients gave a positive correlation (r2 = 0.84, p<0.001) and the salivary CRP concentration capable of distinguishing healthy from diseased patients. Conclusions: The results suggest that this minimally invasive, rapid and sensitive assay will be useful in large patient screening studies for risk assessment of coronary events. (C) 2011 Elsevier B.V. All rights reserved.

A thermal expulsion approach to homogeneous large-volume methacrylate monolith preparation; enabling large-scale rapid purification of biomolecules

Numerous efforts have been dedicated to the synthesis of large-volume methacrylate monoliths for large-scale biomolecules purification but most were obstructed by the enormous release of exotherms during preparation, thereby introducing structural heterogeneity in the monolith pore system. A significant radial temperature gradient develops along the monolith thickness, reaching a terminal temperature that supersedes the maximum temperature required for structurally homogenous monoliths preparation. The enormous heat build-up is perceived to encompass the heat associated with initiator decomposition and the heat released from free radical-monomer and monomer-monomer interactions. The heat resulting from the initiator decomposition was expelled along with some gaseous fumes before commencing polymerization in a gradual addition fashion. Characteristics of 80 mL monolith prepared using this technique was compared with that of a similar monolith synthesized in a bulk polymerization mode. An extra similarity in the radial temperature profiles was observed for the monolith synthesized via the heat expulsion technique. A maximum radial temperature gradient of only 4.3°C was recorded at the center and 2.1°C at the monolith peripheral for the combined heat expulsion and gradual addition technique. The comparable radial temperature distributions obtained birthed identical pore size distributions at different radial points along the monolith thickness.

A homogeneous surface-enhanced Raman scattering platform for ultra-trace detection of trinitrotoluene in the environment

A facile and sensitive surface-enhanced Raman scattering substrate was prepared by controlled potentiostatic deposition of a closely packed single layer of gold nanostructures (AuNS) over a flat gold (pAu) platform. The nanometer scale inter-particle distance between the particles resulted in high population of ‘hot spots’ which enormously enhanced the scattered Raman photons. A renewed methodology was followed to precisely quantify the SERS substrate enhancement factor (SSEF) and it was estimated to be (2.2 ± 0.17) × 105. The reproducibility of the SERS signal acquired by the developed substrate was tested by establishing the relative standard deviation (RSD) of 150 repeated measurements from various locations on the substrate surface. A low RSD of 4.37 confirmed the homogeneity of the developed substrate. The sensitivity of pAu/AuNS was proven by determining 100 fM 2,4,6-trinitrotoluene (TNT) comfortably. As a proof of concept on the potential of the new pAu/AuNS substrate in field analysis, TNT in soil and water matrices was selectively detected after forming a Meisenheimer complex with cysteamine.