Biomechanical wall properties of human intracranial aneurysms resected following surgical clipping (IRRAs Project)

Background and purpose: Individual rupture risk assessment of intracranial aneurysms is a major issue in the clinical management of asymptomatic aneurysms. Aneurysm rupture occurs when wall tension exceeds the strength limit of the wall tissue. At present, aneurysmal wall mechanics are poorly understood and thus, risk assessment involving mechanical properties is inexistent. Aneurysm computational hemodynamics studies make the assumption of rigid walls, an arguable simplification. We therefore aim to assess mechanical properties of ruptured and unruptured intracranial aneurysms in order to provide the foundation for future patient-specific aneurysmal risk assessment. This work also challenges some of the currently held hypotheses in computational flow hemodynamics research. Methods: A specific conservation protocol was applied to aneurysmal tissues following clipping and resection in order to preserve their mechanical properties. Sixteen intracranial aneurysms (11 female, 5 male) underwent mechanical uniaxial stress tests under physiological conditions, temperature, and saline isotonic solution. These represented 11 unruptured and 5 ruptured aneurysms. Stress/strain curves were then obtained for each sample, and a fitting algorithm was applied following a 3-parameter (C(10), C(01), C(11)) Mooney-Rivlin hyperelastic model. Each aneurysm was classified according to its biomechanical properties and (un)rupture status.Results: Tissue testing demonstrated three main tissue classes: Soft, Rigid, and Intermediate. All unruptured aneurysms presented a more Rigid tissue than ruptured or pre-ruptured aneurysms within each gender subgroup. Wall thickness was not correlated to aneurysmal status (ruptured/unruptured). An Intermediate subgroup of unruptured aneurysms with softer tissue characteristic was identified and correlated with multiple documented risk factors of rupture. Conclusion: There is a significant modification in biomechanical properties between ruptured aneurysm, presenting a soft tissue and unruptured aneurysms, presenting a rigid material. This finding strongly supports the idea that a biomechanical risk factor based assessment should be utilized in the to improve the therapeutic decision making.

Dynamic estimation of three-dimensional cerebrovascular deformation from rotational angiography

Purpose: The objective of this study is to investigate the feasibility of detecting and quantifying 3D cerebrovascular wall motion from a single 3D rotational x-ray angiography (3DRA) acquisition within a clinically acceptable time and computing from the estimated motion field for the further biomechanical modeling of the cerebrovascular wall. Methods: The whole motion cycle of the cerebral vasculature is modeled using a 4D B-spline transformation, which is estimated from a 4D to 2D + t image registration framework. The registration is performed by optimizing a single similarity metric between the entire 2D + t measured projection sequence and the corresponding forward projections of the deformed volume at their exact time instants. The joint use of two acceleration strategies, together with their implementation on graphics processing units, is also proposed so as to reach computation times close to clinical requirements. For further characterizing vessel wall properties, an approximation of the wall thickness changes is obtained through a strain calculation. Results: Evaluation on in silico and in vitro pulsating phantom aneurysms demonstrated an accurate estimation of wall motion curves. In general, the error was below 10% of the maximum pulsation, even in the situation when substantial inhomogeneous intensity pattern was present. Experiments on in vivo data provided realistic aneurysm and vessel wall motion estimates, whereas in regions where motion was neither visible nor anatomically possible, no motion was detected. The use of the acceleration strategies enabled completing the estimation process for one entire cycle in 5-10 min without degrading the overall performance. The strain map extracted from our motion estimation provided a realistic deformation measure of the vessel wall. Conclusions: The authors' technique has demonstrated that it can provide accurate and robust 4D estimates of cerebrovascular wall motion within a clinically acceptable time, although it has to be applied to a larger patient population prior to possible wide application to routine endovascular procedures. In particular, for the first time, this feasibility study has shown that in vivo cerebrovascular motion can be obtained intraprocedurally from a 3DRA acquisition. Results have also shown the potential of performing strain analysis using this imaging modality, thus making possible for the future modeling of biomechanical properties of the vascular wall.

Edge detection by selection of pieces of level lines

