A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid-solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently bench-marked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

Neuroscience, quantum indeterminism and the Cartesian soul.

Quantum indeterminism is frequently invoked as a solution to the problem of how a disembodied soul might interact with the brain (as Descartes proposed), and is sometimes invoked in theories of libertarian free will even when they do not involve dualistic assumptions. Taking as example the Eccles-Beck model of interaction between self (or soul) and brain at the level of synaptic exocytosis, I here evaluate the plausibility of these approaches. I conclude that Heisenbergian uncertainty is too small to affect synaptic function, and that amplification by chaos or by other means does not provide a solution to this problem. Furthermore, even if Heisenbergian effects did modify brain functioning, the changes would be swamped by those due to thermal noise. Cells and neural circuits have powerful noise-resistance mechanisms, that are adequate protection against thermal noise and must therefore be more than sufficient to buffer against Heisenbergian effects. Other forms of quantum indeterminism must be considered, because these can be much greater than Heisenbergian uncertainty, but these have not so far been shown to play a role in the brain.

Energetics of RecA-mediated recombination reactions. Without ATP hydrolysis RecA can mediate polar strand exchange but is unable to recycle.

We demonstrate that the step of DNA strand exchange during RecA-mediated recombination reaction can occur equally efficiently in the presence or absence of ATP hydrolysis. The polarity of strand exchange is the same when instead of ATP its non-hydrolyzable analog adenosine-5'-O-(3-thiotriphosphate) is used. We show that the ATP dependence of recombination reaction is limited to the post-exchange stages of the reactions. The low DNA affinity state of RecA protomers, induced after ATP hydrolysis, is necessary for the dissociation of RecA-DNA complexes at the end of the reaction. This dissociation of RecA from DNA is necessary for the release of recombinant DNA molecules from the complexes formed with RecA and for the recycling of RecA protomers for another round of the recombination reaction.

Mendacastrum n. gen. and Domuzdagia n. gen., two Jurassic spherical Spumellaria (Radiolaria) with hagiastrid medullary shell

Two new spumellarian radiolarian genera, Mendacastrum and Domuzdagia, are described from the lower Tithonian and lower Pliensbachian respectively. Both have a spherical cortical shell of actinommid type and a spherical or subspherical double medullary shell with the inner medullary shell of hagiastrid s.l. type. The inner medullary shell of Mendacastrum is of dactyliosphaerid or higumastrid s. sit. type, whereas that of Domuzdagia is of angulobracchiid type. Since they cannot be assigned to any described Mesozoic pyloniacean families, they are considered as type genera of two new families: Mendacastridae and Domuzdagiidae respectively.

Polar gradients of the DYRK-family kinase Pom1 couple cell length with the cell cycle.

Cells normally grow to a certain size before they enter mitosis and divide. Entry into mitosis depends on the activity of Cdk1, which is inhibited by the Wee1 kinase and activated by the Cdc25 phosphatase. However, how cells sense their size for mitotic commitment remains unknown. Here we show that an intracellular gradient of the dual-specificity tyrosine-phosphorylation regulated kinase (DYRK) Pom1, which emanates from the ends of rod-shaped Schizosaccharomyces pombe cells, serves to measure cell length and control mitotic entry. Pom1 provides positional information both for polarized growth and to inhibit cell division at cell ends. We discovered that Pom1 is also a dose-dependent G2-M inhibitor. Genetic analyses indicate that Pom1 negatively regulates Cdr1 and Cdr2, two previously described Wee1 inhibitors of the SAD kinase family. This inhibition may be direct, because in vivo and in vitro evidence suggest that Pom1 phosphorylates Cdr2. Whereas Cdr1 and Cdr2 localize to a medial cortical region, Pom1 forms concentration gradients from cell tips that overlap with Cdr1 and Cdr2 in short cells, but not in long cells. Disturbing these Pom1 gradients leads to Cdr2 phosphorylation and imposes a G2 delay. In short cells, Pom1 prevents precocious M-phase entry, suggesting that the higher medial Pom1 levels inhibit Cdr2 and promote a G2 delay. Thus, gradients of Pom1 from cell ends provide a measure of cell length to regulate M-phase entry.

An algebraic operator approach to the analysis of Gerber-Shiu functions

We introduce an algebraic operator framework to study discounted penalty functions in renewal risk models. For inter-arrival and claim size distributions with rational Laplace transform, the usual integral equation is transformed into a boundary value problem, which is solved by symbolic techniques. The factorization of the differential operator can be lifted to the level of boundary value problems, amounting to iteratively solving first-order problems. This leads to an explicit expression for the Gerber-Shiu function in terms of the penalty function.

Transhepatic ipsilateral right portal vein embolization extended to segment IV: improving hypertrophy and resection outcomes with spherical particles and coils.

