We developed a reliable probabilistic solar flare forecasting model using a deep neural network, named Deep Flare Net-Reliable (DeFN-R). The model can predict the maximum classes of flares that occur in the following 24 h after observing images, along with the event occurrence probability. We detected active regions from 3x10^5 solar images taken during 2010-2015 by Solar Dynamic Observatory and extracted 79 features for each region, which we annotated with flare occurrence labels of X-, M-, and C-classes. The extracted features are the same as used by Nishizuka et al. (2018); for example, line-of-sight/vector magnetograms in the photosphere, brightening in the corona, and the X-ray emissivity 1 and 2 h before an image. We adopted a chronological split of the database into two for training and testing in an operational setting: the dataset in 2010-2014 for training and the one in 2015 for testing. DeFN-R is composed of multilayer perceptrons formed by batch normalizations and skip connections. By tuning optimization methods, DeFN-R was trained to optimize the Brier skill score (BSS). As a result, we achieved BSS = 0.41 for >=C-class flare predictions and 0.30 for >=M-class flare predictions by improving the reliability diagram while keeping the relative operating characteristic curve almost the same. Note that DeFN is optimized for deterministic prediction, which is determined with a normalized threshold of 50%. On the other hand, DeFN-R is optimized for a probability forecast based on the observation event rate, whose probability threshold can be selected according to users' purposes.

The members of asteroid families have various shapes. We investigate the origin of their shapes by high-resolution impact simulations for catastrophic disruptions using a Smoothed Particle Hydrodynamics code. Collisional remnants produced through our simulations of the catastrophic disruptions mainly have spherical or bilobed shapes. However, no flat remnants with the ratio of minor to major axis lengths c/a < 0.5 are formed. The results of the simulations provide various shapes of asteroids and explain most of the shapes in asteroid families that are supposed to be produced through catastrophic disruptions. However, the present simulations do not explain significantly flat asteroids. We suggest that these flat asteroids may be interlopers or formed through low-velocity collisions between member asteroids.

Light curve observations of a recently discovered interstellar object 1I/`Oumuamua suggest that this object has an extremely elongated shape with the axis ratio 0.3 or smaller. Planetesimal collisions can produce irregular shapes including elongated shapes. In this paper, we suggest that the extremely elongated shape of 1I/`Oumuamua may be the result of such an impact. To find detailed impact conditions to form the extremely elongated objects, we conduct numerical simulations of planetesimal collisions using Smoothed Particle Hydrodynamics method for elastic dynamics with self-gravity and interparticle friction. Impacts into strengthless target planetesimals with radius 50 m are conducted with various ratios of impactor mass to target mass q, friction angles phi_d, impact velocities v_imp, and impact angles theta_imp. We find that impacts with q \geq 0.5, phi_d \geq 40 degrees, v_imp \leq 40 degrees, and theta_imp \leq 30 degrees produce remnants with the ratio of intermediate to major axis length less than 0.3. This impact condition suggests that the parent protoplanetary disk in the planetesimal collision stage was weakly turbulent (alpha < 10^{-4} for the inner disk) and composed of planetesimals smaller than ~ 7 km to ensure small impact velocity.

This paper is a sequel to Takahashi et al. (2016), in which the authors investigated the influences of fluctuations in pre-shock accreting matter on the linear stability of the standing accretion shock in core-collapse supernovae (CCSNe). If one understands that this was concerning the effect of the outer boundary condition for the post-shock accretion flows, the present paper should be regarded as an investigation on the possible influences of the inner boundary conditions. More specifically, we impose a time-dependent, oscillating condition at the inner boundary, having in mind the injection of acoustic power by an oscillating proto-neutron star. We also consider possible correlations between the inner and outer boundary conditions as invoked in the argument for Lepton-number Emission Self-sustained Asymmetry, or LESA. In this paper, we conduct the linear stability analysis of the standing accretion shock commonly encountered in CCSNe based on Laplace transform. We find that the acoustic power enhances the standing accretion shock instability, or SASI, especially when the luminosity is low. On the other hand, the correlation between the fluctuations of neutrino luminosity at the neutrino sphere has little influences on the instability, changing the amplitudes of eigenmodes only slightly. We further investigate steady solution of perturbation equations, being motivated by LESA, and conclude that not the difference but the sum of the fluxes of electron-type neutrinos and anti-neutrino is the key ingredient to production of the self-sustained steady perturbed configuration.

More than a half of asteroids in the main belt have irregular shapes with the ratios of the minor to major axis lengths less than 0.6. One of the mechanisms to create such shapes is collisions between asteroids. The relationship between shapes of collisional outcomes and impact conditions such as impact velocities may provide information on the collisional environments and its evolutionary stages when those asteroids are created. In this study, we perform numerical simulations of collisional destruction of asteroids with radii 50 km and subsequent gravitational reaccumulation using Smoothed Particle Hydrodynamics for elastic dynamics with self-gravity, a model of fracture of rock, and a model of friction of completely damaged rock. We systematically vary the impact velocity from 50 m/s to 400 m/s and the impact angle from 5 degrees to 45 degrees. We investigate shapes of the largest remnants resulting from collisional simulations. As a result, various shapes (bilobed, spherical, flat, elongated, and hemispherical shapes) are formed through equal-mass and low-velocity (50 - 400 m/s) impacts. We clarify a range of the impact angle and velocity to form each shape. Our results indicate that irregular shapes, especially flat shapes, of asteroids with diameters larger than 80 km are likely to be formed through similar-mass and low-velocity impacts, which are likely to occur in the primordial environment prior to the formation of Jupiter.

