Strong Dynamical Heterogeneity and Universal Scaling in Driven Granular Fluids

Large-scale simulations of two-dimensional bidisperse granular fluids allow us to determine spatial correlations of slow particles via the four-point structure factor S-4 (q, t). Both cases, elastic (epsilon = 1) and inelastic (epsilon < 1) collisions, are studied. As the fluid approaches structural arrest, i.e., for packing fractions in the range 0.6 <= phi <= 0.805, scaling is shown to hold: S-4 (q, t)/chi(4)(t) = s(q xi(t)). Both the dynamic susceptibility chi(4)(tau(alpha)) and the dynamic correlation length xi(tau(alpha)) evaluated at the alpha relaxation time tau(alpha) can be fitted to a power law divergence at a critical packing fraction. The measured xi(tau(alpha)) widely exceeds the largest one previously observed for three-dimensional (3d) hard sphere fluids. The number of particles in a slow cluster and the correlation length are related by a robust power law, chi(4)(tau(alpha)) approximate to xi(d-p) (tau(alpha)), with an exponent d - p approximate to 1.6. This scaling is remarkably independent of epsilon, even though the strength of the dynamical heterogeneity at constant volume fraction depends strongly on epsilon.

Sensitivity Gains, Linearity, and Spectral Reproducibility in Nonuniformly Sampled Multidimensional MAS NMR Spectra of High Dynamic Range

Recently, we have demonstrated that considerable inherent sensitivity gains are attained in MAS NMR spectra acquired by nonuniform sampling (NUS) and introduced maximum entropy interpolation (MINT) processing that assures the linearity of transformation between the time and frequency domains. In this report, we examine the utility of the NUS/MINT approach in multidimensional datasets possessing high dynamic range, such as homonuclear C-13-C-13 correlation spectra. We demonstrate on model compounds and on 1-73-(U-C-13,N-15)/74-108-(U-N-15) E. coli thioredoxin reassembly, that with appropriately constructed 50 % NUS schedules inherent sensitivity gains of 1.7-2.1-fold are readily reached in such datasets. We show that both linearity and line width are retained under these experimental conditions throughout the entire dynamic range of the signals. Furthermore, we demonstrate that the reproducibility of the peak intensities is excellent in the NUS/MINT approach when experiments are repeated multiple times and identical experimental and processing conditions are employed. Finally, we discuss the principles for design and implementation of random exponentially biased NUS sampling schedules for homonuclear C-13-C-13 MAS correlation experiments that yield high-quality artifact-free datasets.

Theoretical considerations to the verification of dynamic multileaf collimated fields with a SLIC-type portal image detector

The verification possibilities of dynamically collimated treatment beams with a scanning liquid ionization chamber electronic portal image device (SLIC-EPID) are investigated. The ion concentration in the liquid of a SLIC-EPID and therefore the read-out signal is determined by two parameters of a differential equation describing the creation and recombination of the ions. Due to the form of this equation, the portal image detector describes a nonlinear dynamic system with memory. In this work, the parameters of the differential equation were experimentally determined for the particular chamber in use and for an incident open 6 MV photon beam. The mathematical description of the ion concentration was then used to predict portal images of intensity-modulated photon beams produced by a dynamic delivery technique, the sliding window approach. Due to the nature of the differential equation, a mathematical condition for 'reliable leaf motion verification' in the sliding window technique can be formulated. It is shown that the time constants for both formation and decay of the equilibrium concentration in the chamber is in the order of seconds. In order to guarantee reliable leaf motion verification, these time constants impose a constraint on the rapidity of the image-read out for a given maximum leaf speed. For a leaf speed of 2 cm s(-1), a minimum image acquisition frequency of about 2 Hz is required. Current SLIC-EPID systems are usually too slow since they need about a second to acquire a portal image. However, if the condition is fulfilled, the memory property of the system can be used to reconstruct the leaf motion. It is shown that a simple edge detecting algorithm can be employed to determine the leaf positions. The method is also very robust against image noise.

Soft-decision equalization techniques for frequency selective MIMO channels

