Orientational-induced strain hardening of axisymmetric grains.

The rheological response of oriented axisymmetric grains has additional degrees of complexity associated with their microstructure orientation. These additional kinematic degrees of freedom that give rise to complex transient macroscale rheological responses are not well understood. In this Letter, we study the rheology of axisymmetric grains subjected to transient flow. We identify strong coupling between the microstructure rearrangement and strain hardening which, under certain conditions, can yield jamming. We identify the critical conditions corresponding to jamming and the dependency on the shape of the grains. It is shown that this is a particular form of jamming that is directional in nature, since unjamming occurs if the shear direction is reversed.

Kinematic Model of Transient Shape-Induced Anisotropy in Dense Granular Flow.

Nonspherical particles are ubiquitous in nature and industry, yet previous theoretical models of granular media are mostly limited to systems of spherical particles. The problem is that in systems of nonspherical anisotropic particles, dynamic particle alignment critically affects their mechanical response. To study the tendency of such particles to align, we propose a simple kinematic model that relates the flow to the evolution of particle alignment with respect to each other. The validity of the proposed model is supported by comparison with particle-based simulations for various particle shapes ranging from elongated rice-like (prolate) to flattened lentil-like (oblate) particles. The model shows good agreement with the simulations for both steady-state and transient responses, and advances the development of comprehensive constitutive models for shape-anisotropic particles.

Roy's largest root test under rank-one alternatives.

Roy's largest root is a common test statistic in multivariate analysis, statistical signal processing and allied fields. Despite its ubiquity, provision of accurate and tractable approximations to its distribution under the alternative has been a longstanding open problem. Assuming Gaussian observations and a rank-one alternative, or concentrated noncentrality, we derive simple yet accurate approximations for the most common low-dimensional settings. These include signal detection in noise, multiple response regression, multivariate analysis of variance and canonical correlation analysis. A small-noise perturbation approach, perhaps underused in statistics, leads to simple combinations of standard univariate distributions, such as central and noncentral [Formula: see text] and [Formula: see text]. Our results allow approximate power and sample size calculations for Roy's test for rank-one effects, which is precisely where it is most powerful.

Vectorial phase retrieval for linear characterization of attosecond pulses.

The waveforms of attosecond pulses produced by high-harmonic generation carry information on the electronic structure and dynamics in atomic and molecular systems. Current methods for the temporal characterization of such pulses have limited sensitivity and impose significant experimental complexity. We propose a new linear and all-optical method inspired by widely used multidimensional phase retrieval algorithms. Our new scheme is based on the spectral measurement of two attosecond sources and their interference. As an example, we focus on the case of spectral polarization measurements of attosecond pulses, relying on their most fundamental property-being well confined in time. We demonstrate this method numerically by reconstructing the temporal profiles of attosecond pulses generated from aligned CO(2) molecules.

Required length of guardrails before hazards.

One way to protect against impacts during run-off-road accidents with infrastructure is the use of guardrails. However, real-world accidents indicate that vehicles can leave the road and end up behind the guardrail. These vehicles have no possibility of returning to the lane. Vehicles often end up behind the guardrail because the length of the guardrails installed before hazards is too short; this can lead to a collision with a shielded hazard. To identify the basic speed for determining the necessary length of guardrails, we analyzed the speed at which vehicles leave the roadway from the ZEDATU (Zentrale Datenbank Tödlicher Unfälle) real-world accidents database. The required length of guardrail was considered the length that reduces vehicle speed at a maximum theoretically possible deceleration of 0.3g behind the barrier based on real-world road departure speed. To determine the desired length of a guardrail ahead of a hazard, we developed a relationship between guardrail length and the speed at which vehicles depart the roadway. If the initial elements are flared away from the carriageway, the required length will be reduced by up to an additional 30% The ZEDATU database analysis showed that extending the current length of guardrails to the evaluated required length would reduce the number of fatalities among occupants of vehicles striking bridge abutments by approximately eight percent.

Duration and consistency of historical selection are correlated with adaptive trait evolution in the streamside salamander, Ambystoma barbouri.

