Downflows with high velocities occur in the transition region above many sunspots; however, how these signatures evolve over short time-scales in both spatial and spectral terms is still unknown. In this article, we investigate the evolution of downflows detected within the transition region on time-scales of the order minutes and search for clues as to their formation mechanisms in co-temporal imaging data. The high-resolution spectral and imaging data used to identify downflows here were sampled by IRIS on the 20th and 21st May 2015. Imaging data from the Hinode and SDO satellites provided context about the wider solar atmosphere. Four downflows were identified and analysed. The potential super-sonic components of these downflows had widths of around 2" and evolved over time-scales of the order minutes. The measured apparent downflow velocities were structured both in time and space, with the highest apparent velocities occurring above a bright region detected in Si IV 1400 {\AA} images. Downflows with apparent velocities below 50 km s^{-1} were observed to extend a few arcsecs away from the foot-points suggesting the potential super-sonic components are linked to larger-scale flows. The electron density and mass flux for these events were within the ranges 10^{9.6}-10^{10.2} cm^{-3} and 10^{-6.81}-10^{-7.48} g cm^{-2} s^{-1}, respectively. Each downflow formed at the foot-point of thin `fingers' extending out around 3-5" in Si IV 1400 {\AA} data with smaller widths (<1") than the super-sonic downflows. Downflows can appear, disappear, and recur within time-scales of less than one hour in sunspots. As potential super-sonic downflow signatures were detected at the foot-points of both extended fingers in Si IV 1400 {\AA} SJI data and sub-sonic downflows in Si IV 1394 \AA\ spectra, it is likely that these events are linked to flows within structures such as coronal loops.

In this article, we present the application of the weighted horizontal gradient of magnetic field (WGM) flare prediction method to 3-dimensional (3D) extrapolated magnetic configurations of 13 flaring solar active regions (ARs). The main aim is to identify an optimal height range, if any, in the interface region between the photosphere and lower corona, where the flare onset time prediction capability of WGM is best exploited. The optimal height is where flare prediction, by means of the WGM method, is achieved earlier than at the photospheric level. 3D magnetic structures, based on potential and non-linear force-free field extrapolations, are constructed to study a vertical range from the photosphere up to the low corona with a 45 km step size. The WGM method is applied as a function of height to all 13 flaring AR cases that are subject to certain selection criteria. We found that applying the WGM method between 1000 and 1800 km above the solar surface would improve the prediction of the flare onset time by around 2-8 hours.Certain caveats and an outlook for future work along these lines are also discussed.

Super-sonic downflows have been observed in transition region spectra above numerous sunspots; however, little research has been conducted to date into how persistent these signatures are within sunspots on time-scales longer than a few hours. Here, we analyse the lead sunspot of AR 12526 to infer the properties and evolution of super-sonic downflows occurring within it. Sixteen large, dense raster scans sampled by IRIS are analysed. These rasters tracked the lead sunspot of AR 12526 at discrete times between the 27th March 2016 and the 2nd April 2016. One sit-and-stare observation acquired on the 1st April 2016 is also studied in order to analyse the evolution of super-sonic downflows on shorter time-scales. Super-sonic downflows are variable within this sunspot both in terms of spatial structuring and velocities. 13 of the 16 raster scans display some evidence of super-sonic downflows in the Si IV 1394 {\AA} line co-spatial to a sustained bright structure detected in the 1400 {\AA} slit-jaw imaging channel, with a peak velocity of 112 km s^{-1} being recorded on the 29th March 2016. Evidence for super-sonic downflows in the O IV 1401 {\AA} line was found in 14 of these rasters, with the spatial structuring differing from that inferred from the Si IV 1394 {\AA} line. In the sit-and-stare data, no dual flow is initially detected, however, a super-sonic downflow does develop after 60 minutes. This downflow accelerates from 73 km s^{-1} to close to 80 km s^{-1} in both the Si IV 1394 {\AA} and O IV 1401 {\AA} lines over the course of 20 minutes before the end of the observation. The morphology of these downflows evolved over the course of both hours and days and was often different in the Si IV 1394 {\AA} and O IV 1401 {\AA} lines. These events were found co-spatial to a bright region in the core of the Si IV 1394 {\AA} line which appeared to form at the foot-points of coronal fan loops.

Context. Velocity or intensity swirls have now been shown to be widely present throughout the photosphere and chromosphere. It has been suggested that these events could contribute to the heating of the upper solar atmosphere, via exciting Alfv\'en pulses, which could carry significant amounts of energy. However, the conjectured necessary physical conditions for their excitation, that the magnetic field rotates co-spatially and co-temporally with the velocity field, has not been verified. Aims. We aim to understand whether photospheric velocity swirls exist co-spatially and co-temporally with photospheric magnetic swirls, in order to demonstrate the link between swirls and pulses. Methods. The automated swirl detection algorithm (ASDA) is applied to the photospheric horizontal velocity and vertical magnetic fields obtained from a series of realistic numerical simulations using the RMHD code Bifrost. The spatial relationship between the detected velocity and magnetic swirls is further investigated via a well-defined correlation index (CI) study. Results. On average, there are ~63 short-lived photospheric velocity swirls (with lifetimes mostly less than 20 s, and average radius of ~37 km and rotating speeds of ~2.5 km/s) detected in a field of view (FOV) of 6\times6 Mm^{-2}, implying a total population of velocity swirls of ~1.06\times10^7 in the solar photosphere. More than 80% of the detected velocity swirls are found to be accompanied by local magnetic concentrations in intergranular lanes. On average, ~71% of the detected velocity swirls have been found to co-exist with photospheric magnetic swirls with the same rotating direction. Conclusions. The co-temporal and co-spatial rotation in the photospheric velocity and magnetic fields provide evidence that the conjectured condition for the excitation of Alfv\'en pulses by photospheric swirls is fulfilled.

Numerous apparent signatures of magnetic reconnection have been reported in the solar photosphere, including inverted-Y shaped jets. The reconnection at these sites is expected to cause localised bi-directional flows and extended shock waves; however, these signatures are rarely observed as extremely high spatial-resolution data are required. Here, we use H-alpha imaging data sampled by the Swedish Solar Telescope's CRisp Imaging SpectroPolarimeter to investigate whether bi-directional flows can be detected within inverted-Y shaped jets near the solar limb. These jets are apparent in the H-alpha line wings, while no signature of either jet is observed in the H-alpha line core, implying reconnection took place below the chromospheric canopy. Asymmetries in the H-alpha line profiles along the legs of the jets indicate the presence of bi-directional flows, consistent with cartoon models of reconnection in chromospheric anemone jets. These asymmetries are present for over two minutes, longer than the lifetimes of Rapid Blue Excursions, and beyond \pm 1 \AA\ into the wings of the line indicating that flows within the inverted-Y shaped jets are responsible for the imbalance in the profiles, rather than motions in the foreground. Additionally, surges form following the occurrence of the inverted-Y shaped jets. This surge formation is consistent with models which suggest such events could be caused by the propagation of shock waves from reconnection sites in the photosphere to the upper atmosphere. Overall, our results provide evidence that magnetic reconnection in the photosphere can cause bi-directional flows within inverted-Y shaped jets and could be the driver of surges.

Active region NOAA 12673 was extremely volatile in 2017 September, producing many solar flares, including the largest of solar cycle 24, an X9.3 flare of 06 September 2017. It has been reported that this flare produced a number of sunquakes along the flare ribbon (Sharykin & Kosovichev 2018; Zhao & Chen 2018). We have used co-temporal and co-spatial Helioseismic and Magnetic Imager (HMI) line-of-sight (LOS) and Swedish 1-m Solar Telescope observations to show evidence of the chromospheric response to these sunquakes. Analysis of the Ca II 8542 \AA\space line profiles of the wavefronts revealed that the crests produced a strong blue asymmetry, whereas the troughs produced at most a very slight red asymmetry. We used the combined HMI, SST datasets to create time-distance diagrams and derive the apparent transverse velocity and acceleration of the response. These velocities ranged from 4.5 km s$^{-1}$ to 29.5 km s$^{-1}$ with a constant acceleration of 8.6 x 10$^{-3}$ km s$^{-2}$. We employed NICOLE inversions, in addition to the Center-of-Gravity (COG) method to derive LOS velocities ranging 2.4 km s$^{-1}$ to 3.2 km s$^{-1}$. Both techniques show that the crests are created by upflows. We believe that this is the first chromospheric signature of a flare induced sunquake.

Ever since their detection two decades ago, standing kink oscillations in coronal loops have been extensively studied both observationally and theoretically. Almost all driven coronal loop oscillations (e.g., by flares) are observed to damp through time often with Gaussian or exponential profiles. Intriguingly, however, it has been shown theoretically that the amplitudes of some oscillations could be modified from Gaussian or exponential profiles if cooling is present in the coronal loop systems. Indeed, in some cases the oscillation amplitude can even increase through time. In this article, we analyse a flare-driven coronal loop oscillation observed by the Solar Dynamics Observatory's Atmospheric Imaging Assembly (SDO/AIA) in order to investigate whether models of cooling can explain the amplitude profile of the oscillation and whether hints of cooling can be found in the intensity evolution of several SDO/AIA filters. During the oscillation of this loop system, the kink mode amplitude appears to differ from a typical Gaussian or exponential profile with some hints being present that the amplitude increases. The application of cooling coronal loop modelling allowed us to estimate the density ratio between the loop and the background plasma, with a ratio of between 2.05-2.35 being returned. Overall, our results indicate that consideration of the thermal evolution of coronal loop systems can allow us to better describe oscillations in these structures and return more accurate estimates of the physical properties of the loops (e.g., density, scale height, magnetic field strength).

The scaling of reaction yields in light ion fusion to low reaction energies is important for our understanding of stellar fuel chains and the development of future energy technologies. Experiments become progressively more challenging at lower reaction energies due to the exponential drop of fusion cross sections below the Coulomb barrier. We report on experiments where deuterium-deuterium (D-D) fusion reactions are studied in a pulsed plasma in the glow discharge regime using a benchtop apparatus. We model plasma conditions using particle-in-cell codes. Advantages of this approach are relatively high peak ion currents and current densities (0.1 to several A/cm^2) that can be applied to metal wire cathodes for several days. We detect neutrons from D-D reactions with scintillator-based detectors. For palladium targets, we find neutron yields as a function of cathode voltage that are over 100 times higher than yields expected for bare nuclei fusion at ion energies below 2 keV (center of mass frame). A possible explanation is a correction to the ion energy due to an electron screening potential of 1000+/-250 eV, which increases the probability for tunneling through the repulsive Coulomb barrier. Our compact, robust setup enables parametric studies of this effect at relatively low reaction energies.

The term "ultraviolet (UV) burst" is introduced to describe small, intense, transient brightenings in ultraviolet images of solar active regions. We inventorize their properties and provide a definition based on image sequences in transition-region lines. Coronal signatures are rare, and most bursts are associated with small-scale, canceling opposite-polarity fields in the photosphere that occur in emerging flux regions, moving magnetic features in sunspot moats, and sunspot light bridges. We also compare UV bursts with similar transition-region phenomena found previously in solar ultraviolet spectrometry and with similar phenomena at optical wavelengths, in particular Ellerman bombs. Akin to the latter, UV bursts are probably small-scale magnetic reconnection events occurring in the low atmosphere, at photospheric and/or chromospheric heights. Their intense emission in lines with optically thin formation gives unique diagnostic opportunities for studying the physics of magnetic reconnection in the low solar atmosphere. This paper is a review report from an International Space Science Institute team that met in 2016-2017.

Swirling motions in the solar atmosphere have been widely observed in recent years and suggested to play a key role in channeling energy from the photosphere into the corona. Here, we present a newly-developed Automated Swirl Detection Algorithm (ASDA) and discuss its applications. ASDA is found to be very proficient at detecting swirls in a variety of synthetic data with various levels of noise, implying our subsequent scientific results are astute. Applying ASDA to photospheric observations with a spatial resolution of 39.2 km sampled by the Solar Optical Telescope (SOT) on-board Hinode, suggests a total number of $1.62\times10^5$ swirls in the photosphere, with an average radius and rotating speed of $\sim290$ km and $< 1.0$ km s$^{-1}$, respectively. Comparisons between swirls detected in Bifrost numerical MHD simulations and both ground-based and space-borne observations, suggest that: 1) the spatial resolution of data plays a vital role in the total number and radii of swirls detected; and 2) noise introduced by seeing effects could decrease the detection rate of swirls, but has no significant influences in determining their inferred properties. All results have shown that there is no significant difference in the analysed properties between counter-clockwise or clockwise rotating swirls. About 70% of swirls are located in intergranular lanes. Most of the swirls have lifetimes less than twice of the cadences, meaning future research should aim to use data with much higher cadences than 6 s. In the conclusions, we propose some promising future research applications where ASDA may provide useful insights.

