We present the first application of adaptive machine learning to the identification of anomalies in a data set of non-periodic astronomical light curves. The method follows an active learning strategy where highly informative objects are selected to be labelled. This new information is subsequently used to improve the machine learning model, allowing its accuracy to evolve with the addition of every new classification. For the case of anomaly detection, the algorithm aims to maximize the number of real anomalies presented to the expert by slightly modifying the decision boundary of a traditional isolation forest in each iteration. As a proof of concept, we apply the Active Anomaly Discovery (AAD) algorithm to light curves from the Open Supernova Catalog and compare its results to those of a static Isolation Forest (IF). For both methods, we visually inspected objects within 2% highest anomaly scores. We show that AAD was able to identify 80% more true anomalies than IF. This result is the first evidence that AAD algorithms can play a central role in the search for new physics in the era of large scale sky surveys.

In the upcoming decade large astronomical surveys will discover millions of transients raising unprecedented data challenges in the process. Only the use of the machine learning algorithms can process such large data volumes. Most of the discovered transients will belong to the known classes of astronomical objects. However, it is expected that some transients will be rare or completely new events of unknown physical nature. The task of finding them can be framed as an anomaly detection problem. In this work, we perform for the first time an automated anomaly detection analysis in the photometric data of the Open Supernova Catalog (OSC), which serves as a proof of concept for the applicability of these methods to future large scale surveys. The analysis consists of the following steps: 1) data selection from the OSC and approximation of the pre-processed data with Gaussian processes, 2) dimensionality reduction, 3) searching for outliers with the use of the isolation forest algorithm, 4) expert analysis of the identified outliers. The pipeline returned 81 candidate anomalies, 27 (33%) of which were confirmed to be from astrophysically peculiar objects. Found anomalies correspond to a selected sample of 1.4% of the initial automatically identified data sample of ~2000 objects. Among the identified outliers we recognised superluminous supernovae, non-classical Type Ia supernovae, unusual Type II supernovae, one active galactic nucleus and one binary microlensing event. We also found that 16 anomalies classified as supernovae in the literature are likely to be quasars or stars. Our proposed pipeline represents an effective strategy to guarantee we shall not overlook exciting new science hidden in the data we fought so hard to acquire. All code and products of this investigation are made publicly available.

The aim of this work is to construct an algorithm for visualizing a polydisperse phase of solid particles (dust) in an inhomogeneous flow of a two-phase gas-dust mixture that would allow us to see, within one plot, the degree of polydispersity of the dust phase and the difference in the spatial distributions of individual fractions of dust particles in the computational domain. The developed technique allows us to reproduce concentrations from one to three fractions of dust particles in each cell in the computational domain. Each of the three fractions of dust particles is mapped to one of the main channels of the RGB palette. The intensity of the color shade is set to be proportional to the relative concentration of dust particles in this fraction. The final image for a polydisperse mixture is obtained by adding images in each of the three color channels. To visualize the degree of polydispersity, I propose depicting the spatial distribution of the entropy of the dust mixture. The definition of the entropy of a mixture is generalized to take into account the states of a mixture with zero number of particles in the mixture. They correspond to dust-free sections of the computational domain (voids). The proposed method for visualizing the polydispersity of a mixture of particles is demonstrated using the example of dynamic numerical modeling of the spatial features of dust structures formed in turbulent gas-dust flows and in flows with shock waves.

The evolution of a supernova remnant in a cloudy medium as a function of the volume filling factor of the clouds is studied in a three-dimensional axially symmetrical model. The model includes the mixing of heavy elements (metals) ejected by the supernova and their contribution to radiative losses. The interaction of the supernova envelope with the cloudy phase of the interstellar medium leads to nonsimultaneous, and on average earlier, onsets of the radiative phase in different parts of the supernova envelope. Growth in the volume filling factor $f$ leads to a decrease in the time for the transition of the envelope to the radiative phase and a decrease in the envelope's mean radius, due to the increased energy losses by the envelope in the cloudy medium. When the development of hydrodynamical instabilities in the supernova envelope is efficient, the thermal energy falls as $E_t\sim t^{-2.3}$, for the propagation of the supernova remnant through either a homogeneous or a cloudy medium. When the volume filling factor is $f\simgt 0.1$, a layer with excess kinetic energy andmomentumforms far behind the global shock front from the supernova, which traps the hot gas of the cavity in the central part of the supernova remnant. Metals ejected by the supernova are also enclosed in the central region of the remnant, where the initial (high) metallicity is essentially preserved. Thus, the interaction of the supernova envelope with the cloudy interstellar medium appreciably changes the dynamics and structure of the distribution of the gas in the remnant. This affects the observational characteristics of the remnant, in particularly, leading to substantial fluctuations of the emission measure of the gas with $T>10^5$~K and the velocity dispersion of the ionized gas.

The paper describes an explicit multi-dimensional numerical scheme for Special Relativistic Two-Fluid Magnetohydrodynamics of electron-positron plasma and a suit of test problems. The scheme utilizes Cartesian grid and the third order WENO interpolation. The time integration is carried out using the third order TVD method of Runge-Kutta type, thus ensuring overall third order accuracy on smooth solutions. The magnetic field is kept near divergence-free by means of the method of generalized Lagrange multiplier. The test simulations, which include linear and non-linear continuous plasma waves, shock waves, strong explosions and the tearing instability, show that the scheme is sufficiently robust and confirm its accuracy.