We propose an edge detector based on the selection of wellcontrasted pieces of level lines, following the proposal ofDesolneux-Moisan-Morel (DMM) [1]. The DMM edge detectorhas the problem of over-representation, that is, everyedge is detected several times in slightly different positions.In this paper we propose two modifications of the originalDMM edge detector in order to solve this problem. The firstmodification is a post-processing of the output using a generalmethod to select the best representative of a bundle of curves.The second modification is the use of Canny’s edge detectorinstead of the norm of the gradient to build the statistics. Thetwo modifications are independent and can be applied separately.Elementary reasoning and some experiments showthat the best results are obtained when both modifications areapplied together.

Fermat's treatise on quadrature: A new reading

The Treatise on Quadrature of Fermat (c. 1659), besides containing the first known proof of the computation of the area under a higher parabola, R x+m/n dx, or under a higher hyperbola, R x-m/n dx with the appropriate limits of integration in each case , has a second part which was not understood by Fermat s contemporaries. This second part of the Treatise is obscure and difficult to read and even the great Huygens described it as'published with many mistakes and it is so obscure (with proofs redolent of error) that I have been unable to make any sense of it'. Far from the confusion that Huygens attributes to it, in this paper we try to prove that Fermat, in writing the Treatise, had a very clear goal in mind and he managed to attain it by means of a simple and original method. Fermat reduced the quadrature of a great number of algebraic curves to the quadrature of known curves: the higher parabolas and hyperbolas of the first part of the paper. Others, he reduced to the quadrature of the circle. We shall see how the clever use of two procedures, quite novel at the time: the change of variables and a particular case of the formulaof integration by parts, provide Fermat with the necessary tools to square very easily curves as well-known as the folium of Descartes, the cissoid of Diocles or the witch of Agnesi.

Giffen goods and market making

This paper shows that information effects per se are not responsible forthe Giffen goods anomaly affecting competitive traders demands in multi-asset, noisy rational expectations equilibrium models. The role thatinformation plays in traders strategies also matters. In a market withrisk averse, uninformed traders, informed agents havea dual motive for trading: speculation and market making. Whilespeculation entails using prices to assess the effect of private signalerror terms, market making requires employing them to disentangle noisetraders effects in traders aggregate orders. In a correlated environment,this complicates a trader s signal-extraction problem and maygenerate upward-sloping demand curves. Assuming either (i) that competitive,risk neutral market makers price the assets, or that (ii) the risktolerance coefficient of uninformed traders grows without bound, removesthe market making component from informed traders demands, rendering themwell behaved in prices.

Validation procedures in radiological diagnostic models. Neural network and logistic regression

The objective of this paper is to compare the performance of twopredictive radiological models, logistic regression (LR) and neural network (NN), with five different resampling methods. One hundred and sixty-seven patients with proven calvarial lesions as the only known disease were enrolled. Clinical and CT data were used for LR and NN models. Both models were developed with cross validation, leave-one-out and three different bootstrap algorithms. The final results of each model were compared with error rate and the area under receiver operating characteristic curves (Az). The neural network obtained statistically higher Az than LR with cross validation. The remaining resampling validation methods did not reveal statistically significant differences between LR and NN rules. The neural network classifier performs better than the one based on logistic regression. This advantage is well detected by three-fold cross-validation, but remains unnoticed when leave-one-out or bootstrap algorithms are used.

Bounds for the Betti numbers of generalized Cohen-Macaulay ideals

Upper bounds for the Betti numbers of generalized Cohen-Macaulay ideals are given. In particular, for the case of non-degenerate, reduced and ir- reducible projective curves we get an upper bound which only depends on their degree.

Extensions of maps defined on many fibres

Let S be a fibred surface. We prove that the existence of morphisms from non countably many fibres to curves implies, up to base change, the existence of a rational map from S to another surface fibred over the same base reflecting the properties of the original morphisms. Under some conditions of unicity base change is not needed and one recovers exactly the initial maps.

A note on quadrics through an algebraic curve

In this note we describe the intersection of all quadric hypersur- faces containing a given linearly normal smooth projective curve of genus n and degree 2n + 1