PURPOSE: To analyze outcomes after right portal vein embolization extended to segment IV (right PVE + IV) before extended right hepatectomy, including liver hypertrophy, resection rates, and complications after embolization and resection, and to assess differences in outcomes with two different particulate embolic agents. MATERIALS AND METHODS: Between 1998 and 2004, transhepatic ipsilateral right PVE + IV with particles and coils was performed in 44 patients with malignant hepatobiliary disease, including metastases (n = 24), biliary cancer (n = 14), and hepatocellular carcinoma (n = 6). Right PVE + IV was considered if the future liver remnant (FLR; segments II/III with or without I) was less than 25% of the total estimated liver volume (TELV). Tris-acryl microspheres (100-700 microm; n = 21) or polyvinyl alcohol (PVA) particles (355-1,000 microm; n = 23) were administered in a stepwise fashion. Smaller particles were used to occlude distal branches, followed by larger particles to occlude proximal branches until near-complete stasis. Coils were then placed in secondary portal branches. Computed tomographic volumetry was performed before and 3-4 weeks after right PVE + IV to assess FLR hypertrophy. Liver volumes and postembolization and postoperative outcomes were measured. RESULTS: After right PVE + IV with PVA particles, FLR volume increased 45.5% +/- 40.9% and FLR/TELV ratio increased 6.9% +/- 5.6%. After right PVE + IV with tris-acryl microspheres, FLR volume increased 69.0% +/- 30.7% and FLR/TELV ratio increased 9.7% +/- 3.3%. Differences in FLR volume (P = .0011), FLR/TELV ratio (P = .027), and resection rates (P = .02) were statistically significant. Seventy-one percent of patients underwent extended right hepatectomy (86% after receiving tris-acryl microspheres, 57% after receiving PVA). Thirteen patients (29%) did not undergo resection (extrahepatic spread [n = 9], inadequate hypertrophy [n = 3], other reasons [n = 1]). No patient developed postembolization syndrome or progressive liver insufficiency after embolization or resection. One death after resection occurred as a result of sepsis and hemorrhage. Median hospital stays were 1 day after right PVE + IV and 7 days after resection. CONCLUSION: Transhepatic ipsilateral right PVE + IV with use of particles and coils is a safe, effective method for inducing contralateral hypertrophy before extended right hepatectomy. Embolization with small spherical particles provides improved hypertrophy and resection rates compared with larger, nonspherical particles.

Constitutively active mutants of the alpha 1B-adrenergic receptor: role of highly conserved polar amino acids in receptor activation.

Site-directed mutagenesis and molecular dynamics simulations of the alpha 1B-adrenergic receptor (AR) were combined to explore the potential molecular changes correlated with the transition from R (inactive state) to R (active state). Using molecular dynamics analysis we compared the structural/dynamic features of constitutively active mutants with those of the wild type and of an inactive alpha 1B-AR to build a theoretical model which defines the essential features of R and R. The results of site-directed mutagenesis were in striking agreement with the predictions of the model supporting the following hypothesis. (i) The equilibrium between R and R depends on the equilibrium between the deprotonated and protonated forms, respectively, of D142 of the DRY motif. In fact, replacement of D142 with alanine confers high constitutive activity to the alpha 1B-AR. (ii) The shift of R143 of the DRY sequence out of a conserved 'polar pocket' formed by N63, D91, N344 and Y348 is a feature common to all the active structures, suggesting that the role of R143 is fundamental for mediating receptor activation. Disruption of these intramolecular interactions by replacing N63 with alanine constitutively activates the alpha 1B-AR. Our findings might provide interesting generalities about the activation process of G protein-coupled receptors.

Inherently self-calibrating non-Cartesian parallel imaging.

The use of self-calibrating techniques in parallel magnetic resonance imaging eliminates the need for coil sensitivity calibration scans and avoids potential mismatches between calibration scans and subsequent accelerated acquisitions (e.g., as a result of patient motion). Most examples of self-calibrating Cartesian parallel imaging techniques have required the use of modified k-space trajectories that are densely sampled at the center and more sparsely sampled in the periphery. However, spiral and radial trajectories offer inherent self-calibrating characteristics because of their densely sampled center. At no additional cost in acquisition time and with no modification in scanning protocols, in vivo coil sensitivity maps may be extracted from the densely sampled central region of k-space. This work demonstrates the feasibility of self-calibrated spiral and radial parallel imaging using a previously described iterative non-Cartesian sensitivity encoding algorithm.

Direct Fourier Tomographic Reconstruction Image-To-Image Filter

We present an open-source ITK implementation of a directFourier method for tomographic reconstruction, applicableto parallel-beam x-ray images. Direct Fourierreconstruction makes use of the central-slice theorem tobuild a polar 2D Fourier space from the 1D transformedprojections of the scanned object, that is resampled intoa Cartesian grid. Inverse 2D Fourier transform eventuallyyields the reconstructed image. Additionally, we providea complex wrapper to the BSplineInterpolateImageFunctionto overcome ITKâeuro?s current lack for image interpolatorsdealing with complex data types. A sample application ispresented and extensively illustrated on the Shepp-Loganhead phantom. We show that appropriate input zeropaddingand 2D-DFT oversampling rates together with radial cubicb-spline interpolation improve 2D-DFT interpolationquality and are efficient remedies to reducereconstruction artifacts.