Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

We present a detailed method to simulating sensor distortions using a photon and electron Monte Carlo method. We use three dimensional electrostatic simulations to parameterize the perturbed electric field profile for non-ideal sensor details. We follow the conversion of simulated photons, and the subsequent response of the converted electrons to the electric field pattern. These non-ideal sensor details can be implemented efficiently in a Monte Carlo approach. We demonstrate that the non-ideal sensor distortions have a variety of observable consequence including the modification of the astrometric pattern, the distortion of the electron diffusion size and shape, and the distortion of flats. We show analytic validation of the diffusion physics, reproduce two kinds of edge distortion, and show qualitative validation of field-free regions, lithography errors, and fringing. We also demonstrate that there are two related effects of doping variation having different observable consequences. We show that field distortions from accumulated electrons lead to intensity-dependent point-spread-functions and the sub-linear variance in flats. The method is implemented in the Photon Simulator (PhoSim) and the code is publically available.

We extract 18 candidate short gamma-ray bursts (SGRBs) with precursors from 660 SGRBs observed by {\em Fermi} and {\em Swift} satellites, and carry out a comprehensive analysis on their temporal and spectral features. We obtain the following results: (1) For a large fraction of candidates, the main burst durations are longer than their precursor durations, comparable to their quiescent times from the end of precursors to the beginning of their main bursts. (2) The average flux of precursors tends to increase as their main bursts brighten. (3) As seen from the distributions of hardness ratio and spectral fitting, the precursors are slightly spectrally softer with respect to the main bursts. Moreover, a significant portion of precursors and all main bursts favor a non-thermal spectrum. (4) The precursors might be a probe of the progenitor properties of SGRBs such as the magnetic field strength and the crustal equation of state if they arise from some processes before mergers of binary compact objects rather than post-merger processes.

In this paper, we investigate the energy partition of two homologous M1.1 circular-ribbon flares (CRFs) in active region (AR) 12434. They were observed by \textit{SDO}, \textit{GOES}, and \textit{RHESSI} on 2015 October 15 and 16, respectively. The peak thermal energy, nonthermal energy of flare-accelerated electrons, total radiative loss of hot plasma, and radiant energies in 1$-$8 {\AA} and 1$-$70 {\AA} of the flares are calculated. The two flares have similar energetics. The peak thermal energies are (1.94$\pm$0.13)$\times$10$^{30}$ erg. The nonthermal energies in flare-accelerated electrons are (3.9$\pm$0.7)$\times$10$^{30}$ erg. The radiative outputs of the flare loops in 1$-$70 {\AA}, which are $\sim$200 times greater than the outputs in 1$-$8 {\AA}, account for $\sim$62.5\% of the peak thermal energies. The radiative losses of SXR-emitting plasma are one order of magnitude lower than the peak thermal energies. Therefore, the total heating requirements of flare loops including radiative loss are (2.1$\pm$0.1)$\times$10$^{30}$ erg, which could sufficiently be supplied by nonthermal electrons.

We develop a comprehensive approach to simulate the deformation of mirrors and lenses due to thermal and mechanical stresses that couples efficiently to photon-based optics simulations. This expands upon previous work where we demonstrated a comprehensive ab initio approach to simulate astronomical images using a photon Monte Carlo method. We apply elasticity theory and estimate thermal effects by adapting a three-dimensional numerical method. We also consider the effect of active optics control systems and active cooling systems in further correcting distortions in the optics. We validate the approach by showing convergence to analytic estimates, and then apply the methodology to the WIYN 3.5m telescope primary mirror. We demonstrate that changes in the soak temperature result in second order point spread function (PSF) defocusing, the gravitational sag and positioning errors result in highly structured PSF distortions, and large-scale thermal gradients result in an elliptical PSF distortion patterns. All three aspects of the environment are larger than the intrinsic optical aberrations of the design, and further exploration with a variety of telescopes should lead to detailed PSF size and shape, astrometric distortion, and field variation predictions. The simulation capabilities developed in this work is publicly available with the Photon Simulation (PhoSim) package.

