Nonlinear dynamic transfer functions for certain retinal neuronal systems

The applicability of the white-noise method to the identification of a nonlinear system is investigated. Subsequently, the method is applied to certain vertebrate retinal neuronal systems and nonlinear, dynamic transfer functions are derived which describe quantitatively the information transformations starting with the light-pattern stimulus and culminating in the ganglion response which constitutes the visually-derived input to the brain. The retina of the catfish, Ictalurus punctatus, is used for the experiments.

The Wiener formulation of the white-noise theory is shown to be impractical and difficult to apply to a physical system. A different formulation based on crosscorrelation techniques is shown to be applicable to a wide range of physical systems provided certain considerations are taken into account. These considerations include the time-invariancy of the system, an optimum choice of the white-noise input bandwidth, nonlinearities that allow a representation in terms of a small number of characterizing kernels, the memory of the system and the temporal length of the characterizing experiment. Error analysis of the kernel estimates is made taking into account various sources of error such as noise at the input and output, bandwidth of white-noise input and the truncation of the gaussian by the apparatus.

Nonlinear transfer functions are obtained, as sets of kernels, for several neuronal systems: Light → Receptors, Light → Horizontal, Horizontal → Ganglion, Light → Ganglion and Light → ERG. The derived models can predict, with reasonable accuracy, the system response to any input. Comparison of model and physical system performance showed close agreement for a great number of tests, the most stringent of which is comparison of their responses to a white-noise input. Other tests include step and sine responses and power spectra.

Many functional traits are revealed by these models. Some are: (a) the receptor and horizontal cell systems are nearly linear (small signal) with certain "small" nonlinearities, and become faster (latency-wise and frequency-response-wise) at higher intensity levels, (b) all ganglion systems are nonlinear (half-wave rectification), (c) the receptive field center to ganglion system is slower (latency-wise and frequency-response-wise) than the periphery to ganglion system, (d) the lateral (eccentric) ganglion systems are just as fast (latency and frequency response) as the concentric ones, (e) (bipolar response) = (input from receptors) - (input from horizontal cell), (f) receptive field center and periphery exert an antagonistic influence on the ganglion response, (g) implications about the origin of ERG, and many others.

An analytical solution is obtained for the spatial distribution of potential in the S-space, which fits very well experimental data. Different synaptic mechanisms of excitation for the external and internal horizontal cells are implied.

Three-dimensional nonlinear inelastic analysis of steel moment-frame buildings damaged by earthquake excitations

The Northridge earthquake of January 17, 1994, highlighted the two previously known problems of premature fracturing of connections and the damaging capabilities of near-source ground motion pulses. Large ground motions had not been experienced in a city with tall steel moment-frame buildings before. Some steel buildings exhibited fracture of welded connections or other types of structural degradation.

A sophisticated three-dimensional nonlinear inelastic program is developed that can accurately model many nonlinear properties commonly ignored or approximated in other programs. The program can assess and predict severely inelastic response of steel buildings due to strong ground motions, including collapse.

Three-dimensional fiber and segment discretization of elements is presented in this work. This element and its two-dimensional counterpart are capable of modeling various geometric and material nonlinearities such as moment amplification, spread of plasticity and connection fracture. In addition to introducing a three-dimensional element discretization, this work presents three-dimensional constraints that limit the number of equations required to solve various three-dimensional problems consisting of intersecting planar frames.

Two buildings damaged in the Northridge earthquake are investigated to verify the ability of the program to match the level of response and the extent and location of damage measured. The program is used to predict response of larger near-source ground motions using the properties determined from the matched response.

A third building is studied to assess three-dimensional effects on a realistic irregular building in the inelastic range of response considering earthquake directivity. Damage levels are observed to be significantly affected by directivity and torsional response.

Several strong recorded ground motions clearly exceed code-based levels. Properly designed buildings can have drifts exceeding code specified levels due to these ground motions. The strongest ground motions caused collapse if fracture was included in the model. Near-source ground displacement pulses can cause columns to yield prior to weaker-designed beams. Damage in tall buildings correlates better with peak-to-peak displacements than with peak-to-peak accelerations.

Dynamic response of tall buildings shows that higher mode response can cause more damage than first mode response. Leaking of energy between modes in conjunction with damage can cause torsional behavior that is not anticipated.

Various response parameters are used for all three buildings to determine what correlations can be made for inelastic building response. Damage levels can be dramatically different based on the inelastic model used. Damage does not correlate well with several common response parameters.

Realistic modeling of material properties and structural behavior is of great value for understanding the performance of tall buildings due to earthquake excitations.

Radiation damping effects on the ultrashort x-ray pulses generated by nonlinear Thomson backscattering

Nonlinear Thomson backscattering of an intense Gaussian laser pulse by a counterpropagating energetic electron is investigated by numerically solving the electron equation of motion taking into account the radiative damping force. The backscattered radiation characteristics are different for linearly and circularly polarized lasers because of a difference in their ponderomotive forces acting on the electron. The radiative electron energy loss weakens the backscattered power, breaks the symmetry of the backscattered-pulse profile, and prolongs the duration of the backscattered radiation. With the circularly polarized laser, an adjustable double-peaked backscattered pulse can be obtained. Such a profile has potential applications as a subfemtosecond x-ray pump and probe with adjustable time delay and power ratio. (c) 2006 American Institute of Physics.