We report on high-resolution imaging and spectral observations of eruptions of a spiral structure in the transition region, which were taken with the Interface Region Imaging Spectrometer (IRIS), the Atmospheric Imaging Assembly (AIA) and the Helioseismic and Magnetic Imager (HMI). The eruption coincided with the appearance of two series of jets, with velocities comparable to the Alfv\'en speeds in their footpoints. Several pieces of evidence of magnetic braiding in the eruption are revealed, including localized bright knots, multiple well-separated jet threads, transition region explosive events and the fact that all these three are falling into the same locations within the eruptive structures. Through analysis of the extrapolated three-dimensional magnetic field in the region, we found that the eruptive spiral structure corresponded well to locations of twisted magnetic flux tubes with varying curl values along their lengths. The eruption occurred where strong parallel currents, high squashing factors, and large twist numbers were obtained. The electron number density of the eruptive structure is found to be $\sim3\times10^{12}$ cm$^{-3}$, indicating that significant amount of mass could be pumped into the corona by the jets. Following the eruption, the extrapolations revealed a set of seemingly relaxed loops, which were visible in the AIA 94 \AA\ channel indicating temperatures of around 6.3 MK. With these observations, we suggest that magnetic braiding could be part of the mechanisms explaining the formation of solar eruption and the mass and energy supplement to the corona.

Ellerman bombs (EBs) have been widely studied over the past two decades; however, only recently have counterparts of these events been observed in the quiet-Sun. The aim of this article is to further understand small-scale quiet-Sun Ellerman-like brightenings (QSEBs) through research into their spectral signatures, including investigating whether the hot signatures associated with some EBs are also visible co-spatial to any QSEBs. We combine H$\alpha$ and Ca II $8542$ \AA\ line scans at the solar limb with spectral and imaging data sampled by the Interface Region Imaging Spectrograph (IRIS). Twenty one QSEBs were identified with average lifetimes, lengths, and widths measured to be around $120$ s, $0.63$", and $0.35$", respectively. Three of these QSEBs displayed clear repetitive flaring through their lifetimes, comparable to the behaviour of EBs in Active Regions (ARs). Two QSEBs in this sample occurred co-spatial with increased emission in SDO/AIA $1600$ \AA\ and IRIS slit-jaw imager $1400$ \AA\ data, however, these intensity increases were smaller than reported co-spatial to EBs. One QSEB was also sampled by the IRIS slit during its lifetime, displaying increases in intensity in the Si IV $1393$ \AA\ and Si IV $1403$ \AA\ cores as well as the C II and Mg II line wings, analogous to IRIS bursts (IBs). Using RADYN simulations, we are unable to reproduce the observed QSEB H$\alpha$ and Ca II $8542$ \AA\ line profiles leaving the question of the temperature stratification of QSEBs open. Our results imply that some QSEBs could be heated to Transition Region temperatures, suggesting that IB profiles should be observed throughout the quiet-Sun.

Sunspot atmospheres are highly inhomogeneous hosting both quasi-stable and transient features, such as `umbral micro-jets' and dark fibril-like events. We seek to understand the morphological properties and formation mechanisms of small-scale umbral brightenings (SSUBs; analogous to umbral micro-jets) and to understand whether links between these events and short dynamic fibrils, umbral flashes, and umbral dots can be established. An SST filtergram time-series sampling the Ca II H line and a CRISP full-Stokes 15-point Ca II 8542 A line scan dataset were used. The average lifetime and lengths of 54 SSUBs identified in the sunspot umbra are found to be 44.2 seconds (sigma=20 seconds) and 0.56" (sigma=0.14"). The spatial positioning and morphological evolution of these events was investigated finding no evidence of parabolic or ballistic profiles nor a preference for co-spatial formation with umbral flashes. The presence of Stokes V profile reversals provided evidence that these events could form through the development of shocks in the chromosphere. The application of the weak-field approximation indicated that changes in the line-of-sight magnetic field were not responsible for the modifications to the line profile and suggested that thermodynamic effects are the actual cause of the increased emission. Finally, a sub-set of SSUBs were observed to form at the foot-points of short dynamic fibrils. Overall, we found no correlation between the spatial locations where these events were observed and the occurrence of umbral dots and umbral flashes. SSUBs, however, have lifetimes and spectral signatures comparable to umbral flashes and are located at the footpoints of short dynamic fibrils, during or at the end of the red-shifted stage. It is possible, therefore, that these features form due to the shocking of fibrilar material in the lower atmosphere upon its return under gravity.

Burst-like events with signatures in the UV are often observed co-spatial to strong line-of-sight photospheric magnetic fields. Several authors, for example, have noted the spatial relationship between Ellerman bombs (EBs) and Moving Magnetic Features (MMFs), regions of flux which disconnect from a sunspot or pore before propagating away in the moat flow and often displaying evidence of cancellation. In this article, data collected by the Solar Dynamics Observatory's Helioseismic and Magnetic Imager and Atmospheric Imaging Assembly are analysed in an attempt to understand the potential links between such cancellation and UV burst formation. Two MMFs from AR 11579, three bi-poles from AR 11765, and six bi-poles (four of which were co-spatial to IRIS bursts) in AR 11850 were identified for analysis. All of these cancellation features were found to have lifetimes of the order hours and cancellation rates of the order 10^14-10^15 Mx s^-1. H-alpha line wing data from the Dunn Solar Telescope's Interferometric BIdimensional Spectrometer were also available for AR 11579 facilitating a discussion of links between MMFs and EBs. Using an algebraic model of photospheric magnetic reconnection, the measured cancellation rates are then used to ascertain estimates of certain quantities (such as up-flow speeds, jet extents, and potential energy releases) which compared reasonably to the properties of EBs reported within the literature. Our results suggest that cancellation rates of the order measured here are capable of supplying enough energy to drive certain UV bursts (including EBs), however, they are not a guaranteeing condition for burst formation.

Ellerman Bombs (EBs) are often found co-spatial with bipolar photospheric magnetic fields. We use H$\alpha$ imaging spectroscopy along with Fe I 6302.5 \AA\ spectro-polarimetry from the Swedish 1-m Solar Telescope (SST), combined with data from the Solar Dynamic Observatory (SDO) to study EBs and the evolution of the local magnetic fields at EB locations. The EBs are found via an EB detection and tracking algorithm. We find, using NICOLE inversions of the spectro-polarimetric data, that on average (3.43 $\pm$ 0.49) x 10$^{24}$ ergs of stored magnetic energy disappears from the bipolar region during the EBs burning. The inversions also show flux cancellation rates of 10$^{14}$ - 10$^{15}$ Mx s$^{-1}$, and temperature enhancements of 200 K at the detection footpoints. We investigate near-simultaneous flaring of EBs due to co-temporal flux emergence from a sunspot, which shows a decrease in transverse velocity when interacting with an existing, stationary area of opposite polarity magnetic flux and the EBs are formed. We also show that these EBs can get fueled further by additional, faster moving, negative magnetic flux regions.

Chromospheric observations taken at high cadence and high spatial resolution show a range of spicule like features, including Type I, Type II (as well as RBEs and RREs) and those which seem to appear within a few seconds, which if interpreted as flows would imply mass flow velocities in excess of 1000 km/s. This article seeks to quantify and study rapidly appearing spicular type events. We also compare the MOMFBD and speckle reconstruction techniques in order to understand if such spicules are more favourably observed using a particular technique. We use spectral imaging observations taken with the CRISP on the Swedish 1 m Solar Telescope. Data was sampled at multiple positions within the Halpha line profile for both an ondisk and limb location. The data is host to numerous rapidly appearing features which are observed at different locations within the Halpha line profile. The feature's durations vary between 10 and 20 s and lengths around 3500 km. Sometimes, a time delay in their appearance between the blue and red wings of 3 and 5 s is evident, whereas sometimes they are near simultaneous. In some instances features are observed to fade and then re emerge at the same location several tens of seconds later. We provide the first statistical analysis of these spicules and suggest that these observations can be interpreted as the LOS movement of highly dynamic spicules moving in and out of the narrow 60 mA transmission filter used to observe in different parts of the Halpha line profile. The LOS velocity component of the observed fast chromospheric features, manifested as Doppler shifts, are responsible for their appearance in the red and blue wings of Halpha line. Additional work involving data at other wavelengths is required to investigate the nature of their possible wavelike activity.

Ellerman bombs (EBs) have been widely studied in recent years due to their dynamic, explosive nature and apparent links to the underlying photospheric magnetic field implying that they may be formed by magnetic reconnection in the photosphere. Despite a plethora of researches discussing the morphologies of EBs, there has been a limited investigation of how these events appear at the limb, specifically, whether they manifest as vertical extensions away from the disc. In this article, we make use of high-resolution, high-cadence observations of an AR at the solar limb, collected by the CRISP instrument, to identify EBs and infer their physical properties. The upper atmosphere is also probed using the SDO/AIA. We analyse 22 EB events evident within these data, finding that 20 appear to follow a parabolic path away from the solar surface at an average speed of 9 km s^(-1), extending away from their source by 580 km, before retreating back at a similar speed. These results show strong evidence of vertical motions associated with EBs, possibly explaining the dynamical `flaring' (changing in area and intensity) observed in on-disc events. Two in-depth case studies are also presented which highlight the unique dynamical nature of EBs within the lower solar atmosphere. The viewing angle of these observations allows for a direct linkage between these EBs and other small-scale events in the H-alpha line wings, including a potential flux emergence scenario. The findings presented here suggest that EBs could have a wider-reaching influence on the solar atmosphere than previously thought, as we reveal a direct linkage between EBs and an emerging small-scale loop, and other near-by small-scale explosive events. However, as previous research found, these extensions do not appear to impact upon the H-alpha line core, and are not observed by the SDO/AIA EUV filters.