On the binary nature of the γ-ray sources AGL J2241+4454 (= MWC 656) and HESS J0632+057 (= MWC 148)

We present optical spectroscopy of MWC 656 and MWC 148, the proposed optical counterparts of the gamma-ray sources AGL J2241+4454 and HESS J0632+0 57, respectively. The main parameters of the Halpha emission line (EW, FWHM and centroid velocity) in these stars are modulated on the proposed orbital periods of 60.37 and 321 days, respectively. These modulations are likely produced by the resonant interaction of the Be discs with compact stars in eccentric orbits. We also present radial velocity curves of the optical stars folded on the above periods and obtain the first orbital elements of the two gamma-ray sources thus confirming their binary nature. Our orbital solution support eccentricities e~0.4 and 0.83+-0.08 for MWC 656 and MWC 148, respectively. Further, our orbital elements imply that the X-ray outbursts in HESS J0632+057/MWC 148 are delayed ~0.3 orbital phases after periastron passage, similarly to the case of LS I +61 303. In addition, the optical photometric light curve maxima in AGL J2241+4454/MWC 656 occur ~0.25 phases passed periastron, similar to what is seen in LS I +61 303. We also find that the orbital eccentricity is correlated with orbital period for the known gamma-ray binaries. This is explained by the fact that small stellar separations are required for the efficient triggering of VHE radiation. Another correlation between the EW of Halpha and orbital period is also observed, similarly to the case of Be/X-ray binaries. These correlations are useful to provide estimates of the key orbital parameters Porb and e from the Halpha line in future Be gamma-ray binary candidates.

Time-scale and other invariants of integrative mechanical behavior in living cells.

In dealing with systems as complex as the cytoskeleton, we need organizing principles or, short of that, an empirical framework into which these systems fit. We report here unexpected invariants of cytoskeletal behavior that comprise such an empirical framework. We measured elastic and frictional moduli of a variety of cell types over a wide range of time scales and using a variety of biological interventions. In all instances elastic stresses dominated at frequencies below 300 Hz, increased only weakly with frequency, and followed a power law; no characteristic time scale was evident. Frictional stresses paralleled the elastic behavior at frequencies below 10 Hz but approached a Newtonian viscous behavior at higher frequencies. Surprisingly, all data could be collapsed onto master curves, the existence of which implies that elastic and frictional stresses share a common underlying mechanism. Taken together, these findings define an unanticipated integrative framework for studying protein interactions within the complex microenvironment of the cell body, and appear to set limits on what can be predicted about integrated mechanical behavior of the matrix based solely on cytoskeletal constituents considered in isolation. Moreover, these observations are consistent with the hypothesis that the cytoskeleton of the living cell behaves as a soft glassy material, wherein cytoskeletal proteins modulate cell mechanical properties mainly by changing an effective temperature of the cytoskeletal matrix. If so, then the effective temperature becomes an easily quantified determinant of the ability of the cytoskeleton to deform, flow, and reorganize.

Term structure of interest rates. European financial integration

In this paper we estimate, analyze and compare the term structures of interest rates in six different countries over the period 1992-2004. We apply the Nelson-Siegel model to obtain the term structures of interest rates at weekly intervals. A total of 4,038 curves are estimated and analyzed. Four European Monetary Union countries¿Spain, France, Germany and Italy¿are included. The UK is also included as a European non-member of the Monetary Union. Finally the US completes the analysis. The goal is to determine the differences in the shapes of the term structure of interest rates among these countries. Likewise, we can determine the most usual term structure shapes that appear for each country.*****

Phonon softening in Ni-Mn-Ga alloys

The TA2 phonon dispersion curves of Ni-Mn-Ga alloys with different compositions which transform to different martensitic structures have been measured over a broad temperature range covering both paramagnetic and ferromagnetic phases. The branches show an anomaly (dip) at a wave number that depends on the particular martensitic structure, and there is softening of these anomalous phonons with decreasing temperature. This softening is enhanced below the Curie point, as a consequence of spin-phonon coupling. This effect is stronger for systems with higher electronic concentration.