Improvements of the shock arrival times at the Earth model STOA

Prediction of the shocks' arrival times (SATs) at the Earth is very important for space weather forecast. There is a well-known SAT model, STOA, which is widely used in the space weather forecast. However, the shock transit time from STOA model usually has a relative large error compared to the real measurements. In addition, STOA tends to yield too much `yes' prediction, which causes a large number of false alarms. Therefore, in this work, we work on the modification of STOA model. First, we give a new method to calculate the shock transit time by modifying the way to use the solar wind speed in STOA model. Second, we develop new criteria for deciding whether the shock will arrive at the Earth with the help of the sunspot numbers and the angle distances of the flare events. It is shown that our work can improve the SATs prediction significantly, especially the prediction of flare events without shocks arriving at the Earth.Comment: Submitted to JG

Dynamical Properties of a Two-gene Network with Hysteresis

A mathematical model for a two-gene regulatory network is derived and several of their properties analyzed. Due to the presence of mixed continuous/discrete dynamics and hysteresis, we employ a hybrid systems model to capture the dynamics of the system. The proposed model incorporates binary hysteresis with different thresholds capturing the interaction between the genes. We analyze properties of the solutions and asymptotic stability of equilibria in the system as a function of its parameters. Our analysis reveals the presence of limit cycles for a certain range of parameters, behavior that is associated with hysteresis. The set of points defining the limit cycle is characterized and its asymptotic stability properties are studied. Furthermore, the stability property of the limit cycle is robust to small perturbations. Numerical simulations are presented to illustrate the results.Comment: 55 pages, 31 figures.Expanded version of paper in Special Issue on Hybrid Systems and Biology, Elsevier Information and Computation, 201

Model of energy spectrum parameters of ground level enhancement events in solar cycle 23

Mewaldt et al. 2012 fitted the observations of the ground level enhancement (GLE) events during solar cycle 23 to the double power-law equation to obtain the four energy spectra parameters, the normalization parameter $C$, low-energy power-law slope $\gamma_1$, high-energy power-law slope $\gamma_2$, and break energy $E_0$. There are 16 GLEs from which we select $13$ for study by excluding some events with complicated situation. We analyze the four parameters with conditions of the corresponding solar events. According to solar event conditions we divide the GLEs into two groups, one with strong acceleration by interplanetary (IP) shocks and another one without strong acceleration. By fitting the four parameters with solar event conditions we obtain models of the parameters for the two groups of GLEs separately. Therefore, we establish a model of energy spectrum of solar cycle 23 GLEs which may be used in prediction in the future.Comment: 42 pages, 19 figures, 3 table

Modulation of Galactic Cosmic Rays in the Inner Heliosphere over Solar Cycles

The 11-year and 22-year modulation of galactic cosmic rays (GCRs) in the inner heliosphere are studied using a numerical model developed by Qin and Shen in 2017. Based on the numerical solutions of Parker's transport equations, the model incorporates a modified Parker heliospheric magnetic field, a locally static time delayed heliosphere, and a time-dependent diffusion coefficients model in which an analytical expression of the variation of magnetic turbulence magnitude throughout the inner heliosphere is applied. Furthermore, during solar maximum, the solar magnetic polarity is determined randomly with the possibility of $A>0$ decided by the percentage of the north solar polar magnetic field being outward and the south solar polar magnetic field being inward. The computed results are compared with several GCR observations, e.g., IMP 8, SOHO/EPHIN, Ulysses, Voyager 1 \& 2, at various energies and show good agreement. It is shown that our model has successfully reproduced the 11-year and 22-year modulation cycles.Comment: Accepted for publication in The Astrophysical Journa

Neural networks based recognition of 3D freeform surface from 2D sketch

In this paper, the Back Propagation (BP) network and Radial Basis Function (RBF) neural network are employed to recognize and reconstruct 3D freeform surface from 2D freehand sketch. Some tests and comparison experiments have been made to evaluate the performance for the reconstruction of freeform surfaces of both networks using simulation data. The experimental results show that both BP and RBF based freeform surface reconstruction methods are feasible; and the RBF network performed better. The RBF average point error between the reconstructed 3D surface data and the desired 3D surface data is less than 0.05 over all our 75 test sample data

Radiative and Collisional Jet Energy Loss in a Quark-Gluon Plasma

We calculate radiative and collisional energy loss of hard partons traversing the quark-gluon plasma created at RHIC and compare the respective size of these contributions. We employ the AMY formalism for radiative energy loss and include additionally energy loss by elastic collisions. Our treatment of both processes is complete at leading order in the coupling, and accounts for the probabilistic nature of jet energy loss. We find that a solution of the Fokker-Planck equation for the probability density distributions of partons is necessary for a complete calculation of the nuclear modification factor $R_{AA}$ for pion production in heavy ion collisions. It is found that the magnitude of $R_{AA}$ is sensitive to the inclusion of both collisional and radiative energy loss, while the average energy is less affected by the addition of collisional contributions. We present a calculation of $R_{AA}$ for $\pi^0$ at RHIC, combining our energy loss formalism with a relativistic (3+1)-dimensional hydrodynamic description of the thermalized medium.Comment: 4 pages, 4 figures, contributed to Quark Matter 2008, Jaipur, Indi

Optimizing Hartree-Fock orbitals by the density-matrix renormalization group

We have proposed a density-matrix renormalization group (DMRG) scheme to optimize the one-electron basis states of molecules. It improves significantly the accuracy and efficiency of the DMRG in the study of quantum chemistry or other many-fermion system with nonlocal interactions. For a water molecule, we find that the ground state energy obtained by the DMRG with only 61 optimized orbitals already reaches the accuracy of best quantum Monte Carlo calculation with 92 orbitals.Comment: published version, 4 pages, 4 figure