Classification of neocortical interneurons using affinity propagation

In spite of over a century of research on cortical circuits, it is still unknown how many classes of cortical neurons exist. Neuronal classification has been a difficult problem because it is unclear what a neuronal cell class actually is and what are the best characteristics are to define them. Recently, unsupervised classifications using cluster analysis based on morphological, physiological or molecular characteristics, when applied to selected datasets, have provided quantitative and unbiased identification of distinct neuronal subtypes. However, better and more robust classification methods are needed for increasingly complex and larger datasets. We explored the use of affinity propagation, a recently developed unsupervised classification algorithm imported from machine learning, which gives a representative example or exemplar for each cluster. As a case study, we applied affinity propagation to a test dataset of 337 interneurons belonging to four subtypes, previously identified based on morphological and physiological characteristics. We found that affinity propagation correctly classified most of the neurons in a blind, non-supervised manner. In fact, using a combined anatomical/physiological dataset, our algorithm differentiated parvalbumin from somatostatin interneurons in 49 out of 50 cases. Affinity propagation could therefore be used in future studies to validly classify neurons, as a first step to help reverse engineer neural circuits.

Manipulation of domain wall dynamics in amorphous microwires through the magnetoelastic anisotropy

We studied the effect of magnetoelastic anisotropy on domain wall (DW) dynamics and remagnetization process of magnetically bistable Fe-Co-rich microwires with metallic nucleus diameters (from 1.4 to 22 mu m). We manipulated the magnetoelastic anisotropy applying the tensile stresses and changing the magnetostriction constant and strength of the internal stresses. Microwires of the same composition of metallic nucleus but with different geometries exhibit different magnetic field dependence of DW velocity with different slopes. Application of stresses resulted in decrease of the DW velocity, v, and DW mobility, S. Quite fast DW propagation (v until 2,500 m/s at H about 30 A/m) has been observed in low magnetostrictive magnetically bistable Co56Fe8Ni10Si10B16 microwires. Consequently, we observed certain correlation between the magnetoelastic energy and DW dynamics in microwires: decreasing the magnetoelastic energy, K (me), DW velocity increases.

Domain wall dynamics of magnetically bistable microwires

EMM-FM2011 – First Euro Mediterranean Meeting on Functionalized Materials, edited by Cheikhrouhou, A. 1st Euro Mediterranean Meeting on Functionalized Materials (EMM-FM). Sousse, TUNISIA . Sep. 06-10, 2011

Scale-similarity model for Lagrangian velocity correlations in isotropic and stationary turbulence

A scale-similarity model for Lagrangian two-point, two-time velocity correlations LVCs in isotropic turbulence is developed from the Kolmogorov similarity hypothesis. It is a second approximation to the isocontours of LVCs, while the Smith-Hay model is only a first approximation. This model expresses the LVC by its space correlation and a dispersion velocity. We derive the analytical expression for the dispersion velocity from the Navier-Stokes equations using the quasinormality assumption. The dispersion velocity is dependent on enstrophy spectra and shown to be smaller than the sweeping velocity for the Eulerian velocity correlation. Therefore, the Lagrangian decorrelation process is slower than the Eulerian decorrelation process. The data from direct numerical simulation of isotropic turbulence support the scale-similarity model: the LVCs for different space separations collapse into a universal form when plotted against the separation axis defined by the model.

Measurement of velocity profiles in a rectangular microchannel with aspect ratio α = 0.35

In this work, we measured 14 horizontal velocity profiles along the vertical direction of a rectangular microchannel with aspect ratio alpha = h/w = 0.35 (h is the height of the channel and w is the width of the channel) using microPIV at Re = 1.8 and 3.6. The experimental velocity profiles are compared with the full 3D theoretical solution, and also with a Poiseuille parabolic profile. It is shown that the experimental velocity profiles in the horizontal and vertical planes are in agreement with the theoretical profiles, except for the planes close to the wall. The discrepancies between the experimental data and 3D theoretical results in the center vertical plane are less than 3.6%. But the deviations between experimental data and Poiseuille's results approaches 5%. It indicates that 2D Poiseuille profile is no longer a perfect theoretical approximation since a = 0.35. The experiments also reveal that, very near the hydrophilic wall (z = 0.5-1 mu m), the measured velocities are significantly larger than the theoretical velocity based on the no-slip assumption. A proper discussion on some physical effects influencing the near wall velocity measurement is given.

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

