Simulations of X and Y Retinal Ganglion Cell Behavior with a Nonlinear Push-pull Model of Spatiotemporal Retinal Processing

This article describes a nonlinear model of neural processing in the vertebrate retina, comprising model photoreceptors, model push-pull bipolar cells, and model ganglion cells. Previous analyses and simulations have shown that with a choice of parameters that mimics beta cells, the model exhibits X-like linear spatial summation (null response to contrast-reversed gratings) in spite of photoreceptor nonlinearities; on the other hand, a choice of parameters that mimics alpha cells leads to Y-like frequency doubling. This article extends the previous work by showing that the model can replicate qualitatively many of the original findings on X and Y cells with a fixed choice of parameters. The results generally support the hypothesis that X and Y cells can be seen as functional variants of a single neural circuit. The model also suggests that both depolarizing and hyperpolarizing bipolar cells converge onto both ON and OFF ganglion cell types. The push-pull connectivity enables ganglion cells to remain sensitive to deviations about the mean output level of nonlinear photoreceptors. These and other properties of the push-pull model are discussed in the general context of retinal processing of spatiotemporal luminance patterns.

Computational modeling of defects in nanoscale device materials

This PhD thesis concerns the computational modeling of the electronic and atomic structure of point defects in technologically relevant materials. Identifying the atomistic origin of defects observed in the electrical characteristics of electronic devices has been a long-term goal of first-principles methods. First principles simulations are performed in this thesis, consisting of density functional theory (DFT) supplemented with many body perturbation theory (MBPT) methods, of native defects in bulk and slab models of In0.53Ga0.47As. The latter consist of (100) - oriented surfaces passivated with A12O3. Our results indicate that the experimentally extracted midgap interface state density (Dit) peaks are not the result of defects directly at the semiconductor/oxide interface, but originate from defects in a more bulk-like chemical environment. This conclusion is reached by considering the energy of charge transition levels for defects at the interface as a function of distance from the oxide. Our work provides insight into the types of defects responsible for the observed departure from ideal electrical behaviour in III-V metal-oxidesemiconductor (MOS) capacitors. In addition, the formation energetics and electron scattering properties of point defects in carbon nanotubes (CNTs) are studied using DFT in conjunction with Green’s function based techniques. The latter are applied to evaluate the low-temperature, low-bias Landauer conductance spectrum from which mesoscopic transport properties such as the elastic mean free path and localization length of technologically relevant CNT sizes can be estimated from computationally tractable CNT models. Our calculations show that at CNT diameters pertinent to interconnect applications, the 555777 divacancy defect results in increased scattering and hence higher electrical resistance for electron transport near the Fermi level.

Evaluating functional network inference using simulations of complex biological systems.

MOTIVATION: Although many network inference algorithms have been presented in the bioinformatics literature, no suitable approach has been formulated for evaluating their effectiveness at recovering models of complex biological systems from limited data. To overcome this limitation, we propose an approach to evaluate network inference algorithms according to their ability to recover a complex functional network from biologically reasonable simulated data. RESULTS: We designed a simulator to generate data representing a complex biological system at multiple levels of organization: behaviour, neural anatomy, brain electrophysiology, and gene expression of songbirds. About 90% of the simulated variables are unregulated by other variables in the system and are included simply as distracters. We sampled the simulated data at intervals as one would sample from a biological system in practice, and then used the sampled data to evaluate the effectiveness of an algorithm we developed for functional network inference. We found that our algorithm is highly effective at recovering the functional network structure of the simulated system-including the irrelevance of unregulated variables-from sampled data alone. To assess the reproducibility of these results, we tested our inference algorithm on 50 separately simulated sets of data and it consistently recovered almost perfectly the complex functional network structure underlying the simulated data. To our knowledge, this is the first approach for evaluating the effectiveness of functional network inference algorithms at recovering models from limited data. Our simulation approach also enables researchers a priori to design experiments and data-collection protocols that are amenable to functional network inference.

Enhancing the Hodgkin-Huxley Equations: Simulations Based on the First Publication in the Biophysical Journal.