There have been ubiquitous observations of wave-like motions in the solar atmosphere for decades. Recent improvements to space- and ground-based observatories have allowed the focus to shift to smaller magnetic structures on the solar surface. In this paper, high-resolution ground-based data taken using the Swedish 1 m Solar Telescope is combined with co-spatial and co-temporal data from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) satellite to analyze running penumbral waves (RPWs). RPWs have always been thought to be radial wave propagation that occurs within sunspots. Recent research has suggested that they are in fact upwardly propagating field-aligned waves (UPWs). Here, RPWs within a solar pore are observed for the first time and are interpreted as UPWs due to the lack of a penumbra that is required to support RPWs. These UPWs are also observed co-spatially and co-temporally within several SDO/AIA elemental lines that sample the transition region and low corona. The observed UPWs are traveling at a horizontal velocity of around 17 +- 0.5 km s-1 and a minimum vertical velocity of 42 +- 21 km s-1. The estimated energy of the waves is around 150 W m-2, which is on the lower bound required to heat the quiet-Sun corona. This is a new, yet unconsidered source of wave energy within the solar chromosphere and low corona.

We analyse an absorption event within the H$\alpha$ line wings, identified as a surge, and the co-spatial evolution of an EUV brightening, with spatial and temporal scales analogous to a small blinker. We conduct a multi-wavelength, multi-instrument analysis using high-cadence, high-resolution data, collected by the Interferometric BIdimensional Spectrometer on the Dunn Solar Telescope, as well as the space-borne Atmospheric Imaging Assembly and Helioseismic and Magnetic Imager instruments onboard the Solar Dynamics Observatory. One large absorption event situated within the plage region trailing the lead sunspot of AR 11579 is identified within the H-alpha line wings. This event is found to be co-spatially linked to a medium-scale (around 4 arcseconds in diameter) brightening within the transition region and corona. This ejection appears to have a parabolic evolution, first forming in the H-alpha blue wing before fading and reappearing in the H-alpha red wing, and comprises of a number of smaller fibril events. The line-of-sight photospheric magnetic field shows no evidence of cancellation leading to this event. Our research has identified clear evidence that at least a subset of transient brightening events in the transition region are linked to the influx of cooler plasma from the lower solar atmosphere during large eruptive events, such as surges. These observations agree with previous numerical researches on the nature of blinkers and, therefore, suggest that magnetic reconnection is the driver of the analysed surge events; however, further research is required to confirm this.

The presence of photospheric magnetic reconnection has long been thought to give rise to short and impulsive events, such as Ellerman bombs (EBs) and Type II spicules. In this article, we combine high-resolution, high-cadence observations from the Interferometric BIdimensional Spectrometer (IBIS) and Rapid Oscillations in the Solar Atmosphere (ROSA) instruments at the Dunn Solar Telescope, National Solar Observatory, New Mexico with co-aligned Atmospheric Imaging Assembly (SDO/AIA) and Solar Optical Telescope (Hinode/SOT) data to observe small-scale events situated within an active region. These data are then compared with state-of-the-art numerical simulations of the lower atmosphere made using the MURaM code. It is found that brightenings, in both the observations and the simulations, of the wings of the H alpha line profile, interpreted as EBs, are often spatially correlated with increases in the intensity of the FeI 6302.5A line core. Bi-polar regions inferred from Hinode/SOT magnetic field data show evidence of flux cancellation associated, co-spatially, with these EBs, suggesting magnetic reconnection could be a driver of these high-energy events. Through the analysis of similar events in the simulated lower atmosphere, we are able to infer that line profiles analogous to the observations occur co-spatially with regions of strong opposite polarity magnetic flux. These observed events and their simulated counterparts are interpreted as evidence of photospheric magnetic reconnection at scales observable using current observational instrumentation.

Cross-spectral analysis is a mathematical tool for extracting the power spectral density of a correlated signal from two time series in the presence of uncorrelated interfering signals. We demonstrate and explain a set of conditions where the detection of the desired signal using cross-spectral fails partially or entirely in the presence of a second uncorrelated signal. Not understanding when and how this effect occurs can lead to dramatic underreporting of the desired signal. Theoretical, simulated and experimental demonstrations of this effect as well as mitigating methods are presented.

The solar atmosphere is a dynamic environment, constantly evolving to form a wide range of magnetically dominated structures (coronal loops, spicules, prominences, etc.) which cover a significant percentage of the surface at any one time. Oscillations and waves in many of these structures are now widely observed and have led to the new analytic technique of solar magneto-seismology, where inferences of the background conditions of the plasma can be deduced by studying magneto-hydrodynamic (MHD) waves. Here, we generalise a novel magneto-seismological method designed to infer the density distribution of a bounded plasma structure from the relationship of its fundamental and first harmonics. Observations of the solar atmosphere have emphatically shown that stratification, leading to complex density profiles within plasma structures, is common thereby rendering this work instantly accessible to solar physics. We show, in a dynamic waveguide, how the period ratio differs from the idealised harmonic ratios prevalent in homogeneous structures. These ratios show strong agreement with recent observational work. Next, anti-node shifts are also analysed. Using typical scaling parameters for bulk flows within atmospheric waveguides, e.g., coronal loops, it is found that significant anti-node shifts can be predicted, even to the order of 10 Mm. It would be highly encouraged to design specific observations to confirm the predicted anti-node shifts and apply the developed theory of solar magneto-seismology to gain more accurate waveguide diagnostics of the solar atmosphere.

