Variational Bayes and the reduced dependence approximation for the autologistic model on an irregular grid with applications

Discrete Markov random field models provide a natural framework for representing images or spatial datasets. They model the spatial association present while providing a convenient Markovian dependency structure and strong edge-preservation properties. However, parameter estimation for discrete Markov random field models is difficult due to the complex form of the associated normalizing constant for the likelihood function. For large lattices, the reduced dependence approximation to the normalizing constant is based on the concept of performing computationally efficient and feasible forward recursions on smaller sublattices which are then suitably combined to estimate the constant for the whole lattice. We present an efficient computational extension of the forward recursion approach for the autologistic model to lattices that have an irregularly shaped boundary and which may contain regions with no data; these lattices are typical in applications. Consequently, we also extend the reduced dependence approximation to these scenarios enabling us to implement a practical and efficient non-simulation based approach for spatial data analysis within the variational Bayesian framework. The methodology is illustrated through application to simulated data and example images. The supplemental materials include our C++ source code for computing the approximate normalizing constant and simulation studies.

Mean-field descriptions of collective migration with strong adhesion

Random walk models based on an exclusion process with contact effects are often used to represent collective migration where individual agents are affected by agent-to-agent adhesion. Traditional mean field representations of these processes take the form of a nonlinear diffusion equation which, for strong adhesion, does not predict the averaged discrete behavior. We propose an alternative suite of mean-field representations, showing that collective migration with strong adhesion can be accurately represented using a moment closure approach.

Symmetry of in-shoe force signatures in athletes under field conditions : an application for injury prevention

The aetiology behind overuse injuries such as stress fractures is complex and multi-factorial. In sporting events where the loading is likely to be uneven (e.g. hurdling and jumps), research has suggested that the frequency of stress fractures seems to favour the athlete’s dominant limb. The tendency for an individual to have a preferred limb for voluntary motor acts makes limb selection a possible factor behind the development of unilateral overuse injuries, particularly when repeatedly used during high loading activities. The event of sprint hurdling is well suited for the study of loading asymmetry as the hurdling technique is repetitive and the limb movement asymmetrical. Of relevance to this study is the high incidence of Navicular Stress Fractures (NSF) in hurdlers, with suggestions there is a tendency for the fracture to develop in the trail leg foot, although this is not fully accepted. The Ground Reaction Force (GRF) with each foot contact is influenced by the hurdle action, with research finding step-to-step loading variations. However, it is unknown if this loading asymmetry extends to individual forefoot joints, thereby influencing stress fracture development. The first part of the study involved a series of investigations using a commercially available matrix style in-shoe sensor system (FscanTM, Tekscan Inc.). The suitability of insole sensor systems and custom made discrete sensors for use in hurdling-related training activities was assessed. The methodology used to analyse foot loading with each technology was investigated. The insole and discrete sensors systems tested proved to be unsuitable for use during full pace hurdling. Instead, a running barrier task designed to replicate the four repetitive foot contacts present during hurdling was assessed. This involved the clearance of a series of 6 barriers (low training hurdles), place in a straight line, using 4 strides between each. The second part of the study involved the analysis of "inter-limb" and "within foot loading asymmetries" using stance duration as well as vertical GRF under the Hallux (T1), the first metatarsal head (M1) and the central forefoot peak pressure site (M2), during walking, running, and running with barrier clearances. The contribution to loading asymmetry that each of the four repetitive foot contacts made during a series of barrier clearances was also assessed. Inter-limb asymmetry, in forefoot loading, occurred at discrete forefoot sites in a non-uniform manner across the three gait conditions. When the individual barrier foot contacts were compared, the stance duration was asymmetrical and the proportion of total forefoot load at M2 was asymmetrical. There were no significant differences between the proportion of forefoot load at M1, compared to M2; for any of the steps involved in the barrier clearance. A case study testing experimental (discrete) sensors during full pace sprinting and hurdling found that during both gait conditions, the trail limb experienced the greater vertical GRF at M1 and M2. During full pace hurdling, increased stance duration and vertical loading was a characteristic of the trail limb hurdle foot contacts. Commercially available in-shoe systems are not suitable for on field assessment of full pace hurdling. For the use of discrete sensor technology to become commonplace in the field, more robust sensors need to be developed.

