Search

We consider $N$ classical particles interacting via the Coulomb potential in spatial dimension $d$ and in the presence of an external trap, at equilibrium at inverse temperature $\beta$. In the large $N$ limit, the particles are confined within a droplet of finite size. We study smooth linear statistics, i.e. the fluctuations of sums of the form ${\cal L}_N = \sum_{i=1}^N f({\bf x}_i)$, where ${\bf x}_i$'s are the positions of the particles and where $f({\bf x}_i)$ is a sufficiently regular function. There exists at present standard results for the first and second moments of ${\cal L}_N$ in the large $N$ limit, as well as associated Central Limit Theorems in general dimension and for a wide class of confining potentials. Here we obtain explicit expressions for the higher order cumulants of ${\cal L}_N$ at large $N$, when the function $f({\bf x})=f(|{\bf x}|)$ and the confining potential are both rotationnally invariant. A remarkable feature of our results is that these higher cumulants depend only on the value of $f'(|{\bf x}|)$ and its higher order derivatives evaluated exactly at the boundary of the droplet, which in this case is a $d$-dimensional sphere. In the particular two-dimensional case $d=2$ at the special value $\beta=2$, a connection to the Ginibre ensemble allows us to derive these results in an alternative way using the tools of determinantal point processes. Finally we also obtain the large deviation form of the full probability distribution function of ${\cal L}_N$., Comment: 19 pages

This thesis is devoted to the study of extreme value statistics in stochastic processes and their applications. In the first part, we obtain exact analytical results on the extreme value statistics of both discrete-time and continuous-time random walks. In particular, we focus on the gap statistics of random walks and exhibit their asymptotic universality with respect to the jump distribution in the limit of a large number of steps. In addition, we compute the asymptotic behavior of the expected maximum of random walks in the presence of a bridge constraint and reveal a rich behavior in their finite-size correction. Moreover, we compute the expected length of the convex hull of Brownian motion confined in a disk and show that it converges slowly to the perimeter of the disk with a stretched exponential decay. In the second part, we focus on numerically sampling rare trajectories of stochastic processes. We introduce an efficient method to sample bridge discrete-time random walks. We illustrate it and apply it to various examples. We further extend the method to other stochastic processes, both Markovian and non-Markovian. We apply our method to sample surviving particles in the presence of a periodic trapping environment. Finally, we discuss several optimization problems in stochastic processes involving extreme value statistics. In particular, we introduce a new technique to optimally control dynamical systems undergoing a resetting policy., Comment: 134 pages, 54 figures

We consider a one-dimensional Brownian motion with diffusion coefficient $D$ in the presence of $n$ partially absorbing traps with intensity $\beta$, separated by a distance $L$ and evenly spaced around the initial position of the particle. We study the transport properties of the process conditioned to survive up to time $t$. We find that the surviving particle first diffuses normally, before it encounters the traps, then undergoes a period of transient anomalous diffusion, after which it reaches a final diffusive regime. The asymptotic regime is governed by an effective diffusion coefficient $D_\text{eff}$, which is induced by the trapping environment and is typically different from the original one. We show that when the number of traps is \emph{finite}, the environment enhances diffusion and induces an effective diffusion coefficient that is systematically equal to $D_\text{eff}=2D$, independently of the number of the traps, the trapping intensity $\beta$ and the distance $L$. On the contrary, when the number of traps is \emph{infinite}, we find that the environment inhibits diffusion with an effective diffusion coefficient that depends on the traps intensity $\beta$ and the distance $L$ through a non-trivial scaling function $D_\text{eff}=D \mathcal{F}(\beta L/D)$, for which we obtain a closed-form. Moreover, we provide a rejection-free algorithm to generate surviving trajectories by deriving an effective Langevin equation with an effective repulsive potential induced by the traps. Finally, we extend our results to other trapping environments., Comment: 28 pages, 10 figures

