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.

Advancing EDL technologies for future space missions : from ground testing facilities to ablative heatshields

Motivated by recent MSL results where the ablation rate of the PICA heatshield was over-predicted, and staying true to the objectives outlined in the NASA Space Technology Roadmaps and Priorities report, this work focuses on advancing EDL technologies for future space missions.

Due to the difficulties in performing flight tests in the hypervelocity regime, a new ground testing facility called the vertical expansion tunnel is proposed. The adverse effects from secondary diaphragm rupture in an expansion tunnel may be reduced or eliminated by orienting the tunnel vertically, matching the test gas pressure and the accelerator gas pressure, and initially separating the test gas from the accelerator gas by density stratification. If some sacrifice of the reservoir conditions can be made, the VET can be utilized in hypervelocity ground testing, without the problems associated with secondary diaphragm rupture.

The performance of different constraints for the Rate-Controlled Constrained-Equilibrium (RCCE) method is investigated in the context of modeling reacting flows characteristic to ground testing facilities, and re-entry conditions. The effectiveness of different constraints are isolated, and new constraints previously unmentioned in the literature are introduced. Three main benefits from the RCCE method were determined: 1) the reduction in number of equations that need to be solved to model a reacting flow; 2) the reduction in stiffness of the system of equations needed to be solved; and 3) the ability to tabulate chemical properties as a function of a constraint once, prior to running a simulation, along with the ability to use the same table for multiple simulations.

Finally, published physical properties of PICA are compiled, and the composition of the pyrolysis gases that form at high temperatures internal to a heatshield is investigated. A necessary link between the composition of the solid resin, and the composition of the pyrolysis gases created is provided. This link, combined with a detailed investigation into a reacting pyrolysis gas mixture, allows a much needed consistent, and thorough description of many of the physical phenomena occurring in a PICA heatshield, and their implications, to be presented.

Through the use of computational fluid mechanics and computational chemistry methods, significant contributions have been made to advancing ground testing facilities, computational methods for reacting flows, and ablation modeling.

Toward Understanding Astrophysical Phenomena

Fast radio bursts (FRBs), a novel type of radio pulse, whose physics is not yet understood at all. Only a handful of FRBs had been detected when we started this project. Taking account of the scant observations, we put physical constraints on FRBs. We excluded proposals of a galactic origin for their extraordinarily high dispersion measures (DM), in particular stellar coronas and HII regions. Therefore our work supports an extragalactic origin for FRBs. We show that the resolved scattering tail of FRB 110220 is unlikely to be due to propagation through the intergalactic plasma. Instead the scattering is probably caused by the interstellar medium in the FRB's host galaxy, and indicates that this burst sits in the central region of that galaxy. Pulse durations of order $\ms$ constrain source sizes of FRBs implying enormous brightness temperatures and thus coherent emission. Electric fields near FRBs at cosmological distances would be so strong that they could accelerate free electrons from rest to relativistic energies in a single wave period. When we worked on FRBs, it was unclear whether they were genuine astronomical signals as distinct from `perytons', clearly terrestrial radio bursts, sharing some common properties with FRBs. Recently, in April 2015, astronomers discovered that perytons were emitted by microwave ovens. Radio chirps similar to FRBs were emitted when their doors opened while they were still heating. Evidence for the astronomical nature of FRBs has strengthened since our paper was published. Some bursts have been found to show linear and circular polarizations and Faraday rotation of the linear polarization has also been detected. I hope to resume working on FRBs in the near future. But after we completed our FRB paper, I decided to pause this project because of the lack of observational constraints.