Detection of banana bunchy top virus in virus-infected plants using polymerase chain reaction

Polymerase chain reaction (PCR) was developed for the detection of Banana bunchy top virus (BBTV) at maximum after 210 min and at minimum after 90 min using Pc-1 and Pc-2, respectively. PCR detection of BBTV in crude sap indicated that the freezing of banana tissue in liquid nitrogen (LN2) before extraction was more effective than using sand as the extraction technique. BBTV was also detected using PCR assay in 69 healthy and diseased plants using Na-PO4 buffer containing 1 % SDS. PCR detection of BBTV in nucleic acid extracts using seven different extraction buffers to adapt the use of PCR in routine detection in the field was studied. Results proved that BBTV was detected with high sensitivity in nucleic acid extracts more than in infectious sap. The results also suggested the common aetiology for the BBTV by the PCR reactions of BBTV in nucleic acid extracts from Australia, Burundi, Egypt, France, Gabon, Philippines and Taiwan. Results also proved a positive relation between the Egyptian-BBTV isolate and abaca bunchy top isolate from the Philippines, but there no relation was found with the Cucumber mosaic cucumovirus (CMV) isolates from Egypt and Philippines and Banana bract mosaic virus (BBMV) were found.

Effects of mobile phone distraction on drivers’ reaction times

Distraction resulting from mobile phone use whilst driving has been shown to increase the reaction times of drivers, thereby increasing the likelihood of a crash. This study compares the effects of mobile phone conversations on reaction times of drivers responding to traffic events that occur at different points in a driver’s field of view. The CARRS-Q Advanced Driving Simulator was used to test a group of young drivers on various simulated driving tasks including a traffic event that occurred within the driver’s central vision—a lead vehicle braking suddenly—and an event that occurred within the driver’s peripheral—a pedestrian entering a zebra crossing from a footpath. Thirty-two licensed drivers drove the simulator in three phone conditions: baseline (no phone conversation), and while engaged in hands-free and handheld phone conversations. The drivers were aged between 21 to 26 years and split evenly by gender. Differences in reaction times for an event in a driver’s central vision were not statistically significant across phone conditions, probably due to a lower speed selection by the distracted drivers. In contrast, the reaction times to detect an event that originated in a distracted driver’s peripheral vision were more than 50% longer compared to the baseline condition. A further statistical analysis revealed that deterioration of reaction times to an event in the peripheral vision was greatest for distracted drivers holding a provisional licence. Many critical events originate in a driver’s periphery, including vehicles, bicyclists, and pedestrians emerging from side streets. A reduction in the ability to detect these events while distracted presents a significant safety concern that must be addressed.

A finite volume method with linearisation in time for the solution of advection-reaction-diffusion systems

The numerical solution in one space dimension of advection--reaction--diffusion systems with nonlinear source terms may invoke a high computational cost when the presently available methods are used. Numerous examples of finite volume schemes with high order spatial discretisations together with various techniques for the approximation of the advection term can be found in the literature. Almost all such techniques result in a nonlinear system of equations as a consequence of the finite volume discretisation especially when there are nonlinear source terms in the associated partial differential equation models. This work introduces a new technique that avoids having such nonlinear systems of equations generated by the spatial discretisation process when nonlinear source terms in the model equations can be expanded in positive powers of the dependent function of interest. The basis of this method is a new linearisation technique for the temporal integration of the nonlinear source terms as a supplementation of a more typical finite volume method. The resulting linear system of equations is shown to be both accurate and significantly faster than methods that necessitate the use of solvers for nonlinear system of equations.

Structure- and composition-dependent electron field emission from nitrogenated carbon nanotips