Short-pulse laser absorption in very steep plasma density gradients

An analytical fluid model for resonance absorption during the oblique incidence by femtosecond laser pulses on a small-scale-length density plasma [k(0)L is an element of(0.1,10)] is proposed. The physics of resonance absorption is analyzed more clearly as we separate the electric field into an electromagnetic part and an electrostatic part. It is found that the characteristics of the physical quantities (fractional absorption, optimum angle, etc.) in a small-scale-length plasma are quite different from the predictions of classical theory. Absorption processes are generally dependent on the density scale length. For shorter scale length or higher laser intensity, vacuum heating tends to be dominant. It is shown that the electrons being pulled out and then returned to the plasma at the interface layer by the wave field can lead to a phenomenon like wave breaking. This can lead to heating of the plasma at the expanse of the wave energy. It is found that the optimum angle is independent of the laser intensity while the absorption rate increases with the laser intensity, and the absorption rate can reach as high as 25%. (c) 2006 American Institute of Physics.

Short-pulse laser absorption via J x B heating in ultrahigh intensity laser plasma interaction

An analytical fluid model for JxB heating during the normal incidence by a short ultraintense linearly polarized laser on a solid-density plasma is proposed. The steepening of an originally smooth electron density profile as the electrons are pushed inward by the laser is included self-consistently. It is shown that the JxB heating includes two distinct coupling processes depending on the initial laser and plasma conditions: for a moderate intensity (a <= 1), the ponderomotive force of the laser light can drive a large plasma wave at the point n(e)=4 gamma(0)n(c) resonantly. When this plasma wave is damped, the energy is transferred to the plasma. At higher intensity, the electron density is steepened to a high level by the time-independent ponderomotive force, n(e)> 4 gamma(0)n(c), so that no 2 omega resonance will occur, but the longitudinal component of the oscillating ponderomotive field can lead to an absorption mechanism similar to "vacuum heating." (c) 2006 American Institute of Physics.

A probabilistic treatment of uncertainty in nonlinear dynamical systems

In this work, computationally efficient approximate methods are developed for analyzing uncertain dynamical systems. Uncertainties in both the excitation and the modeling are considered and examples are presented illustrating the accuracy of the proposed approximations.

For nonlinear systems under uncertain excitation, methods are developed to approximate the stationary probability density function and statistical quantities of interest. The methods are based on approximating solutions to the Fokker-Planck equation for the system and differ from traditional methods in which approximate solutions to stochastic differential equations are found. The new methods require little computational effort and examples are presented for which the accuracy of the proposed approximations compare favorably to results obtained by existing methods. The most significant improvements are made in approximating quantities related to the extreme values of the response, such as expected outcrossing rates, which are crucial for evaluating the reliability of the system.

Laplace's method of asymptotic approximation is applied to approximate the probability integrals which arise when analyzing systems with modeling uncertainty. The asymptotic approximation reduces the problem of evaluating a multidimensional integral to solving a minimization problem and the results become asymptotically exact as the uncertainty in the modeling goes to zero. The method is found to provide good approximations for the moments and outcrossing rates for systems with uncertain parameters under stochastic excitation, even when there is a large amount of uncertainty in the parameters. The method is also applied to classical reliability integrals, providing approximations in both the transformed (independently, normally distributed) variables and the original variables. In the transformed variables, the asymptotic approximation yields a very simple formula for approximating the value of SORM integrals. In many cases, it may be computationally expensive to transform the variables, and an approximation is also developed in the original variables. Examples are presented illustrating the accuracy of the approximations and results are compared with existing approximations.

Enhancement of Kerr nonlinearity with vanishing absorption in a tripod scheme

The giant enhancement of Kerr nonlinearity in a four-level tripod type system is investigated theoretically. By tuning the value of the Rabi frequency of the coherent control field, owing to the double dark resonances, the giant-enhanced Kerr nonlinearity can be achieved within the right transparency window. The in fluence of Doppler broadening is also discussed.

Upconversion luminescence and optical power limiting effect based on two- and three-photon absorption processes of ZnO crystal

Near-infrared to UV and visible upconversion luminescence was observed in single-crystalline ZnO under an 800 nm infrared femtosecond laser irradiation. The optical properties of the crystal reveal that the UV and VIS emission band are due to the exciton transition (D0X) bound to neutral donors and the deep luminescent centers in ZnO, respectively. The relationship between the upconversion luminescence intensity and the pump power of the femtosecond laser reveals that the UV emission belongs to three-photon sequential band-to-band excitation and the VIS emission belongs to two-photon simultaneous defect-absorption induced luminescence. A saturation phenomenon and polarization-dependent effect are also observed in the upconversion process of ZnO. A very good optical power limiting performance at 800 nm has been demonstrated. The two- and three-photon absorption coefficients of ZnO crystal were measured to be 0.2018 cm GW(-1) and 7.102 x 10(-3) cm(3) GW(-2), respectively. The two- and three-photon cross sections were calculated to be 1.189 x 10(-51) cm(4) s and 1.040 x 10(-80) cm(6) s(2), respectively. (c) 2007 Elsevier B.V. All rights reserved.