Two-episode emission components separated by quiescent gaps in the prompt emission of gamma-ray bursts (GRBs) have been observed in the {\em Swift} era, but there is a lack of spectral information due to the narrow energy band of the {\em Swift}/Burst Alert Telescope. In this paper, a systematic analysis of the spectral and temporal properties of the prompt emission of 101 {\em Fermi}/Gamma-ray Burst Monitor detected long GRBs show the existence of two-episode emission components in the light curves, with quiescent times of up to hundreds of seconds. We focus on investigating the differences of those two emission episodes. We find that the light curves of the two emission components exhibit different behavior, e.g., a soft emission component that either precedes or follows the main prompt emission or that the intensity of the two emission episodes are comparable with each other. No statistically significant correlation in the duration of the two emission episodes can be claimed. We define a new parameter $\varepsilon$ as the ratio of the peak flux of the first and second emission episodes and find that a higher $\varepsilon$ corresponds to a larger fluence. The preferred spectral model in our analysis is a cutoff power-law model for most GRBs. The distribution of $E_p$ for episodes I and II range from tens of keV to 1000 keV with a lognormal fit and there are no significant differences between them. Moreover, we do not find significant relationships between $\varepsilon$ and $E_p$ for the two emission episodes. Those results suggest that these two-episode emission components likely share the same physical origin.

Motivating by the discovery of association between GW 170817 and sGRB 170817A, we present a comprehensive analysis for sGRBs observed with Fermi/GBM in 9 operation years and study the properties of sGRB 170817A -like events. We derive a catalog of 275 typical sGRBs and 48 sGRB 170817A-like weak events from the GBM data of 2217 GRBs. We visibly identify two patterns of their light curve, single episode (Pattern I, 61\% of the SGRBs) and multiple episodes (Pattern II, 39\% of the SGRBs). Their duration distribution shows a tentative bimodal feature. Their spectra can be fitted with a cutoff power-law model, except for 4 sGRBs, and the spectral indices normally distribute at $\Gamma=0.69\pm 0.40$. Their $E_p$ values show a tentative bimodal distribution with peaks at 145 keV and 500 keV. No correlation among $T_{90}$, $E_p$, and $\Gamma$ is found. GRB 170817A is a soft, weak sGRB with $ E_{p}=124\pm 106$ keV, $L_{\rm iso}=(5.67\pm4.65)\times10^{46}\rm ~erg~s^{-1}$, and $E_{\rm iso}=(3.23\pm2.65)\times10^{46}\rm ~erg$. It follows the $E_{\rm iso}-E_{\rm p}$ relation of typical short GRBs. Its lightcurve is of Pattern II. Two lightcurve patterns, together with the potential two components in the $E_{\rm p}$ and $T_{90}$ distributions, we suspect that the current sample may include two distinct types of sGRBs from different progenitors. sGRB 170817A-like events may be from NS-NS mergers and those sGRBs with a Pattern I lightcurve may be from another distinct type of compact binary.

We report observations of a two-stage coronal dimming in an eruptive event of a two-ribbon flare and a fast coronal mass ejection (CME). Weak gradual dimming persists for more than half an hour before the onset of the two-ribbon flare and the fast rise of the CME. It is followed by abrupt rapid dimming. The two-stage dimming occurs in a pair of conjugate dimming regions adjacent to the two flare ribbons, and the flare onset marks the transition between the two stages of dimming. At the onset of the two-ribbon flare, transient brightenings are also observed inside the dimming regions, before rapid dimming occurs at the same places. These observations suggest that the CME structure, most probably anchored at the twin dimming regions, undergoes a slow rise before the flare onset, and its kinematic evolution has significantly changed at the onset of flare reconnection. We explore diagnostics of the CME evolution in the early phase with analysis of the gradual dimming signatures prior to the CME eruption.