The properties of Ellerman bombs (EBs), small-scale brightenings in the H-alpha line wings, have proved difficult to establish due to their size being close to the spatial resolution of even the most advanced telescopes. Here, we aim to infer the size and lifetime of EBs using high-resolution data of an emerging active region collected using the Interferometric BIdimensional Spectrometer (IBIS) and Rapid Oscillations of the Solar Atmosphere (ROSA) instruments as well as the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We develop an algorithm to track EBs through their evolution, finding that EBs can often be much smaller (around 0.3") and shorter lived (less than 1 minute) than previous estimates. A correlation between G-band magnetic bright points and EBs is also found. Combining SDO/HMI and G-band data gives a good proxy of the polarity for the vertical magnetic field. It is found that EBs often occur both over regions of opposite polarity flux and strong unipolar fields, possibly hinting at magnetic reconnection as a driver of these events.The energetics of EB events is found to follow a power-law distribution in the range of "nano-flare" (10^{22-25} ergs).

Brown dwarfs are important objects because they may provide a missing link between stars and planets, two populations that have dramatically different formation history. In this paper, we present the candidate binaries with brown dwarf companions that are found by analyzing binary microlensing events discovered during 2004 - 2011 observation seasons. Based on the low mass ratio criterion of q < 0.2, we found 7 candidate events, including OGLE-2004-BLG-035, OGLE-2004-BLG-039, OGLE-2007-BLG-006, OGLE-2007-BLG-399/MOA-2007-BLG-334, MOA-2011-BLG-104/OGLE-2011-BLG-0172, MOA-2011-BLG-149, and MOA-201-BLG-278/OGLE-2011-BLG-012N. Among them, we are able to confirm that the companions of the lenses of MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149 are brown dwarfs by determining the mass of the lens based on the simultaneous measurement of the Einstein radius and the lens parallax. The measured mass of the brown dwarf companions are (0.02 +/- 0.01) M_Sun and (0.019 +/- 0.002) M_Sun for MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149, respectively, and both companions are orbiting low mass M dwarf host stars. More microlensing brown dwarfs are expected to be detected as the number of lensing events with well covered light curves increases with new generation searches.

We present the analysis of the light curves of 9 high-magnification single-lens gravitational microlensing events with lenses passing over source stars, including OGLE-2004-BLG-254, MOA-2007-BLG-176, MOA-2007-BLG-233/OGLE-2007-BLG-302, MOA-2009-BLG-174, MOA-2010-BLG-436, MOA-2011-BLG-093, MOA-2011-BLG-274, OGLE-2011-BLG-0990/MOA-2011-BLG-300, and OGLE-2011-BLG-1101/MOA-2011-BLG-325. For all events, we measure the linear limb-darkening coefficients of the surface brightness profile of source stars by measuring the deviation of the light curves near the peak affected by the finite-source effect. For 7 events, we measure the Einstein radii and the lens-source relative proper motions. Among them, 5 events are found to have Einstein radii less than 0.2 mas, making the lenses candidates of very low-mass stars or brown dwarfs. For MOA-2011-BLG-274, especially, the small Einstein radius of $\theta_{\rm E}\sim 0.08$ mas combined with the short time scale of $t_{\rm E}\sim 2.7$ days suggests the possibility that the lens is a free-floating planet. For MOA-2009-BLG-174, we measure the lens parallax and thus uniquely determine the physical parameters of the lens. We also find that the measured lens mass of $\sim 0.84\ M_\odot$ is consistent with that of a star blended with the source, suggesting that the blend is likely to be the lens. Although we find planetary signals for none of events, we provide exclusion diagrams showing the confidence levels excluding the existence of a planet as a function of the separation and mass ratio.

Microlensing can provide a useful tool to probe binary distributions down to low-mass limits of binary companions. In this paper, we analyze the light curves of 8 binary lensing events detected through the channel of high-magnification events during the seasons from 2007 to 2010. The perturbations, which are confined near the peak of the light curves, can be easily distinguished from the central perturbations caused by planets. However, the degeneracy between close and wide binary solutions cannot be resolved with a $3\sigma$ confidence level for 3 events, implying that the degeneracy would be an important obstacle in studying binary distributions. The dependence of the degeneracy on the lensing parameters is consistent with a theoretic prediction that the degeneracy becomes severe as the binary separation and the mass ratio deviate from the values of resonant caustics. The measured mass ratio of the event OGLE-2008-BLG-510/MOA-2008-BLG-369 is $q\sim 0.1$, making the companion of the lens a strong brown-dwarf candidate.

We present the first example of binary microlensing for which the parameter measurements can be verified (or contradicted) by future Doppler observations. This test is made possible by a confluence of two relatively unusual circumstances. First, the binary lens is bright enough (I=15.6) to permit Doppler measurements. Second, we measure not only the usual 7 binary-lens parameters, but also the 'microlens parallax' (which yields the binary mass) and two components of the instantaneous orbital velocity. Thus we measure, effectively, 6 'Kepler+1' parameters (two instantaneous positions, two instantaneous velocities, the binary total mass, and the mass ratio). Since Doppler observations of the brighter binary component determine 5 Kepler parameters (period, velocity amplitude, eccentricity, phase, and position of periapsis), while the same spectroscopy yields the mass of the primary, the combined Doppler + microlensing observations would be overconstrained by 6 + (5 + 1) - (7 + 1) = 4 degrees of freedom. This makes possible an extremely strong test of the microlensing solution. We also introduce a uniform microlensing notation for single and binary lenses, we define conventions, summarize all known microlensing degeneracies and extend a set of parameters to describe full Keplerian motion of the binary lenses.

Nuclear regions of galaxies generally host a mixture of components with different exitation, composition, and kinematics. Derivation of emission line ratios and kinematics could then be misleading, if due correction is not made for the limited spatial and spectral resolutions of the observations.The aim of this paper is to demonstrate, with application to a long slit spectrum of the Seyfert 2 galaxy NGC 1358, how line intensities and velocities, together with modelling and knowledge of the point spread function, may be used to resolve the differing structures. In the situation outlined above, the observed kinematics differs for different spectral lines. From the observed intensity and velocity distributions of a number of spectral lines and with some reasonable assumptions to diminish the number of free parameters, the true line ratios and velocity structures may be deduced. A preliminary solution for the nuclear structure of NGC 1358 is obtained, involving a nuclear point source and an emerging outflow of high exitation, as well as a nuclear emission line disk rotating in the potential of a stellar bulge and expressing a radial excitation gradient. The method results in a likely scenario for the nuclear structure of the Seyfert 2 galaxy NGC 1358. For definitive results an extrapolation of the method to two dimensions combined with the use of integral field spectroscopy will generally be necessary.

We present a HST/STIS spectroscopic and optical/radio imaging study of the Seyfert NGC 2110 aiming to measure the dynamics and understand the nature of the nuclear outflow in the galaxy. Previous HST studies have revealed the presence of a linear structure in the Narrow-Line Region (NLR) aligned with the radio jet. We show that this structure is strongly accelerated, probably by the jet, but is unlikely to be entrained in the jet flow. The ionisation properties of this structure are consistent with photoionisation of dusty, dense gas by the active nucleus. We present a plausible geometrical model for the NLR, bringing together various components of the nuclear environment of the galaxy. We highlight the importance of the circum-nuclear disc in determining the appearance of the emission line gas and the morphology of the jet. From the dynamics of the emission line gas, we place constraints on the accelerating mechanism of the outflow and discuss the relative importance of radio source synchrotron pressure, radio jet ram pressure and nuclear radiation pressure in accelerating the gas. While all three mechanisms can account for the energetics of the emission line gas, gravitational arguments support radio jet ram pressure as the most likely source of the outflow.

We report the discovery of kinematic shock signatures associated with a localized radio jet interaction in the merging Seyfert galaxy NGC 5929. We explore the velocity-dependent ionization structure of the gas and find that low ionization gas at the interaction site is significantly more disturbed than high ionization gas, which we attribute to a local enhancement of shock ionization due to the influence of the jet. The characteristic width of the broad low-ionization emission is consistent with shock velocities predicted from the ionization conditions of the gas. We interpret the relative prominence of shocks to the high density of gas in nuclear environment of the galaxy and place some constraints on their importance as a feedback mechanism in Seyferts.

We present Hubble Space Telescope (HST) WFPC2 photometry of 13 microlensed source stars from the 5.7 year Large Magellanic Cloud (LMC) survey conducted by the MACHO Project. The microlensing source stars are identified by deriving accurate centroids in the ground-based MACHO images using difference image analysis (DIA) and then transforming the DIA coordinates to the HST frame. None of these sources is coincident with a background galaxy, which rules out the possibility that the MACHO LMC microlensing sample is contaminated with misidentified supernovae or AGN in galaxies behind the LMC. This supports the conclusion that the MACHO LMC microlensing sample has only a small amount of contamination due to non-microlensing forms of variability. We compare the WFPC2 source star magnitudes with the lensed flux predictions derived from microlensing fits to the light curve data. In most cases the source star brightness is accurately predicted. Finally, we develop a statistic which constrains the location of the Large Magellanic Cloud (LMC) microlensing source stars with respect to the distributions of stars and dust in the LMC and compare this to the predictions of various models of LMC microlensing. This test excludes at > 90% confidence level models where more than 80% of the source stars lie behind the LMC. Exotic models that attempt to explain the excess LMC microlensing optical depth seen by MACHO with a population of background sources are disfavored or excluded by this test. Models in which most of the lenses reside in a halo or spheroid distribution associated with either the Milky Way or the LMC are consistent which these data, but LMC halo or spheroid models are favored by the combined MACHO and EROS microlensing results.

We investigate the origin of the intrinsic scatter in the correlation between black hole mass (MBH) and bulge luminosity [L(bulge)] in a sample of 45 massive, local (z < 0.35) type~1 active galactic nuclei (AGNs). We derive MBH from published optical spectra assuming a spherical broad-line region, and L(bulge) from detailed two-dimensional decomposition of archival optical Hubble Space Telescope images. AGNs follow the MBH-L(bulge) relation of inactive galaxies, but the zero point is shifted by an average of \Delta log MBH ~ -0.3 dex. We show that the magnitude of the zero point offset, which is responsible for the intrinsic scatter in the MBH-L(bulge) relation, is correlated with several AGN and host galaxy properties, all of which are ultimately related to, or directly impact, the BH mass accretion rate. At a given bulge luminosity, sources with higher Eddington ratios have lower MBH. The zero point offset can be explained by a change in the normalization of the virial product used to estimate MBH, in conjunction with modest BH growth (~ 10%--40%) during the AGN phase. Galaxy mergers and tidal interactions appear to play an important role in regulating AGN fueling in low-redshift AGNs.

We present HST spectroscopy of the nucleus of M31 obtained with STIS. Spectra taken around the CaT lines at 8500 see only the red giants in the double bright- ness peaks P1 and P2. In contrast, spectra taken at 3600-5100 A are sensitive to the tiny blue nucleus embedded in P2, the lower surface brightness red nucleus. P2 has a K-type spectrum, but the embedded blue nucleus has an A-type spectrum with strong Balmer absorption lines. Given the small likelihood for stellar collisions, a 200 Myr old starburst appears to be the most plausible origin of the blue nucleus. In stellar population, size, and velocity dispersion, the blue nucleus is so different from P1 and P2 that we call it P3. The line-of-sight velocity distributions of the red stars in P1+P2 strengthen the support for Tremaine s eccentric disk model. The kinematics of P3 is consistent with a circular stellar disk in Keplerian rotation around a super-massive black hole with M_bh = 1.4 x 10^8 M_sun. The P3 and the P1+P2 disks rotate in the same sense and are almost coplanar. The observed velocity dispersion of P3 is due to blurred rotation and has a maximum value of sigma = 1183+-201 km/s. The observed peak rotation velocity of P3 is V = 618+-81 km/s at radius 0.05" = 0.19 pc corresponding to a circular rotation velocity at this radius of ~1700 km/s. Any dark star cluster alternative to a black hole must have a half-mass radius <= 0.03" = 0.11 pc. We show that this excludes clusters of brown dwarfs or dead stars on astrophysical grounds.

Using 7 years of MACHO survey data, we present a new determination of the optical depth to microlensing towards the Galactic bulge. We select the sample of 62 microlensing events (60 unique) on clump giant sources and perform a detailed efficiency analysis. We use only the clump giant sources because these are bright bulge stars and are not as strongly affected by blending as other events. Using a subsample of 42 clump events concentrated in just 4.5 square degrees, we find tau = 2.17^{+0.47}_{-0.38} x 10^{-6} at (l,b) = (1.50, -2.68), somewhat smaller than found in most previous MACHO studies, but in excellent agreement with recent theoretical predictions. We also present the optical depth in each of the 19 fields in which we detected events, and find limits on optical depth for fields with no events. The errors in optical depth in individual fields are dominated by Poisson noise. We measure optical depth gradients (1.06 +/- 0.71) x 10^{-6} deg^{-1} and (0.29 +/- 0.43) x 10^{-6} deg^{-1} in the galactic latitude b and longitude l directions, respectively. Finally, we discuss the possibility of anomalous duration distribution of events in the field 104 centered on (l,b) = (3.11, -3.01) as well as investigate spatial clustering of events in all fields.

Using The Hubble Space Telescope's Space Telescope Imaging Spectrograph HST's STIS, observations of the OIII emission from the narrow-line region (NLR) of NGC 4151 were obtained and radial velocities determined. Five orbits of HST time were used to obtain spectra at five parallel slit configurations, at a position angle of 58 degrees, with spatial resolution 0.2 arcseconds across and 0.1 arcseconds along each slit. A spectral resolving power of ~ 9,000 with the G430M grating gave velocity measurements accurate to ~ 34 km/s. A kinematic model was generated to match the radial velocities, for comparison to previous kinematic models of biconical radial outflow developed for low-dispersion spectra at two slit positions. The new high-resolution spectra permit the measurement of accurate velocity dispersions for each radial-velocity component. The full-width at half-maximum (FWHM) reaches a maximum of 1000 km/s near the nucleus, and generally decreases with increasing distance to about 100 km/s in the extended narrow-line region (ENLR), starting at about 6 arcseconds from the nucleus. In addition to the bright emission knots, which generally fit our model, there are faint high velocity clouds which do not fit the biconical outflow pattern of our kinematic model. These faint clouds occur at the turnover points of the outflowing bright clouds. We suggest possible scenarios that could explain these rogue clouds: (1) backflow resulting from shocks and (2) outflow outside of the bicones, although the latter does not explain how the knots are ionized and accelerated. A comparison of our observations with a high-resolution radio map shows that there is no evidence that the kinematics of the NLR clouds are affected by the radio lobes that comprise the inner jet.

We present a catalog of 450 high signal-to-noise microlensing events observed by the MACHO collaboration between 1993 and 1999. The events are distributed throughout our fields and, as expected, they show clear concentration toward the Galactic center. No optical depth is given for this sample since no blending efficiency calculation has been performed, and we find evidence for substantial blending. In a companion paper we give optical depths for the sub-sample of events on clump giant source stars, where blending is not a significant effect. Several events with sources that may belong to the Sagittarius dwarf galaxy are identified. For these events even relatively low dispersion spectra could suffice to classify these events as either consistent with Sagittarius membership or as non-Sagittarius sources. Several unusual events, such as microlensing of periodic variable source stars, binary lens events, and an event showing extended source effects are identified. We also identify a number of contaminating background events as cataclysmic variable stars.

Black hole masses in active galactic nuclei (AGN) are difficult to measure using conventional dynamical methods, but can be determined using the technique of reverberation mapping. However, it is important to verify that the results of these different methods are equivalent. This can be done indirectly, using scaling relations between the black hole and the host galaxy spheroid. For this purpose, we have obtained new measurements of the bulge stellar velocity dispersion, sigma, in Seyfert 1 galaxies. These are used in conjunction with the M_bh -- sigma relation to validate nuclear black hole masses, M_bh, in active galaxies determined through reverberation mapping. We find that Seyfert galaxies follow the same M_bh -- sigma relation as non-active galaxies, indicating that reverberation mapping measurements of M_bh are consistent with those obtained using other methods. We also reconsider the relationship between bulge absolute magnitude, M_bulge, and black hole mass. We find that Seyfert galaxies are offset from non-active galaxies, but that the deviation can be entirely understood as a difference in bulge luminosity, not black hole mass; Seyfert hosts are brighter than normal galaxies for a given value of their velocity dispersion, perhaps as a result of younger stellar populations.

For the past 10 years there has been an active debate over whether fast shocks play an important role in ionizing emission line regions in Seyfert galaxies. To investigate this claim, we have studied the Seyfert 2 galaxy Mkn 78, using HST UV/optical images and spectroscopy. Since Mkn 78 provides the archetypal jet-driven bipolar velocity field, if shocks are important anywhere they should be important in this object. Having mapped the emission line fluxes and velocity field, we first compare the ionization conditions to standard photoionization and shock models. We find coherent variations of ionization consistent with photoionization model sequences which combine optically thick and thin gas, but are inconsistent with either autoionizing shock models or photoionization models of just optically thick gas. Furthermore, we find absolutely no link between the ionization of the gas and its kinematic state, while we do find a simple decline of ionization degree with radius. We feel this object provides the strongest case to date against the importance of shock related ionization in Seyferts.

Fourier coefficents have been derived for the $V$ and $R$ light curves of 785 overtone RR Lyrae variables in 16 MACHO fields near the bar of the LMC. The $\phi_{31}$ and $R_{21}$ coefficients have been compared with those of the first overtone RR Lyrae variables in the Galactic globular clusters NGC 6441, M107, M5, M3, M2, $\omega$ Centauri and M68. The results indicate that many of the LMC variables have properties similar to the ones in M2, M3, M5 and the Oosterhoff type I variables in $\omega$ Cen, but they are different from the Oosterhoff type II variables in $\omega$ Cen. Equations derived from hydrodynamic pulsation models have been used to calculate the luminosity and temperature for the 330 bona fide first-overtone variables. The results indicate that they have $\log L$ in the range 1.6 to $1.8\lsun$ and $\log T_{eff}$ between 3.85 and 3.87. Based on these temperatures, a mean color excess $E(V-R) =0.08$ mag, equivalent to $E(B-V)=0.14$ mag, has been estimated for these 330 stars. The 80 M5-like variables (selected according to their location in the $\phi_{31}-\log P$ plot) are used to determine a LMC distance. After correcting for the effects of extinction and crowding, a mean apparent magnitude $<V_0>=18.99 \pm 0.02$ (statistical) $\pm 0.16$ (systematic) has been estimated for these 80 stars. Combining this with a mean absolute magnitude $M_V=0.56\pm 0.06$ for M5-like stars derived from Baade-Wesselink analyses, main sequence fitting, Fourier parameters and the trigonometric parallax of RR Lyrae, we derive an LMC distance modulus $\mu=18.43\pm 0.06$ (statistical) $\pm 0.16$ (systematic) mag. The large systematic error arises from the difficulties of correcting for interstellar extinction and for crowding.

We have frequency analyzed 6391 variables classified earlier as fundamental mode RR Lyrae (RR0) stars in the MACHO database on the Large Magellanic Cloud (LMC). The overwhelming majority of these variables have been proved to be indeed RR0 stars, whereas the remaining ones have fallen in one of the following categories: single- and double-mode Cepheids, binaries, first overtone and double-mode RR Lyrae stars and non-classified variables. Special attention has been paid to the Blazhko stars. We found altogether 731 Blazhko variables showing either a doublet or an equidistant triplet pattern at the main pulsation component in their frequency spectra. The incidence rate of the Blazhko variables among the RR0 stars in the LMC is 11.9%, which is three times higher than their rate among the first overtone RR Lyrae stars. From the frequency spectra we found that variables with larger modulation amplitudes at the higher frequency side of the main pulsation component are three times more frequent than the ones showing the opposite amplitude pattern. A search made for a modulation component with the Blazhko period in the average brightness of the individual variables showed the existence of such a modulation with an overall amplitude of ~0.006 mag. On the other hand, a similar search for quadruple modulation patterns around the main pulsation component have failed to clearly detect such components at the ~0.004 mag level. This means that the amplitudes of the quadruple components should be, on the average, at least ten times smaller than those of the triplet components. This finding and the existence of Blazhko variables with highly asymmetric modulation amplitudes not only question the validity of the magnetic oblique rotator model, but also puts stringent constraints on models based on mode coupling theories.

We present near-infrared J and K-band spectra of the z = 6.4 quasar SDSS J114816.64+525150.3 obtained with the NIRSPEC spectrograph at the Keck-II telescope, covering the rest-frame spectral regions surrounding the C IV 1549 and Mg II 2800 emission lines. The iron emission blend at rest wavelength 2900-3000 A is clearly detected and its strength appears nearly indistinguishable from that of typical quasars at lower redshifts. The Fe II / Mg II ratio is also similar to values found for lower-redshift quasars, demonstrating that there is no strong evolution in Fe/alpha broad-line emission ratios even out to z=6.4. In the context of current models for iron enrichment from Type Ia supernovae, this implies that the SN Ia progenitor stars formed at z > 10. We apply the scaling relations of Vestergaard and of McLure & Jarvis to estimate the black hole mass from the widths of the C IV and Mg II emission lines and the ultraviolet continuum luminosity. The derived mass is in the range (2-6)x10^9 solar masses, with an additional uncertainty of a factor of 3 due to the intrinsic scatter in the scaling relations. This result is in agreement with the previous mass estimate of 3x10^9 solar masses by Willott, McLure, & Jarvis, and supports their conclusion that the quasar is radiating close to its Eddington luminosity.

We describe a few recent microlensing results from the MACHO Collaboration. The aim of the MACHO Project was the identification and quantitative description of dark and luminous matter in the Milky Way using microlensing toward the Magellanic Clouds and Galactic bulge. We start with a discussion of the HST follow-up observations of the microlensing events toward the LMC detected in the first 5 years of the experiment. Using color-magnitude diagrams we attempt to distinguish between two possible locations of the microlensing sources: 1) in the LMC or 2) behind the LMC. We conclude that unless the extinction is extremely patchy, it is very unlikely that most of the LMC events have sources behind the LMC. During an examination of the HST images of the 13 LMC events we found a very red object next to the source star of event LMC-5. Astrometry, microlensing parallax fit, and a spectrum suggest that in this case we directly image the lens - a low-mass disk star. Then we focus on the majority of events observed by the MACHO Project, which are detected toward the Galactic bulge. We argue that the microlensing optical depth toward the bulge is best measured using events that have clump giant sources, which are almost unaffected by blending. From this sample we derive a low optical depth toward the Galactic bulge of (1.4 +/- 0.3) x 10^{-6}, in good agreement with other observational constraints and with theoretical models. The presence of many long-duration events among the bulge candidates allows us to investigate the microlensing parallax effect. Events with the strongest parallax signal are probably due to massive remnants. Events MACHO-96-BLG-5 and MACHO-98-BLG-6 might have been caused by the 6-solar-mass black holes.

