Planar radial spots in a three-component FitzHugh-Nagumo system

Localized planar patterns arise in many reaction-diffusion models. Most of the paradigm equations that have been studied so far are two-component models. While stationary localized structures are often found to be stable in such systems, travelling patterns either do not exist or are found to be unstable. In contrast, numerical simulations indicate that localized travelling structures can be stable in three-component systems. As a first step towards explaining this phenomenon, a planar singularly perturbed three-component reaction-diffusion system that arises in the context of gas-discharge systems is analysed in this paper. Using geometric singular perturbation theory, the existence and stability regions of radially symmetric stationary spot solutions are delineated and, in particular, stable spots are shown to exist in appropriate parameter regimes. This result opens up the possibility of identifying and analysing drift and Hopf bifurcations, and their criticality, from the stationary spots described here.

The renormalization group: A perturbation method for the graduate curriculum

In this paper the renormalization group (RG) method of Chen, Goldenfeld, and Oono [Phys. Rev. Lett., 73 (1994), pp.1311-1315; Phys. Rev. E, 54 (1996), pp.376-394] is presented in a pedagogical way to increase its visibility in applied mathematics and to argue favorably for its incorporation into the corresponding graduate curriculum.The method is illustrated by some linear and nonlinear singular perturbation problems. Key word. © 2012 Society for Industrial and Applied Mathematics.

A survey in mathematics for industry: Two-timing and matched asymptotic expansions for singular perturbation problems