The key to photoelectric X-ray polarimetry is the determination of the emission direction of photoelectrons. Because of the low mass of an electron, the ionisation trajectory is not straight and the useful information needed for polarimetry is stored mostly in the initial part of the track where less energy is deposited. We present a new algorithm, based on the shortest path problem in graph theory, to reconstruct the 2D electron track from the measured image that is blurred due to transversal diffusion along drift and multiplication in the gas chamber. Compared with previous methods based on moment analysis, this algorithm allows us to identify the photoelectric interaction point more accurately and precisely for complicated tracks resulting from high energy photons or low pressure chambers. This leads to a better position resolution and a higher degree of modulation toward high energy X-rays. The new algorithm is justified using simulations and measurements with the gas pixel detector (GPD), and it should also work for other polarimetric techniques such as a time projection chamber (TPC). As the improvement is restricted in the high energy band, this new algorithm shows limited improvement for the sensitivity of GPD polarimeters, but it may have a larger potential for low-pressure TPC polarimeters.

The shear measurement is a crucial task in the current and the future weak lensing survey projects. And the reconstruction of the point spread function(PSF) is one of the essential steps. In this work, we present three different methods, including Gaussianlets, Moffatlets and EMPCA to quantify their efficiency on PSF reconstruction using four sets of simulated LSST star images. Gaussianlets and Moffatlets are two different sets of basis functions whose profiles are based on Gaussian and Moffat functions respectively. Expectation Maximization(EM) PCA is a statistical method performing iterative procedure to find principal components of an ensemble of star images. Our tests show that: 1) Moffatlets always perform better than Gaussianlets. 2) EMPCA is more compact and flexible, but the noise existing in the Principal Components (PCs) will contaminate the size and ellipticity of PSF while Moffatlets keeps them very well.

Coronal mass ejections (CMEs) are often accompanied by coronal dimming evident in extreme ultraviolet (EUV) and soft X-ray observations. The locations of dimming are sometimes considered to map footpoints of the erupting flux rope. As the emitting material expands in the corona, the decreased plasma density leads to reduced emission observed in spectral and irradiance measurements. Therefore, signatures of dimming may reflect properties of CMEs in the early phase of its eruption. In this study, we analyze the event of flare, CME, and coronal dimming on December 26, 2011. We use the data from the Atmospheric Imaging Assembly (AIA) on Solar Dynamics Observatories (SDO) for disk observations of the dimming, and analyze images taken by EUVI, COR1, and COR2 onboard the Solar Terrestrial Relations Observatories to obtain the height and velocity of the associated CMEs observed at the limb. We also measure magnetic reconnection rate from flare observations. Dimming occurs in a few locations next to the flare ribbons, and it is observed in multiple EUV passbands. Rapid dimming starts after the onset of fast reconnection and CME acceleration, and its evolution well tracks the CME height and flare reconnection. The spatial distribution of dimming exhibits cores of deep dimming with a rapid growth, and their light curves are approximately linearly scaled with the CME height profile. From the dimming analysis, we infer the process of the CME expansion, and estimate properties of the CME.

Chromospheric evaporation refers to dynamic mass motions in flare loops as a result of rapid energy deposition in the chromosphere. These have been observed as blueshifts in X-ray and extreme-ultraviolet (EUV) spectral lines corresponding to upward motions at a few tens to a few hundreds of km/s. Past spectroscopic observations have also revealed a dominant stationary component, in addition to the blueshifted component, in emission lines formed at high temperatures (~10 MK). This is contradictory to evaporation models predicting predominant blueshifts in hot lines. The recently launched Interface Region Imaging Spectrograph (IRIS) provides high resolution imaging and spectroscopic observations that focus on the chromosphere and transition region in the UV passband. Using the new IRIS observations, combined with coordinated observations from the EUV Imaging Spectrometer, we study the chromospheric evaporation process from the upper chromosphere to corona during an X1.0 flare on 2014 March 29. We find evident evaporation signatures, characterized by Doppler shifts and line broadening, at two flare ribbons separating from each other, suggesting that chromospheric evaporation takes place in successively formed flaring loops throughout the flare. More importantly, we detect dominant blueshifts in the high temperature Fe XXI line (~10 MK), in agreement with theoretical predictions. We also find that, in this flare, gentle evaporation occurs at some locations in the rise phase of the flare, while explosive evaporation is detected at some other locations near the peak of the flare. There is a conversion from gentle to explosive evaporation as the flare evolves.

