Theory and simulation in neuroscience.

Modeling work in neuroscience can be classified using two different criteria. The first one is the complexity of the model, ranging from simplified conceptual models that are amenable to mathematical analysis to detailed models that require simulations in order to understand their properties. The second criterion is that of direction of workflow, which can be from microscopic to macroscopic scales (bottom-up) or from behavioral target functions to properties of components (top-down). We review the interaction of theory and simulation using examples of top-down and bottom-up studies and point to some current developments in the fields of computational and theoretical neuroscience.

Direct numerical simulation of turbulent flows over superhydrophobic surfaces with gas pockets using linearized boundary conditions

Superhydrophobic surfaces are shown to be effective for surface drag reduction under laminar regime by both experiments and simulations (see for example, Ou and Rothstein, Phys. Fluids 17:103606, 2005). However, such drag reduction for fully developed turbulent flow maintaining the Cassie-Baxter state remains an open problem due to high shear rates and flow unsteadiness of turbulent boundary layer. Our work aims to develop an understanding of mechanisms leading to interface breaking and loss of gas pockets due to interactions with turbulent boundary layers. We take advantage of direct numerical simulation of turbulence with slip and no-slip patterned boundary conditions mimicking the superhydrophobic surface. In addition, we capture the dynamics of gas-water interface, by deriving a proper linearized boundary condition taking into account the surface tension of the interface and kinematic matching of interface deformation and normal velocity conditions on the wall. We will show results from our simulations predicting the dynamical behavior of gas pocket interfaces over a wide range of dimensionless surface tensions.

Investigation of gain recovery for InAs/GaAs quantum dot semiconductor optical amplifiers by rate equation simulation

The gain recoveries in quantum dot semiconductor optical amplifiers (QD SOAs) are numerically studied by rate equation simulation. Similar to the optical pump-probe experiment, the injection of double 150 fs optical pulses is used to simulate the gain recovery of a weak continuous signal under different injection levels, inhomogeneous broadenings, detuning wavelengths, and pulse signal energies for the QD SOAs. The obtained gain recoveries are then fitted by a response function with multiple exponential terms to determine the response times. The gain recovery can be described by three exponential terms with the time constants, which can be explained as carrier relaxation from the excited state to the ground state, carrier captured by the excited state from the wetting layer, and the supply of the wetting layer carriers. The fitted lifetimes decrease with the increase of the injection currents under gain unsaturation, slightly decrease with the decrease of inhomogeneous broadening of QDs, and increase with the increase of detuning wavelength between continuous signal and pulse signal and the increase of the pulse energy.

Molecular Dynamics Simulation on Characteristics of Decomposition of Methane Hydrate in the Presence of Hydrogen Peroxide Solution

In this work, the characteristics of the decomposition of methane hydrate Structure I (SI) in the presence of hydrogen peroxide solution is investigated using the molecular dynamics simulation. The mechanism of the transformation process from the solid hydrate to the liquid is analyzed with the effect of hydrogen peroxide (HP) solution. In addition, the effect of ethylene glycol (EG) with the same molar concentration with HP on the methane hydrate dissociation is also studied. The results illustrate that both HP and EG promote well the hydrate dissociation. The work provides the important reference value for the experimental investigation into the promotion effect of HP on the hydrate dissociation.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Numerical simulation of thermal loading produced by shaped high power laser onto engine parts

Recently a new method for simulating the thermal loading on pistons of diesel engines was reported. The spatially shaped high power laser is employed as the heat source, and some preliminary experimental and numerical work was carried out. In this paper, a further effort was made to extend this simulation method to some other important engine parts such as cylinder heads. The incident Gaussian beam was transformed into concentric multi-circular patterns of specific intensity distributions, with the aid of diffractive optical elements (DOEs). By incorporating the appropriate repetitive laser pulses, the designed transient temperature fields and thermal loadings in the engine parts could be simulated. Thermal-structural numerical models for pistons and cylinder heads were built to predict the transient temperature and thermal stress. The models were also employed to find the optimal intensity distributions of the transformed laser beam that could produce the target transient temperature fields. Comparison of experimental and numerical results demonstrated that this systematic approach is effective in simulating the thermal loading on the engine parts. (C) 2009 Elsevier Ltd. All rights reserved.

Design and simulation of charge sensitive preamplifier with CMOS FET implemented as feedback capacitor C-fp

In this paper, to design a new preamplifier for optimum performances with charged-particle or heavy-ion detectors, the CMOS FET is implemented as a feedback capacitor C-fp, so that the entire system should be built only with MOSFET. This work is a revolution design because to realize an ASIC for a preamplifier circuit, the capacitor will also be included. We succeed after a simulation to maintain a rise time less than 3 ns, the output resistance less than 94 Omega and the linearity almost good.

Maxwell's Demon, Rectifiers, and the Second Law: Computer Simulation of Smoluchowski's Trapdoor

