An analytical model probing the internal state of coronal mass ejections based on observations of their expansions and propagations
In this paper, a generic self-similar flux rope model is proposed to probe the internal state of CMEs in order to understand the thermodynamic process and expansion of CMEs in interplanetary space. Using this model, three physical parameters and their variations with heliocentric distance can be inferred based on coronagraph observations of CMEs' propagation and expansion. One is the polytropic index Γ of the CME plasma, and the other two are the average Lorentz force and the thermal pressure force inside CMEs. By applying the model to the 8 October 2007 CME observed by STEREO/SECCHI, we find that (1) the polytropic index of the CME plasma increased from initially 1.24 to more than 1.35 quickly and then slowly decreased to about 1.34; it suggests that there be continuously heat injected/converted into the CME plasma and the value of Γ tends to be , a critical value inferred from the model for a force-free flux rope; (2) the Lorentz force directed inward while the thermal pressure force outward, and both of them decreased rapidly as the CME moved out; the direction of the two forces reveals that the thermal pressure force is the internal driver of the CME expansion, whereas the Lorentz force prevented the CME from expanding. Some limitations of the model and approximations are discussed meanwhile.

title: Statistical study of coronal mass ejection source locations: Understanding CMEs viewed in coronagraphs
abstract: How to properly understand coronal mass ejections (CMEs) viewed in white light coronagraphs is crucial to many relative researches in solar and space physics. The issue is now particularly addressed in this paper through studying the source locations of all the 1078 Large Angle and Spectrometric Coronagraph (LASCO) CMEs listed in Coordinated Data Analysis Workshop (CDAW) CME catalog during 1997–1998 and their correlation with CMEs' apparent parameters. By manually checking LASCO and Extreme Ultraviolet Imaging Telescope (EIT) movies of these CMEs, we find that, except 231 CMEs whose source locations cannot be identified due to poor data, there are 288 CMEs with location identified on the frontside solar disk, 234 CMEs appearing above solar limb, and 325 CMEs without evident eruptive signatures in the field of view of EIT. On the basis of the statistical results of CMEs' source locations, there are four physical issues: (1) the missing rate of CMEs by SOHO LASCO and EIT, (2) the mass of CMEs, (3) the causes of halo CMEs, and (4) the deflections of CMEs in the corona, are exhaustively analyzed. It is found that (1) about 32% frontside CMEs cannot be recognized by SOHO, (2) the brightness of a CME at any heliocentric distance is roughly positively correlated with its speed, and the CME mass derived from the brightness is probably overestimated, (3) both projection effect and violent eruption are the major causes of halo CMEs, and especially for limb halo CMEs the latter is the primary one, and (4) most CMEs deflected toward equator near the solar minimum; these deflections can be classified into three types: the asymmetrical expansion, the nonradial ejection, and the deflected propagation.

title: Full halo coronal mass ejections: Do we need to correct the projection effect in terms of velocity?
abstract: The projection effect is one of the biggest obstacles in learning the real properties of coronal mass ejections (CMEs) and forecasting their geoeffectiveness. To evaluate the projection effect, 86 full halo CMEs (FHCMEs) listed in the Coordinated Data Analysis Workshop CME catalog from 1 March 2007 to 31 May 2012 are investigated. By applying the Graduated Cylindrical Shell model, we obtain the deprojected values of the propagation velocity, direction, and angular width of these FHCMEs and compare them with the projected values measured in the plane-of-sky. Although these CMEs look full halo in the view angle of SOHO, it is found that their propagation directions and angular widths could vary in a large range, implying projection effect is a major reason causing a CME being halo, but not the only one. Furthermore, the comparison of the deprojected and projected velocities reveals that most FHCMEs originating within 45° of the Sun-Earth line with a projected speed slower than 900 km s−1 suffer from large projection effect, while the FHCMEs originating far from the vicinity of solar disk center or moving faster than 900 km s−1have small projection effect. Thus, for the latter class of FHCMEs, it is not necessary to correct the measured velocities.