The near wall velocity measurements in microchannels with different diameter particles by microPIV/PTV. 8th International symposium on PIV

In near wall measurements with microPIV/PTV, whether seeding particles can be effectively used to detect local fluid velocity is a crucial problem. This talk presents our recent measurements in microchannels [1][2]. Based on measured velocity profiles with 200nm and 50nm in pure water, we found that the measured velocity profiles are agreed with the theoretical values in the middle of channel, but large deviations between measured data and theoretical prediction appear close to wall (0.25mm velocity deviation appeared in shear flows.

Experimentally Determined Flame Properties near Flammability Limits Under Gravity and Microgravity Conditions

Flammability limits for flames propagating in a rich propane/air mixture under gravity conditions appeared to be 6.3% C3H8 for downward propagation and 9.2% C3H8 for upward propagation. Different limits might be explained by the action of preferential diffusion of the deficient reactant (Le < 1) on the limit flames, which are in different states of instability. In one of the previous studies, the flammability limits under microgtravity conditions were found to be between the upward and downward limits obtained in a standard flammability tube under normal gravity conditions. It was found in those experiments that there are two limits under microgravity conditions: one indicated by visible flame propagation and another indicated by an increase of pressure without observed flame propagation. These limits were found to be far behind the limit for downward-propagating flame at 1 g (6.3% C3H8) and close to the limit for upward-propagating flame at 1 g (9.2% C3H8). It was decided in the present work to apply a special schlieren system and instant temperature measuring system for drop tower experiments to observe combustion development during propagation of the flame front. A small cubic closed vessel (inner side, 9 cm 9 cm 9 cm) with schlieren quality glass windows were used to study limit flames under gravity and microgravity conditions. Flame development in rich limit mixtures, not visible in previous experiments under microgravity conditions for strait photography, was identified with the use of the schlieren method and instant temperature measuring system. It was found in experiments in a small vessel that there is practically no difference in flammability limits under gravity and microgravity conditions. In this paper, the mechanism of flame propagation under these different conditions is systematically studied and compared and limit burning velocity is estimated.

MODELING EVOLUTION OF INTERNAL WAVES IN SOUTH CHINA SEA

Internal waves are an important factor in the design of drill operations and production in deep water, because the waves have very large amplitude and may induce large horizontal velocity. How the internal waves occur and propagate over benthal terrain is of great concern for ocean engineers. In the present paper, we have formulated a mathematical model of internal wave propagation in a two-layer deep water, which involves the effects of friction, dissipation and shoaling, and is capable of manifesting the variation of the amplitude and the velocity pattern. After calibration by field data measured at the Continental Slope in the Northern South China Sea, we have applied the model to the South China Sea, investigating the westward propagation of internal waves from the Luzon Strait, where internal waves originate due to the interaction of benthal ridge and tides. We find that the internal wave induced velocity profile is obviously characterized by the opposite flow below and above the pycnocline, which results in a strong shear, threatening safety of ocean structures, such as mooring system of oil platform, risers, etc. When internal waves propagate westwards, the amplitude attenuates due to the effects of friction and dissipation. The preliminary results show that the amplitude is likely to become half of its initial value at Luzon Strait when the internal waves propagate about 400 kilometers westwards.

Attenuation tomography. Modelling regional love waves: Imperial Valley to Pasadena

Abstract to Part I

The inverse problem of seismic wave attenuation is solved by an iterative back-projection method. The seismic wave quality factor, Q, can be estimated approximately by inverting the S-to-P amplitude ratios. Effects of various uncertain ties in the method are tested and the attenuation tomography is shown to be useful in solving for the spatial variations in attenuation structure and in estimating the effective seismic quality factor of attenuating anomalies.

Back-projection attenuation tomography is applied to two cases in southern California: Imperial Valley and the Coso-Indian Wells region. In the Coso-Indian Wells region, a highly attenuating body (S-wave quality factor (Q_β ≈ 30) coincides with a slow P-wave anomaly mapped by Walck and Clayton (1987). This coincidence suggests the presence of a magmatic or hydrothermal body 3 to 5 km deep in the Indian Wells region. In the Imperial Valley, slow P-wave travel-time anomalies and highly attenuating S-wave anomalies were found in the Brawley seismic zone at a depth of 8 to 12 km. The effective S-wave quality factor is very low (Q_β ≈ 20) and the P-wave velocity is 10% slower than the surrounding areas. These results suggest either magmatic or hydrothermal intrusions, or fractures at depth, possibly related to active shear in the Brawley seismic zone.

No-block inversion is a generalized tomographic method utilizing the continuous form of an inverse problem. The inverse problem of attenuation can be posed in a continuous form , and the no-block inversion technique is applied to the same data set used in the back-projection tomography. A relatively small data set with little redundancy enables us to apply both techniques to a similar degree of resolution. The results obtained by the two methods are very similar. By applying the two methods to the same data set, formal errors and resolution can be directly computed for the final model, and the objectivity of the final result can be enhanced.

Both methods of attenuation tomography are applied to a data set of local earthquakes in Kilauea, Hawaii, to solve for the attenuation structure under Kilauea and the East Rift Zone. The shallow Kilauea magma chamber, East Rift Zone and the Mauna Loa magma chamber are delineated as attenuating anomalies. Detailed inversion reveals shallow secondary magma reservoirs at Mauna Ulu and Puu Oo, the present sites of volcanic eruptions. The Hilina Fault zone is highly attenuating, dominating the attenuating anomalies at shallow depths. The magma conduit system along the summit and the East Rift Zone of Kilauea shows up as a continuous supply channel extending down to a depth of approximately 6 km. The Southwest Rift Zone, on the other hand, is not delineated by attenuating anomalies, except at a depth of 8-12 km, where an attenuating anomaly is imaged west of Puu Kou. The Ylauna Loa chamber is seated at a deeper level (about 6-10 km) than the Kilauea magma chamber. Resolution in the Mauna Loa area is not as good as in the Kilauea area, and there is a trade-off between the depth extent of the magma chamber imaged under Mauna Loa and the error that is due to poor ray coverage. Kilauea magma chamber, on the other hand, is well resolved, according to a resolution test done at the location of the magma chamber.

Abstract to Part II

Long period seismograms recorded at Pasadena of earthquakes occurring along a profile to Imperial Valley are studied in terms of source phenomena (e.g., source mechanisms and depths) versus path effects. Some of the events have known source parameters, determined by teleseismic or near-field studies, and are used as master events in a forward modeling exercise to derive the Green's functions (SH displacements at Pasadena that are due to a pure strike-slip or dip-slip mechanism) that describe the propagation effects along the profile. Both timing and waveforms of records are matched by synthetics calculated from 2-dimensional velocity models. The best 2-dimensional section begins at Imperial Valley with a thin crust containing the basin structure and thickens towards Pasadena. The detailed nature of the transition zone at the base of the crust controls the early arriving shorter periods (strong motions), while the edge of the basin controls the scattered longer period surface waves. From the waveform characteristics alone, shallow events in the basin are easily distinguished from deep events, and the amount of strike-slip versus dip-slip motion is also easily determined. Those events rupturing the sediments, such as the 1979 Imperial Valley earthquake, can be recognized easily by a late-arriving scattered Love wave that has been delayed by the very slow path across the shallow valley structure.