A fundamental challenge in evolutionary biology concerns estimating the extent to which ecological trade-offs may impose constraints on adaptive evolution. Novel ecological stressors may limit adaptive evolution of naive lineages that have experienced historically different selective regimes. Regarded as recently derived from a pond-breeding ancestor, streamside salamanders face the novel and strong selection pressure of breeding in streams with fish predators. A statistical phylogenetic approach was used to test whether adaptive evolution of antipredator performance phenotypes in streamside salamanders was positively associated with: (1) estimated per-lineage duration of coexistence with predatory fish; and (2) consistency of this predator selective-regime within lineages. Average durations of fish contact were computed for each salamander lineage on a set of chronograms. Selection consistency was determined by estimating the number of ecological transitions between fish and fishless states using stochastic character mapping. Historical selection in streamside salamanders can be generally characterized as unstable, apparently punctuated by the stochastic loss and recolonization of predatory fish in most lineages. We found that the efficacy of antipredator phenotypes in salamanders is strongly related to historical duration, as well as consistency, of selection imposed by predatory fish.

Childhood absence epilepsy requiring more than one medication for seizure control.

INTRODUCTION: Absence epilepsy is the most common primary generalized epilepsy syndrome encountered in pediatric practice. Treatment is pharmacologically specific and usually successful with a single medication. The objective of this study was to identify any clinical or electroencephalographic features at initial presentation in a consecutive cohort of children with absence epilepsy that may be associated with the need for a second medication. METHODS: A computerized pediatric neurology database (1991-2007 inclusive) was retrospectively searched for all patients with typical absence seizures, 3 Hz spike and wave on EEG and no apparent symptomatic etiology who were over the age of two years at seizure onset with at least one year of follow-up. All such children were then divided into two groups; a) those requiring a single medication for seizure control (Group 1), and b) those requiring two medications for seizure control despite optimal management with the initial medication as determined by serum drug monitoring (Group 2). Clinical and electrographic features evident at diagnosis were then contrasted between Group 1 and 2. RESULTS: Seventy-five children with absence seizures were initially identified with 52 meeting the study's inclusion and none of the exclusion criteria. Of these 52 children, 43 required a single medication for seizure control (Group 1), while 9 required two or more medications for seizure control (Group 2). A significant difference (p < 0.05) was apparent between Group 1 and 2 with respect to gender (16/43 males vs 8/9 males) and mean age of diagnosis (8.19 years +/- 3.00 vs 6.06 years +/- 2.22). Age of onset of seizures, interval duration of seizures prior to treatment initiation, duration of seizures, presence of automatisms, family history, presence of co-morbid conditions and EEG findings were not found to be significantly different between the two Groups. CONCLUSIONS: Male gender and an earlier age of diagnosis is associated with the need for two medications for seizure control in children with absence epilepsy. This observation may suggest the need for more intensive early programmatic follow-up for young male children with newly diagnosed absence epilepsy to effect more rapid attainment of seizure control.

Langevin trajectories between fixed concentrations.

We consider the trajectories of particles diffusing between two infinite baths of fixed concentrations connected by a channel, e.g., a protein channel of a biological membrane. The steady state influx and efflux of Langevin trajectories at the boundaries of a finite volume containing the channel and parts of the two baths is replicated by termination of outgoing trajectories and injection according to a residual phase space density. We present a simulation scheme that maintains averaged fixed concentrations without creating spurious boundary layers, consistent with the assumed physics.

Geometric diffusions as a tool for harmonic analysis and structure definition of data: multiscale methods.

In the companion article, a framework for structural multiscale geometric organization of subsets of R(n) and of graphs was introduced. Here, diffusion semigroups are used to generate multiscale analyses in order to organize and represent complex structures. We emphasize the multiscale nature of these problems and build scaling functions of Markov matrices (describing local transitions) that lead to macroscopic descriptions at different scales. The process of iterating or diffusing the Markov matrix is seen as a generalization of some aspects of the Newtonian paradigm, in which local infinitesimal transitions of a system lead to global macroscopic descriptions by integration. This article deals with the construction of fast-order N algorithms for data representation and for homogenization of heterogeneous structures.

Geometric diffusions as a tool for harmonic analysis and structure definition of data: diffusion maps.

We provide a framework for structural multiscale geometric organization of graphs and subsets of R(n). We use diffusion semigroups to generate multiscale geometries in order to organize and represent complex structures. We show that appropriately selected eigenfunctions or scaling functions of Markov matrices, which describe local transitions, lead to macroscopic descriptions at different scales. The process of iterating or diffusing the Markov matrix is seen as a generalization of some aspects of the Newtonian paradigm, in which local infinitesimal transitions of a system lead to global macroscopic descriptions by integration. We provide a unified view of ideas from data analysis, machine learning, and numerical analysis.