Multi-input multi-output (MIMO) technology is an emerging solution for high data rate wireless communications. We develop soft-decision based equalization techniques for frequency selective MIMO channels in the quest for low-complexity equalizers with BER performance competitive to that of ML sequence detection. We first propose soft decision equalization (SDE), and demonstrate that decision feedback equalization (DFE) based on soft-decisions, expressed via the posterior probabilities associated with feedback symbols, is able to outperform hard-decision DFE, with a low computational cost that is polynomial in the number of symbols to be recovered, and linear in the signal constellation size. Building upon the probabilistic data association (PDA) multiuser detector, we present two new MIMO equalization solutions to handle the distinctive channel memory. With their low complexity, simple implementations, and impressive near-optimum performance offered by iterative soft-decision processing, the proposed SDE methods are attractive candidates to deliver efficient reception solutions to practical high-capacity MIMO systems. Motivated by the need for low-complexity receiver processing, we further present an alternative low-complexity soft-decision equalization approach for frequency selective MIMO communication systems. With the help of iterative processing, two detection and estimation schemes based on second-order statistics are harmoniously put together to yield a two-part receiver structure: local multiuser detection (MUD) using soft-decision Probabilistic Data Association (PDA) detection, and dynamic noise-interference tracking using Kalman filtering. The proposed Kalman-PDA detector performs local MUD within a sub-block of the received data instead of over the entire data set, to reduce the computational load. At the same time, all the inter-ference affecting the local sub-block, including both multiple access and inter-symbol interference, is properly modeled as the state vector of a linear system, and dynamically tracked by Kalman filtering. Two types of Kalman filters are designed, both of which are able to track an finite impulse response (FIR) MIMO channel of any memory length. The overall algorithms enjoy low complexity that is only polynomial in the number of information-bearing bits to be detected, regardless of the data block size. Furthermore, we introduce two optional performance-enhancing techniques: cross- layer automatic repeat request (ARQ) for uncoded systems and code-aided method for coded systems. We take Kalman-PDA as an example, and show via simulations that both techniques can render error performance that is better than Kalman-PDA alone and competitive to sphere decoding. At last, we consider the case that channel state information (CSI) is not perfectly known to the receiver, and present an iterative channel estimation algorithm. Simulations show that the performance of SDE with channel estimation approaches that of SDE with perfect CSI.

Synthesis of computations and experiments for obtaining pulsatile gas flow rates from dynamic pressure difference measurements across an orifice-plate meter

The use of conventional orifice-plate meter is typically restricted to measurements of steady flows. This study proposes a new and effective computational-experimental approach for measuring the time-varying (but steady-in-the-mean) nature of turbulent pulsatile gas flows. Low Mach number (effectively constant density) steady-in-the-mean gas flows with large amplitude fluctuations (whose highest significant frequency is characterized by the value fF) are termed pulsatile if the fluctuations have a direct correlation with the time-varying signature of the imposed dynamic pressure difference and, furthermore, they have fluctuation amplitudes that are significantly larger than those associated with turbulence or random acoustic wave signatures. The experimental aspect of the proposed calibration approach is based on use of Coriolis-meters (whose oscillating arm frequency fcoriolis >> fF) which are capable of effectively measuring the mean flow rate of the pulsatile flows. Together with the experimental measurements of the mean mass flow rate of these pulsatile flows, the computational approach presented here is shown to be effective in converting the dynamic pressure difference signal into the desired dynamic flow rate signal. The proposed approach is reliable because the time-varying flow rate predictions obtained for two different orifice-plate meters exhibit the approximately same qualitative, dominant features of the pulsatile flow.

Method of Detecting System Function by Measuring Frequency Response

Real-time battery impedance spectrum is acquired using a one-time record. Fast Summation Transformation (FST) is a parallel method of acquiring a real-time battery impedance spectrum using a one-time record that enables battery diagnostics. An excitation current to a battery is a sum of equal amplitude sine waves of frequencies that are octave harmonics spread over a range of interest. A sample frequency is also octave and harmonically related to all frequencies in the sum. The time profile of this signal has a duration that is a few periods of the lowest frequency. The voltage response of the battery, average deleted, is the impedance of the battery in the time domain. Since the excitation frequencies are known and octave and harmonically related, a simple algorithm, FST, processes the time record by rectifying relative to the sine and cosine of each frequency. Another algorithm yields real and imaginary components for each frequency.

