The Patterns of High-Level Magnetic Activity Occurring on the Surface of V1285 Aql: The OPEA Model of Flares and DFT Models of Stellar Spots

Statistically analyzing Johnson UBVR observations of V1285 Aql during the three observing seasons, both activity level and behavior of the star are discussed in respect to obtained results. We also discuss the out-of-flare variation due to rotational modulation. Eighty-three flares were detected in the U-band observations of season 2006 . First, depending on statistical analyses using the independent samples t-test, the flares were divided into two classes as the fast and the slow flares. According to the results of the test, there is a difference of about 73 s between the flare-equivalent durations of slow and fast flares. The difference should be the difference mentioned in the theoretical models. Second, using the one-phase exponential association function, the distribution of the flare-equivalent durations versus the flare total durations was modeled. Analyzing the model, some parameters such as plateau, half-life values, mean average of the flare-equivalent durations, maximum flare rise, and total duration times are derived. The plateau value, which is an indicator of the saturation level of white-light flares, was derived as 2.421{\pm}0.058 s in this model, while half-life is computed as 201 s. Analyses showed that observed maximum value of flare total duration is 4641 s, while observed maximum flare rise time is 1817 s. According to these results, although computed energies of the flares occurring on the surface of V1285 Aql are generally lower than those of other stars, the length of its flaring loop can be higher than those of more active stars.Comment: 44 pages, 10 figures, 5 tables, 2011PASP..123..659

A detailed study of the rise phase of a long duration X-ray flare in the young star TWA 11B

We analyzed a long duration flare observed in a serendipitous XMM-Newton detection of the M star CD-39 7717B (TWA 11B), member of the young stellar association TW Hya (~ 8 Myr). Only the rise phase (with a duration of ~ 35 ks) and possibly the flare peak were observed. We took advantage of the high count-rate of the X-ray source to carry out a detailed analysis of its spectrum during the whole exposure. After a careful analysis, we interpreted the rise phase as resulting from the ignition of a first group of loops (event A) which triggered a subsequent two-ribbon flare (event B). Event A was analyzed using a single-loop model, while a two-ribbon model was applied for event B. Loop semi-lengths of ~ 4 R* were obtained. Such large structures had been previously observed in very young stellar objects (~ 1 - 4 Myr). This is the first time that they have been inferred in a slightly more evolved star. The fluorescent iron emission line at 6.4 keV was detected during event B. Since TWA 11B seems to have no disk, the most plausible explanation found for its presence in the X-ray spectrum of this star is collisional - or photo- ionization. As far as we are concerned, this is only the third clear detection of Fe photospheric fluorescence in stars other than the Sun.Comment: Accepted for publication in ApJ. 15 pages, 9 figure

Where Are The M Dwarf Disks Older Than 10 Million Years?

We present 11.7-micron observations of nine late-type dwarfs obtained at the Keck I 10-meter telescope in December 2002 and April 2003. Our targets were selected for their youth or apparent IRAS 12-micron excess. For all nine sources, excess infrared emission is not detected. We find that stellar wind drag can dominate the circumstellar grain removal and plausibly explain the dearth of M Dwarf systems older than 10 Myr with currently detected infrared excesses. We predict M dwarfs possess fractional infrared excess on the order of L_{IR}/L_{*}\sim10^{-6} and this may be detectable with future efforts.Comment: 24 pages, 2 figures, accepted to Ap

Bayesian Estimation of Hardness Ratios: Modeling and Computations

A commonly used measure to summarize the nature of a photon spectrum is the so-called Hardness Ratio, which compares the number of counts observed in different passbands. The hardness ratio is especially useful to distinguish between and categorize weak sources as a proxy for detailed spectral fitting. However, in this regime classical methods of error propagation fail, and the estimates of spectral hardness become unreliable. Here we develop a rigorous statistical treatment of hardness ratios that properly deals with detected photons as independent Poisson random variables and correctly deals with the non-Gaussian nature of the error propagation. The method is Bayesian in nature, and thus can be generalized to carry out a multitude of source-population--based analyses. We verify our method with simulation studies, and compare it with the classical method. We apply this method to real world examples, such as the identification of candidate quiescent Low-mass X-ray binaries in globular clusters, and tracking the time evolution of a flare on a low-mass star.Comment: 43 pages, 10 figures, 3 tables; submitted to Ap

Performance optimization and load-balancing modeling for superparametrization by 3D LES

In order to eliminate climate uncertainty w.r.t. cloud and convection parametrizations, superpramaterization (SP) [1] has emerged as one of the possible ways forward. We have implemented (regional) superparametrization of the ECMWF weather model OpenIFS [2] by cloud-resolving, three-dimensional large-eddy simulations. This setup, described in [3], contains a two-way coupling between a global meteorological model that resolves large-scale dynamics, with many local instances of the Dutch Atmospheric Large Eddy Simulation (DALES) [4], resolving cloud and boundary layer physics. The model is currently prohibitively expensive to run over climate or even seasonal time scales, and a global SP requires the allocation of millions of cores. In this paper, we study the performance and scaling behavior of the LES models and the coupling code and present our implemented optimizations. We mimic the observed load imbalance with a simple performance model and present strategies to improve hardware utilization in order to assess the feasibility of a world-covering superparametrization. We conclude that (quasi-)dynamical load-balancing can significantly reduce the runtime for such large-scale systems with wide variability in LES time-stepping speeds