Collective vibrational states with fast iterative QRPA method

An iterative method we previously proposed to compute nuclear strength functions is developed to allow it to accurately calculate properties of individual nuclear states. The approach is based on the quasi-particle-random-phase approximation (QRPA) and uses an iterative non-hermitian Arnoldi diagonalization method where the QRPA matrix does not have to be explicitly calculated and stored. The method gives substantial advantages over conventional QRPA calculations with regards to the computational cost. The method is used to calculate excitation energies and decay rates of the lowest lying 2+ and 3- states in Pb, Sn, Ni and Ca isotopes using three different Skyrme interactions and a separable gaussian pairing force.Comment: 10 pages, 11 figure

Determination of propofol by GC/MS and fast GC/MS-TOF in two cases of poisoning

Two cases of suspected acute and lethal intoxication caused by propofol were delivered by the judicial authority to the Department of Sciences for Health Promotion and Mother-Child Care in Palermo, Sicily. In the first case a female nurse was found in a hotel room, where she lived with her mother; four 10 mg/mL vials and two 20 mg/mL vials of propofol were found near the decedent along with syringes and needles. In the second case a male nurse was found in the operating room of a hospital, along with a used syringe. In both cases a preliminary systematic and toxicological analysis indicated the presence of propofol in the blood and urine. As a result, a method for the quantitative determination of propofol in biological fluids was optimized and validated using a liquid-liquid extraction protocol followed by GC/MS and fast GC/MS-TOF. In the first case, the concentration of propofol in blood was determined to be 8.1 \u3bcg/mL while the concentration of propofol in the second case was calculated at 1.2 \u3bcg/mL. Additionally, the tissue distribution of propofol was determined for both cases. Brain and liver concentrations of propofol were, respectively, 31.1 and 52.2 \u3bcg/g in Case 1 and 4.7 and 49.1 \u3bcg/g in Case 2. Data emerging from the autopsy findings, histopathological exams as well as the toxicological results aided in establishing that the deaths were due to poisoning, however, the manner of death in each were different: homicide in Case 1 and suicide in Case 2

Kinematic dynamo action in a sphere. I. Effects of differential rotation and meridional circulation on solutions with axial dipole symmetry

A sphere containing electrically conducting fluid can generate a magnetic field by dynamo action, provided the flow is sufficiently complicated and vigorous. The dynamo mechanism is thought to sustain magnetic fields in planets and stars. The kinematic dynamo problem tests steady flows for magnetic instability, but rather few dynamos have been found so far because of severe numerical difficulties. Dynamo action might, therefore, be quite unusual, at least for large-scale steady flows. We address this question by testing a two-parameter class of flows for dynamo generation of magnetic fields containing an axial dipole. The class of flows includes two completely different types of known dynamos, one dominated by differential rotation (D) and one with none. We find that 36% of the flows in seven distinct zones in parameter space act as dynamos, while the remaining 64% either fail to generate this type of magnetic field or generate fields that are too small in scale to be resolved by our numerical method. The two previously known dynamo types lie in the same zone, and it is therefore possible to change the flow continuously from one to the other without losing dynamo action. Differential rotation is found to promote large-scale axisymmetric toroidal magnetic fields, while meridional circulation (M) promotes large-scale axisymmetric poloidal fields concentrated at high latitudes near the axis. Magnetic fields resembling that of the Earth are generated by D > 0, corresponding to westward flow at the surface, and M of either sign but not zero. Very few oscillatory solutions are found

Bogoliubov modes of a dipolar condensate in a cylindrical trap

The calculation of properties of Bose-Einstein condensates with dipolar interactions has proven a computationally intensive problem due to the long range nature of the interactions, limiting the scope of applications. In particular, the lowest lying Bogoliubov excitations in three dimensional harmonic trap with cylindrical symmetry were so far computed in an indirect way, by Fourier analysis of time dependent perturbations, or by approximate variational methods. We have developed a very fast and accurate numerical algorithm based on the Hankel transform for calculating properties of dipolar Bose-Einstein condensates in cylindrically symmetric traps. As an application, we are able to compute many excitation modes by directly solving the Bogoliubov-De Gennes equations. We explore the behavior of the excited modes in different trap geometries. We use these results to calculate the quantum depletion of the condensate by a combination of a computation of the exact modes and the use of a local density approximation

Spontaneous breathing pattern as respiratory functional outcome in children with spinal muscular atrophy (SMA)