title: Critical Height for the Destabilization of Solar Prominences: Statistical Results from STEREO Observations
abstract: At which height is a prominence inclined to be unstable, or where is the most probable critical height for the prominence destabilization? This question was statistically studied based on 362 solar limb prominences well recognized by Solar Limb Prominence Catcher and Tracker from 2007 April to the end of 2009. We found that there are about 71% disrupted prominences (DPs), among which about 42% of them did not erupt successfully and about 89% of them experienced a sudden destabilization process. After a comprehensive analysis of the DPs, we

title: Quantitative Analysis of CME Deflections in the Corona
abstract: In this paper, ten CME events viewed by the STEREO twin spacecraft are analyzed to study the deflections of CMEs during their propagation in the corona. Based on the three-dimensional information of the CMEs derived by the graduated cylindrical shell (GCS) model (Thernisien, Howard, and Vourlidas in  2006), it is found that the propagation directions of eight CMEs had changed. By applying the theoretical method proposed by Shen et al. (Solar Phys.269, 389, 2011) to all the CMEs, we found that the deflections are consistent, in strength and direction, with the gradient of the magnetic energy density. There is a positive correlation between the deflection rate and the strength of the magnetic energy density gradient and a weak anti-correlation between the deflection rate and the CME speed. Our results suggest that the deflections of CMEs are mainly controlled by the background magnetic field and can be quantitatively described by the magnetic energy density gradient (MEDG) model.

title: CONTRACTING AND ERUPTING COMPONENTS OF SIGMOIDAL ACTIVE REGIONS
abstract: It has recently been noted that solar eruptions can be associated with the contraction of coronal loops that are not involved in magnetic reconnection processes. In this paper, we investigate five coronal eruptions originating from four sigmoidal active regions, using high-cadence, high-resolution narrowband EUV images obtained by the Solar Dynamic Observatory (SDO). The magnitudes of the flares associated with the eruptions range from GOES class B to class X. Owing to the high-sensitivity and broad temperature coverage of the Atmospheric Imaging Assembly (AIA) on board SDO, we are able to identify both the contracting and erupting components of the eruptions: the former is observed in cold AIA channels as the contracting coronal loops overlying the elbows of the sigmoid, and the latter is preferentially observed in warm/hot AIA channels as an expanding bubble originating from the center of the sigmoid. The initiation of eruption always precedes the contraction, and in the energetically mild events (B- and C-flares), it also precedes the increase in GOES soft X-ray fluxes. In the more energetic events, the eruption is simultaneous with the impulsive phase of the nonthermal hard X-ray emission. These observations confirm that loop contraction is an integrated process in eruptions with partially opened arcades. The consequence of contraction is a new equilibrium with reduced magnetic energy, as the contracting loops never regain their original positions. The contracting process is a direct consequence of flare energy release, as evidenced by the strong correlation of the maximal contracting speed, and strong anti-correlation of the time delay of contraction relative to expansion, with the peak soft X-ray flux. This is also implied by the relationship between contraction and expansion, i.e., their timing and speed.

title: A relativistic mixing-layer model for jets in low-luminosity radio galaxies
abstract: We present an analytical model for jets in Fanaroff & Riley Class I (FR I) radio galaxies, in which an initially laminar, relativistic flow is surrounded by a shear layer. We apply the appropriate conservation laws to constrain the jet parameters, starting the model where the radio emission is observed to brighten abruptly. We assume that the laminar flow fills the jet there and that pressure balance with the surroundings is maintained from that point outwards. Entrainment continuously injects new material into the jet and forms a shear layer, which contains material from both the environment and the laminar core. The shear layer expands rapidly with distance until finally the core disappears, and all of the material is mixed into the shear layer. Beyond this point, the shear layer expands in a cone and decelerates smoothly. We apply our model to the well-observed FR I source 3C 31 and show that there is a self-consistent solution. We derive the jet power, together with the variations of mass flux and entrainment rate with distance from the nucleus. The predicted variation of bulk velocity with distance in the outer parts of the jets is in good agreement with model fits to Very Large Array observations. Our prediction for the shape of the laminar core can be tested with higher-resolution imaging.