The China Dark Matter Experiment reports results on light WIMP dark matter searches at the China Jinping Underground Laboratory with a germanium detector array with a total mass of 20 g. The physics threshold achieved is 177 eVee ("ee" represents electron equivalent energy) at 50% signal efficiency. With 0.784 kg-days of data, exclusion region on spin-independent coupling with the nucleon is derived, improving over our earlier bounds at WIMP mass less than 4.6 GeV.

We report the timing and spatial resolution from the Muon Telescope Detector (MTD) installed in the STAR experiment at RHIC. Cosmic ray muons traversing the STAR detector have an average transverse momentum of 6 GeV/c. Due to their very small multiple scattering, these cosmic muons provide an ideal tool to calibrate the detectors and measure their timing and spatial resolution. The values obtained were ~100 ps and ~1-2 cm, respectively. These values are comparable to those obtained from cosmic-ray bench tests and test beams.

The China Dark matter Experiment collaboration reports the first experimental limit on WIMP dark matter from 14.6 kg-day of data taken with a 994 g p-type point-contact germanium detector at the China Jinping underground Laboratory where the rock overburden is more than 2400 m. The energy threshold achieved was 400 eVee. According to the 14.6 kg-day live data, we placed the limit of N= 1.75 * 10^{-40} cm^{2} at 90% confidence level on the spin-independent cross-section at WIMP mass of 7 GeV before differentiating bulk signals from the surface backgrounds.

China JinPing underground Laboratory (CJPL) is the deepest underground laboratory presently running in the world. In such a deep underground laboratory, the cosmic ray flux is a very important and necessary parameter for rare event experiments. A plastic scintillator telescope system has been set up to measure the cosmic ray flux. The performance of the telescope system has been studied using the cosmic ray on the ground laboratory near CJPL. Based on the underground experimental data taken from November 2010 to December 2011 in CJPL, which has effective live time of 171 days, the cosmic ray muon flux in CJPL is measured to be (2.0+-0.4)*10^(-10)/(cm^2)/(s). The ultra-low cosmic ray background guarantees CJPL's ideal environment for dark matter experiment.

The CDEX Collaboration has been established for direct detection of light dark matter particles, using ultra-low energy threshold p-type point-contact germanium detectors, in China JinPing underground Laboratory (CJPL). The first 1 kg point-contact germanium detector with a sub-keV energy threshold has been tested in a passive shielding system located in CJPL. The outputs from both the point-contact p+ electrode and the outside n+ electrode make it possible to scan the lower energy range of less than 1 keV and at the same time to detect the higher energy range up to 3 MeV. The outputs from both p+ and n+ electrode may also provide a more powerful method for signal discrimination for dark matter experiment. Some key parameters, including energy resolution, dead time, decay times of internal X-rays, and system stability, have been tested and measured. The results show that the 1 kg point-contact germanium detector, together with its shielding system and electronics, can run smoothly with good performances. This detector system will be deployed for dark matter search experiments.

Weakly Interacting Massive Particles (WIMPs) are the candidates of dark matter in our universe. Up to now any direct interaction of WIMP with nuclei has not been observed yet. The exclusion limits of the spin-independent cross section of WIMP-nucleon which have been experimentally obtained is about 10^{-7}pb at high mass region and only 10^{-5}pb} at low mass region. China Jin-Ping underground laboratory CJPL is the deepest underground lab in the world and provides a very promising environment for direct observation of dark matter. The China Dark Matter Experiment (CDEX) experiment is going to directly detect the WIMP flux with high sensitivity in the low mass region. Both CJPL and CDEX have achieved a remarkable progress in recent two years. The CDEX employs a point-contact germanium semi-conductor detector PCGe whose detection threshold is less than 300 eV. We report the measurement results of Muon flux, monitoring of radioactivity and Radon concentration carried out in CJPL, as well describe the structure and performance of the 1 kg PCGe detector CDEX-1 and 10kg detector array CDEX-10 including the detectors, electronics, shielding and cooling systems. Finally we discuss the physics goals of the CDEX-1, CDEX-10 and the future CDEX-1T detectors.