Method of Detecting System Function by Measuring Frequency Response

Real time battery impedance spectrum is acquired using one time record, Compensated Synchronous Detection (CSD). This parallel method enables battery diagnostics. The excitation current to a test battery is a sum of equal amplitude sin waves of a few frequencies spread over range of interest. The time profile of this signal has duration that is a few periods of the lowest frequency. The voltage response of the battery, average deleted, is the impedance of the battery in the time domain. Since the excitation frequencies are known, synchronous detection processes the time record and each component, both magnitude and phase, is obtained. For compensation, the components, except the one of interest, are reassembled in the time domain. The resulting signal is subtracted from the original signal and the component of interest is synchronously detected. This process is repeated for each component.

On the detection of dynamic responses in a drought-perturbed tropical rainforest in Borneo

The dynamics of aseasonal lowland dipterocarp forest in Borneo is influenced by perturbation from droughts. These events might be increasing in frequency and intensity in the future. This paper describes drought-affected dynamics between 1986 and 2001 in Sabah, Malaysia, and considers how it is possible, reliably and accurately, to measure both coarse- and fine-scale responses of the forest. Some fundamental concerns about methodology and data analysis emerge. In two plots forming 8 ha, mortality, recruitment, and stem growth rates of trees ≥10 cm gbh (girth at breast height) were measured in a ‘pre-drought’ period (1986–1996), and in a period (1996–2001) including the 1997–1998 ENSO-drought. For 2.56 ha of subplots, mortality and growth rates of small trees (10–<50 cm gbh) were found also for two sub-periods (1996–1999, 1999–2001). A total of c. 19 K trees were recorded. Mortality rate increased by 25% while both recruitment and relative growth rates increased by 12% for all trees at the coarse scale. For small trees, at the fine scale, mortality increased by 6% and 9% from pre-drought to drought and on to ‘post-drought’ sub-periods. Relative growth rates correspondingly decreased by 38% and increased by 98%. Tree size and topography interacted in a complex manner with between-plot differences. The forest appears to have been sustained by off-setting elevated tree mortality by highly resilient stem growth. This last is seen as the key integrating tree variable which links the external driver (drought causing water stress) and population dynamics recorded as mortality and recruitment. Suitably sound measurements of stem girth, leading to valid growth rates, are needed to understand and model tree dynamic responses to perturbations. The proportion of sound data, however, is in part determined by the drought itself.

Region Specific Response of Intervertebral Disc Cells to Complex Dynamic Loading: An Organ Culture Study Using a Dynamic Torsion-Compression Bioreactor

The spine is routinely subjected to repetitive complex loading consisting of axial compression, torsion, flexion and extension. Mechanical loading is one of the important causes of spinal diseases, including disc herniation and disc degeneration. It is known that static and dynamic compression can lead to progressive disc degeneration, but little is known about the mechanobiology of the disc subjected to combined dynamic compression and torsion. Therefore, the purpose of this study was to compare the mechanobiology of the intervertebral disc when subjected to combined dynamic compression and axial torsion or pure dynamic compression or axial torsion using organ culture. We applied four different loading modalities 1. control: no loading (NL), 2. cyclic compression (CC), 3. cyclic torsion (CT), and 4. combined cyclic compression and torsion (CCT) on bovine caudal disc explants using our custom made dynamic loading bioreactor for disc organ culture. Loads were applied for 8 h/day and continued for 14 days, all at a physiological magnitude and frequency. Our results provided strong evidence that complex loading induced a stronger degree of disc degeneration compared to one degree of freedom loading. In the CCT group, less than 10\% nucleus pulposus (NP) cells survived the 14 days of loading, while cell viabilities were maintained above 70\% in the NP of all the other three groups and in the annulus fibrosus (AF) of all the groups. Gene expression analysis revealed a strong up-regulation in matrix genes and matrix remodeling genes in the AF of the CCT group. Cell apoptotic activity and glycosaminoglycan content were also quantified but there were no statistically significant differences found. Cell morphology in the NP of the CCT was changed, as shown by histological evaluation. Our results stress the importance of complex loading on the initiation and progression of disc degeneration.

