Wavelengths in bioconvection patterns

Unicellular bottom-heavy swimming microorganisms are usually denser than the fluid in which they swim. In shallow suspensions, the bottom heaviness results in a gravitational torque that orients the cells to swim vertically upwards in the absence of fluid flow. Swimming cells thus accumulate at the upper surface to form a concentrated layer of cells. When the cell concentration is high enough, the layer overturns to form bioconvection patterns. Thin concentrated plumes of cells descend rapidly and cells return to the upper surface in wide, slowly moving upwelling plumes. When there is fluid flow, a second viscous torque is exerted on the swimming cells. The balance between the local shear flow viscous and the gravitational torques determines the cells' swimming direction, (gyrotaxis). In this thesis, the wavelengths of bioconvection patterns are studied experimentally as well as theoretically as follow; First, in aquasystem it is rare to find one species lives individually and when they swim they can form complex patterns. Thus, a protocol for controlled experiments to mix two species of swimming algal cells of \emph{C. rienhardtii} and \emph{C. augustae} is systematically described and images of bioconvection patterns are captured. A method for analysing images using wavelets and extracting the local dominant wavelength in spatially varying patterns is developed. The variation of the patterns as a function of the total concentration and the relative concentration between two species is analysed. Second, the linear stability theory of bioconvection for a suspension of two mixed species is studied. The dispersion relationship is computed using Fourier modes in order to calculate the neutral curves as a function of wavenumbers $k$ and $m$. The neutral curves are plotted to compare the instability onset of the suspension of the two mixed species with the instability onset of each species individually. This study could help us to understand which species contributes the most in the process of pattern formation. Finally, predicting the most unstable wavelength was studied previously around a steady state equilibrium situation. Since assuming steady state equilibrium contradicts with reality, the pattern formation in a layer of finite depth of an evolving basic state is studied using the nonnormal modes approach. The nonnormal modes procedure identifies the optimal initial perturbation that can be obtained for a given time $t$ as well as a given set of parameters and wavenumber $k$. Then, we measure the size of the optimal perturbation as it grows with time considering a range of wavenumbers for the same set of parameters to be able to extract the most unstable wavelength.

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.

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.

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.

Novel solutions for a model of wound healing angiogenesis

We prove the existence of novel, shock-fronted travelling wave solutions to a model of wound healing angiogenesis studied in Pettet et al (2000 IMA J. Math. App. Med. 17 395–413) assuming two conjectures hold. In the previous work, the authors showed that for certain parameter values, a heteroclinic orbit in the phase plane representing a smooth travelling wave solution exists. However, upon varying one of the parameters, the heteroclinic orbit was destroyed, or rather cut-off, by a wall of singularities in the phase plane. As a result, they concluded that under this parameter regime no travelling wave solutions existed. Using techniques from geometric singular perturbation theory and canard theory, we show that a travelling wave solution actually still exists for this parameter regime. We construct a heteroclinic orbit passing through the wall of singularities via a folded saddle canard point onto a repelling slow manifold. The orbit leaves this manifold via the fast dynamics and lands on the attracting slow manifold, finally connecting to its end state. This new travelling wave is no longer smooth but exhibits a sharp front or shock. Finally, we identify regions in parameter space where we expect that similar solutions exist. Moreover, we discuss the possibility of more exotic solutions.

Canards in advection-reaction-diffusion systems in one spatial dimension

This thesis contains a mathematical investigation of the existence of travelling wave solutions to singularly perturbed advection-reaction-diffusion models of biological processes. An enhanced mathematical understanding of these solutions and models is gained via the identification of canards (special solutions of fast/slow dynamical systems) and their role in the existence of the most biologically relevant, shock-like solutions. The analysis focuses on two existing models. A new proof of existence of a whole family of travelling waves is provided for a model describing malignant tumour invasion, while new solutions are identified for a model describing wound healing angiogenesis.

Butterfly catastrophe for fronts in a three-component reaction–diffusion system

We study the dynamics of front solutions in a three-component reaction–diffusion system via a combination of geometric singular perturbation theory, Evans function analysis, and center manifold reduction. The reduced system exhibits a surprisingly complicated bifurcation structure including a butterfly catastrophe. Our results shed light on numerically observed accelerations and oscillations and pave the way for the analysis of front interactions in a parameter regime where the essential spectrum of a single front approaches the imaginary axis asymptotically.

