Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTopological data analysis and diagnostics of compressible MHD turbulence
66
0
0.0
(
0
)
Authors
Irina Makarenko
Ask ChatGPT about the research
No Arabic abstract
The predictions of mean-field electrodynamics can now be probed using direct numerical simulations of random flows and magnetic fields. When modelling astrophysical MHD, it is important to verify that such simulations are in agreement with observations. One of the main challenges in this area is to identify robust it{quantitative} measures to compare structures found in simulations with those inferred from astrophysical observations. A similar challenge is to compare quantitatively results from different simulations. Topological data analysis offers a range of techniques, including the Betti numbers and persistence diagrams, that can be used to facilitate such a comparison. After describing these tools, we first apply them to synthetic random fields and demonstrate that, when the data are standardized in a straightforward manner, some topological measures are insensitive to either large-scale trends or the resolution of the data. Focusing upon one particular astrophysical example, we apply topological data analysis to HI observations of the turbulent interstellar medium (ISM) in the Milky Way and to recent MHD simulations of the random, strongly compressible ISM. We stress that these topological techniques are generic and could be applied to any complex, multi-dimensional random field.
rate research
Read More
We discuss magnetic pressure and density fluctuations in strongly turbulent isothermal MHD flows in 1+2/3 dimensions. We first consider simple nonlinear MHD waves, which allow us show that the slow and fast modes have different asymptotic dependences of the magnetic pressure B^2 vs. rho. For the slow mode, B^2 ~= c_1-c_2 rho, while for the fast mode, B^2 ~= rho^2. We also perform a perturbative analysis to investigate Alfven wave pressure, recovering previous results that B^2 ~= rho^gamma_e, with gamma_e ~= 2, 3/2 and 1/2 at respectively large, moderate and low M_a. This variety of scalings implies that a single polytropic description of magnetic pressure is not possible in general, since the relation between B^2 and rho depends on which mode dominates the density fluctuation production, which in turn depends on the angle between the magnetic field and the direction of wave propagation, and on the Alfvenic Mach number M_a. Typically, at small M_a, the slow mode dominates, and B is ANTIcorrelated with rho. At large M_a, both modes contribute to density fluctuation production, and the magnetic pressure decorrelates from density, exhibiting a large scatter, which however decreases towards higher densities. In this case, the unsystematic behavior of the magnetic pressure causes the density PDF to generally maintain the lognormal shape corresponding to non-magnetic isothermal turbulence, except when the slow mode dominates, in which case the PDF develops an excess at low densities. Our results are consistent with the low values and apparent lack of correlation between the magnetic field strength and density in surveys of the lower-density molecular gas, and also with the recorrelation apparently seen at higher densities, if M_a is relatively large there.
We discuss the problem of extending data mining approaches to cases in which data points arise in the form of individual graphs. Being able to find the intrinsic low-dimensionality in ensembles of graphs can be useful in a variety of modeling contexts, especially when coarse-graining the detailed graph information is of interest. One of the main challenges in mining graph data is the definition of a suitable pairwise similarity metric in the space of graphs. We explore two practical solutions to solving this problem: one based on finding subgraph densities, and one using spectral information. The approach is illustrated on three test data sets (ensembles of graphs); two of these are obtained from standard graph generating algorithms, while the graphs in the third example are sampled as dynamic snapshots from an evolving network simulation.
It is generally known that counting statistics is not correctly described by a Gaussian approximation. Nevertheless, in neutron scattering, it is common practice to apply this approximation to the counting statistics; also at low counting numbers. We show that the application of this approximation leads to skewed results not only for low-count features, such as background level estimation, but also for its estimation at double-digit count numbers. In effect, this approximation is shown to be imprecise on all levels of count. Instead, a Multinomial approach is introduced as well as a more standard Poisson method, which we compare with the Gaussian case. These two methods originate from a proper analysis of a multi-detector setup and a standard triple axis instrument.We devise a simple mathematical procedure to produce unbiased fits using the Multinomial distribution and demonstrate this method on synthetic and actual inelastic scattering data. We find that the Multinomial method provide almost unbiased results, and in some cases outperforms the Poisson statistics. Although significantly biased, the Gaussian approach is in general more robust in cases where the fitted model is not a true representation of reality. For this reason, a proper data analysis toolbox for low-count neutron scattering should therefore contain more than one model for counting statistics.
LISA is the upcoming space-based Gravitational Wave telescope. LISA Pathfinder, to be launched in the coming years, will prove and verify the detection principle of the fundamental Doppler link of LISA on a flight hardware identical in design to that of LISA. LISA Pathfinder will collect a picture of all noise disturbances possibly affecting LISA, achieving the unprecedented pureness of geodesic motion necessary for the detection of gravitational waves. The first steps of both missions will crucially depend on a very precise calibration of the key system parameters. Moreover, robust parameters estimation is of fundamental importance in the correct assessment of the residual force noise, an essential part of the data processing for LISA. In this paper we present a maximum likelihood parameter estimation technique in time domain being devised for this calibration and show its proficiency on simulated data and validation through Monte Carlo realizations of independent noise runs. We discuss its robustness to non-standard scenarios possibly arising during the real-life mission, as well as its independence to the initial guess and non-gaussianities. Furthermore, we apply the same technique to data produced in mission-like fashion during operational exercises with a realistic simulator provided by ESA.
In high-energy physics, with the search for ever smaller signals in ever larger data sets, it has become essential to extract a maximum of the available information from the data. Multivariate classification methods based on machine learning techniques have become a fundamental ingredient to most analyses. Also the multivariate classifiers themselves have significantly evolved in recent years. Statisticians have found new ways to tune and to combine classifiers to further gain in performance. Integrated into the analysis framework ROOT, TMVA is a toolkit which hosts a large variety of multivariate classification algorithms. Training, testing, performance evaluation and application of all available classifiers is carried out simultaneously via user-friendly interfaces. With version 4, TMVA has been extended to multivariate regression of a real-valued target vector. Regression is invoked through the same user interfaces as classification. TMVA 4 also features more flexible data handling allowing one to arbitrarily form combined MVA methods. A generalised boosting method is the first realisation benefiting from the new framework.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