We present 47 spectroscopically-confirmed quasars discovered behind the Magellanic Clouds identified via photometric variability in the MACHO database. Thirty-eight quasars lie behind the Large Magellanic Cloud and nine behind the Small Magellanic Cloud, more than tripling the number of quasars previously known in this region. The quasars cover the redshift interval 0.2 < z < 2.8 and apparent mean magnitudes 16.6 < V < 20.1. We discuss the details of quasar candidate selection based on time variability in the MACHO database and present results of spectroscopic follow-up observations. Our follow-up detection efficiency was 20%; the primary contaminants were emission-line Be stars in Magellanic Clouds. For the 47 quasars discovered behind the Magellanic Clouds plus an additional 12 objects previously identified in this region, we present 7.5-year MACHO V- and R-band lightcurves with average sampling times of 2-10 days.

We present an analysis of the longest timescale microlensing events discovered by the MACHO Collaboration during a 7 year survey of the Galactic bulge. We find 6 events that exhibit very strong microlensing parallax signals due, in part, to accurate photometric data from the GMAN and MPS collaborations. The microlensing parallax fit parameters are used in a likelihood analysis, which is able to estimate the distance and masses of the lens objects based upon a standard model of the Galactic velocity distribution. This analysis indicates that the most likely masses of 5 of the 6 lenses are > 1 Msun, which suggests that a substantial fraction of the Galactic lenses may be massive stellar remnants. This could explain the observed excess of long timescale microlensing events. The lenses for events MACHO-96-BLG-5 and MACHO-98-BLG-6 are the most massive, with mass estimates of M/Msun = 6 +10/-3 and M/Msun = 6 +7/-3, respectively. The observed upper limits on the absolute brightness of main sequence stars for these lenses are < 1 Lsun, so both lenses are black hole candidates. The black hole interpretation is also favored by a likelihood analysis with a Bayesian prior using a conventional model for the lens mass function. We consider the possibility that the source stars for some of these 6 events may lie in the foreground or background of the Galactic bulge, but we find that this is unlikely. Future HST observations of these events can either confirm the black hole lens hypothesis or detect the lens stars and provide a direct measurement of their masses. Future observations of similar events by SIM or the Keck or VLTI interferometers will allow direct measurements of the lens masses for stellar remnant lenses as well.

We present a method for solving the lightcurve of an eclipsing binary system which contains a Cepheid variable as one of its components as well as the solutions for three eclipsing Cepheids in the Large Magellanic Cloud (LMC). A geometric model is constructed in which the component stars are assumed to be spherical and on circular orbits. The emergent system flux is computed as a function of time, with the intrinsic variations in temperature and radius of the Cepheid treated self-consistently. Fitting the adopted model to photometric observations, incorporating data from multiple bandpasses, yields a single parameter set best describing the system. This method is applied to three eclipsing Cepheid systems from the MACHO Project LMC database: MACHO ID's 6.6454.5, 78.6338.24 and 81.8997.87. A best-fit value is obtained for each system's orbital period and inclination and for the relative radius, color and limb-darkening coefficients of each star. Pulsation periods and parameterizations of the intrinsic color variations of the Cepheids are also obtained and the amplitude of the radial pulsation of each Cepheid is measured directly. The system 6.6454.5 is found to contain a 4.97-day Cepheid, which cannot be definitely classified as Type I or Type II, with an unexpectedly brighter companion. The system 78.6338.24 consists of a 17.7-day, W Vir Class Type II Cepheid with a smaller, dimmer companion. The system 81.8997.87 contains an intermediate-mass, 2.03-day overtone Cepheid with a dimmer, red giant secondary.

We have performed a search for halo white dwarfs as high proper motion objects in a second epoch WFPC2 image of the Groth-Westphal strip. We identify 24 high proper motion objects with mu > 0.014 ''/yr. Five of these high proper motion objects are identified as strong white dwarf candidates on the basis of their position in a reduced proper motion diagram. We create a model of the Milky Way thin disk, thick disk and stellar halo and find that this sample of white dwarfs is clearly an excess above the < 2 detections expected from these known stellar populations. The origin of the excess signal is less clear. Possibly, the excess cannot be explained without invoking a fourth galactic component: a white dwarf dark halo. We present a statistical separation of our sample into the four components and estimate the corresponding local white dwarf densities using only the directly observable variables, V, V-I, and mu. For all Galactic models explored, our sample separates into about 3 disk white dwarfs and 2 halo white dwarfs. However, the further subdivision into the thin and thick disk and the stellar and dark halo, and the subsequent calculation of the local densities are sensitive to the input parameters of our model for each Galactic component. Using the lowest mean mass model for the dark halo we find a 7% white dwarf halo and six times the canonical value for the thin disk white dwarf density (at marginal statistical significance), but possible systematic errors due to uncertainty in the model parameters likely dominate these statistical error bars. The white dwarf halo can be reduced to around 1.5% of the halo dark matter by changing the initial mass function slightly. The local thin disk white dwarf density in our solution can be made consistent with the canonical value by assuming a larger thin disk scaleheight of 500 pc.

