Yu-Qing Lou

Modeling solar force-free magnetic fields

A class of nonlinear force-free magnetic fields is presented, described in terms of the solutions to a second-order, nonlinear ordinary differential equation. These magnetic fields are three-dimensional, filling the infinite half-space above a plane where the lines of force are anchored. They model the magnetic fields of the sun over active regions with a striking geometric realism. The total energy and the free energy associated with the electric current are finite and can be calculated directly from the magnetic field at the plane boundary using the virial theorem. In the study of solar magnetic fields with data from vector magnetographs, there is a long-standing interest in devising algorithms to extrapolate for the force-free magnetic field in a given domain from prescribed field values at the boundary. The closed-form magnetic fields of this paper open up an opportunity for testing the reliability and accuracy of algorithms that claim the capability of performing this extrapolation. The extrapolation procedure as an ill-posed mathematical problem is discussed. 22 refs.

OPTICAL AND NEAR-INFRARED OBSERVATIONS OF THE HIGHLY REDDENED, RAPIDLY EXPANDING TYPE Ia SUPERNOVA SN 2006X IN M100

We present extensive optical (UBVRI), near-infrared (JK) light curves and optical spectroscopy of the Type Ia supernova SN 2006X in the nearby galaxy NGC 4321 (M100). Our observations suggest that either SN 2006X has an intrinsically peculiar color evolution or it is highly reddened [E(B − V)host = 1.42 ± 0.04 mag ] with RV = 1.48 ± 0.06, much lower than the canonical value of 3.1 for the average Galactic dust. SN 2006X also has one of the highest expansion velocities ever published for an SN Ia. Compared with the other SNe Ia we analyzed, SN 2006X has a broader light curve in the U band, a more prominent bump/shoulder feature in the V and R bands, a more pronounced secondary maximum in the I and NIR bands, and a remarkably smaller late-time decline rate in the B band. The B − V color evolution shows an obvious deviation from the Lira-Phillips relation at 1-3 months after maximum brightness. At early times, optical spectra of SN 2006X displayed strong, high-velocity features of both intermediate-mass elements (Si, Ca, and S) and iron peak elements, while at late times they showed a relatively blue continuum, consistent with the blue U − B and B − V colors at similar epochs. A light echo and/or the interaction of the SN ejecta and its circumstellar material may provide a plausible explanation for its late-time photometric and spectroscopic behavior. Using the Cepheid distance of M100, we derive a Hubble constant of 72.8 ± 8.2 km s−1 Mpc−1 (statistical) from the normalized dereddened luminosity of SN 2006X. We briefly discuss whether abnormal dust is a universal signature for all SNe Ia and whether the most rapidly expanding objects form a subclass with distinct photometric and spectroscopic properties.

The All-Particle Spectrum of Primary Cosmic Rays in the Wide Energy Range from 1014 to 1017 eV Observed with the Tibet-III Air-Shower Array

We present an updated all-particle energy spectrum of primary cosmic rays in a wide range from 10 14 to 10 17 eVusing 5:5 ; 10 7 events collected from 2000 November through 2004 October by the Tibet-III air-shower array located 4300 m abovesealevel(anatmosphericdepthof 606gcm � 2 ).Thesizespectrumexhibitsasharpkneeatacorrespondingprimary energy around 4 PeV. This work uses increased statistics and new simulation calculations for the analysis. We discuss our extensive Monte Carlo calculations and the model dependencies involved in thefinal result, assuming interaction models QGSJET01c and SIBYLL2.1, and heavy dominant (HD) and proton dominant (PD) primary composition models. Pure protonandpureironprimarymodelsarealsoexaminedasextremecases.A detector simulationwasalsoperformedtoimproveouraccuracyindeterminingthesizeof theairshowersandtheenergyof theprimaryparticle.Weconfirmedthatthe all-particle energy spectra obtained under various plausible model parameters are not significantly different from each other, which was the expected result given the characteristics of the experiment at high altitude, where the air showers of the primary energy around the knee reach near-maximum development, with their features dominated by electromagnetic components,leadingtoaweakdependenceontheinteractionmodel or theprimarymass.Thisisthehighest statistical and the best systematics-controlled measurement covering the widest energy range around the knee energy region.

Periodicities in Solar Coronal Mass Ejections