Characterization of alluvial formation by stochastic modelling of paleo-fluvial processes: The concept and method

Modelling fluvial processes is an effective way to reproduce basin evolution and to recreate riverbed morphology. However, due to the complexity of alluvial environments, deterministic modelling of fluvial processes is often impossible. To address the related uncertainties, we derive a stochastic fluvial process model on the basis of the convective Exner equation that uses the statistics (mean and variance) of river velocity as input parameters. These statistics allow for quantifying the uncertainty in riverbed topography, river discharge and position of the river channel. In order to couple the velocity statistics and the fluvial process model, the perturbation method is employed with a non-stationary spectral approach to develop the Exner equation as two separate equations: the first one is the mean equation, which yields the mean sediment thickness, and the second one is the perturbation equation, which yields the variance of sediment thickness. The resulting solutions offer an effective tool to characterize alluvial aquifers resulting from fluvial processes, which allows incorporating the stochasticity of the paleoflow velocity.

Constraining the low-energy pion electromagnetic form factor with space-like and phase of time-like data

The Taylor coefficients c and d of the EM form factor of the pion are constrained using analyticity, knowledge of the phase of the form factor in the time-like region, 4m(pi)(2) <= t <= t(in) and its value at one space-like point, using as input the (g - 2) of the muon. This is achieved using the technique of Lagrange multipliers, which gives a transparent expression for the corresponding bounds. We present a detailed study of the sensitivity of the bounds to the choice of time-like phase and errors present in the space-like data, taken from recent experiments. We find that our results constrain c stringently. We compare our results with those in the literature and find agreement with the chiral perturbation-theory results for c. We obtain d similar to O(10) GeV-6 when c is set to the chiral perturbation-theory values.

Conductivity Scale In Disordered-Systems

A reanalysis of the correction to the Boltzmann conductivity due to maximally crossed graphs for degenerate bands explains why the conductivity scale in many-valley semiconductors is an order of magnitude higher than Mott's "minimum metallic conductivity." With the use of a reasonable assumption for the Boltzmann mean free path, the lowest-order perturbation theory is seen to give a remarkably good, semiquantitative, description of the conductivity variation in both uncompensated doped semiconductors and amorphous alloys.

The onset of warps in Spitzer observations of edge-on spiral galaxies

We analyse warps in the nearby edge-on spiral galaxies observed in the Spitzer/Infrared Array Camera (IRAC)4.5-mu m band. In our sample of 24 galaxies, we find evidence of warp in 14 galaxies. We estimate the observed onset radii for the warps in a subsample of 10 galaxies. The dark matter distribution in each of these galaxies are calculated using the mass distribution derived from the observed light distribution and the observed rotation curves. The theoretical predictions of the onset radii for the warps are then derived by applying a self-consistent linear response theory to the obtained mass models for six galaxies with rotation curves in the literature. By comparing the observed onset radii to the theoretical ones, we find that discs with constant thickness can not explain the observations; moderately flaring discs are needed. The required flaring is consistent with the observations. Our analysis shows that the onset of warp is not symmetric in our sample of galaxies. We define a new quantity called the onset-asymmetry index and study its dependence on galaxy properties. The onset asymmetries in warps tend to be larger in galaxies with smaller dis scalelengths. We also define and quantify the global asymmetry in the stellar light distribution, that we call the edge-on asymmetry in edge-on galaxies. It is shown that in most cases the onset asymmetry in warp is actually anticorrelated with the measured edge-on asymmetry in our sample of edge-on galaxies and this could plausibly indicate that the surrounding dark matter distribution is asymmetric.

Theoretical investigations of the two-bond proton-carbon-13 coupling constants. Angular variations of the couplings involving carboxyl carbon

A formulation has been developed using perturbation theory to evaluate the π-contribution to the nuclear spin coupling constants involving nuclei at least one of which is an unsaturated center. This fromulation accounts for the π-contribution in terms of the core polarization and one-center exchange at the π-center. The formulation developed together with the Dirac vector model and Penney-Dirac bond-order formalisms was employed to calculate the geminal (two-bond) proton coupling constants of carboxyl carbons in α-disubstituted acetic acids. The calculated coupling constants were found to have an orientational dependence. The results of the calculation are in good agreement with the experimental values.