We consider one-dimensional discrete-time random walks (RWs) of $n$ steps, starting from $x_0=0$, with arbitrary symmetric and continuous jump distributions $f(\eta)$, including the important case of L\'evy flights. We study the statistics of the gaps $\Delta_{k,n}$ between the $k^\text{th}$ and $(k+1)^\text{th}$ maximum of the set of positions $\{x_1,\ldots,x_n\}$. We obtain an exact analytical expression for the probability distribution $P_{k,n}(\Delta)$ valid for any $k$ and $n$, and jump distribution $f(\eta)$, which we then analyse in the large $n$ limit. For jump distributions whose Fourier transform behaves, for small $q$, as $\hat f (q) \sim 1 - |q|^\mu$ with a L\'evy index $0< \mu \leq 2$, we find that, the distribution becomes stationary in the limit of $n\to \infty$, i.e. $\lim_{n\to \infty} P_{k,n}(\Delta)=P_k(\Delta)$. We obtain an explicit expression for its first moment $\mathbb{E}[\Delta_{k}]$, valid for any $k$ and jump distribution $f(\eta)$ with $\mu>1$, and show that it exhibits a universal algebraic decay $ \mathbb{E}[\Delta_{k}]\sim k^{1/\mu-1} \Gamma\left(1-1/\mu\right)/\pi$ for large $k$. Furthermore, for $\mu>1$, we show that in the limit of $k\to\infty$ the stationary distribution exhibits a universal scaling form $P_k(\Delta) \sim k^{1-1/\mu} \mathcal{P}_\mu(k^{1-1/\mu}\Delta)$ which depends only on the L\'evy index $\mu$, but not on the details of the jump distribution. We compute explicitly the limiting scaling function $\mathcal{P}_\mu(x)$ in terms of Mittag-Leffler functions. For $1< \mu <2$, we show that, while this scaling function captures the distribution of the typical gaps on the scale $k^{1/\mu-1}$, the atypical large gaps are not described by this scaling function since they occur at a larger scale of order $k^{1/\mu}$., Comment: 40 pages, 11 figures

We investigate a moving boundary problem for a Brownian particle on the semi-infinite line in which the boundary moves by a distance proportional to the time between successive collisions of the particle and the boundary. Phenomenologically rich dynamics arises. In particular, the probability for the particle to first reach the moving boundary for the $n^\text{th}$ time asymptotically scales as $t^{-(1+2^{-n})}$. Because the tail of this distribution becomes progressively fatter, the typical time between successive first passages systematically gets longer. We also find that the number of collisions between the particle and the boundary scales as $\ln\ln t$, while the time dependence of the boundary position varies as $t/\ln t$., Comment: 11 pages, 4 figure, for a special issue of JPA on resetting processes

We introduce a resetting Brownian bridge as a simple model to study search processes where the total search time $t_f$ is finite and the searcher returns to its starting point at $t_f$. This is simply a Brownian motion with a Poissonian resetting rate $r$ to the origin which is constrained to start and end at the origin at time $t_f$. We first provide a rejection-free algorithm to generate such resetting bridges in all dimensions by deriving an effective Langevin equation with an explicit space-time dependent drift $\tilde \mu({\bf x},t)$ and resetting rate $\tilde r({\bf x}, t)$. We also study the efficiency of the search process in one-dimension by computing exactly various observables such as the mean-square displacement, the hitting probability of a fixed target and the expected maximum. Surprisingly, we find that there exists an optimal resetting rate $r^*$ that maximizes the search efficiency, even in the presence of a bridge constraint. We show however that the physical mechanism responsible for this optimal resetting rate for bridges is entirely different from resetting Brownian motions without the bridge constraint., Comment: Main text: 6 pages + 3 figs, Supp. Mat.: 9 pages + 3 figs

``When in a difficult situation, it is sometimes better to give up and start all over again''. While this empirical truth has been regularly observed in a wide range of circumstances, quantifying the effectiveness of such a heuristic strategy remains an open challenge. In this paper, we combine the notions of optimal control and stochastic resetting to address this problem. The emerging analytical framework allows not only to measure the performance of a given restarting policy but also to obtain the optimal strategy for a wide class of dynamical systems. We apply our technique to a system with a final reward and show that the reward value must be larger than a critical threshold for resetting to be effective. Our approach, analogous to the celebrated Hamilton-Jacobi-Bellman paradigm, provides the basis for the investigation of realistic restarting strategies across disciplines. As an application, we show that the framework can be applied to an epidemic model to predict the optimal lockdown policy., Comment: 11 pages, 7 figures