Memoryless control of boundary concentrations of diffusing particles.

Flux between regions of different concentration occurs in nearly every device involving diffusion, whether an electrochemical cell, a bipolar transistor, or a protein channel in a biological membrane. Diffusion theory has calculated that flux since the time of Fick (1855), and the flux has been known to arise from the stochastic behavior of Brownian trajectories since the time of Einstein (1905), yet the mathematical description of the behavior of trajectories corresponding to different types of boundaries is not complete. We consider the trajectories of noninteracting particles diffusing in a finite region connecting two baths of fixed concentrations. Inside the region, the trajectories of diffusing particles are governed by the Langevin equation. To maintain average concentrations at the boundaries of the region at their values in the baths, a control mechanism is needed to set the boundary dynamics of the trajectories. Different control mechanisms are used in Langevin and Brownian simulations of such systems. We analyze models of controllers and derive equations for the time evolution and spatial distribution of particles inside the domain. Our analysis shows a distinct difference between the time evolution and the steady state concentrations. While the time evolution of the density is governed by an integral operator, the spatial distribution is governed by the familiar Fokker-Planck operator. The boundary conditions for the time dependent density depend on the model of the controller; however, this dependence disappears in the steady state, if the controller is of a renewal type. Renewal-type controllers, however, produce spurious boundary layers that can be catastrophic in simulations of charged particles, because even a tiny net charge can have global effects. The design of a nonrenewal controller that maintains concentrations of noninteracting particles without creating spurious boundary layers at the interface requires the solution of the time-dependent Fokker-Planck equation with absorption of outgoing trajectories and a source of ingoing trajectories on the boundary (the so called albedo problem).

Derivation of Poisson and Nernst-Planck equations in a bath and channel from a molecular model.

Permeation of ions from one electrolytic solution to another, through a protein channel, is a biological process of considerable importance. Permeation occurs on a time scale of micro- to milliseconds, far longer than the femtosecond time scales of atomic motion. Direct simulations of atomic dynamics are not yet possible for such long-time scales; thus, averaging is unavoidable. The question is what and how to average. In this paper, we average a Langevin model of ionic motion in a bulk solution and protein channel. The main result is a coupled system of averaged Poisson and Nernst-Planck equations (CPNP) involving conditional and unconditional charge densities and conditional potentials. The resulting NP equations contain the averaged force on a single ion, which is the sum of two components. The first component is the gradient of a conditional electric potential that is the solution of Poisson's equation with conditional and permanent charge densities and boundary conditions of the applied voltage. The second component is the self-induced force on an ion due to surface charges induced only by that ion at dielectric interfaces. The ion induces surface polarization charge that exerts a significant force on the ion itself, not present in earlier PNP equations. The proposed CPNP system is not complete, however, because the electric potential satisfies Poisson's equation with conditional charge densities, conditioned on the location of an ion, while the NP equations contain unconditional densities. The conditional densities are closely related to the well-studied pair-correlation functions of equilibrium statistical mechanics. We examine a specific closure relation, which on the one hand replaces the conditional charge densities by the unconditional ones in the Poisson equation, and on the other hand replaces the self-induced force in the NP equation by an effective self-induced force. This effective self-induced force is nearly zero in the baths but is approximately equal to the self-induced force in and near the channel. The charge densities in the NP equations are interpreted as time averages over long times of the motion of a quasiparticle that diffuses with the same diffusion coefficient as that of a real ion, but is driven by the averaged force. In this way, continuum equations with averaged charge densities and mean-fields can be used to describe permeation through a protein channel.

Oral and oropharyngeal carcinoma: surgical treatment in irradiated and nonirradiated patients.

Surgical treatment was used in 45 cases of oral cavity carcinoma and in 23 cases of oropharyngeal carcinoma. The three-year cure rate was 47% for oral cavity tumors and 9% for orophayngeal lesions. The cure rate was substantially higher in female patients than in males and in white patients than in nonwhites. Results of composite resections in 23 previously irradiated patients and 23 nonirradiated patients with T-2 and T-3 lesions are compared. The irradiated patients with oral cavity carcinoma had a lower cure rate and a much greater incidence of postoperative morbidity than the patients treated with operation alone. Surgical results for oropharyngeal carcinomas were poor in both radiated and nonirradiated patients, although the incidence of postoperative morbidity and mortality was higher in the irradiated group.