We have simulated numerically an automated Maxwell's demon inspired by Smoluchowski's ideas of 1912. Two gas chambers of equal area are connected via an opening that is covered by a trapdoor. The trapdoor can open to the left but not to the right, and is intended to rectify naturally occurring variations in density between the two chambers. Our results confirm that though the trapdoor behaves as a rectifier when large density differences are imposed by external means, it can not extract useful work from the thermal motion of the molecules when left on its own.

Exhibiting the behaviour of time-delayed systems via an extension to qualitative simulation

I. Miguel and Q. Shen. Exhibiting the behaviour of time-delayed systems via an extension to qualitative simulation. IEEE Transactions on Systems, Man and Cybernetics, Part A: Systems and Humans, 35(2):298-305, 2005.

Experimental and numerical simulation of dough kneading in filled geometries

Binding, David; Couch, M.A.; Sujatha, K.S.; Webster, M.F., (2003) 'Experimental and numerical simulation of dough kneading in filled geometries', Journal of Food Engineering 58 pp.111-123 RAE2008

A simulation-based design method to transfer surface mount RF system to flip-chip die implementation

The flip-chip technology is a high chip density solution to meet the demand for very large scale integration design. For wireless sensor node or some similar RF applications, due to the growing requirements for the wearable and implantable implementations, flip-chip appears to be a leading technology to realize the integration and miniaturization. In this paper, flip-chip is considered as part of the whole system to affect the RF performance. A simulation based design is presented to transfer the surface mount PCB board to the flip-chip die package for the RF applications. Models are built by Q3D Extractor to extract the equivalent circuit based on the parasitic parameters of the interconnections, for both bare die and wire-bonding technologies. All the parameters and the PCB layout and stack-up are then modeled in the essential parts' design of the flip-chip RF circuit. By implementing simulation and optimization, a flip-chip package is re-designed by the parameters given by simulation sweep. Experimental results fit the simulation well for the comparison between pre-optimization and post-optimization of the bare die package's return loss performance. This design method could generally be used to transfer any surface mount PCB to flip-chip package for the RF systems or to predict the RF specifications of a RF system using the flip-chip technology.

Simulation chamber studies of the atmospheric degradation of naphthalene, 1-nitronaphthalene and phthaldialdehyde

A detailed series of simulation chamber experiments has been performed on the atmospheric degradation pathways of the primary air pollutant naphthalene and two of its photooxidation products, phthaldialdehyde and 1-nitronaphthalene. The measured yields of secondary organic aerosol (SOA) arising from the photooxidation of naphthalene varied from 6-20%, depending on the concentrations of naphthalene and nitrogen oxides as well as relative humidity. A range of carbonyls, nitro-compounds, phenols and carboxylic acids were identified among the gas- and particle-phase products. On-line analysis of the chemical composition of naphthalene SOA was performed using aerosol time-of-flight mass spectrometry (ATOFMS) for the first time. The results indicate that enhanced formation of carboxylic acids may contribute to the observed increase in SOA yields at higher relative humidity. The photolysis of phthaldialdehyde and 1-nitronaphthalene was investigated using natural light at the European Photoreactor (EUPHORE) in Valencia, Spain. The photolysis rate coefficients were measured directly and used to confirm that photolysis is the major atmospheric loss process for these compounds. For phthaldialdehyde, the main gas-phase products were phthalide and phthalic anhydride. SOA yields in the range 2-11% were observed, with phthalic acid and dihydroxyphthalic acid identified among the particle phase products. The photolysis of 1-nitronaphthalene yielded nitric oxide and a naphthoxy radical which reacted to form several products. SOA yields in the range 57-71% were observed, with 1,4-naphthoquinone, 1-naphthol and 1,4-naphthalenediol identified in the particle phase. On-line analysis of the SOA generated in an indoor chamber using ATOFMS provided evidence for the formation of high-molecular-weight products. Further investigations revealed that these products are oxygenated polycyclic compounds most likely produced from the dimerization of naphthoxy radicals. These results of this work indicate that naphthalene is a potentially large source of SOA in urban areas and should be included in atmospheric models. The kinetic and mechanistic information could be combined with existing literature data to produce an overall degradation mechanism for naphthalene suitable for inclusion in photochemical models that are used to predict the effect of emissions on air quality.

Simulation of multigate SOI transistors with silicon, germanium and III-V channels