The pulsar triple system, J0733+1715, has its orbital parameters fitted to high accuracy owing to the precise timing of the central $\ms$ pulsar. The two orbits are highly hierarchical, namely $P_{\mathrm{orb,1}}\ll P_{\mathrm{orb,2}}$, where 1 and 2 label the inner and outer white dwarf (WD) companions respectively. Moreover, their orbital planes almost coincide, providing a unique opportunity to study secular interaction associated purely with eccentricity beyond the solar system. Secular interaction only involves effect averaged over many orbits. Thus each companion can be represented by an elliptical wire with its mass distributed inversely proportional to its local orbital speed. Generally there exists a mutual torque, which vanishes only when their apsidal lines are parallel or anti-parallel. To maintain either mode, the eccentricity ratio, $e_1/e_2$, must be of the proper value, so that both apsidal lines precess together. For J0733+1715, $e_1\ll e_2$ for the parallel mode, while $e_1\gg e_2$ for the anti-parallel one. We show that the former precesses $\sim 10$ times slower than the latter. Currently the system is dominated by the parallel mode. Although only a little anti-parallel mode survives, both eccentricities especially $e_1$ oscillate on $\sim 10^3\yr$ timescale. Detectable changes would occur within $\sim 1\yr$. We demonstrate that the anti-parallel mode gets damped $\sim 10^4$ times faster than its parallel brother by any dissipative process diminishing $e_1$. If it is the tidal damping in the inner WD, we proceed to estimate its tidal quantity parameter ($Q$) to be $\sim 10^6$, which was poorly constrained by observations. However, tidal damping may also happen during the preceding low-mass X-ray binary (LMXB) phase or hydrogen thermal nuclear flashes. But, in both cases, the inner companion fills its Roche lobe and probably suffers mass/angular momentum loss, which might cause $e_1$ to grow rather than decay.

Several pairs of solar system satellites occupy mean motion resonances (MMRs). We divide these into two groups according to their proximity to exact resonance. Proximity is measured by the existence of a separatrix in phase space. MMRs between Io-Europa, Europa-Ganymede and Enceladus-Dione are too distant from exact resonance for a separatrix to appear. A separatrix is present only in the phase spaces of the Mimas-Tethys and Titan-Hyperion MMRs and their resonant arguments are the only ones to exhibit substantial librations. When a separatrix is present, tidal damping of eccentricity or inclination excites overstable librations that can lead to passage through resonance on the damping timescale. However, after investigation, we conclude that the librations in the Mimas-Tethys and Titan-Hyperion MMRs are fossils and do not result from overstability.

Rubble piles are common in the solar system. Monolithic elements touch their neighbors in small localized areas. Voids occupy a significant fraction of the volume. In a fluid-free environment, heat cannot conduct through voids; only radiation can transfer energy across them. We model the effective thermal conductivity of a rubble pile and show that it is proportional the square root of the pressure, $P$, for $P\leq \epsy^3\mu$ where $\epsy$ is the material's yield strain and $\mu$ its shear modulus. Our model provides an excellent fit to the depth dependence of the thermal conductivity in the top $140\,\mathrm{cm}$ of the lunar regolith. It also offers an explanation for the low thermal inertias of rocky asteroids and icy satellites. Lastly, we discuss how rubble piles slow down the cooling of small bodies such as asteroids.

Electromagnetic (EM) follow-up observations of gravitational wave (GW) events will help shed light on the nature of the sources, and more can be learned if the EM follow-ups can start as soon as the GW event becomes observable. In this paper, we propose a computationally efficient time-domain algorithm capable of detecting gravitational waves (GWs) from coalescing binaries of compact objects with nearly zero time delay. In case when the signal is strong enough, our algorithm also has the flexibility to trigger EM observation {\it before} the merger. The key to the efficiency of our algorithm arises from the use of chains of so-called Infinite Impulse Response (IIR) filters, which filter time-series data recursively. Computational cost is further reduced by a template interpolation technique that requires filtering to be done only for a much coarser template bank than otherwise required to sufficiently recover optimal signal-to-noise ratio. Towards future detectors with sensitivity extending to lower frequencies, our algorithm's computational cost is shown to increase rather insignificantly compared to the conventional time-domain correlation method. Moreover, at latencies of less than hundreds to thousands of seconds, this method is expected to be computationally more efficient than the straightforward frequency-domain method.

Phase-space analysis of wavefront coding imaging systems

We explore the use of the Radon-Wigner transform, which is associated with the fractional Fourier transform of the pupil function, for determining the point spread function (PSF) of an incoherant defocused optical system. Then we introduce these phase-space tools to analyse the wavefront coding imaging system. It is shown that the shape of the PSF for such a system is highly invarient to the defocous-related aberrations except for a lateral shift. The optical transfer function of this system is also investigated briefly from a new understanding of ambiguity function.

Free space adaptive optical interconnect, using a ferroelectric liquid crystal SLM for beam steering

A free-space, board-to-board, adaptive optical interconnect demonstrator has been developed. Binary phase gratings displayed on a Ferroelectric Liquid Crystal Spatial Light Modulator are used to maintain data transfer at 1.25Gbps, given varying optical misalignment © 2005 Optical Society of America.

Multi-objective optimization for multiphase orbital rendezvous missions

