Pinsker estimators for local helioseismology

A major goal of helioseismology is the three-dimensional reconstruction of the three velocity components of convective flows in the solar interior from sets of wave travel-time measurements. For small amplitude flows, the forward problem is described in good approximation by a large system of convolution equations. The input observations are highly noisy random vectors with a known dense covariance matrix. This leads to a large statistical linear inverse problem. Whereas for deterministic linear inverse problems several computationally efficient minimax optimal regularization methods exist, only one minimax-optimal linear estimator exists for statistical linear inverse problems: the Pinsker estimator. However, it is often computationally inefficient because it requires a singular value decomposition of the forward operator or it is not applicable because of an unknown noise covariance matrix, so it is rarely used for real-world problems. These limitations do not apply in helioseismology. We present a simplified proof of the optimality properties of the Pinsker estimator and show that it yields significantly better reconstructions than traditional inversion methods used in helioseismology, i.e.\ Regularized Least Squares (Tikhonov regularization) and SOLA (approximate inverse) methods. Moreover, we discuss the incorporation of the mass conservation constraint in the Pinsker scheme using staggered grids. With this improvement we can reconstruct not only horizontal, but also vertical velocity components that are much smaller in amplitude

Self-similar solutions with fat tails for Smoluchowski's coagulation equation with locally bounded kernels

The existence of self-similar solutions with fat tails for Smoluchowski's coagulation equation has so far only been established for the solvable and the diagonal kernel. In this paper we prove the existence of such self-similar solutions for continuous kernels $K$ that are homogeneous of degree $\gamma \in [0,1)$ and satisfy $K(x,y) \leq C (x^{\gamma} + y^{\gamma})$. More precisely, for any $\rho \in (\gamma,1)$ we establish the existence of a continuous weak self-similar profile with decay $x^{-(1{+}\rho)}$ as $x \to \infty$

Bi-defects of Nematic Surfactant Bilayers

We consider the effects of the coupling between the orientational order of the two monolayers in flat nematic bilayers. We show that the presence of a topological defect on one bilayer generates a nontrivial orientational texture on both monolayers. Therefore, one cannot consider isolated defects on one monolayer, but rather associated pairs of defects on either monolayer, which we call bi-defects. Bi-defects generally produce walls, such that the textures of the two monolayers are identical outside the walls, and different in their interior. We suggest some experimental conditions in which these structures could be observed.Comment: RevTeX, 4 pages, 3 figure

Universal analytic properties of noise. Introducing the J-Matrix formalism

We propose a new method in the spectral analysis of noisy time-series data for damped oscillators. From the Jacobi three terms recursive relation for the denominators of the Pad\'e Approximations built on the well-known Z-transform of an infinite time-series, we build an Hilbert space operator, a J-Operator, where each bound state (inside the unit circle in the complex plane) is simply associated to one damped oscillator while the continuous spectrum of the J-Operator, which lies on the unit circle itself, is shown to represent the noise. Signal and noise are thus clearly separated in the complex plane. For a finite time series of length 2N, the J-operator is replaced by a finite order J-Matrix J_N, having N eigenvalues which are time reversal covariant. Different classes of input noise, such as blank (white and uniform), Gaussian and pink, are discussed in detail, the J-Matrix formalism allowing us to efficiently calculate hundreds of poles of the Z-transform. Evidence of a universal behaviour in the final statistical distribution of the associated poles and zeros of the Z-transform is shown. In particular the poles and zeros tend, when the length of the time series goes to infinity, to a uniform angular distribution on the unit circle. Therefore at finite order, the roots of unity in the complex plane appear to be noise attractors. We show that the Z-transform presents the exceptional feature of allowing lossless undersampling and how to make use of this property. A few basic examples are given to suggest the power of the proposed method.Comment: 14 pages, 8 figure

Effect of latitudinal differential rotation on solar Rossby waves: Critical layers, eigenfunctions, and momentum fluxes in the equatorial β\beta plane