We consider a Brownian particle with diffusion coefficient $D$ in a $d$-dimensional ball of radius $R$ with reflecting boundaries. We study the maximum $M_x(t)$ of the trajectory of the particle along the $x$-direction at time $t$. In the long time limit, the maximum converges to the radius of the ball $M_x(t) \to R$ for $t\to \infty$. We investigate how this limit is approached and obtain an exact analytical expression for the distribution of the fluctuations $\Delta(t) = [R-M_x(t)]/R$ in the limit of large $t$ in all dimensions. We find that the distribution of $\Delta(t)$ exhibits a rich variety of behaviors depending on the dimension $d$. These results are obtained by establishing a connection between this problem and the narrow escape time problem. We apply our results in $d=2$ to study the convex hull of the trajectory of the particle in a disk of radius $R$ with reflecting boundaries. We find that the mean perimeter $\langle L(t)\rangle$ of the convex hull exhibits a slow convergence towards the perimeter of the circle $2\pi R$ with a stretched exponential decay $2\pi R-\langle L(t)\rangle \propto \sqrt{R}(Dt)^{1/4} \,e^{-2\sqrt{2Dt}/R}$. Finally, we generalise our results to other confining geometries, such as the ellipse with reflecting boundaries. Our results are corroborated by thorough numerical simulations., Comment: 19 pages, 7 figures

We propose a method to exactly generate Brownian paths $x_c(t)$ that are constrained to return to the origin at some future time $t_f$, with a given fixed area $A_f = \int_0^{t_f}dt\, x_c(t)$ under their trajectory. We derive an exact effective Langevin equation with an effective force that accounts for the constraint. In addition, we develop the corresponding approach for discrete-time random walks, with arbitrary jump distributions including L\'evy flights, for which we obtain an effective jump distribution that encodes the constraint. Finally, we generalise our method to other types of dynamical constraints such as a fixed occupation time on the positive axis $T_f=\int_0^{t_f}dt\, \Theta\left[x_c(t)\right]$ or a fixed generalised quadratic area $\mathcal{A}_f=\int_0^{t_f}dt \,x_c^2(t)$., Comment: 32 pages, 7 figures

We consider one-dimensional discrete-time random walks (RWs) in the presence of finite size traps of length $\ell$ over which the RWs can jump. We study the survival probability of such RWs when the traps are periodically distributed and separated by a distance $L$. We obtain exact results for the mean first-passage time and the survival probability in the special case of a double-sided exponential jump distribution. While such RWs typically survive longer than if they could not leap over traps, their survival probability still decreases exponentially with the number of steps. The decay rate of the survival probability depends in a non-trivial way on the trap length $\ell$ and exhibits an interesting regime when $\ell\rightarrow 0$ as it tends to the ratio $\ell/L$, which is reminiscent of strongly chaotic deterministic systems. We generalize our model to continuous-time RWs, where we introduce a power-law distributed waiting time before each jump. In this case, we find that the survival probability decays algebraically with an exponent that is independent of the trap length. Finally, we derive the diffusive limit of our model and show that, depending on the chosen scaling, we obtain either diffusion with uniform absorption, or diffusion with periodically distributed point absorbers., Comment: 29 pages, 5 figures

We introduce a minimalist dynamical model of wealth evolution and wealth sharing among $N$ agents as a platform to compare the relative merits of altruism and individualism. In our model, the wealth of each agent independently evolves by diffusion. For a population of altruists, whenever any agent reaches zero wealth (that is, the agent goes bankrupt), the remaining wealth of the other $N-1$ agents is equally shared among all. The population is collectively defined to be bankrupt when its total wealth falls below a specified small threshold value. For individualists, each time an agent goes bankrupt (s)he is considered to be "dead" and no wealth redistribution occurs. We determine the evolution of wealth in these two societies. Altruism leads to more global median wealth at early times; eventually, however, the longest-lived individualists accumulate most of the wealth and are richer and more long lived than the altruists., Comment: 16 pages, 6 figures, IOP format

We propose a method to exactly generate bridge run-and-tumble trajectories that are constrained to start at the origin with a given velocity and to return to the origin after a fixed time with another given velocity. The method extends the concept of effective Langevin equations, valid for Markovian stochastic processes such as Brownian motion, to a non-Markovian stochastic process driven by a telegraphic noise, with exponentially decaying correlations. We obtain effective space-time dependent tumbling rates that implicitly accounts for the bridge constraint. We extend the method to other types of constrained run-and-tumble particles such as excursions and meanders. The method is implemented numerically and is shown to be very efficient., Comment: 21 pages, 7 figures