In this work by employing numerical three-dimensional simulations we study the electrical performance and short channel behavior of several multi-gate transistors based on advanced SOI technology. These include FinFETs, triple-gate and gate-all-around nanowire FETs with different channel material, namely Si, Ge, and III-V compound semiconductors, all most promising candidates for future nanoscale CMOS technologies. Also, a new type of transistor called “junctionless nanowire transistor” is presented and extensive simulations are carried out to study its electrical characteristics and compare with the conventional inversion- and accumulation-mode transistors. We study the influence of device properties such as different channel material and orientation, dimensions, and doping concentration as well as quantum effects on the performance of multi-gate SOI transistors. For the modeled n-channel nanowire devices we found that at very small cross sections the nanowires with silicon channel are more immune to short channel effects. Interestingly, the mobility of the channel material is not as significant in determining the device performance in ultrashort channels as other material properties such as the dielectric constant and the effective mass. Better electrostatic control is achieved in materials with smaller dielectric constant and smaller source-to-drain tunneling currents are observed in channels with higher transport effective mass. This explains our results on Si-based devices. In addition to using the commercial TCAD software (Silvaco and Synopsys TCAD), we have developed a three-dimensional Schrödinger-Poisson solver based on the non-equilibrium Green’s functions formalism and in the framework of effective mass approximation. This allows studying the influence of quantum effects on electrical performance of ultra-scaled devices. We have implemented different mode-space methodologies in our 3D quantum-mechanical simulator and moreover introduced a new method to deal with discontinuities in the device structures which is much faster than the coupled-mode-space approach.

First principles simulation of amorphous silicon bulk, interfaces, and nanowires for photovoltaics

Amorphous silicon has become the material of choice for many technologies, with major applications in large area electronics: displays, image sensing and thin film photovoltaic cells. This technology development has occurred because amorphous silicon is a thin film semiconductor that can be deposited on large, low cost substrates using low temperature. In this thesis, classical molecular dynamics and first principles DFT calculations have been performed to generate structural models of amorphous and hydrogenated amorphous silicon and interfaces of amorphous and crystalline silicon, with the ultimate aim of understanding the photovoltaic properties of core-shell crystalline amorphous Si nanowire structures. We have shown, unexpectedly, from the simulations, that our understanding of hydrogenated bulk a-Si needs to be revisited, with our robust finding that when fully saturated with hydrogen, bulk a-Si exhibits a constant optical energy gap, irrespective of the hydrogen concentration in the sample. Unsaturated a-Si:H, with a lower than optimum hydrogen content, shows a smaller optical gap, that increases with hydrogen content until saturation is reached. The mobility gaps obtained from an analysis of the electronic states show similar behavior. We also obtained that the optical and mobility gaps show a volcano curve as the H content is varied from 7% (undersaturation) to 18% (mild oversaturation). In the case of mild over saturation, the mid-gap states arise exclusively from an increase in the density of strained Si-Si bonds. Analysis of our structures shows the extra H atoms in this case form a bridge between neighboring silicon atoms which increases the corresponding Si-Si distance and promotes bond length disorder in the sample. That has the potential to enhance the Staebler-Wronski effect. Planar interface models of amorphous-crystalline silicon have been generated in Si (100), (110) and (111) surfaces. The interface models are characterized by structure, RDF, electronic density of states and optical absorption spectrum. We find that the least stable (100) surface will result in the formation of the thickest amorphous silicon layer, while the most stable (110) surface forms the smallest amorphous region. We calculated for the first time band offsets of a-Si:H/c-Si heterojunctions from first principles and examined the influence of different surface orientations and amorphous layer thickness on the offsets and implications for device performance. The band offsets depend on the amorphous layer thickness and increase with thickness. By controlling the amorphous layer thickness we can potentially optimise the solar cell parameters. Finally, we have successfully generated different amorphous layer thickness of the a-Si/c-Si and a-Si:H/c-Si 5 nm nanowires from heat and quench. We perform structural analysis of the a-Si-/c-Si nanowires. The RDF, Si-Si bond length distributions, and the coordination number distributions of amorphous regions of the nanowires reproduce similar behaviour compared to bulk amorphous silicon. In the final part of this thesis we examine different surface terminating chemical groups, -H, - OH and –NH2 in (001) GeNW. Our work shows that the diameter of Ge nanowires and the nature of surface terminating groups both play a significant role in both the magnitude and the nature of the nanowire band gaps, allowing tuning of the band gap by up to 1.1 eV. We also show for the first time how the nanowire diameter and surface termination shifts the absorption edge in the Ge nanowires to longer wavelengths. Thus, the combination of nanowire diameter and surface chemistry can be effectively utilised to tune the band gaps and thus light absorption properties of small diameter Ge nanowires.

The numerical simulation of enclosure fires using a Cfd fire field model coupled with a pyrolysis based solid fuel combustion submodel—a first approximation

In this paper, we present some early work concerned with the development of a simple solid fuel combustion model incorporated within a Computational Fluid Dynamics (CFD) framework. The model is intended for use in engineering applications of fire field modeling and represents an extension of this technique to situations involving the combustion of solid cellulosic fuels. A simple solid fuel combustion model consisting of a thermal pyrolysis model, a six flux radiation model and an eddy-dissipation model for gaseous combustion have been developed and implemented within the CFD code CFDS-FLOW3D. The model is briefly described and demonstrated through two applications involving fire spread in a compartment with a plywood lined ceiling. The two scenarios considered involve a fire in an open and closed compartment. The model is shown to be able to qualitatively predict behaviors similar to "flashover"—in the case of the open room—and "backdraft"— in the case of the initially closed room.