title: A new model for the Milky Way bar
abstract: We use Schwarzschild's orbit-superposition technique to construct self-consistent models of the Galactic bar. Using χ2 minimization, we find that the best-fitting Galactic bar model has a pattern speed Ωp = 60 km s−1 kpc−1, disc mass Md = 1.0 × 1011 M⊙ and bar angle θbar = 20° for an adopted bar mass Mbar = 2 × 1010 M⊙. The model can reproduce not only the three-dimensional and projected density distributions but also velocity and velocity dispersion data from the Bulge Radial Velocity Assay (BRAVA) survey. We also predict the proper motions in the range l = [−12°, 12°], b = [−10°, 10°], which appear to be higher than observations in the longitudinal direction. The model is stable within a time-scale of 0.5 Gyr, but appears to deviate from steady state on longer time-scales. Our model can be further tested by future observations such as those from GAIA.

title: Is the energy generation rate of nuclear reactions in hot accretion flows important
abstract: The temperature of hot accretion flows around black holes is sufficiently high for the ignition of nuclear reactions. This is potentially an important nucleosynthesis mechanism in the universe. As the first step in studying this problem, we need to measure physical quantities such as density and temperature of the accretion flow. In usual studies of the hot accretion flow, viscous dissipation is considered to be the only heating mechanism, while the heating caused by nuclear reactions is not considered. In this paper, we investigate whether the energy generation rate of nuclear reaction is important compared to the viscous heating. Our calculation indicates that the former is at most one percent of the latter and thus is not important. The dynamics of accretion flow can be therefore calculated in the usual way, without the need to consider heating due to nuclear reactions.

title: EXPLORATION OF PARALLEL OPTIMIZATION TECHNIQUES FOR ACCELERATOR DESIGN
abstract: Optimization through simulation is one of the most timeconsuming tasks in accelerator design, especially for highdimensional problems. We explored several parallel optimization techniques, including parallel genetic algorithm (PGA), parallel hybrid simplex (PHS), and parallel particle swarm optimization (PPSO), to solve some real-world optimization problems. The serial version of the simplex method in elegant [1] is used as a benchmark to compare with the newly developed parallel optimization algorithms in Pelegant . Our experimental results show almost all of the tested parallel optimization methods converge to an equal or better optimization result compared with the serial simplex method, which indicates the serial version of the simplex method could be easily trapped to a local optimum for some applications. The PPSO method is reliable for global optimization and is well suited for parallel computing. The optimization time for a coupling minimization problem was reduced significantly from more than six hours with serial simplex to about twenty minutes on a BlueGene/P supercomputer at the Argonne Leadership Computer Facility (ALCF). The HPS method requires much less time than other optimization algorithms for an experiment where the optimal solution is very close to the initial point for the optimization. These parallel optimization methods not only make the optimization result more reliable but also provide a feasible approach to on-demand accelerator optimization.

title: Projectile fragmentation reactions of Ar-40 at 57 MeV/nucleon
abstract: The projectile fragmentation reactions of Ar-40 at 57 MeV/nucleon on Be-9 and Ta-181 targets have been studied by the Radioactive Ion Beam Line in Lanzhou (RIBLL) at the Heavy-Ion Research Facility in Lanzhou (HIRFL). The momentum distributions of fragments have been measured and the competition between different mechanisms are observed. The production cross sections have been obtained and compared with the empirical parametrization of fragmentation cross section (EPAX), abrasion-ablation (AA), and heavy-ion phase space exploration (HIPSE) models, and the target dependence of fragment cross sections has also been discussed