We consider one-dimensional discrete-time random walks (RWs) with arbitrary symmetric and continuous jump distributions $f(\eta)$, including the case of L\'evy flights. We study the expected maximum ${\mathbb E}[M_n]$ of bridge RWs, i.e., RWs starting and ending at the origin after $n$ steps. We obtain an exact analytical expression for ${\mathbb E}[M_n]$ valid for any $n$ and jump distribution $f(\eta)$, which we then analyze in the large $n$ limit up to second leading order term. For jump distributions whose Fourier transform behaves, for small $k$, as $\hat f(k) \sim 1 - |a\, k|^\mu$ with a L\'evy index $0<\mu \leq 2$ and an arbitrary length scale $a>0$, we find that, at leading order for large $n$, ${\mathbb E}[M_n]\sim a\, h_1(\mu)\, n^{1/\mu}$. We obtain an explicit expression for the amplitude $h_1(\mu)$ and find that it carries the signature of the bridge condition, being different from its counterpart for the free random walk. For $\mu=2$, we find that the second leading order term is a constant, which, quite remarkably, is the same as its counterpart for the free RW. For generic $0< \mu < 2$, this second leading order term is a growing function of $n$, which depends non-trivially on further details of $\hat f (k)$, beyond the L\'evy index $\mu$. Finally, we apply our results to compute the mean perimeter of the convex hull of the $2d$ Rouse polymer chain and of the $2d$ run-and-tumble particle, as well as to the computation of the survival probability in a bridge version of the well-known "lamb-lion" capture problem., Comment: 31 pages, 6 figures

We introduce a method to exactly generate bridge trajectories for discrete-time random walks, with arbitrary jump distributions, that are constrained to initially start at the origin and return to the origin after a fixed time. The method is based on an effective jump distribution that implicitly accounts for the bridge constraint. It is illustrated on various jump distributions and is shown to be very efficient in practice. In addition, we show how to generalize the method to other types of constrained random walks such as generalized bridges, excursions, and meanders., Comment: 18 pages, 11 figures