The electron field emission (EFE) properties of nitrogenated carbon nanotips (NCNTPs) were studied under high-vacuum conditions. The NCNTPs were prepared in a plasma-assisted hot filament chemical vapor deposition system using CH4 and N2 as the carbon and nitrogen sources, respectively. The work functions of NCNTPs were measured using x-ray photoelectron spectroscopy. The morphological and structural properties of NCNTPs were studied by field emission scanning electron microscopy, micro-Raman spectroscopy, and x-ray photoelectron spectroscopy. The field enhancement factors of NCNTPs were calculated using relevant EFE models based on the Fowler-Nordheim approximation. Analytical characterization and modeling results were used to establish the relations between the EFE properties of NCNTPs and their morphology, structure, and composition. It is shown that the EFE properties of NCNTPs can be enhanced by the reduction of oxygen termination on the surface as well as by increasing the ratio of the NCNTP height to the radius of curvature at its top. These results also suggest that a significant amount of electrons is emitted from other surface areas besides the NCNTP tops, contrary to the common belief. The outcomes of this study advance our knowledge on the electron emission properties of carbonnanomaterials and contribute to the development of the next-generation of advanced applications in the fields of micro- and opto-electronics.

Can estimates of antimalarial efficacy from field studies be improved?

Monitoring therapeutic efficacy of antimalarial drugs is important because treatment failure rates are the primary basis for changing antimalarial treatment policy. An important aspect of efficacy studies is the use of PCR genotyping to distinguish recrudescent from new infections. The conclusions reached using this technique might be misleading if there is insufficient parasite diversity or a non-uniform haplotype frequency distribution in the study area. Statistical techniques can be used to overcome this problem, but only when data describing the haplotype frequency distribution are available. Therefore, assessing haplotype frequency and distribution should form an integral part of all studies investigating the therapeutic efficacy of antimalarial treatment regimes.

A fluorenone based low band gap solution processable copolymer for air stable and high mobility organic field effect transistors

A fluorenone based alternating copolymer (PFN-DPPF) with a furan based fused aromatic moiety has been designed and synthesized. PFN-DPPF exhibits a small band gap with a lower HOMO value. Testing this polymer semiconductor as the active layer in organic thin-film transistors results in hole mobilities as high as 0.15 cm2 V-1 s-1 in air.

Across the divide - Criminology as other: Observations on the construction of the field

In 1978 Donald Cressey commented on an emerging division in the study of crime with some scholars concentrating on the development of a “crime fi ghting coalition” and others concerned with the processes associated with “making laws, breaking laws, and the reaction to the breaking of laws” (1978: 175). Since Cressey’s paper, many others have refl ected on the distinction between criminology and the sociology of crime and deviance (Akers, 1992; Garland, 1999; Garland & Sparks, 2000; Konty, 2007). But does such a distinction actually exist? Adopting a pragmatic position, the immediate answer is yes, if we assume that these categories have substance on the basis that they are grounded in everyday beliefs, institutional preferences and research practice (Konty, 2007). Moreover, these are viable categories in that some people studying crime label themselves criminologists (or are given this label by others) while others prefer or are given the label sociologist. Of course, there are further labels that may apply to persons studying crime, which include psychologist, penologist, biologist, chemist, and so on. One could argue that such labels are unimportant, however, it remains that these categories have a practical character. For criminology and the sociology of crime in particular, scholarly discourse frames these categories as oppositional (Bader et al., 1996.; Bendle, 1989; Laub & Sampson, 1991; Sibley, 2002) and to the extent that this has occurred, the categories have social relevance.

A preconditioned numerical solver for stiff nonlinear reaction-diffusion equations with fractional Laplacians that avoids dense matrices

The numerical solution of fractional partial differential equations poses significant computational challenges in regard to efficiency as a result of the spatial nonlocality of the fractional differential operators. The dense coefficient matrices that arise from spatial discretisation of these operators mean that even one-dimensional problems can be difficult to solve using standard methods on grids comprising thousands of nodes or more. In this work we address this issue of efficiency for one-dimensional, nonlinear space-fractional reaction–diffusion equations with fractional Laplacian operators. We apply variable-order, variable-stepsize backward differentiation formulas in a Jacobian-free Newton–Krylov framework to advance the solution in time. A key advantage of this approach is the elimination of any requirement to form the dense matrix representation of the fractional Laplacian operator. We show how a banded approximation to this matrix, which can be formed and factorised efficiently, can be used as part of an effective preconditioner that accelerates convergence of the Krylov subspace iterative solver. Our approach also captures the full contribution from the nonlinear reaction term in the preconditioner, which is crucial for problems that exhibit stiff reactions. Numerical examples are presented to illustrate the overall effectiveness of the solver.