Context. Hard X-ray (HXR) spikes refer to fine time structures on timescales of seconds to milliseconds in high-energy HXR emission profiles during solar flare eruptions. Aims. We present a preliminary statistical investigation of temporal and spectral properties of HXR spikes. Methods. Using a three-sigma spike selection rule, we detected 184 spikes in 94 out of 322 flares with significant counts at given photon energies, which were detected from demodulated HXR light curves obtained by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). About one fifth of these spikes are also detected at photon energies higher than 100 keV. Results. The statistical properties of the spikes are as follows. (1) HXR spikes are produced in both impulsive flares and long-duration flares with nearly the same occurrence rates. Ninety percent of the spikes occur during the rise phase of the flares, and about 70% occur around the peak times of the flares. (2) The time durations of the spikes vary from 0.2 to 2 s, with the mean being 1.0 s, which is not dependent on photon energies. The spikes exhibit symmetric time profiles with no significant difference between rise and decay times. (3) Among the most energetic spikes, nearly all of them have harder count spectra than their underlying slow-varying components. There is also a weak indication that spikes exhibiting time lags in high-energy emissions tend to have harder spectra than spikes with time lags in low-energy emissions.

In this paper, we analyze hard X-ray spikes observed by RHESSI to understand their temporal, spectral, and spatial properties. A recently developed demodulation code was applied to hard X-ray light curves in several energy bands observed by RHESSI. Hard X-ray spikes were selected from the demodulated flare light curves. We measured the spike duration, the energy-dependent time delay, and count spectral index of these spikes. We also located the hard X-ray source emitting these spikes from RHESSI mapping that was coordinated with imaging observations in visible and UV wavelengths. We identify quickly varying structures of <1 s during the rise of hard X-rays in five flares. These hard X-ray spikes can be observed at photon energies over 100 keV. They exhibit sharp rise and decay with a duration (FWHM) of less than 1 s. Energy-dependent time lags are present in some spikes. It is seen that the spikes exhibit harder spectra than underlying components, typically by 0.5 in the spectral index when they are fitted to power-law distributions. RHESSI clean maps at 25-100 keV with an integration of 2 s centered on the peak of the spikes suggest that hard X-ray spikes are primarily emitted by double foot-point sources in magnetic fields of opposite polarities. With the RHESSI mapping resolution of ~ 4 arsec, the hard X-ray spike maps do not exhibit detectable difference in the spatial structure from sources emitting underlying components. Coordinated high-resolution imaging UV and infrared observations confirm that hard X-ray spikes are produced in magnetic structures embedded in the same magnetic environment of the underlying components. The coordinated high-cadence TRACE UV observations of one event possibly reveal new structures on spatial scales <1-2 arsec at the time of the spike superposed on the underlying component. They are probably sources of hard X-ray spikes.

We examine the \alpha-element abundance ratio, [\alpha/Fe], of 5620 stars, observed by the Sloan Extension for Galactic Understanding and Exploration survey in the region 6 kpc < R < 16 kpc, 0.15 kpc < |Z| < 1.5 kpc, as a function of Galactocentric radius R and distance from the Galactic plane |Z|. Our results show that the high-\alpha\ thick disk population has a short scale length (L_thick ~ 1.8 kpc) compared to the low-\alpha population, which is typically associated with the thin disk. We find that the fraction of high-\alpha\ stars in the inner disk increases at large |Z|, and that high-\alpha\ stars lag in rotation compared to low-\alpha\ stars. In contrast, the fraction of high-\alpha\ stars in the outer disk is low at all |Z|, and high- and low-\alpha\ stars have similar rotational velocities up to 1.5 kpc from the plane. We interpret these results to indicate that different processes were responsible for the high-\alpha\ populations in the inner and outer disk. The high-\alpha\ population in the inner disk has a short scale length and large scale height, consistent with a scenario in which the thick disk forms during an early gas-rich accretion phase. Stars far from the plane in the outer disk may have reached their current locations through heating by minor mergers. The lack of high-\alpha\ stars at large R and |Z| also places strict constraints on the strength of radial migration via transient spiral structure.