We present preliminary results of the analysis of 5 years of MACHO data on the Galactic bulge microlensing events with clump giants as sources. This class of events allows one to obtain robust conclusions because relatively bright clump stars are not strongly affected by blending. We discuss: 1) the selection of `giant' events, 2) the distribution of event durations, 3) the anomalous character of event durations and optical depth in the MACHO field 104 centered on (l,b) = (3.1, -3.0). We report the preliminary average optical depth of (2.0 +/- 0.4) x 10^{-6} (internal) at (l,b) = (3.9, -3.8), and present a map of the spatial distribution of the optical depth. When field 104 is removed from the sample, the optical depth drops to (1.4 +/- 0.3) x 10^{-6}, which is in excellent agreement with infrared-based models of the central Galactic region.

Galaxies that contain bulges appear to contain central black holes whose masses correlate with the velocity dispersion of the bulge. We show that no corresponding relationship applies in the pure disk galaxy M33. Three-integral dynamical models fit Hubble Space Telescope WFPC2 photometry and STIS spectroscopy best if the central black hole mass is zero. The upper limit is 1500 M_sun. This is significantly below the mass expected from the velocity dispersion of the nucleus and far below any mass predicted from the disk kinematics. Our results suggest that supermassive black holes are associated only with galaxy bulges and not with their disks.

We report the discovery of eight new R Coronae Borealis (RCB) stars in the Large Magellanic Cloud (LMC) using the MACHO project photometry database. The discovery of these new stars increases the number of known RCB stars in the LMC to thirteen. We have also discovered four stars similar to the Galactic variable DY Per. These stars decline much more slowly and are cooler than the RCB stars. The absolute luminosities of the Galactic RCB stars are unknown since there is no direct measurement of the distance to any Galactic RCB star. Hence, the importance of the LMC RCB stars. We find a much larger range of absolute magnitudes (M(V) = -2.5 to -5 mag) than inferred from the small pre-MACHO sample of LMC RCB stars. It is likely that there is a temperature - M(V)relationship with the cooler stars being intrinsically fainter. Cool (~5000 K) RCB stars are much more common than previously thought based on the Galactic RCB star sample. Using the fairly complete sample of RCB stars discovered in the MACHO fields, we have estimated the likely number of RCB stars in the Galaxy to be ~3,200. The SMC MACHO fields were also searched for RCB stars but none were found.

By cross-correlating the results of two recent large-scale surveys, the general properties of a well defined sample of semi-regular variable stars have been determined. ISOGAL mid-infrared photometry and MACHO lightcurves are assembled for approximately 300 stars in the Baade's Windows of low extinction towards the Galactic bulge. These stars are mainly giants of late M spectral type, evolving along the asymptotic giant branch (AGB). They are found to possess a wide and continuous distribution of pulsation periods and to obey an approximate log~period -- bolometric magnitude relation or set of such relations. Approximate mass-loss rates in the range of 1e-8 to 5e-7 M_sun per year are derived from ISOGAL mid-infrared photometry and models of stellar spectra adjusted for the presence of optically-thin circumstellar silicate dust. Mass-loss rates depend on luminosity and pulsation period. Some stars lose mass as rapidly as short-period Miras but do not show Mira-like amplitudes. A period of 70 days or longer is a necessary but not a sufficient condition for mass loss to occur. For AGB stars in the mass-loss ranges that we observe, the functional dependence of mass-loss rate on temperature and luminosity is found to be in agreement with recent theoretical predictions. If we include our mass-loss rates with a sample of extreme mass-losing AGB stars in the Large Magellanic Cloud, we get the general result for AGB stars that mass-loss rate is proportional to luminosity^{2.7}, valid for AGB stars with 10^{-8} to 10^{-4} M_sun per year (Abridged).

We attempt to determine whether the MACHO microlensing source stars are drawn from the average population of the LMC or from a population behind the LMC by examining the HST color-magnitude diagram (CMD) of microlensing source stars. We present WFPC2 HST photometry of eight MACHO microlensing source stars and the surrounding fields in the LMC. The microlensing source stars are identified by deriving accurate centroids in the ground-based MACHO images using difference image analysis (DIA) and then transforming the DIA coordinates to the HST frame. We consider in detail a model for the background population of source stars based on that presented by Zhao, Graff & Guhathakurta. In this model, the source stars have an additional reddening <E(B-V)> = 0.13 mag and a slightly larger distance modulus <Delta mu> ~ 0.3 mag than the average LMC population. We also investigate a series of source star models, varying the relative fraction of source stars drawn from the average and background populations and the displacement of the background population from the LMC. Due to the small number of analyzed events the distribution of probabilities of different models is rather flat. A shallow maximum occurs at a fraction s_LMC ~ 0.8 of the source stars in the LMC. This is consistent with the interpretation that a significant fraction of observed microlensing events are due to lenses in the Milky Way halo, but does not definitively exclude other models.

We have examined the physical conditions in intrinsic UV-absorbing gas in the Seyfert galaxy NGC 4151, using echelle spectra obtained with the Space Telescope Imaging Spectrograph (STIS). We confirm the presence of the kinematic components detected in earlier GHRS observations as well as a new broad absorption feature at a radial velocity of -1680 km/s. The UV continuum of NGC 4151 decreased by a factor of 4 over the previous two years, and we argue the changes in the column density of the low ionization absorption lines associated with the broad component at -490 km/s reflect the decrease in the ionizing flux. Most of the strong absorption lines (e.g., N V, C IV, Si IV) from this component are saturated, but show substantial residual flux in their cores, indicating that the absorber does not fully cover the source of emission. Our interpretation is that the unocculted light is due to scattering by free electrons from an extended region, which reflects continuum, emission lines, and absorption lines. We have been able to constrain the densities for the kinematic components based on absorption lines from metastable states of C III and Fe II, and/or the ratios of ground and fine structure lines of O I,C II, and Si II. We have generated a set of photoionization models which match the ionic column densities for each component during the present low flux state and those seen in previous high flux states with the GHRS and STIS, confirming that the absorbers are photoionized and respond to the changes in the continuum flux. We have been able to map the relative radial positions of the absorbers, and find that the gas decreases in density with distance. None of the UV absorbers is of sufficiently large column density or high enough ionization state to account for the X-ray absorption.

We present photometry and analysis of the microlensing alert MACHO 96-LMC-2. The ~3% photometry provided by the Global Microlensing Alert Network follow--up effort reveals a periodic modulation in the lightcurve. We attribute this to binarity of the lensed source. Microlensing fits to a rotating binary source magnified by a single lens converge on two minima, separated by delta chi^2 ~ 1. The most significant fit X1 predicts a primary which contributes ~100% of the light, a dark secondary, and an orbital period (T) of 9.2 days. The second fit X2 yields a binary source with two stars of roughly equal mass and luminosity, and T = 21.2 days. The lensed object appears to lie on the upper LMC main sequence. We estimate the mass of the primary component of the binary system, M ~2 M_sun. For the preferred model X1, we explore the range of dark companions by assuming 0.1 M_sun and 1.4 M_sun objects in models X1a and X1b, respectively. We find lens velocities projected to the LMC in these models of v^hat_X1a = 18.3 +/- 3.1 km/s and v^hat_X1b = 188 +/- 32 k/ms. In both these cases, a likelihood analysis suggests an LMC lens is preferred over a Galactic halo lens, although only marginally so in model X1b. We also find v^hat_X2 = 39.6 +/- 6.1 k/ms, where the likelihood for the lens location is strongly dominated by the LMC disk. In all cases, the lens mass is consistent with that of an M-dwarf. The LMC self-lensing rate contributed by 96-LMC-2 is consistent with model self-lensing rates. (Abridged)

We report on a search for long duration microlensing events towards the Large Magellanic Cloud. We find none, and therefore put limits on the contribution of high mass objects to the Galactic dark matter. At 95% confidence level we exclude objects in the mass range 0.3 solar masses to 30.0 solar masses from contributing more than 4 times 10^11 solar masses to the Galactic halo. Combined with earlier results, this means that objects with masses under 30 solar masses cannot make up the entire dark matter halo if the halo is of typical size. For a typical dark halo, objects with masses under 10 solar masses contribute less than 40% of the dark matter.

We analyze the nuclear stellar dynamics of the SB0 galaxy NGC 1023, utilizing observational data both from the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope and from the ground. The stellar kinematics measured from these long-slit spectra show rapid rotation (V = 70 km/s at a distance of 0.1 arcsec = 4.9 pc from the nucleus) and increasing velocity dispersion toward the nucleus (where sigma = 295 +/- 30 km/s). We model the observed stellar kinematics assuming an axisymmetric mass distribution with both two and three integrals of motion. Both modeling techniques point to the presence of a central dark compact mass (which presumably is a supermassive black hole) with confidence > 99%. The isotropic two-integral models yield a best-fitting black hole mass of (6.0 +/- 1.4) x 10^7 M_sun and mass-to-light ratio (M/L_V) of 5.38 +/- 0.08, and the goodness-of-fit (chi^2) is insensitive to reasonable values for the galaxy's inclination. The three-integral models, which non-parametrically fit the observed line-of-sight velocity distribution as a function of position in the galaxy, suggest a black hole mass of (3.9 +/- 0.4) x 10^7 M_sun and M/L_V of 5.56 +/- 0.02 (internal errors), and the edge-on models are vastly superior fits over models at other inclinations. The internal dynamics in NGC 1023 as suggested by our best-fit three-integral model shows that the velocity distribution function at the nucleus is tangentially anisotropic, suggesting the presence of a nuclear stellar disk. The nuclear line of sight velocity distribution has enhanced wings at velocities >= 600 km/s from systemic, suggesting that perhaps we have detected a group of stars very close to the central dark mass.

We present a ~ 5-yr optical light curve of the recurrent Be/X-ray transient A0538-66 obtained as a by-product of the MACHO Project. These data reveal both a long-term modulation at P = 420.8 +/- 0.8 d and a short-term modulation at 16.6510 +/- 0.0022 d which, within errors, confirms the previously found orbital period. Furthermore, the orbital activity is only seen at certain phases of the 421 d cycle suggesting that the long-term modulation is related to variations in the Be star envelope.

We present echelle observations of the intrinsic UV absorption lines in the Seyfert galaxy NGC 4151, which were obtained with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope on 1999 July 19. The UV continuum flux at 1450 Angstroms decreased by factor of about four over the previous two years and there was a corresponding dramatic increase in the column densities of the low-ionization absorption lines (e.g., Si II, Fe II, and Al II), presumably as a result of a decrease in the ionizing continuum. In addition to the absorption lines seen in previous low states, we identify a large number of Fe II absorption lines that arise from metastable levels as high as 4.1 eV above the ground state, indicating high densities (> 10^6 cm^-3). We find that the transient absorption feature in the blue wing of the broad C IV emission, seen in a Goddard High Resolution Spectrograph spectrum and thought to be a high-velocity C~IV component, is actually a Si II fine-structure absorption line at a radial velocity of -560 km/s (relative to systemic). We also demonstrate that the ``satellite'' emission lines of C~IV found in International Ultraviolet Explorer spectra are actually regions of unabsorbed continuum plus broad emission that become prominent when the UV continuum of NGC 4151 is in a low state.

We undertake a comparison of observed Algol-type binaries with a library of computed Case A binary evolution tracks. The library consists of 5500 binary tracks with various values of initial primary mass M_{10}, mass ratio q_0, and period P_0, designed to sample the phase-space of Case A binaries in the range -0.10 < log M_{10} < 1.7. Each binary is evolved using a standard code with the assumption that both total mass and orbital angular momentum are conserved. This code follows the evolution of both stars until the point where contact or reverse mass transfer occurs. The resulting binary tracks show a rich variety of behavior which we sort into several subclasses of Case A and Case B. We present the results of this classification, the final mass ratio and the fraction of time spent in Roche Lobe overflow for each binary system. The conservative assumption under which we created this library is expected to hold for a broad range of binaries, where both components have spectra in the range G0 to B1 and luminosity class III - V. We gather a list of relatively well-determined observed hot Algol-type binaries meeting this criterion, as well as a list of cooler Algol-type binaries where we expect significant dynamo-driven mass loss and angular momentum loss. We fit each observed binary to our library of tracks using a chi^2-minimizing procedure. We find that the hot Algols display overall acceptable chi^2, confirming the conservative assumption, while the cool Algols show much less acceptable chi^2 suggesting the need for more free parameters, such as mass and angular momentum loss.

The recent discovery of a correlation between nuclear black hole mass, M_bh, and the stellar velocity dispersion (Gebhardt et al. 2000, Ferrarese and Merritt 2000), in elliptical galaxies and spiral bulges, has raised the question whether such a relationship exists for AGN. Estimates of M_bh for many AGN, made using reverberation mapping techniques, allow exploration of the relationship between black hole mass, the host galaxy and the energetics of nuclear emission. However, since only a few AGN have both M_bh and velocity dispersion measurements, we use the [OIII] 5007 emission line widths on the assumption that for most AGN the forbidden line kinematics are dominated by virial motion in the host galaxy bulge. We find that a relation does exist between M_bh and [OIII] line width for AGN which is similar to the one found by Gebhardt et al. 2000, although with more scatter as expected if secondary influences on the gas kinematics are also present. Our conclusion is that both active and inactive galaxies follow the same relationship between black hole mass and bulge gravitational potential. We find no compelling evidence for systematic differences in the mass estimates from reverberation mapping and stellar dynamics. We also find that for radio quiet AGN the radio power and black hole mass are highly correlated linking emission on scales of kiloparsecs with the nuclear energy source.

We present a simple kinematic model for the narrow-line region (NLR) of the Seyfert 1 galaxy NGC 4151, based on our previous observations of extended [O III] emission with the Space Telescope Imaging Spectrograph (STIS). The model is similar to a biconical radial outflow model developed for the Seyfert 2 galaxy NGC 1068, except that the bicone axis is tilted much more into our line of sight (40 degrees out of the plane of the sky instead of 5 degrees), and the maximum space velocities are lower (750 km/s instead of 1300 km/s. We find evidence for radial acceleration of the emission-line knots to a distance of 160 pc, followed by deceleration that approaches the systemic velocity at a distance of 290 pc (for a distance to NGC 4151 of 13.3 Mpc). Other similarities to the kinematics of NGC 1068 are: 1) there are a number of high-velocity clouds that are not decelerated, suggesting that the medium responsible for the deceleration is patchy, and 2) the bicone in NGC 4151 is at least partially evacuated along its axis. Together, these two Seyfert galaxies provide strong evidence for radial outflow (e.g., due to radiation and/or wind pressure) and against gravitational motion or expansion away from the radio jets as the principal kinematic component in the NLR.

The rotation of the disk of the Large Magellanic Cloud (LMC) is derived from the radial velocities of 422 carbon stars. New aspects of this analysis include the propagation of uncertainties in the LMC proper motion with a Monte Carlo, and a self-consistent modeling of the rotation curve and disk kinematics. The rotation curve is well fit by a truncated, finite-thickness exponential disk model with no dark halo. The velocity dispersion in radial bins decreases from 22 to 15 km/s, then increases to 20 km/s at larger radii. Constant-thickness models in virial equilibrium cannot be reconciled with the data even if the effects of LMC or Galactic dark halos are included. If the disk is virialized, the LMC disk is flared. We model the velocity dispersion at large radii (R > 6 kpc) as a maximal flared disk under the influence of the Galactic dark halo, which favors a mean density for the latter of 0.00025 M_sol/pc^3 at the LMC distance. LMC stellar kinematics play an important role in elucidating the nature of MACHOs, a dark population inferred from LMC microlensing. We have constructed a truncated and flared maximal disk model for the LMC which is kinematically based. Our model does not include a nonvirialized component such as tidal debris. Our upper limit on the self-lensing optical depth is in good agreement with that obtained from less sophisticated models, and is an order of magnitude too small to account for the MACHO microlensing signal [Abridged].

We present preliminary results of the analysis of 5 years of MACHO data on the Galactic bulge microlensing events with clump giants as sources. In particular, we discuss: 1) the selection of `giant' events, 2) distribution of impact parameters, 3) distribution of event durations, 4) the concentration of long duration events in MACHO field 104 centered on (l,b) = (3.1,-3.0). We report the preliminary average optical depth of (2.0 +/- 0.4) x 10^{-6} (internal) at (l,b) = (3.9,-3.8). We discuss future work and prospects for building a coherent model of the Galaxy.