Dynamic Chemical Shift Imaging for Image-Guided Thermal Therapy

Magnetic resonance temperature imaging (MRTI) is recognized as a noninvasive means to provide temperature imaging for guidance in thermal therapies. The most common method of estimating temperature changes in the body using MR is by measuring the water proton resonant frequency (PRF) shift. Calculation of the complex phase difference (CPD) is the method of choice for measuring the PRF indirectly since it facilitates temperature mapping with high spatiotemporal resolution. Chemical shift imaging (CSI) techniques can provide the PRF directly with high sensitivity to temperature changes while minimizing artifacts commonly seen in CPD techniques. However, CSI techniques are currently limited by poor spatiotemporal resolution. This research intends to develop and validate a CSI-based MRTI technique with intentional spectral undersampling which allows relaxed parameters to improve spatiotemporal resolution. An algorithm based on autoregressive moving average (ARMA) modeling is developed and validated to help overcome limitations of Fourier-based analysis allowing highly accurate and precise PRF estimates. From the determined acquisition parameters and ARMA modeling, robust maps of temperature using the k-means algorithm are generated and validated in laser treatments in ex vivo tissue. The use of non-PRF based measurements provided by the technique is also investigated to aid in the validation of thermal damage predicted by an Arrhenius rate dose model.

Dynamic directed interictal connectivity in left and right temporal lobe epilepsy.

OBJECTIVE There is increasing evidence that epileptic activity involves widespread brain networks rather than single sources and that these networks contribute to interictal brain dysfunction. We investigated the fast-varying behavior of epileptic networks during interictal spikes in right and left temporal lobe epilepsy (RTLE and LTLE) at a whole-brain scale using directed connectivity. METHODS In 16 patients, 8 with LTLE and 8 with RTLE, we estimated the electrical source activity in 82 cortical regions of interest (ROIs) using high-density electroencephalography (EEG), individual head models, and a distributed linear inverse solution. A multivariate, time-varying, and frequency-resolved Granger-causal modeling (weighted Partial Directed Coherence) was applied to the source signal of all ROIs. A nonparametric statistical test assessed differences between spike and baseline epochs. Connectivity results between RTLE and LTLE were compared between RTLE and LTLE and with neuropsychological impairments. RESULTS Ipsilateral anterior temporal structures were identified as key drivers for both groups, concordant with the epileptogenic zone estimated invasively. We observed an increase in outflow from the key driver already before the spike. There were also important temporal and extratemporal ipsilateral drivers in both conditions, and contralateral only in RTLE. A different network pattern between LTLE and RTLE was found: in RTLE there was a much more prominent ipsilateral to contralateral pattern than in LTLE. Half of the RTLE patients but none of the LTLE patients had neuropsychological deficits consistent with contralateral temporal lobe dysfunction, suggesting a relationship between connectivity changes and cognitive deficits. SIGNIFICANCE The different patterns of time-varying connectivity in LTLE and RTLE suggest that they are not symmetrical entities, in line with our neuropsychological results. The highest outflow region was concordant with invasive validation of the epileptogenic zone. This enhanced characterization of dynamic connectivity patterns could better explain cognitive deficits and help the management of epilepsy surgery candidates.

The scaling of genetic diversity in a changing and fragmented world

Most species do not live in a constant environment over space or time. Their environment is often heterogeneous with a huge variability in resource availability and exposure to pathogens or predators, which may affect the local densities of the species. Moreover, the habitat might be fragmented, preventing free and isotropic migrations between local sub-populations (demes) of a species, making some demes more isolated than others. For example, during the last ice age populations of many species migrated towards refuge areas from which re-colonization originated when conditions improved. However, populations that could not move fast enough or could not adapt to the new environmental conditions faced extinctions. Populations living in these types of dynamic environments are often referred to as metapopulations and modeled as an array of subdivisions (or demes) that exchange migrants with their neighbors. Several studies have focused on the description of their demography, probability of extinction and expected patterns of diversity at different scales. Importantly, all these evolutionary processes may affect genetic diversity, which can affect the chance of populations to persist. In this chapter we provide an overview on the consequences of fragmentation, long-distance dispersal, range contractions and range shifts on genetic diversity. In addition, we describe new methods to detect and quantify underlying evolutionary processes from sampled genetic data.