The observed radial and vertical metallicity distribution of old stars in the Milky Way disk provides a powerful constraint on the chemical enrichment and dynamical history of the disk. We present the radial metallicity gradient, \Delta[Fe/H]/\Delta R, as a function of height above the plane, |Z|, using 7010 main sequence turnoff stars observed by the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey. The sample consists of mostly old thin and thick disk stars, with a minimal contribution from the stellar halo, in the region 6 < R < 16 kpc, 0.15 < |Z| < 1.5 kpc. The data reveal that the radial metallicity gradient becomes flat at heights |Z| > 1 kpc. The median metallicity at large |Z| is consistent with the metallicities seen in outer disk open clusters, which exhibit a flat radial gradient at [Fe/H] ~ -0.5. We note that the outer disk clusters are also located at large |Z|; because the flat gradient extends to small R for our sample, there is some ambiguity in whether the observed trends for clusters are due to a change in R or |Z|. We therefore stress the importance of considering both the radial and vertical directions when measuring spatial abundance trends in the disk. The flattening of the gradient at high |Z| also has implications on thick disk formation scenarios, which predict different metallicity patterns in the thick disk. A flat gradient, such as we observe, is predicted by a turbulent disk at high redshift, but may also be consistent with radial migration, as long as mixing is strong. We test our analysis methods using a mock catalog based on the model of Sch\"onrich & Binney, and we estimate our distance errors to be ~25%. We also show that we can properly correct for selection biases by assigning weights to our targets.

In this paper, we present comprehensive analysis of a two-ribbon flare observed in UV 1600{\AA} by Transition Region and Coronal Explorer and in HXRs by Reuven Ramaty High Energy Solar Spectroscopic Imager. HXR (25-100 keV) imaging observations show two kernels of size (FWHM) 15?? moving along the two UV ribbons. We find the following results. (1) UV brightening is substantially enhanced wherever and whenever the compact HXR kernel is passing, and during the HXR transit across a certain region, the UV count light curve in that region is temporally correlated with the HXR total flux light curve. After the passage of the HXR kernel, the UV light curve exhibits smooth monotonical decay. (2)We measure the apparent motion speed of the HXR sources and UV ribbon fronts, and decompose the motion into parallel and perpendicular motions with respect to the magnetic polarity inversion line (PIL). It is found that HXR kernels and UV fronts exhibit similar apparent motion patterns and speeds. The parallel motion dominates during the rise of the HXR emission, and the perpendicular motion starts and dominates at the HXR peak, the apparent motion speed being 10-40 km s-1. (3) We also find that UV emission is characterized by a rapid rise correlated with HXRs, followed by a long decay on timescales of 15-30 minutes. The above analysis provides evidence that UV brightening is primarily caused by beam heating, which also produces thick-target HXR emission. The thermal origin of UV emission cannot be excluded, but would produce weaker heating by one order of magnitude. The extended UV ribbons in this event are most likely a result of sequential reconnection along the PIL, which produces individual flux tubes (post-flare loops), subsequent non-thermal energy release and heating in these flux tubes, and then the very long cooling time of the transition region at the feet of these flux tubes.