A multi-objective optimization approach was proposed for multiphase orbital rendezvous missions and validated by application to a representative numerical problem. By comparing the Pareto fronts obtained using the proposed method, the relationships between the three objectives considered are revealed, and the influences of other mission parameters, such as the sensors' field of view, can also be analyzed effectively. For multiphase orbital rendezvous missions, the tradeoff relationships between the total velocity increment and the trajectory robustness index as well as between the total velocity increment and the time of flight are obvious and clear. However, the tradeoff relationship between the time of flight and the trajectory robustness index is weak, especially for the four- and five-phase missions examined. The proposed approach could be used to reorganize a stable rendezvous profile for an engineering rendezvous mission, when there is a failure that prevents the completion of the nominal mission.

Development of a time and space resolved sampling probe diagnostic for engine exhaust hydrocarbons

In order to understand how unburned hydrocarbons emerge from SI engines and, in particular, how non-fuel hydrocarbons are formed and oxidized, a new gas sampling technique has been developed. A sampling unit, based on a combination of techniques used in the Fast Flame Ionization Detector (FFID) and wall-mounted sampling valves, was designed and built to capture a sample of exhaust gas during a specific period of the exhaust process and from a specific location within the exhaust port. The sampling unit consists of a transfer tube with one end in the exhaust port and the other connected to a three-way valve that leads, on one side, to a FFID and, on the other, to a vacuum chamber with a high-speed solenoid valve. Exhaust gas, drawn by the pressure drop into the vacuum chamber, impinges on the face of the solenoid valve and flows radially outward. Once per cycle during a specified crank angle interval, the solenoid valve opens and traps exhaust gas in a storage unit, from which gas chromatography (GC) measurements are made. The port end of the transfer tube can be moved to different locations longitudinally or radially, thus allowing spatial resolution and capturing any concentration differences between port walls and the center of the flow stream. Further, the solenoid valve's opening and closing times can be adjusted to allow sampling over a window as small as 0.6 ms during any portion of the cycle, allowing resolution of a crank angle interval as small as 15°CA. Cycle averaged total HC concentration measured by the FFID and that measured by the sampling unit are in good agreement, while the sampling unit goes one step further than the FFID by providing species concentrations. Comparison with previous measurements using wall-mounted sampling valves suggests that this sampling unit is fully capable of providing species concentration information as a function of air/fuel ratio, load, and engine speed at specific crank angles. © Copyright 1996 Society of Automotive Engineers, Inc.

Free space adaptive optical interconnect, using a ferroelectric liquid crystal SLM for beam steering

A free-space, board-to-board, adaptive optical interconnect demonstrator has been developed. Binary phase gratings displayed on a Ferroelectric Liquid Crystal Spatial Light Modulator are used to maintain data transfer at 1.25Gbps, given varying optical misalignment.© 2005 Optical Society of America.

Effects of phase-breaking on long-range charge transfer in DNA: Partially-coherent-tunneling model study

The mechanism of hole charge transfer in DNA of various lengths and sequences is investigated based on a partially coherent tunneling theory (Zhang et al., J Chem Phys 117:4578, 2002), where the effects of phase-breaking in adenine-thymine and guanine-cytosine base pairs are treated on equal foot. This work aims at providing a self-consistent microscopic interpretation for rate experiments on various DNA systems. We will also clarify the condition under which the simple superexchange-mediated-hopping picture is valid, and make some comments on the further development of present theory.

Scattering matrix approach to electronic dephasing in long-range electron transfer

Based on the Buttiker dephasing model, we propose an analytical scattering matrix approach to the long-range electron transfer phenomena. The present efficient scheme smoothly interpolates between the superexchange and the sequential hopping mechanisms. Various properties such as the drastic dephasing-assisted enhancement and turnover behaviors are demonstrated in good agreement with those obtained via the dynamical reduced density-matrix methods. These properties are further elucidated as results of the interplay among the dephasing strength, the tunneling parameter, and the bridge length of the electron transfer system. (C) 2001 American Institute of Physics.

Optimized logarithmic phase masks used to generate defocus invariant modulation transfer function for wavefront coding system

In a previous Letter [Opt. Lett. 33, 1171 (2008)], we proposed an improved logarithmic phase mask by making modifications to the original one designed by Sherif. However, further studies in another paper [Appl. Opt. 49, 229 (2010)] show that even when the Sherif mask and the improved one are optimized, their corresponding defocused modulation transfer functions (MTFs) are still not stable with respect to focus errors. So, by further modifying their phase profiles, we design another two logarithmic phase masks that exhibit more stable defocused MTF. However, with the defocus-induced phase effect considered, we find that the performance of the two masks proposed in this Letter is better than the Sherif mask, but worse than our previously proposed phase mask, according to the Hilbert space angle. (C) 2010 Optical Society of America