Using the exceptional long-term monitoring capabilities of the MACHO project, we present here the optical history of LMC X-2 for a continuous 6-yr period. These data were used to investigate the previously claimed periodicities for this source of 8.15 h and 12.54 d : we find upper amplitude limits of 0.10 mag and 0.09 mag, respectively.

More than 1300 variables classified provisionally as first overtone RR Lyrae pulsators in the MACHO variable star database of the Large Magellanic Cloud (LMC) have been subjected to standard frequency analysis. Based on the remnant power in the prewhitened spectra, we found 70% of the total population to be monoperiodic. The remaining 30% (411 stars) are classified as one of 9 types according to their frequency spectra. Several types of RR Lyrae pulsational behavior are clearly identified here for the first time. Together with the earlier discovered double-mode (fundamental & first overtone) variables this study increased the number of the known double-mode stars in the LMC to 181. During the total 6.5yr time span of the data, 10% of the stars show strong period changes. We also discovered two additional types of multifrequency pulsators with low occurrence rates of 2% for each. In the first type there remains one closely spaced component after prewhitening by the main pulsation frequency. In the second type the number of remnant components is two, they are also closely spaced, and, in addition, they are symmetric in their frequency spacing relative to the central component. This latter type of variables is associated with their relatives among the fundamental pulsators, known as Blazhko variables. Their high frequency (~20%) among the fundamental mode variables versus the low occurrence rate of their first overtone counterparts makes it more difficult to explain Blazhko phenomenon by any theory depending mainly on the role of aspect angle or magnetic field. Current theoretical models invoke nonradial pulsation components in these stars.

We present the microlensing optical depth towards the Galactic bulge based on the detection of 99 events found in our Difference Image Analysis (DIA) survey. This analysis encompasses three years of data, covering ~ 17 million stars in ~ 4 deg^2, to a source star baseline magnitude limit of V = 23. The DIA technique improves the quality of photometry in crowded fields, and allows us to detect more microlensing events with faint source stars. We find this method increases the number of detection events by 85% compared with the standard analysis technique. DIA light curves of the events are presented and the microlensing fit parameters are given. The total microlensing optical depth is estimated to be tau_(total)= 2.43^(+0.39/-0.38) x 10^(-6) averaged over 8 fields centered at l=2.68 and b=-3.35. For the bulge component we find tau_(bulge)=3.23^(+0.52/-0.50) x 10^(-6) assuming a 25% stellar contribution from disk sources. These optical depths are in good agreement with the past determinations of the MACHO Alcock et al. (1997) and OGLE Udalski et al. (1994) groups, and are higher than predicted by contemporary Galactic models. We show that our observed event timescale distribution is consistent with the distribution expected from normal mass stars, if we adopt the stellar mass function of Scalo (1986) as our lens mass function. However, we note that as there is still disagreement about the exact form of the stellar mass function, there is uncertainty in this conclusion. Based on our event timescale distribution we find no evidence for the existence of a large population of brown dwarfs in the direction of the Galactic bulge.

The MACHO project is a search for dark matter in the form of massive compact halo objects (MACHOs). The project has photometrically monitored tens of millions of stars in the Large Magellanic Cloud (LMC), Small Magellanic Cloud (SMC), and Galactic bulge in search of rare gravitational microlensing events caused by these otherwise invisible objects. In 5.7 years of observations toward the LMC some 13-17 microlensing events have been observed by the MACHO survey, allowing powerful statements to be made about the nature of the dark population in the halo of our Galaxy. A critical component of these statements is an accurate determination of the survey's detection efficiency. The detection efficiency is a complicated function of temporal sampling, stellar crowding (the luminosity function), image quality, photometry, time-series analysis, and criteria used to select the microlensing candidates. Such a complex interdependence is most naturally solved using a Monte Carlo approach. Here we describe the details of the Monte Carlo used to calculate the efficiency presented in the MACHO 5.7-year LMC results. Here we correct several shortcomings of past determinations, including (1) adding fainter source stars (2.5 magnitudes below our faintest detected "stars"), (2) an up-to-date luminosity function for the LMC, (3) better sampling of real images in both stellar density and observing conditions, (4) an improved scheme for adding artificial microlensing onto a random sample of real lightcurves, and many other improvements. [Abridged]

We present a 9 million star color-magnitude diagram (9M CMD) of the LMC bar. The 9M CMD reveals a complex superposition of different age and metallicity stellar populations, with important stellar evolutionary phases occurring over 3 orders of magnitude in number density. First, we count the non-variable supergiants, the associated Cepheids, and measure the effective temperatures defining the instability strip. Lifetime predictions of stellar evolution theory are tested, with implications for the origin of low-luminosity Cepheids. The highly-evolved AGB stars have a bimodal distribution in brightness, which we interpret as discrete old populations (>1 Gyr). The faint AGB may be metal-poor and very old. We identify the clusters NGC 411 and M3 as templates for the admixture of old stellar populations. However, there are indications that the old and metal-poor field population has a red HB morphology: the RR Lyraes lie on the red edge of the instability strip, the AGB-bump is very red, and the ratio of AGB-bump stars to RR Lyraes is quite large. If the HB second parameter is age, the old and metal-poor field population likely formed after the oldest clusters. Lifetime predictions of stellar evolution theory lead us to associate a significant fraction of the red HB clump giants with the same old and metal-poor population producing the RR Lyraes and the AGB-bump. In this case, compared to the age-dependent luminosity predictions of stellar evolution theory, the red HB clump is too bright relative to the RR Lyraes and AGB-bump. Last, the surface density profile of RR Lyraes is fit by an exponential, favoring a disk-like rather than spheroidal distribution. We conclude that the age of the LMC disk is probably similar to the age of the Galactic disk. (ABRIDGED)