MUSTANG is a bolometer camera for the Green Bank Telescope (GBT) working at a frequency of 90 GHz. The detector has a field of view of 40 arcseconds. To cancel out random emission change from atmosphere and other sources, requires a fast scanning reflecting system with a few arcminute ranges. In this paper, the aberrations of an off-axis system are reviewed. The condition for an optimized system is provided. In an optimized system, as additional image transfer mirrors are introduced, new aberrations of the off-axis system may be reintroduced, resulting in a limited field of view. In this paper, different scanning mirror arrangements for the GBT system are analyzed through the ray tracing analysis. These include using the subreflector as the scanning mirror, chopping a flat mirror and transferring image with an ellipse mirror, and chopping a flat mirror and transferring image with a pair of face-to-face paraboloid mirrors. The system analysis shows that chopping a flat mirror and using a well aligned pair of paraboloids can generate the required field of view for the MUSTUNG detector system, while other systems all suffer from larger off-axis aberrations added by the system modification. The spot diagrams of the well aligned pair of paraboloids produced is only about one Airy disk size within a scanning angle of about 3 arcmin.

We describe a three hour inquiry activity involving converging lenses and telescopes as part of a semester-long astronomy lab course for non-science majors at Hartnell Community College in Salinas, CA. Students were shown several short demonstrations and given the chance to experiment with the materials, after which there was a class discussion about the phenomena they observed. Students worked in groups of 2-4 to design their own experiments to address a particular question of interest to them and then presented their findings to the class. An instructor-led presentation highlighted the students' discoveries and the lab's content goals, followed by a short worksheet-based activity that guided them in applying their new knowledge to build a simple telescope using two converging lenses. The activity was successful in emphasizing communication skills and giving students opportunities to engage in the process of science in different ways. One of the biggest challenges in designing this activity was covering all of the content given the short amount of time available. Future implementations may have more success by splitting the lab into two sessions, one focusing on converging lenses and the other focusing on telescopes.

(abridged) In the last decade, the advent of enormous galaxy surveys has motivated the development of automated morphological classification schemes to deal with large data volumes. Existing automated schemes can successfully distinguish between early and late type galaxies and identify merger candidates, but are inadequate for studying detailed morphologies of red sequence galaxies. To fill this need, we present a new automated classification scheme that focuses on making finer distinctions between early types roughly corresponding to Hubble types E, S0, and Sa. We visually classify a sample of 984 non-starforming SDSS galaxies with apparent sizes >14". We then develop an automated method to closely reproduce the visual classifications, which both provides a check on the visual results and makes it possible to extend morphological analysis to much larger samples. We visually classify the galaxies into three bulge classes (BC) by the shape of the light profile in the outer regions: discs have sharp edges and bulges do not, while some galaxies are intermediate. We separately identify galaxies with features: spiral arms, bars, clumps, rings, and dust. We find general agreement between BC and the bulge fraction B/T measured by the galaxy modeling package GIM2D, but many visual discs have B/T>0.5. Three additional automated parameters -- smoothness, axis ratio, and concentration -- can identify many of these high-B/T discs to yield automated classifications that agree ~70% with the visual classifications (>90% within one BC). Both methods are used to study the bulge vs. disc frequency as a function of four measures of galaxy 'size': luminosity, stellar mass, velocity dispersion, and radius. All size indicators show a fall in disc fraction and a rise in bulge fraction among larger galaxies.

We have calculated the H$\alpha$ and Ca {\sc ii} 8542 {\AA} line profiles based on four different atmospheric models, including the effects of nonthermal electron beams with various energy fluxes. These two lines have different responses to thermal and nonthermal effects, and can be used to diagnose the thermal and nonthermal heating processes. We apply our method to an X-class flare that occurred on 2001 October 19. We are able to identify quantitatively the heating effects during the flare eruption. We find that the nonthermal effects at the outer edge of the flare ribbon are more notable than that at the inner edge, while the temperature at the inner edge seems higher. On the other hand, the results show that nonthermal effects increase rapidly in the rise phase and decrease quickly in the decay phase, but the atmospheric temperature can still keep relatively high for some time after getting to its maximum. For the two kernels that we analyze, the maximum energy fluxes of the electron beams are $\sim$ 10$^{10}$ and 10$^{11}$ ergs cm$^{-2}$ s$^{-1}$, respectively. However, the atmospheric temperatures are not so high, i.e., lower than or slightly higher than that of the weak flare model F1 at the two kernels. We discuss the implications of the results for two-ribbon flare models.