ABSTRACT Mid-term quasi-periodicities in solar coronal mass ejections (CMEs) during the mostrecent solar maximum cycle 23 are reported here for the first time using the four-year data (February 5, 1999 to February 10, 2003) of the Large Angle SpectrometricCoronagraph (LASCO) onboard the Solar and Heliospheric Observatory (SOHO). Inparallel, mid-term quasi-periodicities in solar X-ray flares (class >M5.0) from theGeosynchronous Operational Environment Satellites (GOES) and in daily averages ofAp index for geomagnetic disturbances from the World Data Center (WDC) at theInternationalAssociation for Geomagnetism andAeronomy (IAGA)are alsoexaminedfor the same four-year time span. By Fourier power spectral analyses, the CME dataappears to contain significant power peaks at periods of ∼ 358 ± 38, ∼ 272 ± 26,∼ 196±13 days and so forth, while except for the ∼ 259±24-day period, X-ray solarflares of class >∼ M5.0 show the familiar Rieger-type quasi-periods at ∼ 157 ± 11,∼ 122± 5, ∼ 98±3 days and shorter ones until ∼ 34±0.5 days. In the data of dailyaverages of Ap index, the two significant peaks at periods ∼ 273±26 and ∼ 187±12days (the latter is most prominent) could imply that CMEs (periods at ∼ 272±26and∼ 196 ± 13 days) may be proportionally correlated with quasi-periodic geomagneticstorm disturbances; at the speculative level, the ∼ 138 ± 6-day period might implythat X-ray flares of class >∼ M5.0 (period at ∼ 157±11 days) may drive certain typesof geomagnetic disturbances; and the ∼ 28±0.2-day periodicity is most likely causedby recurrent high-speed solar winds at the Earth’s magnetosphere. For the same threedata sets, we further perform Morlet wavelet analysis to derive period-time contoursand identify wavelet power peaks and timescales at the 99 percent confidence level forcomparisons. Several conceptual aspects of possible equatorially trapped Rossby-typewaves at and beneath the solar photosphere are discussed.Key words: oscillations — space weather — Sun: activities — corona — coronalmass ejections — magnetic fields

New self-similar solutions of polytropic gas dynamics

We explore semicomplete self-similar solutions for the polytropic gas dynamics involving self-gravity under spherical symmetry, examine behaviours of the sonic critical curve and present new asymptotic collapse solutions that describe ‘quasi-static’ asymptotic behaviours at small radii and large times. These new ‘quasi-static’ solutions with divergent mass density approaching the core can have self-similar oscillations. Earlier known solutions are summarized. Various semicomplete self-similar solutions involving such novel asymptotic solutions are constructed, either with or without a shock. In contexts of stellar core collapse and supernova explosion, a hydrodynamic model of a rebound shock initiated around the stellar degenerate core of a massive progenitor star is presented. With this dynamic model framework, we attempt to relate progenitor stars and the corresponding remnant compact stars: neutron stars, black holes and white dwarfs.

A novel color parameter as a luminosity calibrator for type Ia supernovae

Type Ia supernovae (SNe Ia) provide us with a unique tool for measuring extragalactic distances and determining cosmological parameters. As a result, the precise and effective calibration for peak luminosities of SNe Ia becomes extremely crucial and thus is critically scrutinized for cosmological explorations. In this Letter, we reveal clear evidence for a tight linear correlation between peak luminosities of SNe Ia and their B - V colors ~12 days after the B maximum denoted by ΔC12. By introducing such a novel color parameter, ΔC12, this empirical correlation allows us to uniformly standardize SNe Ia with decline rates Δm15 in the range of 0.8 < Δm15 < 2.0 and to reduce scatters in estimating their peak luminosities from ~0.5 mag to the levels of 0.18 and 0.12 mag in the V and I bands, respectively. For a sample of SNe Ia with insignificant reddenings of host galaxies [e.g., E(B - V)host 0.06 mag], the scatter drops further to only 0.07 mag (or 3%-4% in distance), which is comparable to observational accuracies and is better than other calibrations for SNe Ia. This would impact observational and theoretical studies of SNe Ia and cosmological scales and parameters.

Multi-TeV Gamma-Ray Observation from the Crab Nebula Using the Tibet-III Air Shower Array Finely Tuned by the Cosmic Ray Moon's Shadow

The Tibet-III air shower array, consisting of 533 scintillation detectors, has been operating successfully at Yangbajing in Tibet, China since 1999. Using the data set collected by this array from 1999 November through 2005 November, we obtained the energy spectrum of γ-rays from the Crab Nebula, expressed by a power law as (dJ/dE) = (2.09 ± 0.32) × 10–12(E/3 TeV)–2.96±0.14 cm–2 s–1 TeV–1 in the energy range of 1.7-40 TeV. This result is consistent with other independent γ-ray observations by imaging air Cherenkov telescopes. In this paper, we carefully checked and tuned the performance of the Tibet-III array using data on the Moons shadow in comparison with a detailed Monte Carlo (MC) simulation. The shadow is shifted to the west of the Moons apparent position as an effect of the geomagnetic field, although the extent of this displacement depends on the primary energy of positively charged cosmic rays. This finding enables us to estimate the systematic error in determining the primary energy from its shower size. This error is estimated to be less than ±12% in our experiment. This energy scale estimation is the first attempt among cosmic ray experiments at ground level. The systematic pointing error is also estimated to be smaller than 0011. The deficit rate and the position of the Moons shadow are shown to be very stable within a statistical error of ±6% year by year. This guarantees the long-term stability of pointlike source observation with the Tibet-III array. These systematic errors are adequately taken into account in our study of the Crab Nebula.