Quantum-mechanical calculations of vibrationally resolved cross sections for non-dissociative charge transfer of O3+ with H-2

Charge transfer due to collisions of ground state O3+ (2s(2)2p P-2) ions with molecular hydrogen is investigated using the quantum-mechanical molecular-orbital close-coupling (MOCC) method, and electronic and vibrational state-selective cross sections along with the corresponding differential cross sections are calculated for projectile energies of 100, 500, 1000 and 5000 eV/u at the orientation angles of 25 degrees,45 degrees and 89 degrees. The adiabatic potentials and radial coupling matrix elements utilized in the QMOCC calculations were obtained with the spin-coupled valence-bond approach. The infinite order sudden approximation (IOSA) and the vibrational sudden approximation (VSA) are utilized to deal with the rotation of H-2 and the coupling between the electron and the vibration of H-2. It is found that the distribution of vibrationally resolved cross sections with the vibrational quantum number upsilon' of H-2(+) (upsilon') varies with the increment of the projectile energy; and the electronic and vibrational stateselective differential cross sections show similar behaviors: there is a highest platform within a very small scattering angle, beyond which the differential cross sections decrease as the scattering angle increases and lots of oscillating structures appear, where the scattering angle of the first structure decreases as E-P(-1/2) with the increment of the projectile energy E-P; and the structure and amplitude of the differential cross sections are sensitive to the orientation of molecule H-2, which provides a possibility to identify the orientations of molecule H-2 by the vibrational state-selective differential scattering processes.

Ab initio calculations of differential cross sections for single charge transfer in He-3(2+)+He-4 collisions

The single charge transfer process in He-3(2+)+He-4 collisions is investigated using the quantum-mechanical molecular-orbital close-coupling method, in which the adiabatic potentials and radial couplings are calculated by using the ab initio multireference single- and double-excitation configuration interaction methods. The differential cross sections for the single charge transfer are presented at the laboratorial energies E = 6 keV and 10 keV for the projectile He-3(2+). Comparison with the existing data shows that the present results are better in agreement with the experimental measurements than other calculations in the dominant small angle scattering, which is attributed to the accurate calculations of the adiabatic potentials and the radial couplings.

Medium-spin states in Ge-70 and the role of the g(9/2) orbital

Medium-spin states of Ge-70 have been studied via the Ni-60(C-12,2p gamma)Ge-70 reaction at 45 MeV. The ground-state band and the second 0(+) band have been extended to the 12(+) and 8(+) states, respectively. Two negative-parity bands, one of which has a coupled structure and the other has a decoupled structure, have been observed additionally. Although the latter decoupled structure was known up to the (21(-)) state from a previous experiment, the part of the level scheme up to the 15(-) state has been largely modified by the present experiment. Backbendings observed in the positive- and negative-parity yrast bands have been compared with those of the neighboring even Ge isotopes. The experimental level structure has been compared with the shell-model calculations in the model space (2p(3/2), 1f(5/2), 2(p1/2), 1g(9/2)) employing two kinds of effective interactions, one of which is an extended P + QQ interaction with monopole interactions and the other is developed from a renormalized G matrix. Microscopic structures of the observed bands have been discussed with the help of the shell-model calculations.

Proton transfer reactions within the pyridazine-methanol cluster ion

The multiphoton ionization of the hydrogen-bonding cluster pyridazine-methanol (C4H4N2-CH3OH) was studied using a time-of-flight mass spectrometer at the wavelengths of 355 and 532 nm. At both wavelengths, a series of protonated C4H4N2-(CH3OH)(n)-H+ cluster ions were obtained. Relevant ab initio calculations were performed with HF and B3LYP methods. Equilibrium geometries of both neutral and ionic C4H4N2-CH3OH clusters, and dissociation channels and dissociation energies of ionic clusters, are presented. The results show that when C4H4N2-CH3OH is vertically ionized, C4H4N2H+ and CH3O are the dominant products via proton transfer reaction. A high energy barrier makes another channel corresponding to the production of C4H4N2H+ and CH2OH disfavored. (C) 2002 Elsevier Science B.V. All rights reserved.