Retrograde-propagating waves of vertical vorticity with longitudinal wavenumbers between 3 and 15 have been observed on the Sun with a dispersion relation close to that of classical sectoral Rossby waves. The observed vorticity eigenfunctions are symmetric in latitude, peak at the equator, switch sign near $20^\circ$-$30^\circ$, and decrease at higher latitudes. We search for an explanation that takes into account solar latitudinal differential rotation. In the equatorial $\beta$ plane, we study the propagation of linear Rossby waves (phase speed $c <0$) in a parabolic zonal shear flow, $U = - \overline{U}\ \xi^2<0$, where $\overline{U} = 244$ m/s and $\xi$ is the sine of latitude. In the inviscid case, the eigenvalue spectrum is real and continuous and the velocity stream functions are singular at the critical latitudes where $U = c$. We add eddy viscosity in the problem to account for wave attenuation. In the viscous case, the stream functions are solution of a fourth-order modified Orr-Sommerfeld equation. Eigenvalues are complex and discrete. For reasonable values of the eddy viscosity corresponding to supergranular scales and above (Reynolds number $100 \le Re \le 700$), all modes are stable. At fixed longitudinal wavenumber, the least damped mode is a symmetric mode with a real frequency close to that of the classical Rossby mode, which we call the R mode. For $Re \approx 300$, the attenuation and the real part of the eigenfunction is in qualitative agreement with the observations (unlike the imaginary part of the eigenfunction, which has a larger amplitude in the model. Conclusion: Each longitudinal wavenumber is associated with a latitudinally symmetric R mode trapped at low latitudes by solar differential rotation. In the viscous model, R modes transport significant angular momentum from the dissipation layers towards the equator.Comment: Submitted to Astron. Astrophys. on 29 May 202

Asymptotics of self-similar solutions to coagulation equations with product kernel

We consider mass-conserving self-similar solutions for Smoluchowski's coagulation equation with kernel $K(\xi,\eta)= (\xi \eta)^{\lambda}$ with $\lambda \in (0,1/2)$. It is known that such self-similar solutions $g(x)$ satisfy that $x^{-1+2\lambda} g(x)$ is bounded above and below as $x \to 0$. In this paper we describe in detail via formal asymptotics the qualitative behavior of a suitably rescaled function $h(x)=h_{\lambda} x^{-1+2\lambda} g(x)$ in the limit $\lambda \to 0$. It turns out that $h \sim 1+ C x^{\lambda/2} \cos(\sqrt{\lambda} \log x)$ as $x \to 0$. As $x$ becomes larger $h$ develops peaks of height $1/\lambda$ that are separated by large regions where $h$ is small. Finally, $h$ converges to zero exponentially fast as $x \to \infty$. Our analysis is based on different approximations of a nonlocal operator, that reduces the original equation in certain regimes to a system of ODE

Improving the Acoustic Performance of Linear Multi-Element Transducers

The electro-acoustic performance of transducers has a direct impact on the performance of ultrasound inspections. The signal/noise ratio and the resolution (both axial and lateral) are key factors for detecting and/or proportioning the indications being sought. The signal/noise ratio partly depends on the sensitivity and the signal/noise ratio of the transducer itself. The axial resolution depends on the length of the signal and therefore, for a given maximum frequency, on the damping of the transducer. Sensitivity and damping are often considered antagonistic, as damping traditionally reduces resonance and therefore sensitivity. Earlier studies have demonstrated the advantages gained through using piezocomposite technology to improve this compromise. These two parameters also depend on the acoustic adaptation to the coupling medium (water, plexiglass, rexolite, steel, etc.), and according to the design used, performance deteriorates more or less as one moves further from the nominal use. In addition to sensitivity and the signal/noise ratio, other parameters such as the angular acceptance and resistance to abrasion are sometimes to be integrated in the expected performances. This article presents the recent developments undertaken and tested in the context of improving the acoustic performance of multi-element probes: - Identification of the components that influence performance; - Simulations; - Selection of the configurations that meet the needs of various applications; - The experimental results obtained; - Comparison with the simulations. These studies have led to the development of a design expertise for responding to requests for custom-made, industrial, multi-element probes with improved performance, for production runs from a single item to dozens, even hundreds. The detailed results will be presented, as well as the possibilities for future development

The scaling attractor and ultimate dynamics for Smoluchowski's coagulation equations

We describe a basic framework for studying dynamic scaling that has roots in dynamical systems and probability theory. Within this framework, we study Smoluchowski's coagulation equation for the three simplest rate kernels $K(x,y)=2$, $x+y$ and $xy$. In another work, we classified all self-similar solutions and all universality classes (domains of attraction) for scaling limits under weak convergence (Comm. Pure Appl. Math 57 (2004)1197-1232). Here we add to this a complete description of the set of all limit points of solutions modulo scaling (the scaling attractor) and the dynamics on this limit set (the ultimate dynamics). The main tool is Bertoin's L\'{e}vy-Khintchine representation formula for eternal solutions of Smoluchowski's equation (Adv. Appl. Prob. 12 (2002) 547--64). This representation linearizes the dynamics on the scaling attractor, revealing these dynamics to be conjugate to a continuous dilation, and chaotic in a classical sense. Furthermore, our study of scaling limits explains how Smoluchowski dynamics ``compactifies'' in a natural way that accounts for clusters of zero and infinite size (dust and gel)

Control of mucocutaneous leishmaniasis, a neglected disease: results of a control programme in Satipo Province, Peru.

Mucocutaneous leishmaniasis (MCL) is an important health problem in many rural areas of Latin America, but there are few data on the results of programmatic approaches to control the disease. We report the results of a control programme in San Martin de Pangoa District, which reports one of the highest prevalences of MCL in Peru. For 2 years (2001--2002), the technicians at the health post were trained in patient case management, received medical support and were supplied with antimonials. An evaluation after 2 years showed the following main achievements: better diagnosis of patients, who were confirmed by microscopy in 34% (82/240) of the cases in 2001 and 60% of the cases (153/254) in 2002; improved follow-up during treatment: 237 of 263 (90%) patients who initiated an antimonial therapy ended the full treatment course; improved follow-up after treatment: 143 of 237 (60%) patients who ended their full treatment were correctly monitored during the required period of 6 (cutaneous cases) or 12 (mucosal cases) months after the end of treatment. These achievements were largely due to the human and logistical resources made available, the constant availability of medications and the close collaboration between the Ministry of Health, a national research institute and an international non-governmental organization. At the end of this period, the health authorities decided to register a generic brand of sodium stibogluconate, which is now in use. This should allow the treatment of a significant number of additional patients, while saving money to invest in other facets of the case management