CLASSICAL NOVAE IN ANDROMEDA: LIGHT CURVES FROM THE PALOMAR TRANSIENT FACTORY AND GALEX

We present optical light curves of 29 novae in M31 during the 2009 and 2010 observing seasons of the Palomar Transient Factory (PTF). The dynamic and rapid cadences in PTF monitoring of M31, from one day to ten minutes, provide excellent temporal coverage of nova light curves, enabling us to record the photometric evolution of M31 novae in unprecedented detail. We also detect eight of these novae in the near-ultraviolet (UV) band with the Galaxy Evolution Explorer (GALEX) satellite. Novae M31N 2009-10b and M31N 2010-11a show prominent UV emission peaking a few days prior to their optical maxima, possibly implying aspherical outbursts. Additionally, our blueshifted spectrum of the recent outburst of PT And (M31N 2010-12a) indicates that it is a recurrent nova in M31 and not a dwarf nova in the Milky Way as was previously assumed. Finally, we systematically searched for novae in all confirmed globular clusters (GCs) of M31 and found only M31N 2010-10f associated with Bol 126. The specific nova rate in the M31 GC system is thus about one per year, which is not enhanced relative to the rate outside the GC system.

Spherical isothermal self-similar shock flows

We explore self-similar dynamical processes in a spherical isothermal self-gravitational fluid with an emphasis on shocks and outline astrophysical applications of such shock solutions. The previous similarity shock solutions of Tsai & Hsu and of Shu et al. may be classified into two types: Class I solutions with downstream being free-fall collapses; and Class II solutions with downstream being Larson‐Penston (LP) type solutions. By the analyses of Lou & Shen and Shen & Lou, we further construct similarity shock solutions in the ‘semicomplete space’. These general shock solutions can accommodate and model dynamical processes of radial outflows (wind), inflows (accretion or contraction), subsonic oscillations and free-fall core collapses all with shocks in various settings such as star-forming molecular clouds, ‘champagne flows’ in H II regions around luminous massive OB stars or surrounding distant quasars, dynamical connections between the asymptotic giant branch phase to the protoplanetary nebula phase with a central hot white dwarf as well as accretion shocks around compact objects such as white dwarfs, neutron stars and supermassive black holes. By a systematic exploration, we are able to construct families of infinitely many discrete Class I and Class II solutions matching asymptotically with a static outer envelope of a singular isothermal sphere; the shock solutions of Tsai & Hsu form special subsets. These similarity shocks travel at either subsonic or supersonic constant speeds. We also construct twin shocks as well as an ‘isothermal shock’ separating two fluid regions of two different yet constant temperatures. Ke yw ords: hydrodynamics ‐ stars: AGB and post AGB ‐ stars: formation ‐ H II regions ‐ planetary nebulae: general ‐ quasars: general.

THE PROGENITOR OF SN 2004dj IN A STAR CLUSTER

The progenitor of Type II-P (P p plateau) supernova SN 2004dj is identified with a supergiant in a compact star cluster known as “Sandage’s star 96” (S96) in the nearby spiral galaxy NGC 2403, which was fortuitously imaged as part of the Beijing-Arizona-Taiwan-Connecticut (BATC) Multicolor Sky Survey from 1995 February to 2003 December prior to SN 2004dj. The superior photometry of BATC images for S96, taken with 14 intermediate-band filters covering 3000–10000 , unambiguously establishes the star cluster nature of S96, with u A an age of ∼20 Myr, a reddening of mag, and a total mass of ∼96,000 M,. The compact star E(B V) ∼ 0.35 cluster nature of S96 is also consistent with the lack of light variations in the past decade. The SN progenitor is estimated to have a main-sequence mass of ∼12 M,. The comparison of our intermediate-band data of S96 with the postoutburst photometry obtained as the SN significantly dims with time may hopefully conclusively establish the nature of the progenitor. Subject headings: galaxies: individual (NGC 2403) — galaxies: star clusters — stars: evolution — supergiants — supernovae: general — supernovae: individual (SN 2004dj)