Introduction: SMA is characterised by progressive motor and respiratory muscle weakness. We aimed to verify if in SMA children 1)each form is characterized by specific ventilatory and thoracoabdominal pattern(VTAp) during quiet breathing(QB); 2)VTAp is affected by salbutamol therapy, currently suggested as standard treatment, or by the natural history(NH) of SMA; 3)the severity of global motor impairment linearly correlates with VTAp. Materials and methods: VTAp was analysed on 32 SMA type I (SMA1, the most severe form), 51 type II (SMA2, the moderate), 8 type III (SMA3, the mildest) and 20 healthy (HC) using opto-electronic plethysmography. Spirometry, cough and motor function were measured in a subgroup of patients. Results: In SMA1, a normal ventilation is obtained in supine position by rapid and shallow breathing with paradoxical ribcage motion. In SMA2, ventilation is within a normal range in seated position due to an increased respiratory rate(p0.05) while tachypnea occurred in type I NH. A linear correlation(p<0.001) was found between motor function scales and VTAp. Conclusion: A negative or reduced %ΔVRC,P, indicative of ribcage muscle weakness, is a distinctive feature of SMA1 and SMA2 since infancy. Its quantitative assessment represents a non-invasive, non-volitional index that can be obtained in all children, even uncollaborative, and provides useful information on the action of ribcage muscles that are known to be affected by the disease. Low values of motor function scales indicate impairment of motor but also of respiratory function

Kinematic dynamo action in a sphere: Effects of periodic time-dependent flows on solutions with axial dipole symmetry

Choosing a simple class of flows, with characteristics that may be present in the Earth's core, we study the ability to generate a magnetic field when the flow is permitted to oscillate periodically in time. The flow characteristics are parameterised by D, representing a differential rotation, M, a meridional circulation, and C, a component characterising convective rolls. Dynamo action is sensitive to these flow parameters and fails spectacularly for much of the parameter space where magnetic flux is concentrated into small regions. Oscillations of the flow are introduced by varying the flow parameters in time, defining a closed orbit in the space (D,M). Time-dependence appears to smooth out flux concentrations, often enhancing dynamo action. Dynamo action can be impaired, however, when flux concentrations of opposite signs occur close together as smoothing destroys the flux by cancellation. It is possible to produce geomagnetic-type reversals by making the orbit stray into a region where the steady flows generate oscillatory fields. In this case, however, dynamo action was not found to be enhanced by the time-dependence. A novel approach is taken to solving the time-dependent eigenvalue problem, where by combining Floquet theory with a matrix-free Krylov-subspace method we avoid large memory requirements for storing the matrix required by the standard approach.Comment: 22 pages, 12 figures. Geophys. Astrophys. Fluid Dynam., as accepted (2004

Stability and Complexity of Minimising Probabilistic Automata

We consider the state-minimisation problem for weighted and probabilistic automata. We provide a numerically stable polynomial-time minimisation algorithm for weighted automata, with guaranteed bounds on the numerical error when run with floating-point arithmetic. Our algorithm can also be used for "lossy" minimisation with bounded error. We show an application in image compression. In the second part of the paper we study the complexity of the minimisation problem for probabilistic automata. We prove that the problem is NP-hard and in PSPACE, improving a recent EXPTIME-result.Comment: This is the full version of an ICALP'14 pape

Diffusion with critically correlated traps and the slow relaxation of the longest wavelength mode

We study diffusion on a substrate with permanent traps distributed with critical positional correlation, modeled by their placement on the perimeters of a critical percolation cluster. We perform a numerical analysis of the vibrational density of states and the largest eigenvalue of the equivalent scalar elasticity problem using the method of Arnoldi and Saad. We show that the critical trap correlation increases the exponent appearing in the stretched exponential behavior of the low frequency density of states by approximately a factor of two as compared to the case of no correlations. A finite size scaling hypothesis of the largest eigenvalue is proposed and its relation to the density of states is given. The numerical analysis of this scaling postulate leads to the estimation of the stretch exponent in good agreement with the density of states result.Comment: 15 pages, LaTeX (RevTeX

Markov chain analysis of random walks on disordered medium

We study the dynamical exponents $d_{w}$ and $d_{s}$ for a particle diffusing in a disordered medium (modeled by a percolation cluster), from the regime of extreme disorder (i.e., when the percolation cluster is a fractal at $p=p_{c}$) to the Lorentz gas regime when the cluster has weak disorder at $p>p_{c}$ and the leading behavior is standard diffusion. A new technique of relating the velocity autocorrelation function and the return to the starting point probability to the asymptotic spectral properties of the hopping transition probability matrix of the diffusing particle is used, and the latter is numerically analyzed using the Arnoldi-Saad algorithm. We also present evidence for a new scaling relation for the second largest eigenvalue in terms of the size of the cluster, $|\ln{\lambda}_{max}|\sim S^{-d_w/d_f}$, which provides a very efficient and accurate method of extracting the spectral dimension $d_s$ where $d_s=2d_f/d_w$.Comment: 34 pages, REVTEX 3.