Duration-dependent influence of dynamic torsion on the intervertebral disc: an intact disc organ culture study.

PURPOSE Mechanical loading is an important parameter that alters the homeostasis of the intervertebral disc (IVD). Studies have demonstrated the role of compression in altering the cellular metabolism, anabolic and catabolic events of the disc, but little is known how complex loading such as torsion-compression affects the IVD cell metabolism and matrix homeostasis. Studying how the duration of torsion affects disc matrix turnover could provide guidelines to prevent overuse injury to the disc and suggest possible beneficial effect of torsion. The aim of the study was to evaluate the biological response of the IVD to different durations of torsional loading. METHODS Intact bovine caudal IVD were isolated for organ culture in a bioreactor. Different daily durations of torsion were applied over 7 days at a physiological magnitude (±2°) in combination with 0.2 MPa compression, at a frequency of 1 Hz. RESULTS Nucleus pulpous (NP) cell viability and total disc volume decreased with 8 h of torsion-compression per day. Gene expression analysis suggested a down-regulated MMP13 with increased time of torsion. 1 and 4 h per day torsion-compression tended to increase the glycosaminoglycans/hydroxyproline ratio in the NP tissue group. CONCLUSIONS Our result suggests that load duration thresholds exist in both torsion and compression with an optimal load duration capable of promoting matrix synthesis and overloading can be harmful to disc cells. Future research is required to evaluate the specific mechanisms for these observed effects.

Higher order parametric excitation modes for spaceborne quadrupole mass spectrometers

This paper describes a technique to significantly improve upon the mass peak shape and mass resolution of spaceborne quadrupolemass spectrometers (QMSs) through higher order auxiliary excitation of the quadrupole field. Using a novel multiresonant tank circuit, additional frequency components can be used to drive modulating voltages on the quadrupole rods in a practical manner, suitable for both improved commercial applications and spaceflight instruments. Auxiliary excitation at frequencies near twice that of the fundamental quadrupole RF frequency provides the advantages of previously studied parametric excitation techniques, but with the added benefit of increased sensed excitation amplitude dynamic range and the ability to operate voltage scan lines through the center of upper stability islands. Using a field programmable gate array, the amplitudes and frequencies of all QMS signals are digitally generated and managed, providing a robust and stable voltage control system. These techniques are experimentally verified through an interface with a commercial Pfeiffer QMG422 quadrupole rod system. When operating through the center of a stability island formed from higher order auxiliary excitation, approximately 50% and 400% improvements in 1% mass resolution and peak stability were measured, respectively, when compared with traditional QMS operation. Although tested with a circular rod system, the presented techniques have the potential to improve the performance of both circular and hyperbolic rod geometry QMS sensors.

A compact lightweight multipurpose ground-penetrating radar for glaciological applications

We describe a compact lightweight impulse radar for radio-echo sounding of subsurface structures designed specifically for glaciological applications. The radar operates at frequencies between 10 and 75 MHz. Its main advantages are that it has a high signal-to-noise ratio and a corresponding wide dynamic range of 132 dB due mainly to its ability to perform real-time stacking (up to 4096 traces) as well as to the high transmitted power (peak voltage 2800 V). The maximum recording time window, 40 ?s at 100 MHz sampling frequency, results in possible radar returns from as deep as 3300 m. It is a versatile radar, suitable for different geophysical measurements (common-offset profiling, common midpoint, transillumination, etc.) and for different profiling set-ups, such as a snowmobile and sledge convoy or carried in a backpack and operated by a single person. Its low power consumption (6.6 W for the transmitter and 7.5 W for the receiver) allows the system to operate under battery power for mayor que7 hours with a total weight of menor que9 kg for all equipment, antennas and batteries.