On the mathematical modeling of wound healing angiogenesis in skin as a reaction-transport process

Over the last 30 years, numerous research groups have attempted to provide mathematical descriptions of the skin wound healing process. The development of theoretical models of the interlinked processes that underlie the healing mechanism has yielded considerable insight into aspects of this critical phenomenon that remain difficult to investigate empirically. In particular, the mathematical modeling of angiogenesis, i.e., capillary sprout growth, has offered new paradigms for the understanding of this highly complex and crucial step in the healing pathway. With the recent advances in imaging and cell tracking, the time is now ripe for an appraisal of the utility and importance of mathematical modeling in wound healing angiogenesis research. The purpose of this review is to pedagogically elucidate the conceptual principles that have underpinned the development of mathematical descriptions of wound healing angiogenesis, specifically those that have utilized a continuum reaction-transport framework, and highlight the contribution that such models have made toward the advancement of research in this field. We aim to draw attention to the common assumptions made when developing models of this nature, thereby bringing into focus the advantages and limitations of this approach. A deeper integration of mathematical modeling techniques into the practice of wound healing angiogenesis research promises new perspectives for advancing our knowledge in this area. To this end we detail several open problems related to the understanding of wound healing angiogenesis, and outline how these issues could be addressed through closer cross-disciplinary collaboration.

Fourier spectral methods for fractional-in-space reaction-diffusion equations

Fractional differential equations are becoming increasingly used as a powerful modelling approach for understanding the many aspects of nonlocality and spatial heterogeneity. However, the numerical approximation of these models is demanding and imposes a number of computational constraints. In this paper, we introduce Fourier spectral methods as an attractive and easy-to-code alternative for the integration of fractional-in-space reaction-diffusion equations described by the fractional Laplacian in bounded rectangular domains ofRn. The main advantages of the proposed schemes is that they yield a fully diagonal representation of the fractional operator, with increased accuracy and efficiency when compared to low-order counterparts, and a completely straightforward extension to two and three spatial dimensions. Our approach is illustrated by solving several problems of practical interest, including the fractional Allen–Cahn, FitzHugh–Nagumo and Gray–Scott models, together with an analysis of the properties of these systems in terms of the fractional power of the underlying Laplacian operator.

Double diffusive Marangoni convection flow of electrically conducting fluid in a square cavity with chemical reaction

Double diffusive Marangoni convection flow of viscous incompressible electrically conducting fluid in a square cavity is studied in this paper by taking into consideration of the effect of applied magnetic field in arbitrary direction and the chemical reaction. The governing equations are solved numerically by using alternate direct implicit (ADI) method together with the successive over relaxation (SOR) technique. The flow pattern with the effect of governing parameters, namely the buoyancy ratio W, diffusocapillary ratio w, and the Hartmann number Ha, is investigated. It is revealed from the numerical simulations that the average Nusselt number decreases; whereas the average Sherwood number increases as the orientation of magnetic field is shifted from horizontal to vertical. Moreover, the effect of buoyancy due to species concentration on the flow is stronger than the one due to thermal buoyancy. The increase in diffusocapillary parameter, w caus

Theoretical description of time-resolved pump/probe photoemission in TaS2: a single-band DFT+DMFT(NRG) study within the quasiequilibrium approximation

In this work, we theoretically examine recent pump/probe photoemission experiments on the strongly correlated charge-density-wave insulator TaS2.We describe the general nonequilibrium many-body formulation of time-resolved photoemission in the sudden approximation, and then solve the problem using dynamical mean-field theory with the numerical renormalization group and a bare density of states calculated from density functional theory including the charge-density-wave distortion of the ion cores and spin-orbit coupling. We find a number of interesting results: (i) the bare band structure actually has more dispersion in the perpendicular direction than in the two-dimensional planes; (ii) the DMFT approach can produce upper and lower Hubbard bands that resemble those in the experiment, but the upper bands will overlap in energy with other higher energy bands; (iii) the effect of the finite width of the probe pulse is minimal on the shape of the photoemission spectra; and (iv) the quasiequilibrium approximation does not fully describe the behavior in this system.