The Wigner function $W_N({\bf x}, {\bf p})$ is a useful quantity to characterize the quantum fluctuations of an $N$-body system in its phase space. Here we study $W_N({\bf x}, {\bf p})$ for $N$ noninteracting spinless fermions in a $d$-dimensional spherical hard box of radius $R$ at temperature $T=0$. In the large $N$ limit, the local density approximation (LDA) predicts that $W_N({\bf x}, {\bf p}) \approx 1/(2 \pi \hbar)^d$ inside a finite region of the $({\bf x}, {\bf p})$ plane, namely for $|{\bf x}| < R$ and $|{\bf p}| < k_F$ where $k_F$ is the Fermi momentum, while $W_N({\bf x}, {\bf p})$ vanishes outside this region, or "droplet", on a scale determined by quantum fluctuations. In this paper we investigate systematically, in this quantum region, the structure of the Wigner function along the edge of this droplet, called the Fermi surf. In one dimension, we find that there are three distinct edge regions along the Fermi surf and we compute exactly the associated nontrivial scaling functions in each regime. We also study the momentum distribution $\hat \rho_N(p)$ and find a striking algebraic tail for very large momenta $\hat \rho_N(p) \propto 1/p^4$, well beyond $k_F$, reminiscent of a similar tail found in interacting quantum systems (discussed in the context of Tan's relation). We then generalize these results to higher $d$ and find, remarkably, that the scaling function close to the edge of the box is universal, i.e., independent of the dimension~$d$., Comment: 31 pages, 14 figures

We consider a one-dimensional run-and-tumble particle, or persistent random walk, in the presence of an absorbing boundary located at the origin. After each tumbling event, which occurs at a constant rate $\gamma$, the (new) velocity of the particle is drawn randomly from a distribution $W(v)$. We study the survival probability $S(x,t)$ of a particle starting from $x \geq 0$ up to time $t$ and obtain an explicit expression for its double Laplace transform (with respect to both $x$ and $t$) for an arbitrary velocity distribution $W(v)$, not necessarily symmetric. This result is obtained as a consequence of Spitzer's formula, which is well known in the theory of random walks and can be viewed as a generalization of the Sparre Andersen theorem. We then apply this general result to the specific case of a two-state particle with velocity $\pm v_0$, the so-called persistent random walk (PRW), and in the presence of a constant drift $\mu$ and obtain an explicit expression for $S(x,t)$, for which we present more detailed results. Depending on the drift $\mu$, we find a rich variety of behaviours for $S(x,t)$, leading to three distinct cases: (i) subcritical drift $-v_0\!<\!\mu\!<\! v_0$, (ii) supercritical drift $\mu < -v_0$ and (iii) critical drift $\mu=-v_0$. In these three cases, we obtain exact analytical expressions for the survival probability $S(x,t)$ and establish connections with existing formulae in the mathematics literature. Finally, we discuss some applications of these results to record statistics and to the statistics of last-passage times., Comment: 58 pages, 14 figures

We combine the processes of resetting and first-passage to define \emph{first-passage resetting}, where the resetting of a random walk to a fixed position is triggered by a first-passage event of the walk itself. In an infinite domain, first-passage resetting of isotropic diffusion is non-stationary, with the number of resetting events growing with time as $\sqrt{t}$. We calculate the resulting spatial probability distribution of the particle analytically, and also obtain this distribution by a geometric path decomposition. In a finite interval, we define an optimization problem that is controlled by first-passage resetting; this scenario is motivated by reliability theory. The goal is to operate a system close to its maximum capacity without experiencing too many breakdowns. However, when a breakdown occurs the system is reset to its minimal operating point. We define and optimize an objective function that maximizes the reward (being close to maximum operation) minus a penalty for each breakdown. We also investigate extensions of this basic model to include delay after each reset and to two dimensions. Finally, we study the growth dynamics of a domain in which the domain boundary recedes by a specified amount whenever the diffusing particle reaches the boundary after which a resetting event occurs. We determine the growth rate of the domain for the semi-infinite line and the finite interval and find a wide range of behaviors that depend on how much the recession occurs when the particle hits the boundary., Comment: 31 pages in IOP format, 7 figures. Version 2: various additions and corrections; for publication in JSTAT

We investigate classic diffusion with the added feature that a diffusing particle is reset to its starting point each time the particle reaches a specified threshold. In an infinite domain, this process is non-stationary and its probability distribution exhibits rich features. In a finite domain, we define a non-trivial optimization in which a cost is incurred whenever the particle is reset and a reward is obtained while the particle stays near the reset point. We derive the condition to optimize the net gain in this system, namely, the reward minus the cost., Comment: 4 pages, 3 figures, revtex 4-1 format. Version 1 contains changes in response to referee comments. Version 2: A missing factor of 2 in an inline formula has been corrected

Emotion lexica are commonly used resources to combat data poverty in automatic emotion detection. However, vocabulary coverage issues, differences in construction method and discrepancies in emotion framework and representation result in a heterogeneous landscape of emotion detection resources, calling for a unified approach to utilising them. To combat this, we present an extended emotion lexicon of 30,273 unique entries, which is a result of merging eight existing emotion lexica by means of a multi-view variational autoencoder (VAE). We showed that a VAE is a valid approach for combining lexica with different label spaces into a joint emotion label space with a chosen number of dimensions, and that these dimensions are still interpretable. We tested the utility of the unified VAE lexicon by employing the lexicon values as features in an emotion detection model. We found that the VAE lexicon outperformed individual lexica, but contrary to our expectations, it did not outperform a naive concatenation of lexica, although it did contribute to the naive concatenation when added as an extra lexicon. Furthermore, using lexicon information as additional features on top of state-of-the-art language models usually resulted in a better performance than when no lexicon information was used., Comment: Computer Speech and Language journal, to appear

Constructing a quantum description of crystals from scattering experiments is of paramount importance to explain their macroscopic properties and to evaluate the pertinence of theoretical ab-initio models. While reconstruction methods of the one-electron reduced density matrix have already been proposed, they are usually tied to strong assumptions that limit and may introduce bias in the model. The goal of this paper is to infer a one-electron reduced density matrix (1-RDM) with minimal assumptions. We have found that the mathematical framework of Semidefinite Programming can achieve this goal. Additionally, it conveniently addresses the nontrivial constraints on the 1-RDM which were major hindrances for the existing models. The framework established in this work can be used as a reference to interpret experimental results. This method has been applied to the crystal of dry ice and provides very satisfactory results when compared with periodic ab-initio calculations., Comment: 7 pages, 6 figures

In this paper a decentralized control algorithm for systems composed of $N$ dynamically decoupled agents, coupled by feasibility constraints, is presented. The control problem is divided into $N$ optimal control sub-problems and a communication scheme is proposed to decouple computations. The derivative of the solution of each sub-problem is used to approximate the evolution of the system allowing the algorithm to decentralize and parallelize computations. The effectiveness of the proposed algorithm is shown through simulations in a cooperative driving scenario.

In this paper we present a method to analyse absorption line spectra of a galaxy designed to determine the stellar dynamics and the stellar populations by a direct fit to the spectra. This paper is the first one to report on the application of the method to data. The modelling results in the knowledge of distribution functions that are sums of basis functions. The practical implementation of the method is discussed and a new type of basis functions is introduced. With this method, a dynamical model for NGC 3258 is constructed. This galaxy can be successfully modelled with a potential containing 30% dark matter within 1r_e with a mass of 1.6x10^11 M_o. The total mass within 2r_e is estimated as 5x10^11 M_o, containing 63% dark matter. The model is isotropic in the centre, is radially anisotropic between 0.2 and 2 kpc (0.88 r_e) and becomes tangentially anisotropic further on. The photometry reveals the presence of a dust disk near the centre.

In this paper we adopt a method to analyse absorption line spectra from elliptical galaxies that determines the dynamics of the galaxy and at the same time offers a way to study the stellar populations in that galaxy by a direct fit to the spectra. The result of the modelling is a distribution function for each stellar population that is considered in the fit. The method is described in detail in an accompanying paper. This paper reports on a dynamical stellar population study in NGC3258, based on spectra in two different wavelength regions, the near-IR Ca II triplet around 8600A and the Ca H and K lines around 3900A. These absorption lines have discriminating power toward various stellar types. The dynamical model shows an increase in dwarfs (represented by G2V stars) toward the centre. Most of the rotation in the model is delivered by the giants (represented by M1III stars). Moreover, the different models that were considered indicate that establishing a potential for a galaxy is dependent on the wavelength range used for the modelling., Comment: 16 pages, 18 figures, accepted for publication by MNRAS

Current error estimates on kinematic parameters are based on the assumption that the data points in the spectra follow a Poisson distribution. For realistic data that have undergone several steps in a reduction process, this is generally not the case. Neither is the noise distribution independent in adjacent pixels. Hence, the error estimates on the derived kinematic parameters will (in most cases) be smaller than the real errors. In this paper we propose a method that makes a diagnosis of the characteristics of the observed noise The method also offers the possibility to calculate more realistic error estimates on kinematic parameters. The method was tested on spectroscopic observations of NGC3258. In this particular case, the realistic errors are almost a factor of 2 larger than the errors based on least squares statistics., Comment: 11 pages, 11 figures, accepted for publication by MNRAS

We use a set of Stackel potentials to obtain a local approximation for an effective third integral in axisymmetric systems. We present a study on the feasibility and effectiveness of this approach. We have applied it to three trial potentials of various flattenings, corresponding to nearly ellipsoidal, disky and boxy density isophotes. In all three cases, a good fit to the potential requires only a small set of Stackel potentials, and the associated Stackel third integral provides a very satisfactory, yet analytically simple, approximation to the trial potentials effective third integral., Comment: 11 pages, 12 figures, accepted for publication in MNRAS

Using the photometry and the kinematic data, part of which acquired with the NTT, we constructed a 2-integral and a 3-integral distribution function for NGC 4697. We detected a nuclear dust lane at $3.4''$ or 0.4 kpc from the centre. A comparison of the Lucy-deprojection method and the multi-gaussian expansion method showed that the latter offers a number of advantages, among which a qualitatively better fit and better controllable convergence. Futhermore, we show that a 2-integral model cannot be excluded for as far as the Satoh parameter $k$ is allowed to vary. The 3-integral model is based on a St\"ackel third integral. The original potential is retained where appropriate, in order to minimise errors due to the St\"ackel approximation. On the whole, the 3-integral model produces more satisfying results than the 2-integral model, and is better in handling some peculiarities in the kinematic data. For both models we show detailed contour plots for all deprojected moments up to order 2. The mass-to-light ratio appears to be a fairly well determined $(M/L)_B=4.8h_{50}$, and is much better constrained than would be the case for spherical anisotropic models with comparable data., Comment: 18 pages, uuencoded and Z-compressed postscript file without the 25 figures, complete version in Z-compressed postscript format available at ftp://naos.rug.ac.be/pub/N4697paper.ps.Z, accepted for publication in A&A