Extracting material response from simple mechanical tests on hardening-softening-hardening viscoplastic solids

Compliant foams are usually characterized by a wide range of desirable mechanical properties. These properties include viscoelasticity at different temperatures, energy absorption, recoverability under cyclic loading, impact resistance, and thermal, electrical, acoustic and radiation-resistance. Some foams contain nano-sized features and are used in small-scale devices. This implies that the characteristic dimensions of foams span multiple length scales, rendering modeling their mechanical properties difficult. Continuum mechanics-based models capture some salient experimental features like the linear elastic regime, followed by non-linear plateau stress regime. However, they lack mesostructural physical details. This makes them incapable of accurately predicting local peaks in stress and strain distributions, which significantly affect the deformation paths. Atomistic methods are capable of capturing the physical origins of deformation at smaller scales, but suffer from impractical computational intensity. Capturing deformation at the so-called meso-scale, which is capable of describing the phenomenon at a continuum level, but with some physical insights, requires developing new theoretical approaches.

A fundamental question that motivates the modeling of foams is ‘how to extract the intrinsic material response from simple mechanical test data, such as stress vs. strain response?’ A 3D model was developed to simulate the mechanical response of foam-type materials. The novelty of this model includes unique features such as the hardening-softening-hardening material response, strain rate-dependence, and plastically compressible solids with plastic non-normality. Suggestive links from atomistic simulations of foams were borrowed to formulate a physically informed hardening material input function. Motivated by a model that qualitatively captured the response of foam-type vertically aligned carbon nanotube (VACNT) pillars under uniaxial compression [2011,“Analysis of Uniaxial Compression of Vertically Aligned Carbon Nanotubes,” J. Mech.Phys. Solids, 59, pp. 2227–2237, Erratum 60, 1753–1756 (2012)], the property space exploration was advanced to three types of simple mechanical tests: 1) uniaxial compression, 2) uniaxial tension, and 3) nanoindentation with a conical and a flat-punch tip. The simulations attempt to explain some of the salient features in experimental data, like
1) The initial linear elastic response.
2) One or more nonlinear instabilities, yielding, and hardening.

The model-inherent relationships between the material properties and the overall stress-strain behavior were validated against the available experimental data. The material properties include the gradient in stiffness along the height, plastic and elastic compressibility, and hardening. Each of these tests was evaluated in terms of their efficiency in extracting material properties. The uniaxial simulation results proved to be a combination of structural and material influences. Out of all deformation paths, flat-punch indentation proved to be superior since it is the most sensitive in capturing the material properties.

Random propagation in complex systems : nonlinear matrix recursions and epidemic spread

This dissertation studies long-term behavior of random Riccati recursions and mathematical epidemic model. Riccati recursions are derived from Kalman filtering. The error covariance matrix of Kalman filtering satisfies Riccati recursions. Convergence condition of time-invariant Riccati recursions are well-studied by researchers. We focus on time-varying case, and assume that regressor matrix is random and identical and independently distributed according to given distribution whose probability distribution function is continuous, supported on whole space, and decaying faster than any polynomial. We study the geometric convergence of the probability distribution. We also study the global dynamics of the epidemic spread over complex networks for various models. For instance, in the discrete-time Markov chain model, each node is either healthy or infected at any given time. In this setting, the number of the state increases exponentially as the size of the network increases. The Markov chain has a unique stationary distribution where all the nodes are healthy with probability 1. Since the probability distribution of Markov chain defined on finite state converges to the stationary distribution, this Markov chain model concludes that epidemic disease dies out after long enough time. To analyze the Markov chain model, we study nonlinear epidemic model whose state at any given time is the vector obtained from the marginal probability of infection of each node in the network at that time. Convergence to the origin in the epidemic map implies the extinction of epidemics. The nonlinear model is upper-bounded by linearizing the model at the origin. As a result, the origin is the globally stable unique fixed point of the nonlinear model if the linear upper bound is stable. The nonlinear model has a second fixed point when the linear upper bound is unstable. We work on stability analysis of the second fixed point for both discrete-time and continuous-time models. Returning back to the Markov chain model, we claim that the stability of linear upper bound for nonlinear model is strongly related with the extinction time of the Markov chain. We show that stable linear upper bound is sufficient condition of fast extinction and the probability of survival is bounded by nonlinear epidemic map.

Polaron effects on the optical absorptions in asymmetrical semi-parabolic quantum wells

The linear and nonlinear optical absorptions considering the weak-coupling electron-LO-phonon interaction in asymmetrical semiparabolic quantum wells are theoretically investigated. The numerical results for the typical GaAs/AlxGa1-xAs material show that the factors of Al content x, the relaxation time and the photon energy have great influence on the optical absorption coefficients. Moreover, the theoretical values of the optical absorptions are more than a factor of 2-3 higher than the one in the structure without considering the electron-LO-phonon interaction by calculating. (C) 2007 Elsevier B.V. All rights reserved.