We developed a solar flare prediction model using a deep neural network (DNN), named Deep Flare Net (DeFN). The model can calculate the probability of flares occurring in the following 24 h in each active region, which is used to determine the most likely maximum classes of flares via a binary classification (e.g., >=M class versus <M class or >=C class versus <C class). From 3x10^5 observation images taken during 2010-2015 by Solar Dynamic Observatory, we automatically detected sunspots and calculated 79 features for each region, to which flare occurrence labels of X-, M-, and C-class were attached. We adopted the features used in Nishizuka et al. (2017) and added some features for operational prediction: coronal hot brightening at 131 A (T>=10^7 K) and the X-ray and 131 A intensity data 1 and 2 h before an image. For operational evaluation, we divided the database into two for training and testing: the dataset in 2010-2014 for training and the one in 2015 for testing. The DeFN model consists of deep multilayer neural networks, formed by adapting skip connections and batch normalizations. To statistically predict flares, the DeFN model was trained to optimize the skill score, i.e., the true skill statistic (TSS). As a result, we succeeded in predicting flares with TSS=0.80 for >=M-class flares and TSS=0.63 for >=C-class flares. Note that in usual DNN models, the prediction process is a black box. However, in the DeFN model, the features are manually selected, and it is possible to analyze which features are effective for prediction after evaluation.

Godunov Smoothed Particle Hydrodynamics (Godunov SPH) method is a computational fluid dynamics method that utilizes a Riemann solver and achieves the second-order accuracy in space. In this paper, we extend the Godunov SPH method to elastic dynamics by incorporating deviatoric stress tensor that represents the stress for shear deformation or anisotropic compression. Analogously to the formulation of the original Godunov SPH method, we formulate the equation of motion, the equation of energy, and the time evolution equation of deviatoric stress tensor so that the resulting discretized system achieves the second-order accuracy in space. The standard SPH method tends to suffer from the tensile instability that results in unphysical clustering of particles especially in tension-dominated region. We find that the tensile instability can be suppressed by selecting appropriate interpolation for density distribution in the equation of motion for the Godunov SPH method even in the case of elastic dynamics. Several test calculations for elastic dynamics are performed, and the accuracy and versatility of the present method are shown.

We developed a flare prediction model using machine learning, which is optimized to predict the maximum class of flares occurring in the following 24 h. Machine learning is used to devise algorithms that can learn from and make decisions on a huge amount of data. We used solar observation data during the period 2010-2015, such as vector magnetogram, ultraviolet (UV) emission, and soft X-ray emission taken by the Solar Dynamics Observatory and the Geostationary Operational Environmental Satellite. We detected active regions from the full-disk magnetogram, from which 60 features were extracted with their time differentials, including magnetic neutral lines, the current helicity, the UV brightening, and the flare history. After standardizing the feature database, we fully shuffled and randomly separated it into two for training and testing. To investigate which algorithm is best for flare prediction, we compared three machine learning algorithms: the support vector machine (SVM), k-nearest neighbors (k-NN), and extremely randomized trees (ERT). The prediction score, the true skill statistic (TSS), was higher than 0.9 with a fully shuffled dataset, which is higher than that for human forecasts. It was found that k-NN has the highest performance among the three algorithms. The ranking of the feature importance showed that the previous flare activity is most effective, followed by the length of magnetic neutral lines, the unsigned magnetic flux, the area of UV brightening, and the time differentials of features over 24 h, all of which are strongly correlated with the flux emergence dynamics in an active region.

The modification of Smoothed Particle Hydrodynamics (SPH) method with Riemann Solver is called Godunov SPH. We further extend the Godunov SPH to the description of a medium with negative pressure. Under certain circumstances, the SPH method shows an unphysical instability that results in particle clustering. This instability is called the tensile instability. The tensile instability occurs in positive pressure regions in a regular fluid if a very large number of neighbor particles are used with certain shapes of kernel functions, and it is significant in negative pressure regions that emerge in stretched elastic bodies. We must suppress the tensile instability in SPH for calculations of elastic bodies. In this study, we develop a new technique to remove the tensile instability by extending the Godunov SPH method and conducting a linear stability analysis of the equation of motion for the extended method. We find that the tensile instability can be suppressed by choosing an appropriate order of interpolation in the equation of motion of the Godunov SPH method. We also derive an analytic solution for a Riemann solver for a simple equation of state of an elastic body, and construct a Godunov SPH method for the equation of state that allows negative pressure.