The experiments in the Cole and Moore article in the first issue of the Biophysical Journal provided the first independent experimental confirmation of the Hodgkin-Huxley (HH) equations. A log-log plot of the K current versus time showed that raising the HH variable n to the sixth power provided the best fit to the data. Subsequent simulations using n(6) and setting the resting potential at the in vivo value simplifies the HH equations by eliminating the leakage term. Our article also reported that the K current in response to a depolarizing step to ENa was delayed if the step was preceded by a hyperpolarization. While the interpretation of this phenomenon in the article was flawed, subsequent simulations show that the effect completely arises from the original HH equations.

Numerical Simulations and Coherence Properties of Supercontinuum Generation in Photonic Crystal and Tapered Optical Fibers

Numerical simulations have been used to study broad-band supercontinuum generation in optical fibers with dispersion and nonlinearity characteristics typical and photonic crystal or tapered fibers structures. The simulations include optical shock and Raman nonlinearity terms, with quantum noise taken into account phenomenologically by including in the input field a noise seed of one photon per mode with random phase. For input pulses of 150-fs duration injected in the anomalous dispersion regime, the effect of modulational instability is shown to lead to severe temporal jitter in the output, and associated fluctuations in the spectral amplitude and phase across the generated supercontinuum. The spectral phase fluctuations are quantified by performing multiple simulations and calculating both the standard deviation of the phase and, more rigorously, the degree of first-order coherence as a function of wavelength across the spectrum. By performing simulations over a range of input pulse durations and wavelengths, we can identify the conditions under which coherent supercontinua with a well-defined spectral phase are generated.

Cross-correlation frequency resolved optical gating analysis of broadband continuum generation in photonic crystal fiber: Simulations and experiments

Numerical simulations are used to study the temporal and spectral characteristics of broadband supercontinua generated in photonic crystal fiber. In particular, the simulations are used to follow the evolution with propagation distance of the temporal intensity, the spectrum, and the cross-correlation frequency resolved optical gating (XFROG) trace. The simulations allow several important physical processes responsible for supercontinuum generation to be identified and, moreover, illustrate how the XFROG trace provides an intuitive means of interpreting correlated temporal and spectral features of the supercontinuum. Good qualitative agreement with preliminary XFROG measurements is observed. © 2002 Optical Society of America.

Measurements and simulations of noise imposed on supercontinuum generated in microstructure fiber

We present experimental and theoretical investigations of the highly nonlinear and broadband noise that exists on supercontinuum spectra generated from launching femtosecond Ti:Sapphire pulses into microstructure fiber.

Amplitude noise on supercontinuum generated in microstructure fiber: Measurements and simulations

Supercontinua generated in microstructure fiber can exhibit significant excess amplitude noise. We present experimental and numerical studies of the origins of this excess noise and its dependence on the input laser pulse parameters.

A multi-scale atomistic-continuum modelling of crack propagation in a two-dimensional macroscopic plate

A novel multi-scale seamless model of brittle-crack propagation is proposed and applied to the simulation of fracture growth in a two-dimensional Ag plate with macroscopic dimensions. The model represents the crack propagation at the macroscopic scale as the drift-diffusion motion of the crack tip alone. The diffusive motion is associated with the crack-tip coordinates in the position space, and reflects the oscillations observed in the crack velocity following its critical value. The model couples the crack dynamics at the macroscales and nanoscales via an intermediate mesoscale continuum. The finite-element method is employed to make the transition from the macroscale to the nanoscale by computing the continuum-based displacements of the atoms at the boundary of an atomic lattice embedded within the plate and surrounding the tip. Molecular dynamics (MD) simulation then drives the crack tip forward, producing the tip critical velocity and its diffusion constant. These are then used in the Ito stochastic calculus to make the reverse transition from the nanoscale back to the macroscale. The MD-level modelling is based on the use of a many-body potential. The model successfully reproduces the crack-velocity oscillations, roughening transitions of the crack surfaces, as well as the macroscopic crack trajectory. The implications for a 3-D modelling are discussed.

Computational issues in the modelling of materials-based manufacturing processes

