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

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

Jet energy loss and high pTp_T photon production in hot quark-gluon plasma

Jet-quenching and photon production at high transverse momentum are studied at RHIC energies, together with the correlation between jets and photons. The energy loss of hard partons traversing the hot QGP is evaluated in the AMY formalism, consistently taking into account both induced gluon emission and elastic collisions. The production of high $p_T$ photons in Au+Au collisions is calculated, incorporating a complete set of photon-production channels. Putting all these ingredients together with a (3+1)-dimensional ideal relativistic hydrodynamical description of the thermal medium, we achieve a good description of the current experimental data. Our results illustrate that the interaction between hard jets and the soft medium is important for a complete understanding of jet quenching, photon production, and photon-hadron correlations in relativistic nuclear collisions.Comment: 4 pages, 4 figures - To appear in the conference proceedings for Quark Matter 2009, March 30 - April 4, Knoxville, Tennesse

Modulation of galactic cosmic rays during the unusual solar minimum between cycles 23 and 24

During the recent solar minimum between cycles 23 and 24 (solar minimum $P_{23/24}$) the intensity of Galactic Cosmic Rays (GCRs) measured at the Earth was the highest ever recorded since space age. It is the purpose of this paper to resolve the most plausible mechanism for this unusually high intensity. A GCR transport model in three-dimensional heliosphere based on a simulation of Markov stochastic process is used to find the relation of cosmic ray modulation to various transport parameters, including solar wind (SW) speed, distance of heliospheric boundary, magnitude of interplanetary magnetic field (IMF) at the Earth, tilt angle of heliospheric current sheet (HCS), values of parallel and perpendicular diffusion coefficients. We calculate GCR proton energy spectra at the Earth for the last three solar minima $P_{21/22}$, $P_{22/23}$, and $P_{23/24}$, with the transport parameters obtained from observations. Besides weak IMF magnitude and slow SW speed, we find that a possible low magnetic turbulence, which increases the parallel diffusion and reduces the perpendicular diffusion in the polar direction, might be an additional possible mechanism for the high GCR intensity in the solar minimum $P_{23/24}$.Comment: Accepted for publication in JGR space physic

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

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