Following the derivation of amplitude equations through a new two-time-scale method [O'Malley, R. E., Jr. & Kirkinis, E (2010) A combined renormalization group-multiple scale method for singularly perturbed problems. Stud. Appl. Math. 124, 383-410], we show that a multi-scale method may often be preferable for solving singularly perturbed problems than the method of matched asymptotic expansions. We illustrate this approach with 10 singularly perturbed ordinary and partial differential equations. © 2011 Cambridge University Press.

From robots to animals : Virtual fences for controlling cattle

We consider the problem of monitoring and controlling the position of herd animals, and view animals as networked agents with natural mobility but not strictly controllable. By exploiting knowledge of individual and herd behavior we would like to apply a vast body of theory in robotics and motion planning to achieving the constrained motion of a herd. In this paper we describe the concept of a virtual fence which applies a stimulus to an animal as a function of its pose with respect to the fenceline. Multiple fence lines can define a region, and the fences can be static or dynamic. The fence algorithm is implemented by a small position-aware computer device worn by the animal, which we refer to as a Smart Collar.We describe a herd-animal simulator, the Smart Collar hardware and algorithms for tracking and controlling animals as well as the results of on-farm experiments with up to ten Smart Collars.

Dynamic virtual fences for controlling cows

A virtual fence is created by applying an aversive stimulus to an animal when it approaches a predefined boundary. It is implemented by a small animal-borne computer system with a GPS receiver. This approach allows the implementation of virtual paddocks inside a normal physically-fenced paddock. Since the fence lines are virtual they can be moved by programming to meet the needs of animal or land management. This approach enables us to consider animals as agents with natural mobility that are controllable and to apply a vast body of theory in motion planning. In this paper we describe a herd-animal simulator and physical experiments conducted on a small herd of 10 animals using a Smart Collar. The Smart Collar consists of a GPS, PDA, wireless networking and a sound amplifier. We describe a motion planning algorithm that can move a virtual paddock subject to landscape constraints which is suitable for mustering cows. We present simulation results and data from experiments with 8 cows equipped with Smart Collars.

Improving road safety through deterrence-based initiatives : a review of research

The efficacy of road safety countermeasures to deter motorists from engaging in illegal behaviours is extremely important when considering the personal and economic impact of road accidents on the community. Within many countries, deterrence theory has remained a cornerstone to criminology and criminal justice policy, particularly within the field of road safety, as policy makers and enforcement agencies attempt to increase perceptions regarding the certainty, severity and swiftness of sanctions for those who engage in illegal motoring behaviours. Using the Australian experience (particularly the tremendous amount of research into drink driving), the current paper reviews the principles underpinning deterrence theory, the utilisation of the approach within some contemporary road safety initiatives (e.g., Random Breath Testing) as well as highlights some methods to enhance a deterrent effect. The paper also provides direction for future deterrence-based research, in particular, considering the powerful impact of non-legal sanctions, punishment avoidance as well as creating culturally embedded behavioural change.

Pulse dynamics in a three-component system : existence analysis

In this article, we analyze the three-component reaction-diffusion system originally developed by Schenk et al. (PRL 78:3781–3784, 1997). The system consists of bistable activator-inhibitor equations with an additional inhibitor that diffuses more rapidly than the standard inhibitor (or recovery variable). It has been used by several authors as a prototype three-component system that generates rich pulse dynamics and interactions, and this richness is the main motivation for the analysis we present. We demonstrate the existence of stationary one-pulse and two-pulse solutions, and travelling one-pulse solutions, on the real line, and we determine the parameter regimes in which they exist. Also, for one-pulse solutions, we analyze various bifurcations, including the saddle-node bifurcation in which they are created, as well as the bifurcation from a stationary to a travelling pulse, which we show can be either subcritical or supercritical. For two-pulse solutions, we show that the third component is essential, since the reduced bistable two-component system does not support them. We also analyze the saddle-node bifurcation in which two-pulse solutions are created. The analytical method used to construct all of these pulse solutions is geometric singular perturbation theory, which allows us to show that these solutions lie in the transverse intersections of invariant manifolds in the phase space of the associated six-dimensional travelling wave system. Finally, as we illustrate with numerical simulations, these solutions form the backbone of the rich pulse dynamics this system exhibits, including pulse replication, pulse annihilation, breathing pulses, and pulse scattering, among others.

Computational study of methyl derivatives of ammonia borane for hydrogen storage

The structures and thermodynamic properties of methyl derivatives of ammonia–borane (BH3NH3, AB) have been studied with the frameworks of density functional theory and second-order Møller–Plesset perturbation theory. It is found that, with respect to pure AB, methyl ammonia–boranes show higher complexation energies and lower reaction enthalpies for the release of H2, together with a slight increment of the activation barrier. These results indicate that the methyl substitution can enhance the reversibility of the system and prevent the formation of BH3/NH3, but no enhancement of the release rate of H2 can be expected.

Deposition of hydrocarbon molecules on diamond (001) surfaces : atomic scale modeling

The impact induced chemisorption of hydrocarbon molecules (CH3 and CH2) on H-terminated diamond (001)-(2x1) surface was investigated by molecular dynamics simulation using the many-body Brenner potential. The deposition dynamics of the CH3 radical at impact energies of 0.1-50 eV per molecule was studied and the energy threshold for chemisorption was calculated. The impact-induced decomposition of hydrogen atoms and the dimer opening mechanism on the surface was investigated. Furthermore, the probability for dimer opening event induced by chemisorption of CH, was simulated by randomly varying the impact position as well as the orientation of the molecule relative to the surface. Finally, the energetic hydrocarbons were modeled, slowing down one after the other to simulate the initial fabrication of diamond-like carbon (DLC) films. The structure characteristic in synthesized films with different hydrogen flux was studied. Our results indicate that CH3, CH2 and H are highly reactive and important species in diamond growth. Especially, the fraction of C-atoms in the film having sp(3) hybridization will be enhanced in the presence of H atoms, which is in good agreement with experimental observations. (C) 2002 Elsevier Science B.V. All rights reserved.

Memory effect in the deposition of C20 fullerenes on a diamond surface

In this paper, the deposition of C-20 fullerenes on a diamond (001)-(2x1) surface and the fabrication of C-20 thin film at 100 K were investigated by a molecular dynamics (MD) simulation using the many-body Brenner bond order potential. First, we found that the collision dynamic of a single C-20 fullerene on a diamond surface was strongly dependent on its impact energy. Within the energy range 10-45 eV, the C-20 fullerene chemisorbed on the surface retained its free cage structure. This is consistent with the experimental observation, where it was called the memory effect in "C-20-type" films [P. Melion , Int. J. Mod. B 9, 339 (1995); P. Milani , Cluster Beam Synthesis of Nanostructured Materials (Springer, Berlin, 1999)]. Next, more than one hundred C-20 (10-25 eV) were deposited one after the other onto the surface. The initial growth stage of C-20 thin film was observed to be in the three-dimensional island mode. The randomly deposited C-20 fullerenes stacked on diamond surface and acted as building blocks forming a polymerlike structure. The assembled film was also highly porous due to cluster-cluster interaction. The bond angle distribution and the neighbor-atom-number distribution of the film presented a well-defined local order, which is of sp(3) hybridization character, the same as that of a free C-20 cage. These simulation results are again in good agreement with the experimental observation. Finally, the deposited C-20 film showed high stability even when the temperature was raised up to 1500 K.

Energy dependence of methyl-radical adsorption on diamond (001)-(2×1) surface

The deposition of hyperthermal CH3 on diamond (001)-(2×1) surface at room temperature has been studied by means of molecular dynamics simulation using the many-body hydrocarbon potential. The energy threshold effect has been observed. That is, with fixed collision geometry, chemisorption can occur only when the incident energy of CH3 is above a critical value (Eth). Increasing the incident energy, dissociation of hydrogen atoms from the incident molecule was observed. The chemisorption probability of CH3 as a function of its incident energy was calculated and compared with that of C2H2. We found that below 10 eV, the chemisorption probability of C2H2 is much lower than that of CH3 on the same surface. The interesting thing is that it is even lower than that of CH3 on a hydrogen covered surface at the same impact energy. It indicates that the reactive CH3 molecule is the more important species than C2H2 in diamond synthesis at low energy, which is in good agreement with the experimental observation.

Existence of travelling wave solutions for a model of tumour invasion

The existence of travelling wave solutions to a haptotaxis dominated model is analysed. A version of this model has been derived in Perumpanani et al. (1999) to describe tumour invasion, where diffusion is neglected as it is assumed to play only a small role in the cell migration. By instead allowing diffusion to be small, we reformulate the model as a singular perturbation problem, which can then be analysed using geometric singular perturbation theory. We prove the existence of three types of physically realistic travelling wave solutions in the case of small diffusion. These solutions reduce to the no diffusion solutions in the singular limit as diffusion as is taken to zero. A fourth travelling wave solution is also shown to exist, but that is physically unrealistic as it has a component with negative cell population. The numerical stability, in particular the wavespeed of the travelling wave solutions is also discussed.

A geometric construction of travelling wave solutions to the Keller–Segel model

We study a version of the Keller–Segel model for bacterial chemotaxis, for which exact travelling wave solutions are explicitly known in the zero attractant diffusion limit. Using geometric singular perturbation theory, we construct travelling wave solutions in the small diffusion case that converge to these exact solutions in the singular limit.

Photoelectron spectroscopic study of the oxyallyl diradical

The photoelectron spectrum of the oxyallyl (OXA) radical anion has been measured. The radical anion has been generated in the reaction of the atomic oxygen radical anion (O center dot-) with acetone. Three low-lying electronic states of OXA have been observed in the spectrum. Electronic structure calculations have been performed for the triplet states (B-3(2) and B-3(1)) of OXA and the ground doublet state ((2)A(2)) of the radical anion using density, functional theory (DFT). Spectral simulations have been carried out for the triplet statics based on the results of the DFT calculations. The simulation identifies a vibrational progression of the CCC bending mode of the B-3(2) state of OXA in the lower electron binding energy (eBE) portion of the spectrum. On top of the B-3(2) feature, however, the experimental spectrum exhibits additional photoelectron peaks whose angular distribution is distinct from that for the vibronic peaks of the B-3(2) state. Complete active space self-consistent field (CASSCF) method and second-order perturbation theory based on the CASSCF wave function (CASPT2) have been employed to study the lowest singlet state ((1)A(1)) of OXA. The simulation based on the results of these electronic structure calculations establishes that the overlapping peaks represent the vibrational ground level of the (1)A(1) state and its vibrational progression of the CO stretching mode. The A, state is the lowest electronic state of,OXA, and the electron affinity (EA) of OXA is 1.940 +/- 0.010 eV. The B-3(2) state is the first excited state with an electronic term energy of 55 +/- 2 meV. The widths of the vibronic peaks of the (X) over tilde (1)A(1) state are much broader than those of the (a) over tilde B-3(2) state, implying that the (1)A(1) state is indeed a transition state. The CASSCF and CASPT2 calculations suggest that the (1)A(1) state is at a potential maximum along the nuclear coordinate representing disrotatory motion of the two methylene groups, which leads to three-membered-ring formation, i.e., cydopropanone. The simulation of (b) over tilde B-3(1) OXA reproduces the higher eBE portion of the spectrum very well. The term energy of the B-3(1) state is 0.883 +/- 0.012 eV. Photoelectron spectroscopic measurements have also been conducted for the other ion products of the O center dot- reaction with acetone. The photoelectron imaging spectrum of the acetylcarbene (AC) radical anion exhibits a broad, structureless feature, which is assigned to the (X) over tilde (3)A '' state of AC. The ground ((2)A '') and first excited ((2)A') states of the 1-methylvinoxy (1-MVO) radical have been observed in the photoelectron spectrum of the 1-MVO ion, and their vibronic structure has been analyzed.

Bistable regime of surface magnetoplasma waves excitation at a semiconductor-metal interface

The theoretical analysis of the bistability associated with the excitation of surface magnetoplasma waves (SWs) propagating across an external magnetic field at the semiconductor-metal interface by the attenuated total reflection (ATR) method is presented. The Kretschmann-Raether configuration of the ATR method is considered, i.e. a plane electromagnetic wave is incident onto a metal surface through a coupling prism. The third-order nonlinearity of the semiconductor medium is considered in the general form using the formalism of the third-order nonlinear susceptibilities and of the perturbation theory. The examples of the nonlinear mechanisms which influence the SW propagation are given. The analytical and numerical analyses show that the realization of bistable regimes of the SW excitation is possible. The SW amplitude values providing bistability in the structure are evaluated and are reasonably low to provide the experimental observation.