The manufacture of materials products involves the control of a range of interacting physical phenomena. The material to be used is synthesised and then manipulated into some component form. The structure and properties of the final component are influenced by both interactions of continuum-scale phenomena and those at an atomistic-scale level. Moreover, during the processing phase there are some properties that cannot be measured (typically the liquid-solid phase change). However, it seems there is a potential to derive properties and other features from atomistic-scale simulations that are of key importance at the continuum scale. Some of the issues that need to be resolved in this context focus upon computational techniques and software tools facilitating: (i) the multiphysics modeling at continuum scale; (ii) the interaction and appropriate degrees of coupling between the atomistic through microstructure to continuum scale; and (iii) the exploitation of high-performance parallel computing power delivering simulation results in a practical time period. This paper discusses some of the attempts to address each of the above issues, particularly in the context of materials processing for manufacture.

Surface physics at the nano-scale via scanning probe microscopy and molecular dynamics simulations

Probe-based scanning microscopes, such as the STM and the AFM, are used to obtain the topographical and electronic structure maps of material surfaces, and to modify their morphologies on nanoscopic scales. They have generated new areas of research in condensed matter physics and materials science. We will review some examples from the fields of experimental nano-mechanics, nano-electronics and nano-magnetism. These now form the basis of the emerging field of Nano-technology. A parallel development has been brought about in the field of Computational Nano-science, using quantum-mechanical techniques and computer-based numerical modelling, such as the Molecular Dynamics (MD) simulation method. We will report on the simulation of nucleation and growth of nano-phase films on supporting substrates. Furthermore, a theoretical modelling of the formation of STM images of metallic clusters on metallic substrates will also be discussed within the non-equilibrium Keldysh Green function method to study the effects of coherent tunnelling through different atomic orbitals in a tip-sample geometry.

Magnetic modelling and magneto-hydro-dynamic simulations of an aluminium production electrolytic cell

A simulation of the motion of molten aluminium inside an electrolytic cell is presented. Since the driving term of the aluminium motion is the Lorentz (j × B) body force acting within the fluid,this problem involves the solution of the magneto-hydro-dynamic equations. Different solver modules for the magnetic field computation and for the fluid motion simulation are coupled together. The interactions of all these are presented and discussed.

Correlation of solder paste rheology with computational simulations of the stencil printing process

Soldering technologies continue to evolve to meet the demands of the continuous miniaturisation of electronic products, particularly in the area of solder paste formulations used in the reflow soldering of surface mount devices. Stencil printing continues to be a leading process used for the deposition of solder paste onto printed circuit boards (PCBs) in the volume production of electronic assemblies, despite problems in achieving a consistent print quality at an ultra-fine pitch. In order to eliminate these defects a good understanding of the processes involved in printing is important. Computational simulations may complement experimental print trials and paste characterisation studies, and provide an extra dimension to the understanding of the process. The characteristics and flow properties of solder pastes depend primarily on their chemical and physical composition and good material property data is essential for meaningful results to be obtained by computational simulation.This paper describes paste characterisation and computational simulation studies that have been undertaken through the collaboration of the School of Aeronautical, Mechanical and Manufacturing Engineering at Salford University and the Centre for Numerical Modelling and Process Analysis at the University of Greenwich. The rheological profile of two different paste formulations (lead and lead-free) for sub 100 micron flip-chip devices are tested and applied to computational simulations of their flow behaviour during the printing process.

Low Reynolds number turbulence models for accurate thermal simulations of electronic components

The electronics industry and the problems associated with the cooling of microelectronic equipment are developing rapidly. Thermal engineers now find it necessary to consider the complex area of equipment cooling at some level. This continually growing industry also faces heightened pressure from consumers to provide electronic product miniaturization, which in itself increases the demand for accurate thermal management predictions to assure product reliability. Computational fluid dynamics (CFD) is considered a powerful and almost essential tool for the design, development and optimization of engineering applications. CFD is now widely used within the electronics packaging design community to thermally characterize the performance of both the electronic component and system environment. This paper discusses CFD results for a large variety of investigated turbulence models. Comparison against experimental data illustrates the predictive accuracy of currently used models and highlights the growing demand for greater mathematical modelling accuracy with regards to thermal characterization. Also a newly formulated low Reynolds number (i.e. transitional) turbulence model is proposed with emphasis on hybrid techniques.