We report on our search for microlensing towards the Large Magellanic Cloud (LMC). Analysis of 5.7 years of photometry on 11.9 million stars in the LMC reveals 13 - 17 microlensing events. This is significantly more than the $\sim$ 2 to 4 events expected from lensing by known stellar populations. The timescales ($\that$) of the events range from 34 to 230 days. We estimate the microlensing optical depth towards the LMC from events with $2 < \that < 400$ days to be 1.2 ^{+0.4}_ {-0.3} \ten{-7}$, with an additional 20% to 30% of systematic error. The spatial distribution of events is mildly inconsistent with LMC/LMC disk self-lensing, but is consistent with an extended lens distribution such as a Milky Way or LMC halo. Interpreted in the context of a Galactic dark matter halo, consisting partially of compact objects, a maximum likelihood analysis gives a MACHO halo fraction of 20% for a typical halo model with a 95% confidence interval of 8% to 50%. A 100% MACHO halo is ruled out at the 95% C.L. for all except our most extreme halo model. Interpreted as a Galactic halo population, the most likely MACHO mass is between $ 0.15 \msun$ and $ 0.9 \msun$, depending on the halo model, and the total mass in MACHOs out to 50 kpc is found to be 9+4-3 10^{10} msun, independent of the halo model. These results are marginally consistent with our previous results, but are lower by about a factor of two. Besides a larger data set, this work also includes an improved efficiency determination, improved likelihood analysis, and more thorough testing of systematic errors, especially with respect to the treatment of potential backgrounds to microlensing, such as supernovae in galaxies behind the LMC. [Abridged]

We have obtained new Hubble Space Telescope (HST) observations of M84, a nearby massive elliptical galaxy whose nucleus contains a 1.5 X 10^9 Msun dark compact object, which presumably is a supermassive black hole. Our Space Telescope Imaging Spectrograph (STIS) spectrum provides the first clear detection of emission lines in the blue (e.g., [O II] 3727, Hbeta, and [O III] 4959, 5007), which arise from a compact region 0.28" across centered on the nucleus. Our Near Infrared Camera and Multi-Object Spectrometer (NICMOS) images exhibit the best view through the prominent dust lanes evident at optical wavelengths and provide a more accurate correction for the internal extinction. The relative fluxes of the emission lines we have detected in the blue together with those detected in the wavelength range 6295 - 6867 AA by Bower et al. (1998, ApJ, 492, L111) indicate that the gas at the nucleus is photoionized by a nonstellar process, instead of hot stars. Stellar absorption features from cool stars at the nucleus are very weak. We update the spectral energy distribution of the nuclear point source and find that although it is roughly flat in most bands, the optical to UV continuum is very red, similar to the spectral energy distribution of BL Lac. Thus, the nuclear point source seen in high-resolution optical images (Bower et al. 1997, ApJ, 483, L33) is not a star cluster but is instead a nonstellar source. Assuming isotropic emission from this source, we estimate that the ratio of bolometric luminosity to Eddington luminosity is 5 X 10^(-7). However, this could be underestimated if this source is a misaligned BL Lac object, which is a possibility suggested by the spectral energy distribution and the evidence of optical variability we describe.

We present UBV photometry of the eclipsing binary Harvard Variable 2274 in the Large Magellanic Cloud (LMC). The stellar parameters of this binary system were calculated by Guinan et al. (1998a) who gave both a reddening towards HV 2274 of E(B-V) = 0.083 +/- 0.006 and a distance modulus to the LMC, mu_{LMC} = 18.42 +/- 0.007. The reddening of this system was also determined by Udalski et al. (1998) who found E(B-V) = 0.149 +\- 0.015. With Udalski et al. (1998) B and V photometry Guinan et al.(1998b) obtained E(B-V) = 0.120 +\- 0.009 and mu_{LMC} = 18.30 +\- 0.07. Using our UBV photometry, we derive a reddening of E(B-V)= 0.088 +\- 0.025, consistent with the original value of Guinan et al. (1998a) and supporting a longer distance modulus to the LMC of about mu_{LMC} = 18.40 +\- 0.07. We stress the uncertainties inherent in ground-based UBV photometry and the concomitant uncertainties in determining distances based upon such photometry.

We have examined the physical conditions in the narrow-line region of the well-studied Seyfert galaxy NGC 4151, using long-slit spectra obtained with the Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS). The data were taken along a position angle of 221 degrees, centered on the optical nucleus. We have generated photoionization models for a contiguous set of radial zones, out to 2.3 arcsec in projected position to the southwest of the nucleus, and 2.7 arcsec to the northeast. Given the uncertainties in the reddening correction, the calculated line ratios successfully matched nearly all the dereddened ratios. We find that the narrow-line region consists of dusty atomic gas photoionized by a power-law continuum that has been modified by transmission through a mix of low and high ionization gas, specifically UV and X-ray absorbing components. The physical characteristics of the absorbers resemble those observed along our line of sight to the nucleus, although the column density of the X-ray absorber is a factor of ten less than observed. The large inferred covering factor of the absorbing gas is in agreement with the results of our previous study of UV absorption in Seyfert 1 galaxies. We find evidence, specifically the suppression of L-alpha, that we are observing the back-end of dusty ionized clouds in the region southwest of the nucleus. Since these clouds are blueshifted, this supports the interpretation of the cloud kinematics as being due to radial outflow from the nucleus. We find that the narrow-line gas at each radial position is inhomogeneous, and can be modeled as consisting of a radiation-bounded component and a more tenuous, matter-bounded component. The density of the narrow-line gas dropswith increasing radial distance, which confirms our earlier results.

Longslit spectra of the Seyfert galaxy NGC 4151 from the UV to near infrared have been obtained with STIS to study the kinematics and physical conditions in the NLR. The kinematics show evidence for three components, a low velocity system in normal disk rotation, a high velocity system in radial outflow at a few hundred km/s relative to the systemic velocity and an additional high velocity system also in outflow with velocities up to 1400 km/s, in agreement with results from STIS slitless spectroscopy (Hutchings et al., 1998, Kaiser et al., 1999, Hutchings et al., 1999) We have explored two simple kinematic models and suggest that radial outflow in the form of a wind is the most likely explanation. We also present evidence indicating that the wind may be decelerating with distance from the nucleus. We find that the emission line ratios along our slits are all entirely consistent with photoionization from the nuclear continuum source. A decrease in the [OIII]5007/H-beta and [OIII]5007/[OII]3727 ratios suggests that the density decreases with distance from the nucleus. This trend is borne out by the [SII] ratios as well. We find no strong evidence for interaction between the radio jet and the NLR gas in either the kinematics or the emission line ratios in agreement with the results of Kaiser et al. (1999) who find no spatial coincidence of NLR clouds and knots in the radio jet. These results are in contrast to other recent studies of nearby AGN which find evidence for significant interaction between the radio source and the NLR gas.

The MACHO Project is a microlensing survey that monitors the brightnesses of 60 million stars in the Large Magellanic Cloud (LMC), Small Magellanic Cloud, and Galactic bulge. Our database presently contains about 80 billion photometric measurements, a significant fraction of all astronomical photometry. We describe the calibration of MACHO two-color photometry and transformation to the standard Kron-Cousins V and R system. Calibrated MACHO photometry may be properly compared with all other observations on the Kron-Cousins standard system, enhancing the astrophysical value of these data. For 9 million stars in the LMC bar, independent photometric measurements of 20,000 stars with V < 18 mag in field-overlap regions demonstrate an internal precision sigV = 0.021, sigR = 0.019, and sig(V-R) = 0.028 mag. The accuracy of the zero-point in this calibration is estimated to be +-0.035 mag for stars with colors in the range -0.1 < (V-R) < 1.2 mag. A comparison of calibrated MACHO photometry with published photometric sequences and new Hubble Space Telescope observations shows agreement. The current calibration zero-point uncertainty for the remainder of the MACHO photometry database is estimated to be +- 0.10 mag in V or R and +-0.04 mag in (V-R). We describe the first application of calibrated MACHO photometry data: the construction of a color-magnitude diagram used to calculate our experimental sensitivity to detect microlensing in the LMC.

Narrow-band imaging of the nuclear region of NGC 4151 with the Hubble Space Telescope is presented. The filter bandpasses isolate line emission in various high velocity ranges in several ions. Slitless and long-slit spectra of the region with the Space Telescope Imaging Spectrograph also indicate the locations of high velocity gas. These emission regions are faint and are interspersed among the bright emission clouds seen in direct images. They have radial velocities up to 1400 km/s relative to the nucleus, and are found in both approach and recession on both sides of the nucleus. This contrasts strongly with the bright emission line clouds which have been discussed previously as showing bidirectional outflow with velocities within 400 km/s of the nucleus. We discuss the possible connections of the high velocity material with the radio jet and the nuclear radiation.

We present slitless spectra of the Narrow Line Region (NLR) in NGC4151 from the Space Telescope Imaging Spectrograph (STIS) on HST, and investigate the kinematics and physical conditions of the emission line clouds in this region. Using medium resolution (~0.5 Angstrom) slitless spectra at two roll angles and narrow band undispersed images, we have mapped the NLR velocity field from 1.2 kpc to within 13 pc (H_o=75 km/s/Mpc) of the nucleus. The inner biconical cloud distribution exhibits recessional velocities relative to the nucleus to the NE and approaching velocities to the SW of the nucleus. We find evidence for at least two kinematic components in the NLR. One kinematic component is characterized by Low Velocities and Low Velocity Dispersions (LVLVD clouds: |v| < 400 km/s, and Delta_v < 130 km/s). This population extends through the NLR and their observed kinematics may be gravitationally associated with the host galaxy. Another component is characterized by High Velocities and High Velocity Dispersions (HVHVD clouds: 400 < |v| < ~1700 km/s, Delta_v > 130 km/s). This set of clouds is located within 1.1 arcsec (~70pc) of the nucleus and has radial velocities which are too high to be gravitational in origin, but show no strong correlation between velocity or velocity dispersion and the position of the radio knots. Outflow scenarios will be discussed as the driving mechanism for these HVHVD clouds.

We have measured the strengths of Ca II Triplet and Mgb stellar absorption lines in the nuclear and off-nuclear spectra of Seyfert galaxies. These features are diluted to varying degrees by continuum emission from the active nucleus and from young stars. Ca II Triplet strengths can be enhanced if late-type supergiant stars dominate the near-IR light. Thus, objects with strong Ca II Triplet and weak Mgb lines may be objects with strong bursts of star formation. We find that for most of our sample the line strengths are at least consistent with dilution of a normal galaxy spectrum by a power law continuum, in accord with the standard model for AGN. However, for several Seyferts in our sample, it appears that dilution by a power law continuum cannot simultaneously explain strong Ca II Triplet and relatively weak Mgb. Also, these objects occupy the region of the IRAS color-color diagram characteristic of starburst galaxies. In these objects it appears that the optical to near-IR emission is dominated by late-type supergiants produced in a circumnuclear burst of star formation.

We present and discuss visible-wavelength long-slit spectra of four low redshift 3C galaxies obtained with the STIS instrument on the Hubble Space Telescope. The slit was aligned with near-nuclear jet-like structure seen in HST images of the galaxies, to give unprecedented spatial resolution of the galaxy inner regions. In 3C 135 and 3C 171, the spectra reveal clumpy emission line structures that indicate outward motions of a few hundred km s$^{-1}$ within a centrally illuminated and ionised biconical region. There may also be some low-ionisation high-velocity material associated with 3C 135. In 3C 264 and 3C 78, the jets have blue featureless spectra consistent with their proposed synchrotron origin. There is weak associated line emission in the innermost part of the jets with mild outflow velocity. These jets are bright and highly collimated only within a circumnuclear region of lower galaxy luminosity, which is not dusty. We discuss the origins of these central regions and their connection with relativistic jets.