We investigate the spin dynamics of a millisecond pulsar (MSP) in compact orbit around a Kerr-like massive black hole with an general mass quadrupole. We use the Mathisson-Papetrou-Dixon formulation to compute the orbital and spin evolution of the MSP, accounting for the non-linear interaction of the pulsar's energy-momentum tensor on the background spacetime metric. We investigate how the MSP spin and BH quadrupole moment manifest in the pulsar spin-orbital dynamics. We discuss the astrophysical observational implications of these spin and orbital dynamics on the timing of a radio pulsar in an Extreme Mass Ratio Binary, e.g. a Galactic Centre pulsar. In particular, notable timing variations in the Einstein delay and Roemer delay are observed, along with modifications to the pulsar pulse profile.

This paper presents a class of novel high-order accurate discontinuous Galerkin (DG) schemes for the compressible Euler equations under gravitational fields. A notable feature of these schemes is that they are well-balanced for a general hydrostatic equilibrium state, and at the same time, provably preserve the positivity of density and pressure. In order to achieve the well-balanced and positivity-preserving properties simultaneously, a novel DG spatial discretization is carefully designed with suitable source term reformulation and a properly modified Harten-Lax-van Leer contact (HLLC) flux. Based on some technical decompositions as well as several key properties of the admissible states and HLLC flux, rigorous positivity-preserving analyses are carried out. It is proven that the resulting well-balanced DG schemes, coupled with strong stability preserving time discretizations, satisfy a weak positivity property, which implies that one can apply a simple existing limiter to effectively enforce the positivity-preserving property, without losing high-order accuracy and conservation. The proposed methods and analyses are applicable to the Euler system with general equation of state. Extensive one- and two-dimensional numerical tests demonstrate the desired properties of these schemes, including the exact preservation of the equilibrium state, the ability to capture small perturbation of such state, the robustness for solving problems involving low density and/or low pressure, and good resolution for smooth and discontinuous solutions.

Pulsars (PSRs) orbiting intermediate or supermassive black holes at the centre of galaxies and globular clusters are known as Extreme Mass Ratio Binaries (EMRBs) and have been identified as precision probes of strong-field GR. For appropriate orbital parameters, some of these systems may also emit gravitational radiation in a `burst-like' pattern. The observation of this burst radiation in conjunction with the electromagnetic radio timing signal would allow for multimessenger astronomy in strong-field gravitational regimes. In this work we investigate gravitational radiation from these PSR-EMRBs, calculating the waveforms and SNRs and explore the influence of this GW on the pulsar radio signal. We find that for typical PSR-EMRBs, gravitational burst radiation should be detectable from both the Galactic centre and the centre of stellar clusters, and that this radiation will not meaningfully affect the pulsar timing signal, allowing PSR-EMRB to remain `clean' test-beds of strong-field GR.

We propose and analyze a class of robust, uniformly high-order accurate discontinuous Galerkin (DG) schemes for multidimensional relativistic magnetohydrodynamics (RMHD) on general meshes. A distinct feature of the schemes is their physical-constraint-preserving (PCP) property, i.e., they are proven to preserve the subluminal constraint on the fluid velocity and the positivity of density, pressure, and specific internal energy. Developing PCP high-order schemes for RMHD is highly desirable but remains a challenging task, especially in the multidimensional cases, due to the inherent strong nonlinearity in the constraints and the effect of the magnetic divergence-free condition. Inspired by some crucial observations at the PDE level, we construct the provably PCP schemes by using the locally divergence-free DG schemes of the recently proposed symmetrizable RMHD equations as the base schemes, a limiting technique to enforce the PCP property of the DG solutions, and the strong-stability-preserving methods for time discretization. We rigorously prove the PCP property by using a novel "quasi-linearization" approach to handle the highly nonlinear physical constraints, technical splitting to offset the influence of divergence error, and sophisticated estimates to analyze the beneficial effect of the additional source term in the symmetrizable RMHD system. Several two-dimensional numerical examples are provided to confirm the PCP property and to demonstrate the accuracy, effectiveness and robustness of the proposed PCP schemes.

Outflows from starburst galaxies can be driven by thermal pressure, radiation and cosmic rays. We present an analytic phenomenological model that accounts for these contributions simultaneously to investigate their effects on the hydrodynamical properties of outflows. We assess the impact of energy injection, wind opacity, magnetic field strength and the mass of the host galaxy on flow velocity, temperature, density and pressure profiles. For an M82-like wind, a thermally-dominated driving mechanism is found to deliver the fastest and hottest wind. Radiation-driven winds in typical starburst-galaxy configurations are unable to attain the higher flow velocities and temperatures associated with thermal and cosmic ray-driven systems, leading to higher wind densities which would be more susceptible to cooling and fragmentation at lower altitudes. High opacity winds are more sensitive to radiative driving, but terminal flow velocities are still lower than those achieved by other driving mechanisms at realistic opacities. We demonstrate that variations in the outflow magnetic field can influence its coupling with cosmic rays, where stronger fields enable greater streaming but less driving near the base of the flow, instead with cosmic rays redirecting their driving impact to higher altitudes. The gravitational potential is less important in M82-like wind configurations, and substantial variations in the flow profiles only emerge at high altitude in massive haloes. This model offers a more generalised approach to examine the large scale hydrodynamical properties for a wide variety of starburst galaxies.

The detection of a pulsar (PSR) in a tight, relativistic orbit around a supermassive or intermediate mass black hole - such as those in the Galactic centre or in the centre of Globular clusters - would allow for precision tests of general relativity (GR) in the strong-field, non-linear regime. We present a framework for calculating the theoretical time-frequency signal from a PSR in such an Extreme Mass Ratio Binary (EMRB). This framework is entirely relativistic with no weak-field approximations and so able to account for all higher-order strong-field gravitational effects, relativistic spin dynamics, the convolution with astrophysical effects and the combined impact on the PSR timing signal. Specifically we calculate both the spacetime path of the pulsar radio signal and the complex orbital and spin dynamics of a spinning pulsar around a Kerr black hole, accounting for spacetime curvature and frame dragging, relativistic and gravitational time delay, gravitational light bending, temporal and spatial dispersion induced by the presence of plasma along the line of sight and relativistic aberration. This then allows for a consistent time-frequency solution to be generated. Such a framework is key for assessing the use of PSR as probes of strong field GR, helping to inform the detection of an EMRB system hosting a PSR and, most essentially, for providing an accurate theoretical basis to then compare with observations to test fundamental physics.

Large-scale outflows from starburst galaxies are multi-phase, multi-component fluids. Charge-exchange lines which originate from the interfacing surface between the neutral and ionised components are a useful diagnostic of the cold dense structures in the galactic outflow. From the charge-exchange lines observed in the nearby starburst galaxy M82, we conduct surface-to-volume analyses and deduce that the cold dense clumps in its galactic outflow have flattened shapes, resembling a hamburger or a pancake morphology rather than elongated shapes. The observed filamentary H$\alpha$ features are therefore not prime charge-exchange line emitters. They are stripped material torn from the slow moving dense clumps by the faster moving ionised fluid which are subsequently warmed and stretched into elongated shapes. Our findings are consistent with numerical simulations which have shown that cold dense clumps in galactic outflows can be compressed by ram pressure, and also progressively ablated and stripped before complete disintegration. We have shown that some clumps could survive their passage along a galactic outflow. These are advected into the circumgalactic environment, where their remnants would seed condensation of the circumgalactic medium to form new clumps. The infall of these new clumps back into the galaxy and their subsequent re-entrainment into the galactic outflow form a loop process of galactic material recycling.

Faraday rotation measure at radio wavelengths is commonly used to diagnose large-scale magnetic fields. It is argued that the length-scales on which magnetic fields vary in large-scale diffuse astrophysical media can be inferred from correlations in the observed RM. RM is a variable which can be derived from the polarised radiative transfer equations in restrictive conditions. This paper assesses the usage of RMF (rotation measure fluctuation) analyses for magnetic field diagnostics in the framework of polarised radiative transfer. We use models of various magnetic field configurations and electron density distributions to show how density fluctuations could affect the correlation length of the magnetic fields inferred from the conventional RMF analyses. We caution against interpretations of RMF analyses when a characteristic density is ill-defined, e.g. in cases of log-normal distributed and fractal-like density structures. As the spatial correlations are generally not the same in the line-of-sight longitudinal direction and the sky plane direction, one also needs to clarify the context of RMF when inferring from observational data. In complex situations, a covariant polarised radiative transfer calculation is essential to capture all aspects of radiative and transport processes, which would otherwise ambiguate the interpretations of magnetism in galaxy clusters and larger-scale cosmological structures.

Detailed studies of the Sun have shown that sunspots and solar flares are closely correlated. Photometric data from Kepler/K2 has allowed similar studies to be carried out on other stars. Here, we utilise TESS photometric 2-min cadence of 167 low mass stars from Sectors 1 - 3 to investigate the relationship between starspots and stellar flares. From our sample, 90 percent show clear rotational modulation likely due to the presence of a large, dominant starspot and we use this to determine a rotational period for each star. Additionally, each low mass star shows one or more flares in its lightcurve and using Gaia DR2 parallaxes and SkyMapper magnitudes we can estimate the energy of the flares in the TESS band-pass. Overall, we have 1834 flares from the 167 low mass stars with energies from $6.0\times 10^{29}$ - $2.4\times 10^{35}$~erg. We find none of the stars in our sample show any preference for rotational phase suggesting the lack of a correlation between the large, dominant star spot and flare number. We discuss this finding in greater detail and present further scenarios to account for the origin of flares on these low mass stars.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

This paper presents entropy symmetrization and high-order accurate entropy stable schemes for the relativistic magnetohydrodynamic (RMHD) equations. It is shown that the conservative RMHD equations are not symmetrizable and do not possess an entropy pair. To address this issue, a symmetrizable RMHD system, which admits a convex entropy pair, is proposed by adding a source term into the equations. Arbitrarily high-order accurate entropy stable finite difference schemes are then developed on Cartesian meshes based on the symmetrizable RMHD system. The crucial ingredients of these schemes include (i) affordable explicit entropy conservative fluxes which are technically derived through carefully selected parameter variables, (ii) a special high-order discretization of the source term in the symmetrizable RMHD system, and (iii) suitable high-order dissipative operators based on essentially non-oscillatory reconstruction to ensure the entropy stability. Several benchmark numerical tests demonstrate the accuracy and robustness of the proposed entropy stable schemes of the symmetrizable RMHD equations.

Quenching of star-formation has been identified in many starburst and post-starburst galaxies, indicating burst-like star-formation histories (SFH) in the primordial Universe. We have investigated the role of high energy cosmic rays (CRs) in such environments, particularly how they could contribute to this burst-like SFH via quenching and feedback. These high energy particles interact with the baryon and radiation fields of their host via hadronic processes to produce secondary leptons. The secondary particles then also interact with ambient radiation fields to generate X-rays through inverse-Compton scattering. In addition, they can thermalise directly with the semi-ionised medium via Coulomb processes. Heating at a rate of $\sim 10^{-25} \; \text{erg}~\text{cm}^{-3}~\text{s}^{-1}$ can be attained by Coulomb processes in a star-forming galaxy with one core-collapse SN event per decade, and this is sufficient to cause quenching of star-formation. At high-redshift, a substantial amount of CR secondary electron energy can be diverted into inverse-Compton X-ray emission. This yields an X-ray luminosity of above $10^{41}~\text{erg}~\text{s}^{-1}$ by redshift $z=7$ which drives a further heating effect, operating over larger scales. This would be able to halt inflowing cold gas filaments, strangulating subsequent star-formation. We selected a sample of 16 starburst and post-starburst galaxies at $7\lesssim z \lesssim 9$ and determine the star-formation rates they could have sustained. We applied a model with CR injection, propagation and heating to calculate energy deposition rates in these 16 sources. Our calculations show that CR feedback cannot be neglected as it has the strength to suppress star-formation in these systems. We also show that their currently observed quiescence is consistent with the suffocation of cold inflows, probably by a combination of X-ray and CR heating.

We investigate the spin dynamics of a millisecond pulsar (MSP) in a tightly bounded orbit around a massive black hole. These binaries are progenitors of the extreme-mass-ratio-inspirals (EMRIs) and intermediate-mass-ratio-inspirals (IMRIs) gravitational wave events. The Mathisson-Papapetrou-Dixon (MPD) formulation is used to determine the orbital motion and spin modulation and evolution. We show that the MSP will not be confined in a planar Keplerian orbit and its spin will exhibit precession and nutation induced by spin-orbit coupling and spin-curvature interaction. These spin and orbital behaviours will manifest observationally in the temporal variations in the MSP's pulsed emission and, with certain geometries, in the self-occultation of the pulsar's emitting poles. Radio pulsar timing observations will be able to detect such signatures. These extreme-mass-ratio binaries (EMRBs) and intermediate-mass-ratio binaries (IMRBs) are also strong gravitational wave sources. Combining radio pulsar timing and gravitational wave observations will allow us to determine the dynamics of these systems in high precision and hence the subtle behaviours of spinning masses in strong gravity.

We investigate the propagation of light through a plasma on a background Kerr spacetime via a Hamiltonian formulation. The behaviour of light when propagating through a vacuum and through a plasma is not the same; the convolution of gravitational and plasma effects gives rise to a dispersion in both space and time. The magnitude of the dispersion is a strong function of both the ray frequency and impact parameter. We discuss implications for the detection of gravitationally bent pulsar beams near the Galactic centre.

We investigate the interactions of energetic hadronic particles with the media in outflows from star-forming protogalaxies. These particles undergo pion-producing interactions which can drive a heating effect in the outflow, while those advected by the outflow also transport energy beyond the galaxy, heating the circumgalactic medium. We investigate how this process evolves over the length of the outflow and calculate the corresponding heating rates in advection-dominated and diffusion-dominated cosmic ray transport regimes. In a purely diffusive transport scenario, we find the peak heating rate reaches $10^{-26}\;\! {\rm erg~cm}^{-3}\;\! {\rm s}^{-1}$ at the base of the outflow where the wind is driven by core-collapse supernovae at an event rate of 0.1 $\text{yr}^{-1}$, but does not extend beyond 2 kpc. In the advection limit, the peak heating rate is reduced to $10^{-28}\;\! {\rm erg~cm}^{-3}\;\! {\rm s}^{-1}$, but its extent can reach to tens of kpc. Around 10% of the cosmic rays injected into the system can escape by advection with the outflow wind, while the remaining cosmic rays deliver an important interstellar heating effect. We apply our cosmic ray heating model to the recent observation of the high-redshift galaxy MACS1149-JD1 and show that it could account for the quenching of a previous starburst inferred from spectroscopic observations. Re-ignition of later star-formation may be caused by the presence of filamentary circumgalactic inflows which are reinstated after cosmic ray heating has subsided.

Polarization of radiation is a powerful tool to study cosmic magnetism and analysis of polarization can be used as a diagnostic tool for large-scale structures. In this paper, we present a solid theoretical foundation for using polarized light to investigate large-scale magnetic field structures: the cosmological polarized radiative transfer (CPRT) formulation. The CPRT formulation is fully covariant. It accounts for cosmological and relativistic effects in a self-consistent manner and explicitly treats Faraday rotation, as well as Faraday conversion, emission, and absorption processes. The formulation is derived from the first principles of conservation of phase-space volume and photon number. Without loss of generality, we consider a flat Friedmann-Robertson-Walker (FRW) space-time metric and construct the corresponding polarized radiative transfer equations. We propose an all-sky CPRT calculation algorithm, based on a ray-tracing method, which allows cosmological simulation results to be incorporated and, thereby, model templates of polarization maps to be constructed. Such maps will be crucial in our interpretation of polarized data, such as those to be collected by the Square Kilometer Array (SKA). We describe several tests which are used for verifying the code and demonstrate applications in the study of the polarization signatures in different distributions of electron number density and magnetic fields. We present a pencil-beam CPRT calculation and an all-sky calculation, using a simulated galaxy cluster or a model magnetized universe obtained from GCMHD+ simulations as the respective input structures. The implications on large-scale magnetic field studies are discussed; remarks on the standard methods using rotation measure are highlighted.

The bright transient AT2018cow has been unlike any other known type of transient. Its high brightness, rapid rise and decay and initially nearly featureless spectrum are unprecedented and difficult to explain using models for similar burst sources. We present evidence for faint gamma-ray emission continuing for at least 8 days, and featureless spectra in the ultraviolet bands -- both unusual for eruptive sources. The X-ray variability of the source has a burst-like character. The UV-optical spectrum does not show any CNO line but is well described by a blackbody. We demonstrate that a model invoking the tidal disruption of a 0.1 - 0.4 Msun Helium White Dwarf (WD) by a 100,000 to one million solar mass Black Hole (BH) located in the outskirts of galaxy Z~137-068 could provide an explanation for most of the characteristics shown in the multi-wavelength observations. A blackbody-like emission is emitted from an opaque photosphere, formed by the debris of the WD disruption. Broad features showing up in the optical/infrared spectra in the early stage are probably velocity broadened lines produced in a transient high-velocity outward moving cocoon. The asymmetric optical/infrared lines that appeared at a later stage are emission from an atmospheric layer when it detached from thermal equilibrium with the photosphere, which undergoes more rapid cooling. The photosphere shrinks when its temperature drops, and the subsequent infall of the atmosphere produced asymmetric line profiles. Additionally, a non-thermal jet might be present, emitting X-rays in the 10-150 keV band.

Radial velocity is one of key measurements in understanding the fundamental properties of stars, stellar clusters and the Galaxy. A plate of stars in the Kepler field were observed in May of 2018 with the medium-resolution spectrographs of LAMOST, aiming to test the performance of this new system which is the upgraded equipment of LAMOST after the first five-year regular survey.We present our analysis on the radial velocity measurements (RVs) derived from these data. The results show that slight and significant systematic errors exist among the RVs obtained from the spectra collected by different spectrographs and exposures, respectively. After correcting the systematic errors with different techniques, the precision of RVs reaches ~1.3, ~1.0, ~0.5 and ~0.3 km/s at S/Nr = 10, 20, 50, and 100, respectively. Comparing with the RVs of the standard stars of the APOGEE survey, our RVs are calibrated with a zero-point shift of ~7 km/s. The results indicate that the LAMOST medium-resolution spectroscopic system may provide RVs in a reasonable accuracy and precision for the selected targets.

The intracluster medium of galaxy clusters is permeated by {\mu}G magnetic fields. Observations with current and future facilities have the potential to illuminate the role of these magnetic fields play in the astrophysical processes of galaxy clusters. To obtain a greater understanding of how the initial seed fields evolve to the magnetic fields in the intracluster medium requires magnetohydrodynamic simulations. We critically assess the current Smoothed Particle Magneto-Hydrodynamics (SPMHD) schemes, especially highlighting the impact of a hyperbolic divergence cleaning scheme and artificial resistivity switch on the magnetic field evolution in cosmological simulations of the formation of a galaxy cluster using the N-body/SPMHD code gcmhd++. The impact and performance of the cleaning scheme and two different schemes for the artificial resistivity switch is demonstrated via idealized test cases and cosmological simulations. We demonstrate that the hyperbolic divergence cleaning scheme is effective at suppressing the growth of the numerical divergence error of the magnetic field and should be applied to any SPMHD simulation. Although the artificial resistivity is important in the strong field regime, it can suppress the growth of the magnetic field in the weak field regime, such as galaxy clusters. With sufficient resolution, simulations with divergence cleaning can reproduce observed magnetic fields. We conclude that the cleaning scheme alone is sufficient for galaxy cluster simulations, but our results indicate that the SPMHD scheme must be carefully chosen depending on the regime of the magnetic field.

The magnetosphere of strongly magnetized neutron stars, such as magnetars, can sustain large electric currents. The charged particles return to the surface with large Lorentz factors, producing a particle bombardment. We investigate the transport of radiation in the atmosphere of strongly magnetized neutron stars, in the presence of particle bombardment heating. We solve the radiative transfer equations for a gray atmosphere in the Eddington approximation, accounting for the polarization induced by a strong magnetic field. The solutions show the formation of a hot external layer and a low (uniform) temperature atmospheric interior. This suggests that the emergent spectrum may be described by a single blackbody with the likely formation of a optical/infrared excess (below ~ 1 eV). We also found that the polarization is strongly dependent on both the luminosity and penetration length of the particle bombardment. Therefore, the thermal emission from active sources, such as transient magnetars, in which the luminosity decreases by orders of magnitude, may show a substantial variation in the polarization pattern during the outburst decline. Our results may be relevant in view of future X-ray polarimetric missions such as IXPE and eXTP.

We investigate emission signatures of binary compact star gravitational wave sources consisting of strongly magnetized neutron stars (NSs) and/or white dwarfs (WDs) in their late-time inspiral phase. Because of electromagnetic interactions between the magnetospheres of the two compact stars, a substantial amount of energy will be extracted, and the resultant power is expected to be $\sim 10^{38} - 10^{44}$ erg/s in the last few seconds before the two stars merge, when the binary system contains a NS with a surface magnetic field $10^{12}$ G. The induced electric field in the process can accelerate charged particles up to the EeV energy range. Synchrotron radiation is emitted from energetic electrons, with radiative energies reaching the GeV energy for binary NSs and the MeV energy for NS - WD or double WD binaries. In addition, a blackbody component is also presented and it peaks at several to hundreds keV for binary NSs and at several keV for NS - WD or double WD binaries. The strong angular dependence of the synchrotron radiation and the isotropic nature of the blackbody radiation lead to distinguishable modulation patterns between the two emission components. If coherent curvature radiation is presented, fast radio bursts could be produced. These components provide unique simultaneous electromagnetic signatures as precursors of gravitational wave events associated with magnetized compact star mergers and short gamma ray bursts (e.g., GRB 100717).

We present an empirical stellar spectra library created using spectra from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) DR5. This library represents a uniform data set ranging from 3750 through 8500K in effective temperature (T_eff), from -2.5 through +1.0 dex in metallicity ([Fe/H]), and from 0 to 5.0 dex in gravity (log g). The spectra in the library have resolution R~1800, with well-calibrated fluxes and rest-framed wavelengths. Using a large number of red stars observed by LAMOST, we generated denser K type templates to fill in data missing from current empirical spectral libraries, particularly the late K type. For K giants, we calibrated the spectroscopic surface gravities against the asteroseismic surface gravities. To verify the reliability of the parameters labeled for this library, we performed an internal cross-validation by using a chi^2 minimization method to interpolate parameters of each individual spectrum using the remaining spectra in the library. We obtained precisions of 41 K, 0.11 dex, and 0.05 dex for T_eff, log g, and [Fe/H], respectively, which means the templates are labeled with correct stellar parameters. Through external comparisons, we confirmed that measurements of the stellar parameters through this library can achieve accuracies of approximately 125K in T_eff, 0.1 dex in [Fe/H] and 0.20 dex in log g without systematic offset. This empirical library is useful for stellar parameter measurements because it has large parameter coverage and full wavelength coverage from 3800 to 8900A.

We investigate the interactions of energetic hadronic particles (cosmic ray protons) with photons and baryons in protogalactic environments, where the target photons are supplied by the first generations of stars to form in the galaxy and the cosmological microwave background, while the target baryons are the interstellar and circumgalactic medium. We show that pair-production and photo-pion processes are the dominant interactions at particle energies above $10^{19}\;\! {\rm eV}$, while ${\rm pp}$-interaction pion-production dominates at the lower energies in line with expectations from, for example $\gamma$-ray observations of star-forming galaxies and dense regions of our own galaxy's interstellar medium. We calculate the path lengths for the interaction channels and determine the corresponding rates of energy deposition. We have found that protogalactic magnetic fields and their evolution can significantly affect the energy transport and energy deposition processes of cosmic rays. Within a Myr after the onset of star-formation the magnetic field in a protogalaxy could attain a strength sufficient to confine all but the highest energy particles within the galaxy. This enhances the cosmic ray driven self-heating of the protogalaxy to a rate of around $10^{-24}\;\! {\rm erg} \;\!{\rm cm}^{-3}\;\! {\rm s}^{-1}$ for a galaxy with strong star-forming activity that yields 1 core collapse SN event per year. This heating power exceeds even that due to radiative emission from the protogalaxy's stellar populations. However, in a short window before the protogalaxy is fully magnetised, energetic particles could stream across the galaxy freely, delivering energy into the circumgalactic and intergalactic medium.

This paper proposes and analyzes arbitrarily high-order discontinuous Galerkin (DG) and finite volume methods which provably preserve the positivity of density and pressure for the ideal MHD on general meshes. Unified auxiliary theories are built for rigorously analyzing the positivity-preserving (PP) property of MHD schemes with a HLL type flux on polytopal meshes in any space dimension. The main challenges overcome here include establishing relation between the PP property and discrete divergence of magnetic field on general meshes, and estimating proper wave speeds in the HLL flux to ensure the PP property. In 1D case, we prove that the standard DG and finite volume methods with the proposed HLL flux are PP, under condition accessible by a PP limiter. For multidimensional conservative MHD system, standard DG methods with a PP limiter are not PP in general, due to the effect of unavoidable divergence-error. We construct provably PP high-order DG and finite volume schemes by proper discretization of symmetrizable MHD system, with two divergence-controlling techniques: locally divergence-free elements and a penalty term. The former leads to zero divergence within each cell, while the latter controls the divergence error across cell interfaces. Our analysis reveals that a coupling of them is important for positivity preservation, as they exactly contribute the discrete divergence-terms absent in standard DG schemes but crucial for ensuring the PP property. Numerical tests confirm the PP property and the effectiveness of proposed PP schemes. Unlike conservative MHD system, the exact smooth solutions of symmetrizable MHD system are proved to retain the positivity even if the divergence-free condition is not satisfied. Our analysis and findings further the understanding, at both discrete and continuous levels, of the relation between the PP property and the divergence-free constraint.

We report a detailed analysis of the orbital properties of binary millisecond pulsar (MSP) with a white dwarf (WD) companion. Positive correlations between the orbital period $P_{\rm b}$ and eccentricity $\epsilon$ are found in two classes of MSP binaries with a He WD and with a CO/ONeMg WD, though their trends are different. The distribution of $P_{\rm b}$ is not uniform. Deficiency of sources at $P_{\rm b}\sim35-50$~days (Gap 1) have been mentioned in previous studies. On the other hand, another gap at $P_{\rm b}\sim2.5-4.5$~days (Gap 2) is identified for the first time. Inspection of the relation between $P_{\rm b}$ and the companion masses $M_{\rm c}$ revealed the subpopulations of MSP binaries with a He WD separated by Gap 1, above which $P_{\rm b}$ is independent of $M_{\rm c}$ (horizontal branch) but below which $P_{\rm b}$ correlates strongly with $M_{\rm c}$ (lower branch). Distinctive horizontal branch and lower branch separated by Gap 2 were identified for the MSP binaries with a CO/ONeMg WD at shorter $P_{\rm b}$ and higher $M_{\rm c}$. Generally, $M_{\rm c}$ are higher in the horizontal branch than in the lower branch for the MSP binaries with a He WD. These properties can be explained in terms of a binary orbital evolution scenario in which the WD companion was ablated by a pulsar wind in the post mass-transfer phase.

We present an analysis of K2 short cadence data of 34 M dwarfs which have spectral types in the range M0 - L1. Of these stars, 31 showed flares with a duration between $\sim$10-90 min. Using distances obtained from Gaia DR2 parallaxes, we determined the energy of the flares to be in the range $\sim1.2\times10^{29}-6\times10^{34}$ erg. In agreement with previous studies we find rapidly rotating stars tend to show more flares, with evidence for a decline in activity in stars with rotation periods longer than $\sim$10 days. The rotational modulation seen in M dwarf stars is widely considered to result from a starspot which rotates in and out of view. Flux minimum is therefore the rotation phase where we view the main starspot close to the stellar disk center. Surprisingly, having determined the rotational phase of each flare in our study we find none show any preference for rotational phase. We outline three scenarios which could account for this unexpected finding. The relationship between rotation phase and flare rate will be explored further using data from wide surveys such as NGTS and TESS.

The density and pressure are positive physical quantities in magnetohydrodynamics (MHD). Design of provably positivity-preserving (PP) numerical schemes for ideal compressible MHD is highly desirable, but remains a challenge especially in the multidimensional cases. In this paper, we first develop uniformly high-order discontinuous Galerkin (DG) schemes which provably preserve the positivity of density and pressure for multidimensional ideal MHD. The schemes are constructed by using the locally divergence-free DG schemes for the symmetrizable ideal MHD equations as the base schemes, a PP limiter to enforce the positivity of the DG solutions, and the strong stability preserving methods for time discretization. The significant innovation is that we discover and rigorously prove the PP property of the proposed DG schemes by using a novel equivalent form of the admissible state set and very technical estimates. Several two-dimensional numerical examples further confirm the PP property, and demonstrate the accuracy, effectiveness and robustness of the proposed PP methods.

Isolated Neutron Stars are some of the most exciting stellar objects known to astronomers: they have the most extreme magnetic fields, with values up to $10^{15}$ G, and, with the exception of stellar-mass black holes, they are the most dense stars, with densities of $\approx 10^{14}$ g cm$^{-3}$. As such, they are perfect laboratories to test theories of electromagnetism and nuclear physics under conditions of magnetic field and density unattainable on Earth. In particular, the interaction of radiation with strong magnetic fields is the cause of the {\em vacuum birefringence}, an effect predicted by quantum electrodynamics in 1936 but that lacked an observational evidence until now. Here, we show how the study of the polarisation of the optical radiation from the surface of an isolated neutron star yielded such an observational evidence, opening exciting perspectives for similar studies at other wavelengths.

In this work, we present a catalog of 2651 carbon stars from the fourth Data Release (DR4) of the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST). Using an efficient machine-learning algorithm, we find out these stars from more than seven million spectra. As a by-product, 17 carbon-enhanced metal-poor (CEMP) turnoff star candidates are also reported in this paper, and they are preliminarily identified by their atmospheric parameters. Except for 176 stars that could not be given spectral types, we classify the other 2475 carbon stars into five subtypes including 864 C-H, 226 C-R, 400 C-J, 266 C-N, and 719 barium stars based on a series of spectral features. Furthermore, we divide the C-J stars into three subtypes of CJ( H), C-J(R), C-J(N), and about 90% of them are cool N-type stars as expected from previous literature. Beside spectroscopic classification, we also match these carbon stars to multiple broadband photometries. Using ultraviolet photometry data, we find that 25 carbon stars have FUV detections and they are likely to be in binary systems with compact white dwarf companions.

The routinely flaring events from \sgras\ trace dynamic, high-energy processes in the immediate vicinity of the supermassive black hole. We statistically study temporal and spectral properties, as well as fluence and duration distributions, of the flares detected by the \chandra\ X-ray Observatory from 1999 to 2012. The detection incompleteness and bias are carefully accounted for in determining these distributions. We find that the fluence distribution can be well characterized by a power-law with a slope of $1.73^{+0.20}_{-0.19}$, while the durations ($\tau$ in seconds) by a log-normal function with a mean $\log(\tau)=3.39^{+0.27}_{-0.24}$ and an intrinsic dispersion $\sigma=0.28^{+0.08}_{-0.06}$. No significant correlation between the fluence and duration is detected. The apparent positive correlation, as reported previously, is mainly due to the detection bias (i.e., weak flares can be detected only when their durations are short). These results indicate that the simple self-organized criticality model has difficulties in explaining these flares. We further find that bright flares usually have asymmetric lightcurves with no statistically evident difference/preference between the rising and decaying phases in terms of their spectral/timing properties. Our spectral analysis shows that although a power-law model with a photon index of $2.0\pm0.4$ gives a satisfactory fit to the joint spectra of strong and weak flares, there is weak evidence for a softer spectrum of weaker flares. This work demonstrates the potential to use statistical properties of X-ray flares to probe their trigger and emission mechanisms, as well as the radiation propagation around the black hole.

Galaxies at high redshifts with strong star formation are sources of high-energy cosmic rays. These cosmic rays interact with the baryon and radiation fields of the galactic environment via photo-pair, photo-pion and proton-proton processes to produce charged and neutral pions, neutrons and protons. The cosmic rays thereby deposit energy into the interstellar medium (ISM) as they propagate. We show how energy transport and deposition by ultra high-energy cosmic rays is regulated by the evolution of the galaxy, in particular by the development of the galactic magnetic field. We show how the particle-driven energy deposition can influence the thermal evolution of the host and its surroundings. Using a parametric protogalaxy model, we calculate the heating effect on the ISM as the cosmic rays are increasingly confined by the magnetic evolution of the galaxy.

The general-relativistic magnetohydrodynamical (GRMHD) formulation for black hole-powered jets naturally gives rise to a stagnation surface, wherefrom inflows and outflows along magnetic field lines that thread the black hole event horizon originate. We derive a conservative formulation for the transport of energetic electrons which are initially injected at the stagnation surface and subsequently transported along flow streamlines. With this formulation the energy spectra evolution of the electrons along the flow in the presence of radiative and adiabatic cooling is determined. For flows regulated by synchrotron radiative losses and adiabatic cooling, the effective radio emission region is found to be finite, and geometrically it is more extended along the jet central axis. Moreover, the emission from regions adjacent to the stagnation surface is expected to be the most luminous as this is where the freshly injected energetic electrons concentrate. An observable stagnation surface is thus a strong prediction of the GRMHD jet model with the prescribed non-thermal electron injection. Future millimeter/sub-millimeter (mm/sub-mm) very-long-baseline interferometric (VLBI) observations of supermassive black hole candidates, such as the one at the center of M87, can verify this GRMHD jet model and its associated non-thermal electron injection mechanism.

The recent detection of a large polarization degree in the optical emission of an isolated neutron star led to the suggestion that this has been the first evidence of vacuum polarization in a strong magnetic field, an effect predicted by quantum electrodynamics but never observed before. This claim was challanged in a paper by Capparelli, Maiani & Polosa (2017), according to whom a much higher polarization degree would be necessary to positively identify vacuum polarization. Here we show that their conclusions are biased by several inadequate assumptions and have no impact on the original claim.

We model the intermediate-mass black hole HLX-1, using the Hubble Space Telescope, XMM-Newton and Swift. We quantify the relative contributions of a bluer component, function of X-ray irradiation, and a redder component, constant and likely coming from an old stellar population. We estimate a black hole mass of about (2^{+2}_{-1}) x 10^4 M_{sun}, a spin parameter a/M ~ 0.9 for moderately face-on view, and a peak outburst luminosity of about 0.3 times the Eddington luminosity. We discuss the discrepancy between the characteristic sizes inferred from the short X-ray timescale (R ~ a few 10^{11} cm) and from the optical emitter (R sqrt[cos theta] ~ 2.2 x 10^{13} cm). One possibility is that the optical emitter is a circumbinary disk; however, we disfavour this scenario because it would require a very small donor star. A more plausible scenario is that the disk is large but only the inner annuli are involved in the X-ray outburst. We propose that the recurrent outbursts are caused by an accretion-rate oscillation driven by wind instability in the inner disk. We argue that the system has a long-term-average accretion rate of a few percent Eddington, just below the upper limit of the low/hard state; a wind-driven oscillation can trigger transitions to the high/soft state, with a recurrence period of ~1 year (much longer than the binary period, which we estimate as ~10 days). The oscillation that dominated the system in the last decade is now damped such that the accretion rate no longer reaches the level required to trigger a transition. Finally, we highlight similarities between disk winds in HLX-1 and in the Galactic black hole V404 Cyg.

The task of flux calibration for LAMOST (Large sky Area Multi-Object Spectroscopic Telescope) spectra is difficult due to many factors. For example, the lack of standard stars, flat fielding for large field of view, and variation of reddening between different stars especially at low galactic latitudes etc. Poor selection, bad spectral quality, or extinction uncertainty of standard stars not only might induce errors to the calculated spectral response curve (SRC), but also might lead to failures in producing final 1D spectra. In this paper, we inspected spectra with Galactic latitude |b|>=60 degree and reliable stellar parameters, determined through the LAMOST Stellar Parameter Pipeline (LASP), to study the stability of the spectrograph. To guarantee the selected stars had been observed by each fiber, we selected 37,931 high quality exposures of 29,000 stars from LAMOST DR2, and more than 7 exposures for each fiber. We calculated the SRCs for each fiber for each exposure, and calculated the statistics of SRCs for spectrographs with both the fiber variations and time variations. The result shows that the average response curve of each spectrograph (henceforth ASPSRC) is relatively stable with statistical errors <= 10%. From the comparison between each ASPSRC and the SRCs for the same spectrograph obtained by 2D pipeline, we find that the ASPSRCs are good enough to use for the calibration. The ASPSRCs have been applied to spectra which were abandoned by LAMOST 2D pipeline due to the lack of standard stars, increasing the number of LAMOST spectra by 52,181 in DR2. Comparing those same targets with SDSS, the relative flux differences between SDSS spectra and that of LAMOST spectra with the ASPSRC method are less than 10%, which underlines that the ASPSRC method is feasible for LAMOST flux calibration.

The "Magnificent Seven" (M7) are a group of radio-quiet Isolated Neutron Stars (INSs) discovered in the soft X-rays through their purely thermal surface emission. Owing to the large inferred magnetic fields ($B\approx 10^{13}$ G), radiation from these sources is expected to be substantially polarised, independently on the mechanism actually responsible for the thermal emission. A large observed polarisation degree is, however, expected only if quantum-electrodynamics (QED) polarisation effects are present in the magnetised vacuum around the star. The detection of a strongly linearly polarised signal would therefore provide the first observational evidence of QED effects in the strong-field regime. While polarisation measurements in the soft X-rays are not feasible yet, optical polarisation measurements are within reach also for quite faint targets, like the M7 which have optical counterparts with magnitudes $\approx 26$--$28$. Here, we report on the measurement of optical linear polarisation for the prototype, and brightest member, of the class, RX\, J1856.5$-$3754 ($V\sim 25.5$), the first ever for one of the M7, obtained with the Very Large Telescope. We measured a polarisation degree $\mathrm{P.D.} =16.43\% \pm5.26\%$ and a polarisation position angle $\mathrm{P.A.}=145\fdg39\pm9\fdg44$, computed east of the North Celestial Meridian. The $\mathrm{P.D.}$ that we derive is large enough to support the presence of vacuum birefringence, as predicted by QED.

We present the results of a series of short radio observations of six ultracool dwarfs made using the upgraded VLA in S (2--4GHz) and C (4--7GHz) bands. LSR J1835+3259 exhibits a 100 percent right-hand circularly polarised burst which shows intense narrowband features with a fast negative frequency drift of about $-30$ MHz $\textrm{s}^{-1}$. They are superimposed on a fainter broadband emission feature with a total duration of about 20 minutes, bandwidth of about 1 GHz, centred at about 3.5 GHz, and a slow positive frequency drift of about 1 MHz $\textrm{s}^{-1}$. This makes it the first such event detected below 4 GHz and the first one exhibiting both positive and negative frequency drifts. Polarised radio emission is also seen in 2MASS J00361617+1821104 and NLTT 33370, while LP 349-25 and TVLM 513-46546 have unpolarised emission and BRI B0021-0214 was not detected. We can reproduce the main characteristics of the burst from LSR J1835+3259 using a model describing the magnetic field of the dwarf as a tilted dipole. We also analyse the origins of the quiescent radio emission and estimate the required parameters of the magnetic field and energetic electrons. Although our results are non-unique, we find a set of models which agree well with the observations.

We calculate the radial profiles of galaxies where the nuclear region is self-gravitating, consisting of self-interacting dark matter (SIDM) with $F$ degrees of freedom. For sufficiently high density this dark matter becomes collisional, regardless of its behaviour on galaxy scales. Our calculations show a spike in the central density profile, with properties determined by the dark matter microphysics, and the densities can reach the `mean density' of a black hole (from dividing the black-hole mass by the volume enclosed by the Schwarzschild radius). For a galaxy halo of given compactness ($\chi=2GM/Rc^2$), certain values for the dark matter entropy yield a dense central object lacking an event horizon. For some soft equations of state of the SIDM (e.g. $F\ge6$), there are multiple horizonless solutions at given compactness. Although light propagates around and through a sphere composed of dark matter, it is gravitationally lensed and redshifted. While some calculations give non-singular solutions, others yield solutions with a central singularity. In all cases the density transitions smoothly from the central body to the dark-matter envelope around it, and to the galaxy's dark matter halo. We propose that pulsar timing observations will be able to distinguish between systems with a centrally dense dark matter sphere (for different equations of state) and conventional galactic nuclei that harbour a supermassive black hole.

The observed polarization properties of thermal radiation from isolated, cooling neutron stars depend on both the emission processes at the surface and the effects of the magnetized vacuum which surrounds the star. Here we investigate the polarized thermal emission from X-ray Dim Isolated Neutron Stars, taking RX J1856.5-3754 as a representative case. The physical conditions of the star outermost layers in these sources is still debated, and so we consider emission from a magnetized atmosphere and a condensed surface, accounting for the effects of vacuum polarization as the radiation propagates in the star magnetosphere. We have found that, for a significant range of viewing geometries, measurement of the phase-averaged polarization fraction and phase-averaged polarization angle at both optical and X-ray wavelengths allow us to determine whether this neutron star has an atmosphere or a condensed surface. Our results may therefore be relevant in view of future developments of soft X-ray polarimeters.

The X-ray and radio flares observed in X-ray binaries and active galactic nuclei (AGN) are attributed to energetic electrons in the plasma ejecta from the accretion flows near the black hole in these systems. It is argued that magnetic reconnection could occur in the coronae above the accretion disk around the black hole, and that this drives plasmoid outflows resembling the solar coronal mass ejection (CME) phenomenon. The X-ray and radio flares are emission from energetic electrons produced in the process. As the emission region is located near the black hole event horizon, the flare emission would be subject to special- and general-relativistic effects. We present calculations of the flaring emission from plasmoids orbiting around a black hole and plasmoid ejecta launched from the inner accretion disk when general-relativistic effects are crucial in determining the observed time-dependent properties of the emission. We consider fully general-relativistic radiative transfer calculations of the emission from evolving ejecta from black hole systems, with proper accounting for differential arrival times of photons emitted from the plasmoids, and determine the emission lightcurves of plasmoids when they are in orbit and when they break free from their magnetic confinement. The implications for interpreting time-dependent spectroscopic observations of flaring emission from accreting black holes are discussed.

High-energy neutrinos and photons are complementary messengers, probing violent astrophysical processes and structural evolution of the Universe. X-ray and neutrino observations jointly constrain conditions in active galactic nuclei (AGN) jets: their baryonic and leptonic contents, and particle production efficiency. Testing two standard neutrino production models for local source Cen A \citep{KT2008,BB2009}, we calculate the high-energy neutrino spectra of single AGN sources and derive the flux of high-energy neutrinos expected for the current epoch. Assuming that accretion determines both X-rays and particle creation, our parametric scaling relations predict neutrino yield in various AGN classes. We derive redshift-dependent number densities of each class, from {\it Chandra} and {\it Swift}/BAT X-ray luminosity functions \citep{SGB2008,ACS2009}. We integrate the neutrino spectrum expected from the cumulative history of AGN (correcting for cosmological and source effects, e.g. jet orientation and beaming). Both emission scenarios yield neutrino fluxes well above limits set by {\it IceCube} (by $\sim 4$--$10^6 \times$ at 1 PeV, depending on the assumed jet models for neutrino production). This implies that: (i) Cen A might not be a typical neutrino source as commonly assumed; (ii) both neutrino production models overestimate the efficiency; (iii) neutrino luminosity scales with accretion power differently among AGN classes and hence does not follow X-ray luminosity universally; (iv) some AGN are neutrino-quiet (e.g. below a power threshold for neutrino production); (v) neutrino and X-ray emission have different duty cycles (e.g. jets alternate between baryonic and leptonic flows); or (vi) some combination of the above.

The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) General Survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects both in the pilot survey and the first year general survey are included in the LAMOST First Data Release (DR1). The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The general survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, of which 1,790,879 spectra have observed signal-to-noise S/N >10. All data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not only spectra, but also three stellar catalogues with measured parameters: AFGK-type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 entries), and M stars (121,522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. Description of the FITS structure of spectral files and parameter catalogues is also provided.

General relativistic magnetohydrodynamic (GRMHD) flows along magnetic fields threading a black hole can be divided into inflow and outflow parts, according to the result of the competition between the black hole gravity and magneto-centrifugal forces along the field line. Here we present the first self-consistent, semi-analytical solution for a cold, Poynting flux-dominated (PFD) GRMHD flow, which passes all four critical (inner and outer, Alfven and fast magnetosonic) points along a parabolic streamline. By assuming that the dominating (electromagnetic) component of the energy flux per flux tube is conserved at the surface where the inflow and outflow are separated, the outflow part of the solution can be constrained by the inflow part. The semi-analytical method can provide fiducial and complementary solutions for GRMHD simulations around the rotating black hole, given that the black hole spin, global streamline, and magnetizaion (i.e., a mass loading at the inflow/outflow separation) are prescribed. For reference, we demonstrate a self-consistent result with the work by McKinney in a quantitative level.

Over the past few years two of the largest and highest fidelity experiments conceived have been approved for construction: Euclid is an ESA M-Class mission that will map three-quarters of the extra galactic sky with Hubble Space Telescope resolution optical and NIR imaging, and NIR spectroscopy, its scientific aims (amongst others) are to create a map of the dark Universe and to determine the nature of dark energy. The Square Kilometre Array (SKA) has similar scientific aims (and others) using radio wavelength observations. The two experiments are synergistic in several respects, both through the scientific objectives and through the control of systematic effects. SKA Phase-1 and Euclid will be commissioned on similar timescales offering an exciting opportunity to exploit synergies between these facilities.

We investigate the black hole (BH) scaling relation in galaxies using a model in which the galaxy halo and central BH are a self-gravitating sphere of dark matter (DM) with an isotropic, adiabatic equation of state. The equipotential where the escape velocity approaches the speed of light defines the horizon of the BH. We find that the BH mass ($m_\bullet$) depends on the DM entropy, when the effective thermal degrees of freedom ($F$) are specified. Relations between BH and galaxy properties arise naturally, with the BH mass and DM velocity dispersion following $m_\bullet\propto\sigma^{F/2}$ (for global mean density set by external cosmogony). Imposing observationally derived constraints on $F$ provides insight into the microphysics of DM. Given that DM velocities and stellar velocities are comparable, the empirical correlation between $m_\bullet$ and stellar velocity dispersions $\sigma_\star$ implies that $7<F<10$. A link between $m_\bullet$ and globular cluster properties also arises because the halo potential binds the globular cluster swarm at large radii. Interestingly, for $F>6$ the dense dark envelope surrounding the BH approaches the mean density of the BH itself, while the outer halo can show a nearly uniform kpc-scale core resembling those observed in galaxies.

Dwarf galaxies and globular clusters may contain intermediate mass black holes ($10^{3}$ to $10^{5}$ solar masses) in their cores. Estimates of ~$10^{3}$ neutron stars in the central parsec of the Galaxy and similar numbers in small elliptical galaxies and globular clusters along with an estimated high probability of ms-pulsar formation in those environments has led many workers to propose the use of ms-pulsar timing to measure the mass and spin of intermediate mass black holes. Models of pulsar motion around a rotating black hole generally assume geodesic motion of a "test" particle in the Kerr metric. These approaches account for well-known effects like de Sitter precession and the Lense-Thirring effect but they do not account for the non-linear effect of the pulsar's stress-energy tensor on the space-time metric. Here we model the motion of a pulsar near a black hole with the Mathisson-Papapetrou-Dixon (MPD) equations. Numerical integration of the MPD equations for black holes of mass 2 X $10^{6}$, $10^{5}$ and $10^{3}$ solar masses shows that the pulsar will not remain in an orbital plane with motion vertical to the plane being largest relative to the orbit's radial dimensions for the lower mass black holes. The pulsar's out of plane motion will lead to timing variations that are up to ~10 microseconds different from those predicted by planar orbit models. Such variations might be detectable in long term observations of millisecond pulsars. If pulsar signals are used to measure the mass and spin of intermediate mass black holes on the basis of dynamical models of the received pulsar signal then the out of plane motion of the pulsar should be part of that model.

(Abridged) We combine NUV, optical and IR imaging of the nearby starburst galaxy M82 to explore the properties of the dust both in the interstellar medium of the galaxy and the dust entrained in the superwind. The three NUV filters of Swift/UVOT enable us to probe in detail the properties of the extinction curve in the region around the 2175A bump. The NUV colour-colour diagram strongly rules out a Calzetti-type law, which can either reflect intrinsic changes in the dust properties or in the star formation history compared to starbursts well represented by such an attenuation law. We emphasize that it is mainly in the NUV region where a standard Milky-Way-type law is preferred over a Calzetti law. The age and dust distribution of the stellar populations is consistent with the scenario of an encounter with M81 in the recent 400 Myr. The radial gradients of the NUV and optical colours in the superwind region support the hypothesis that the emission in the wind cone is driven by scattering from dust grains entrained in the ejecta. The observed wavelength dependence reveals either a grain size distribution $n(a)\propto a^{-2.5}$, where $a$ is the size of the grain, or a flatter distribution with a maximum size cutoff, suggesting that only small grains are entrained in the supernovae-driven wind.

We investigate the effects of a cD galaxy's gravity and AGN heating of the host galaxy cluster. We consider a standard prescription for the hydrodynamics, with the structures determined by mass continuity, momentum and energy conservation equations in spherical symmetry. The cluster comprises a dark matter halo (DM) and ionised X-ray emitting intracluster gas (ICM), which jointly determine the gravitational potential. The cD galaxy is an additive gravitational potential component. The DM assumes a polytropic equation of state (determined by its microphysics), which could be non-radiative self-interacting particles or more exotically interacting particles. The AGN provides distributed heating, counteracting radiative cooling. Stationary density and velocity dispersion profiles are obtained by numerically integrating the hydrodynamic equations with appropriate boundary conditions. The minimum gas temperature in the cluster core is higher when a cD galaxy is present than when it is absent. The solutions also yield a point-like mass concentration exceeding a minimum mass: presumably the AGN's supermassive black hole (SMBH). Consistency with observed SMBH masses constrains the possible DM equations of state. The constraints are looser when a cD galaxy is present. Distributed (AGN) heating alters cluster global properties, and also reduces the lower limits for the central point-mass, for the preferred DM models in which the dark particles have greater heat capacity than point particles. Eluding these constraints would require dominant non-spherical or anisotropic effects (e.g. bulk rotation, non-radial streaming, asymmetric lumps, or a strong magnetic field).

X-ray polarimetry, sometimes alone, and sometimes coupled to spectral and temporal variability measurements and to imaging, allows a wealth of physical phenomena in astrophysics to be studied. X-ray polarimetry investigates the acceleration process, for example, including those typical of magnetic reconnection in solar flares, but also emission in the strong magnetic fields of neutron stars and white dwarfs. It detects scattering in asymmetric structures such as accretion disks and columns, and in the so-called molecular torus and ionization cones. In addition, it allows fundamental physics in regimes of gravity and of magnetic field intensity not accessible to experiments on the Earth to be probed. Finally, models that describe fundamental interactions (e.g. quantum gravity and the extension of the Standard Model) can be tested. We describe in this paper the X-ray Imaging Polarimetry Explorer (XIPE), proposed in June 2012 to the first ESA call for a small mission with a launch in 2017 but not selected. XIPE is composed of two out of the three existing JET-X telescopes with two Gas Pixel Detectors (GPD) filled with a He-DME mixture at their focus and two additional GPDs filled with pressurized Ar-DME facing the sun. The Minimum Detectable Polarization is 14 % at 1 mCrab in 10E5 s (2-10 keV) and 0.6 % for an X10 class flare. The Half Energy Width, measured at PANTER X-ray test facility (MPE, Germany) with JET-X optics is 24 arcsec. XIPE takes advantage of a low-earth equatorial orbit with Malindi as down-link station and of a Mission Operation Center (MOC) at INPE (Brazil).

A covariant scattering kernel is a core component in any self-consistent general relativistic radiative transfer formulation in scattering media. An explicit closed-form expression for a covariant Compton scattering kernel with a good dynamical energy range has unfortunately not been available thus far. Such an expression is essential to obtain numerical solutions to the general relativistic radiative transfer equations in complicated astrophysical settings where strong scattering effects are coupled with highly relativistic flows and steep gravitational gradients. Moreover, this must be performed in an efficient manner. With a self-consistent covariant approach, we have derived a closed-form expression for the Compton scattering kernel for arbitrary energy range. The scattering kernel and its angular moments are expressed in terms of hypergeometric functions, and their derivations are shown explicitly in this paper. We also evaluate the kernel and its moments numerically, assessing various techniques for their calculation. Finally, we demonstrate that our closed-form expression produces the same results as previous calculations, which employ fully numerical computation methods and are applicable only in more restrictive settings.

We investigate the orbital dynamics of hierarchical three-body systems containing a double neutron star system orbiting around a massive black hole. These systems show complex dynamical behaviour because of relativistic coupling between orbits of the neutron stars in the double neutron star system and the orbit of the double neutron star system around the black hole. The orbital motion of the neutron stars around each other drives a loop mass current, which gives rise to gravito-magnetism. Generally, gravito-magnetism involves a rotating black hole. The hierarchical three-body system that we consider is an unusual situation in which black hole rotation is not required. Using a gravito-electromagnetic formulation, we calculate the orbital precession and nutation of the double neutron star system. These precession and nutation effects are observable, thus providing probes to the spacetime around black holes as well as tests of gravito-electromagnetism in the framework of general relativity.

We construct a general relativistic radiative transfer (RT) formulation, applicable to particles with or without mass in astrophysical settings. Derived from first principles, the formulation is manifestly covariant. Absorption and emission, as well as relativistic, geometrical and optical depth effects are treated self-consistently. The RT formulation can handle 3D geometrical settings and structured objects with variations and gradients in the optical depths across the objects and along the line-of-sight. The presence of mass causes the intensity variation along the particle bundle ray to be reduced by an aberration factor. We apply the formulation and demonstrate RT calculations for emission from accretion tori around rotating black holes, considering two cases: idealised optically thick tori that have a sharply defined emission boundary surface, and structured tori that allow variations in the absorption coefficient and emissivity within the tori. Intensity images and emission spectra of these tori are calculated. Geometrical effects, such as lensing-induced self-occulation and multiple-image contribution are far more significant in accretion tori than geometrically thin accretion disks. Optically thin accretion tori emission line profiles are distinguishable from the profiles of lines from optically thick accretion tori and optically thick geometrically thin accretion disks. Line profiles of optically thin accretion tori have a weaker dependence on viewing inclination angle than those of the optically thick accretion tori or accretion disks, especially at high viewing inclination angles. Limb effects are present in accretion tori with finite optical depths. Finally, in accretion flows onto relativistic compact objects, gravitationally induced line resonance can occur. This resonance occurs easily in 3D flows, but not in 2D flows, such as a thin accretion disk around a black hole.

We studied the X-ray timing and spectral variability of the X-ray source Sw J1644+57, a candidate for a tidal disruption event. We have separated the long-term trend (an initial decline followed by a plateau) from the short-term dips in the Swift light-curve. Power spectra and Lomb-Scargle periodograms hint at possible periodic modulation. By using structure function analysis, we have shown that the dips were not random but occurred preferentially at time intervals ~ [2.3, 4.5, 9] x 10^5 s and their higher-order multiples. After the plateau epoch, dipping resumed at ~ [0.7, 1.4] x 10^6 s and their multiples. We have also found that the X-ray spectrum became much softer during each of the early dip, while the spectrum outside the dips became mildly harder in its long-term evolution. We propose that the jet in the system undergoes precession and nutation, which causes the collimated core of the jet briefly to go out of our line of sight. The combined effects of precession and nutation provide a natural explanation for the peculiar patterns of the dips. We interpret the slow hardening of the baseline flux as a transition from an extended, optically thin emission region to a compact, more opaque emission core at the base of the jet.

Radio observations of galaxy clusters show that the intra cluster medium is permeated by \mu G magnetic fields. The origin and evolution of these cosmological magnetic fields is currently not well understood and so their impact on the dynamics of structure formation is not known. Numerical simulations are required to gain a greater understanding and produce predictions for the next generation of radio telescopes. We present the galactic chemodynamics smoothed particle magnetohydrodynamic (SPMHD) code (GCMHD+), which is an MHD implementation for the cosmological smoothed particle hydrodynamic code GCD+. The results of 1, 2 and 3 dimensional tests are presented and the performance of the code is shown relative to the ATHENA grid code. GCMHD+ shows good agreement with the reference solutions produced by ATHENA. The code is then used to simulate the formation of a galaxy cluster with a simple primordial magnetic field embedded in the gas. A homogeneous seed field of 10^-11 G is amplified by a factor of 10^3 during the formation of the cluster. The results show good agreement with the profiles found in other magnetic cluster simulations of similar resolution.

We searched for long period variation in V-band, Ic-band and RXTE X-ray light curves of the High Mass X-ray Binaries (HMXBs) LS 1698 / RX J1037.5-5647, HD 110432 / 1H 1249-637 and HD 161103 / RX J1744.7-2713 in an attempt to discover orbitally induced variation. Data were obtained primarily from the ASAS database and were supplemented by shorter term observations made with the 24- and 40-inch ANU telescopes and one of the robotic PROMPT telescopes. Fourier periodograms suggested the existence of long period variation in the V-band light curves of all three HMXBs, however folding the data at those periods did not reveal convincing periodic variation. At this point we cannot rule out the existence of long term V-band variation for these three sources and hints of longer term variation may be seen in the higher precision PROMPT data. Long term V-band observations, on the order of several years, taken at a frequency of at least once per week and with a precision of 0.01 mag, therefore still have a chance of revealing long term variation in these three HMXBs.

The results of a coordinated space-based photometric and ground-based spectroscopic observing campaign on the enigmatic gamma-ray binary LS 5039 are reported. Sixteen days of observations from the MOST satellite have been combined with high-resolution optical echelle spectroscopy from the 2.3m ANU Telescope in Siding Spring, Australia. These observations were used to measure the orbital parameters of the binary and to study the properties of stellar wind from the O primary. We found that any broad-band optical photometric variability at the orbital period is below the 2 mmag level, supporting the scenario that the orbital eccentricity of the system is near the 0.24 +/- 0.08 value implied by our spectroscopy, which is lower than values previously obtained by other workers. The low amplitude optical variability also implies the component masses are at the higher end of estimates based on the primary's O6.5V((f)) spectral type with a primary mass of ~26 solar masses and a mass for the compact star of at least 1.8 solar masses. The mass loss rate from the O primary was determined to be 3.7E-7 to 4.8E-7 solar masses per year.

We re-assess the XMM-Newton and Swift observations of HLX1, to examine the evidence for its identification as an intermediate-mass black hole. We show that the X-ray spectral and timing properties are equally consistent with an intermediate-mass black hole in a high state, or with a foreground neutron star with a luminosity of about a few times 10^{32} erg/s ~ 10^{-6} L_{Edd}, located at a distance of about 1.5 to 3 kpc. Contrary to previously published results, we find that the X-ray spectral change between the two XMM-Newton observations of 2004 and 2008 (going from power-law dominated to thermal dominated) is not associated with a change in the X-ray luminosity. The thermal component becomes more dominant (and hotter) during the 2009 outburst seen by Swift, but in a way that is consistent with either scenario.

We investigated the time-dependent radiative and dynamical properties of light supersonic jets launched into an external medium, using hydrodynamic simulations and numerical radiative transfer calculations. These involved various structural models for the ambient media, with density profiles appropriate for galactic and extragalactic systems. The radiative transfer formulation took full account of emission, absorption, re-emission, Faraday rotation and Faraday conversion explicitly. High time-resolution intensity maps were generated, frame-by-frame, to track the spatial hydrodynamical and radiative properties of the evolving jets. Intensity light curves were computed via integrating spatially over the emission maps. We apply the models to jets in active galactic nuclei (AGN). From the jet simulations and the time-dependent emission calculations we derived empirical relations for the emission intensity and size for jets at various evolutionary stages. The temporal properties of jet emission are not solely consequences of intrinsic variations in the hydrodynamics and thermal properties of the jet. They also depend on the interaction between the jet and the ambient medium. The interpretation of radio jet morphology therefore needs to take account of environmental factors. Our calculations have also shown that the environmental interactions can affect specific emitting features, such as internal shocks and hotspots. Quantification of the temporal evolution and spatial distribution of these bright features, together with the derived relations between jet size and emission, would enable us to set constraints on the hydrodynamics of AGN and the structure of the ambient medium.

The existence of black holes of masses ~ 10^2-10^5 Msun has important implications for the formation and evolution of star clusters and supermassive black holes. One of the strongest candidates to date is the hyperluminous X-ray source HLX1, possibly located in the S0-a galaxy ESO243-49, but the lack of an identifiable optical counterpart had hampered its interpretation. Using the Magellan telescope, we have discovered an unresolved optical source with R = (23.80 +/- 0.25) mag and V = (24.5 +/- 0.3) mag within HLX1's positional error circle. This implies an average X-ray/optical flux ratio ~ 500. Taking the same distance as ESO243-49, we obtain an intrinsic brightness M_R = (-11.0 +/- 0.3) mag, comparable to that of a massive globular cluster. Alternatively, the optical source is consistent with a main-sequence M star in the Galactic halo (for example an M4.4 star at ~ 2.5 kpc). We also examined the properties of ESO243-49 by combining Swift/UVOT observations with stellar population modelling. We found that the overall emission is dominated by a ~5 Gyr old stellar population, but the UV emission at ~2000 Ang is mostly due to ongoing star-formation at a rate of ~ 0.03 Msun/yr. The UV emission is more intense (at least a 9-sigma enhancement above the mean) North East of the nucleus, in the same quadrant as HLX1. With the combined optical and X-ray measurements, we put constraints on the nature of HLX1. We rule out a foreground star and a background AGN. Two alternative scenarios are still viable. HLX1 could be an accreting intermediate-mass black hole in a star cluster, which may itself be the stripped nucleus of a dwarf galaxy that passed through ESO243-49, an event which might have caused the current episode of star formation. Or, it could be a neutron star in the Galactic halo, accreting from an M4-M5 donor star.

Emission line profiles of tracer molecule H$_2$CO 140 GHz transition from gravitational core collapsing clouds in the dynamic process of forming protostars are calculated, using a simple ray-tracing radiative transfer model. Three self-similar dynamic inside-out core collapse models -- the conventional polytropic model, the empirical hybrid model and the isothermal model -- for star-forming molecular clouds are examined and compared. The isothermal model cannot produce observed asymmetric double-peak molecular line profiles. The conventional polytropic model, which gives flow velocity, mass density and temperature profiles self-consistently, can produce asymmetric double-peak line profiles for a core collapsing cloud. In particular, the blue peak is stronger than the red peak, consistent with a broad class of molecular line profile observations. The relative strengths of the blue and red peaks within a molecular line profile are determined by the cloud temperature gradient. The conventional polytropic model can be utilized to produce molecular line-profile templates, for extracting dynamical information from line spectra of molecular globules undergoing a gravitational core collapse. We show a sample fit using the 140 GHz H$_2$CO emission line from the central region of the molecular globule B335 by our model with $\gamma=1.2$. The calculation of line profiles and fitting processes also offer a scenario to estimate the protostellar mass, the kernel mass accretion rate, and the evolution time scale of a core collapsing cloud. Our model can be readily adapted to other tracer molecules with more or less constant abundances in star-forming clouds.

We investigate the polarization of Compton scattered X-rays from relativistic jets in active galactic nuclei using Monte Carlo simulations. We consider three scenarios: scattering of photons from an accretion disk, scattering of cosmic microwave background (CMB) photons, and synchrotron self-Comptonization (SSC) within the jet. For Comptonization of thermal disk photons or CMB photons the maximum linear polarization attained is slightly over 20% at viewing angles close to 90 degrees. The value decreases with the viewing inclination. For SSC, the maximum value may exceed 80%. The angle dependence is complicated, and it varies with the photon injection sites. Our study demonstrates that X-ray polarization, in addition to multi-wavelength spectra, can distinguish certain models for emission and particle acceleration in relativistic jets.

This article reviews the current works on ultra-compact double-degenerate binaries in the presence of magnetic interaction, in particular, unipolar induction. The orbital dynamics and evolution of compact white-dwarf pairs are discussed in detail. Models and predictions of electron cyclotron masers from unipolar-inductor compact binaries and unipolar-inductor white-dwarf planetary systems are presented. Einstein-Laub effects in compact binaries are briefly discussed.

We present our results of Monte-Carlo simulations of polarized Compton X-rays from magnetic cataclysmic variables, with realistic density, temperature and velocity structures in the accretion flow. Our study has shown that the X-ray linear polarization may reach about 8% for systems with high accretion rates viewed at a high viewing inclination angle. This value is roughly twice the maximum value obtained by previous studies which assumed a cold, static emission region with a uniform density. We also investigate the X-ray polarization properties of ultra-compact double-degenerate binaries for the unipolar-inductor and direct-impact accretor models. Our study has shown negligible X-ray polarization for the unipolar-induction model. However, the direct-impact accretor model may give X-ray polarization levels similar to that predicted for the magnetic cataclysmic variables.

We revisit the discovery outburst of the X-ray transient XTE J1550-564 during which relativistic jets were observed in 1998 September, and review the radio images obtained with the Australian Long Baseline Array, and lightcurves obtained with the Molonglo Observatory Synthesis Telescope and the Australia Telescope Compact Array. Based on HI spectra, we constrain the source distance to between 3.3 and 4.9 kpc. The radio images, taken some two days apart, show the evolution of an ejection event. The apparent separation velocity of the two outermost ejecta is at least 1.3c and may be as large as 1.9c; when relativistic effects are taken into account, the inferred true velocity is >0.8c. The flux densities appear to peak simultaneously during the outburst, with a rather flat (although still optically thin) spectral index of -0.2.

We have calculated the relativistic reflection component of the X-ray spectra of accretion disks in active galactic nuclei (AGN). Our calculations have shown that the spectra can be significantly modified by the motion of the accretion flow and the gravity and rotation of the central black hole. The absorption edges in the spectra suffer severe energy shifts and smearing, and the degree of distortion depends on the system parameters, in particular, the inner radius of the accretion disk and the disk viewing inclination angles. The effects are significant. Fluorescent X-ray emission lines from the inner accretion disk could be powerful diagnostic of space-time distortion and dynamical relativistic effects near the event horizons of accreting black holes. However, improper treatment of the reflection component in fitting the X-ray continuum could give rise to spurious line-like features. These features mimic the true fluorescent emission lines and may mask their relativistic signatures. Fully relativistic models for reflection continua together with the emission lines are needed in order to extract black-hole parameters from the AGN X-ray spectra.

Episodic ejection of plasma blobs have been observed in many black hole systems. While steady, continuous jets are believed to be associated with large-scale open magnetic fields, what causes the episodic ejection of blobs remains unclear. Here by analogy with the coronal mass ejection on the Sun, we propose a magnetohydrodynamical model for episodic ejections from black holes associated with the closed magnetic fields in an accretion flow. Shear and turbulence of the accretion flow deform the field and result in the formation of a flux rope in the disk corona. Energy and helicity are accumulated and stored until a threshold is reached. The system then loses its equilibrium and the flux rope is thrust outward by the magnetic compression force in a catastrophic way. Our calculations show that for parameters appropriate for the black hole in our Galactic center, the plasmoid can attain relativistic speeds in about 35 minutes.

This white paper briefly describes the astrophysics of ultra-compact binaries, with emphasis of the challenges and opportunities in the next decade.

We investigate the polarization properties of Comptonized X-rays from relativistic jets in Active Galactic Nuclei (AGN) using Monte Carlo simulations. We consider three scenarios commonly proposed for the observed X-ray emission in AGN: Compton scattering of blackbody photons emitted from an accretion disk; scattering of cosmic microwave background (CMB) photons; and self-Comptonization of intrinsically polarized synchrotron photons emitted by jet electrons. Our simulations show that for Comptonization of disk and CMB photons, the degree of polarization of the scattered photons increases with the viewing inclination angle with respect to the jet axis. In both cases the maximum linear polarization is approximately 20%. In the case of synchrotron self-Comptonization (SSC), we find that the resulting X-ray polarization depends strongly on the seed synchrotron photon injection site, with typical fractional polarizations of approximately P = 10 - 20% when synchrotron emission is localized near the jet base, while P = 20 - 70% for the case of uniform emission throughout the jet. These results indicate that X-ray polarimetry may be capable of providing unique clues to identify the location of particle acceleration sites in relativistic jets. In particular, if synchrotron photons are emitted quasi-uniformly throughout a jet, then the observed degree of X-ray polarization may be sufficiently different for each of the competing X-ray emission mechanisms (synchrotron, SSC or external Comptonization) to determine which is the dominant process. However, X-ray polarimetry alone is unlikely to be able to distinguish between disk and CMB Comptonization.

We present a convariant formulation for radiative transfer in curved space time and demonstrate some applications in the black-hole systems. We calculate the emission from semi-transparent accretion tori around black holes, for opacity provided by the Fe K lines and for opacity dominated by electron scattering. We also calculate the emission from radiative inefficient accretion flow in black holes with opacity provided by electron-positron annihilation lines. Finally we show shadows cast by accreting black holes with different spins and with different distribution of warm material around them.

We investigate the structure of dynamics of large self-gravitating astrophysical systems using a self-interacting two-component model. We consider two cases, galaxy clusters and cosmic walls, for illustrations. In both cases stability analyses are conducted using perturbative expansion. We have found that waves and solitons are easily generated in these systems. Our analysis shows that dark matter can be Jeans unstable in the very inner regions of galaxy clusters if it has a large internal degree of freedom. The dark matter core may collapse under external perturbations. We also discuss dark-matter oscillations in galaxy clusters and how mode growth and decay lead to heating of intracluster medium. Our analysis shows that dark-matter solitons with both positive and negative amplitudes can be excited in cosmic walls. Resonances in soliton interaction could enhance gas condensation. The co-existence of the two types of dark-matter solitons implies that bright filaments can arise in dark voids.

The high-mass X-ray binary RX J0146.9+6121, with optical counterpart LS I+61 235 (V831 Cas), is an intriguing system on the outskirts of the open cluster NGC 663. It contains the slowest X-ray pulsar known with a pulse period of around 1400s and, primarily from the study of variation in the emission line profile of H alpha, it is known to have a Be decretion disk with a one-armed density wave period of approximately 1240d. Here we present the results of an extensive photometric campaign, supplemented with optical spectroscopy, aimed at measuring short time-scale periodicities. We find three significant periodicities in the photometric data at, in order of statistical significance, 0.34d, 0.67d and 0.10d. We give arguments to support the interpretation that the 0.34d and 0.10d periods could be due to stellar oscillations of the B type primary star and that the 0.67d period is the spin period of the Be star with a spin axis inclination of 23 +10 -8 degrees. We measured a systemic velocity of -37.0 +- 4.3 km/s confirming that LS I+61 235 has a high probability of membership in the young cluster NGC 663 from which the system's age can be estimated as 20-25 Myr. From archival RXTE ASM data we further find "super" X-ray outbursts roughly every 450d. If these super outbursts are caused by the alignment of the compact star with the one-armed decretion disk enhancement, then the orbital period is approximately 330d.

We analyse the radial structure of self-gravitating spheres consisting of multiple interpenetrating fluids, such as the X-ray emitting gas and the dark halo of a galaxy cluster. In these dipolytropic models, adiabatic dark matter sits in equilibrium, while the gas develops a gradual, smooth, quasi-stationary cooling flow. Both affect and respond to the collective gravitational field. We find that all subsonic, radially continuous, steady solutions require a non-zero minimum central point mass. For Mpc-sized halos with 7 to 10 effective degrees of freedom (F2), the minimum central mass is compatible with observations of supermassive black holes. Smaller gas mass influxes enable smaller central masses for wider ranges of F2. The halo comprises a sharp spike around the central mass, embedded within a core of nearly constant density (at 10-10^2.5kpc scales), with outskirts that attenuate and naturally truncate at finite radius (several Mpc). The gas density resembles a broken power law in radius, but the temperature dips and peaks within the dark core. A finite minimum temperature occurs due to gravitational self-warming, without cold mass dropout nor needing regulatory heating. X-ray emission from the intracluster medium mimics a beta-model plus bright compact nucleus. Near-sonic points in the gas flow are bottlenecks to allowed steady solutions; the outermost are at kpc scales. These sites may preferentially develop cold mass dropout during strong perturbations off equilibrium. Within the sonic point, the profile of gas specific entropy is flatter than s r^{1/2}, but this is a shallow ramp and not an isentropic core. When F2 is large, the inner halo spike is only marginally Jeans stable in the central pc, suggesting that a large non-linear disturbance could trigger local dark collapse onto the central object.

The temperatures of electrons and ions in the post-shock accretion region of a magnetic cataclysmic variable (mCV) will be equal at sufficiently high mass flow rates or for sufficiently weak magnetic fields. At lower mass flow rates or in stronger magnetic fields, efficient cyclotron cooling will cool the electrons faster than the electrons can cool the ions and a two-temperature flow will result. Here we investigate the differences in polarized radiation expected from mCV post-shock accretion columns modeled with one- and two-temperature hydrodynamics. In an mCV model with one accretion region, a magnetic field >~30 MG and a specific mass flow rate of ~0.5 g/cm/cm/s, along with a relatively generic geometric orientation of the system, we find that in the ultraviolet either a single linear polarization pulse per binary orbit or two pulses per binary orbit can be expected, depending on the accretion column hydrodynamic structure (one- or two-temperature) modeled. Under conditions where the physical flow is two-temperature, one pulse per orbit is predicted from a single accretion region where a one-temperature model predicts two pulses. The intensity light curves show similar pulse behavior but there is very little difference between the circular polarization predictions of one- and two-temperature models. Such discrepancies indicate that it is important to model some aspect of two-temperature flow in indirect imaging procedures, like Stokes imaging, especially at the edges of extended accretion regions, were the specific mass flow is low, and especially for ultraviolet data.

Compton scattering within the accretion column of magnetic cataclysmic variables (mCVs) can induce a net polarization in the X-ray emission. We investigate this process using Monte Carlo simulations and find that significant polarization can arise as a result of the stratified flow structure in the shock-ionized column. We find that the degree of linear polarization can reach levels up to ~8% for systems with high accretion rates and low white-dwarf masses, when viewed at large inclination angles with respect to the accretion column axis. These levels are substantially higher than previously predicted estimates using an accretion column model with uniform density and temperature. We also find that for systems with a relatively low-mass white dwarf accreting at a high accretion rate, the polarization properties may be insensitive to the magnetic field, since most of the scattering occurs at the base of the accretion column where the density structure is determined mainly by bremsstrahlung cooling instead of cyclotron cooling.

Context: We investigate the distribution of optically high and low states of the system AM Herculis (AM Her). Aims: We determine the state duty cycles, and their relationships with the mass transfer process and binary orbital evolution of the system. Methods: We make use of the photographic plate archive of the Harvard College Observatory between 1890 and 1953 and visual observations collected by the American Association of Variable Star Observers between 1978 and 2005. We determine the statistical probability of the two states, their distribution and recurrence behaviors. Results: We find that the fractional high state duty cycle of the system AM Her is 63%. The data show no preference of timescales on which high or low states occur. However, there appears to be a pattern of long and short duty cycle alternation, suggesting that the state transitions retain memories. We assess models for the high/low states for polars (AM Her type systems). We propose that the white-dwarf magnetic field plays a key role in regulating the mass transfer rate and hence the high/low brightness states, due to variations in the magnetic-field configuration in the system.

A standard disk around an accreting black hole may become effectively optically-thin and scattering dominated in the inner region, for high accretion rates (as already predicted by the Shakura-Sunyaev model). Radiative emission from that region is less efficient than blackbody emission, leading to an increase of the colour temperature in the inner region, by an order of magnitude above the effective temperature. We show that the integrated spectrum has a power-law-like shape in the ~ 1-5 keV band, with a soft excess at lower energies and a downward curvature or break at higher energies, in agreement with the observed spectra of many ultraluminous X-ray sources.

We calculate line emission from relativistic accretion tori around Kerr black holes and investigate how the line profiles depend on the viewing inclination, spin of the central black hole, parameters describing the shape of the tori, and spatial distribution of line emissivity on the torus surface. We also compare the lines with those from thin accretion disks. Our calculations show that lines from tori and lines from thin disks share several common features. In particular, at low and moderate viewing inclination angles they both have asymmetric double-peaked profiles with a tall, sharp blue peak and a shorter red peak which has an extensive red wing. At high viewing inclination angles they both have very broad, asymmetric lines which can be roughly considered single-peaked. Torus and disk lines may show very different red and blue line wings, but the differences are due to the models for relativistic tori and disks having differing inner boundary radii. Self-eclipse and lensing play some role in shaping the torus lines, but they are effective only at high inclination angles. If inner and outer radii of an accretion torus are the same as those of an accretion disk, their line profiles show substantial differences only when inclination angles are close to 90 degrees, and those differences manifest mostly at the central regions of the lines instead of the wings.

Unipolar induction (UI) is a fundamental physical process, which occurs when a conducting body transverses a magnetic field. It has been suggested that UI is operating in RX J0806+15 and RX J1914+24, which are believed to be ultra-compact binaries with orbital periods of 5.4 min and 9.6 min respectively. The UI model predicts that those two sources may be electron cyclotron maser sources at radio wavelengths. Other systems in which UI has been predicted to occur are short period extra-solar terrestrial planets with conducting cores. If UI is present, circularly polarised radio emission is predicted to be emitted. We have searched for this predicted radio emission from short period binaries using the VLA and ATCA. In one epoch we find evidence for a radio source, coincident in position with the optical position of RX J0806+15. Although we cannot completely exclude that this is a chance alignment between the position of RX J0806+15 and an artifact in the data reduction process, the fact that it was detected at a significance level of 5.8 sigma and found to be transient, suggests that it is more likely that RX J0806+15 is a transient radio source. We find an upper limit on the degree of circular polarisation to be ~50%. The inferred brightness temperature exceeds 10^18 K, which is too high for any known incoherent process, but is consistent with maser emission and UI being the driving mechanism. We did not detect radio emission from ES Cet, RX J1914+24 or Gliese 876.

We investigate the hydrodynamics of accretion channelled by a dipolar magnetic field (funnel flows). We consider situations in which the electrons and ions in the flow cannot maintain thermal equilibrium (two-temperature effects) due to strong radiative loss, and determine the effects on the keV X-ray properties of the systems. We apply this model to investigate the accretion shocks of white dwarfs in magnetic cataclysmic variables. We have found that the incorporation of two-temperature effects could harden the keV X-rays. Also, the dipolar model yields harder X-ray spectra than the standard planar model if white dwarf is sufficiently massive (>~1M_sun). When fitting observed keV X-ray spectra of magnetic cataclysmic variables, the inclusion of two-temperature hydrodynamics and a dipolar accretion geometry lowers estimates for white-dwarf masses when compared with masses inferred from models excluding these effects. We find mass reductions <~9% in the most massive cases.

We present a model for the prediction of the optical/infra-red emission from ULXs. In the model, ULXs are binary systems with accretion taking place through Roche lobe overflow. We show that irradiation effects and presence of an accretion disk significantly modify the optical/infrared flux compared to single stars, and also that the system orientation is important. We include additional constraints from the mass transfer rate to constrain the parameters of the donor star, and to a lesser extent the mass of the BH. We apply the model to fit photometric data for several ULX counterparts. We find that most donor stars are of spectral type B and are older and less massive than reported elsewhere, but that no late-type donors are admissable. The degeneracy of the acceptable parameter space will be significantly reduced with observations over a wider spectral range, and if time-resolved data become available.

Compton scattering of X-rays in the bulk flow of the accretion column in magnetic cataclysmic variables (mCVs) can significantly shift photon energies. We present Monte Carlo simulations based on a nonlinear algorithm demonstrating the effects of Compton scattering on the H-like, He-like and neutral Fe K alpha lines produced in the post-shock region of the accretion column. The peak line emissivities of the photons in the post-shock flow are taken into consideration and frequency shifts due to Doppler effects are also included. We find that line profiles are most distorted by Compton scattering effects in strongly magnetized mCVs with a low white dwarf mass and high mass accretion rate and which are viewed at an oblique angle with respect to the accretion column. The resulting line profiles are most sensitive to the inclination angle. We have also explored the effects of modifying the accretion column width and using a realistic emissivity profile. We find that these do not have a significant overall effect on the resulting line profiles. A comparison of our simulated line spectra with high resolution Chandra/HETGS observations of the mCV GK Per indicates that a wing feature redward of the 6.4 keV line may result from Compton recoil near the base of the accretion column.

This is a call for amateur astronomers who have the equipment and experience for producing high quality photometry to contribute to a program of finding periods in the optical light curves of high mass X-ray binaries (HMXB). HMXBs are binary stars in which the lighter star is a neutron star or a black hole and the more massive star is an O type supergiant or a Be type main sequence star. Matter is transferred from the ordinary star to the compact object and X-rays are produced as the the gravitational energy of the accreting gas is converted into light. HMXBs are very bright, many are brighter than 10th magnitude, and so make perfect targets for experienced amateur astronomers with photometry capable CCD equipment coupled with almost any size telescope.

We have constructed a model which describes the optical emission from ultraluminous X-ray sources (ULXs), and have used it to constrain the parameters of seven ULX systems. Our model assumes a binary nature for ULXs, and accounts for optical emission from an X-ray irradiated companion star and accretion disk. We apply our model to six different ULX optical counterparts observed with HST, and one observed with the ESO VLT, and determine the mass, radius and age of the donor stars in these systems. In addition, we obtained constraints for the black hole (BH) mass in some cases. We use the mass accretion rate implied by the X-ray luminosity of these sources as an additional constraint on the donor star, by assuming the mass transfer is driven by the stellar nuclear evolution. We find that in general the donors are older and less massive than previously thought, and are consistent with being of spectral type B. We discuss how these results affect our understanding of the evolution and history of ULXs. Where we can constrain the BH masses, we find them to be consistent with stellar mass BHs or intermediate mass BHs of order 100 solar masses. We make predictions for future observations of optical/infrared ULX counterparts, calculating binary periods for different BH masses in each of the seven sources.

We investigate jet formation in black-hole systems using 3-D General Relativistic Particle-In-Cell (GRPIC) and 3-D GRMHD simulations. GRPIC simulations, which allow charge separations in a collisionless plasma, do not need to invoke the frozen condition as in GRMHD simulations. 3-D GRPIC simulations show that jets are launched from Kerr black holes as in 3-D GRMHD simulations, but jet formation in the two cases may not be identical. Comparative study of black hole systems with GRPIC and GRMHD simulations with the inclusion of radiate transfer will further clarify the mechanisms that drive the evolution of disk-jet systems.

We present the results of some recent research on AM CVn systems. We present: X-ray/UV observations made using XMM-Newton; the X-ray grating spectrum of RX J1914+24; preliminary results of a search for radio emission from AM CVn binaries, and discuss the strategy and first results of the RATS project, whose main aim is to discover AM CVn systems.

We report observations with the Chandra X-ray Observatory of the single, cool, magnetic white dwarf GD 356. For consistent comparison with other X-ray observations of single white dwarfs, we also re-analyzed archival ROSAT data for GD 356 (GJ 1205), G 99-47 (GR 290 = V1201 Ori), GD 90, G 195-19 (EG250 = GJ 339.1), and WD 2316+123 and archival Chandra data for LHS 1038 (GJ 1004) and GD 358 (V777 Her). Our Chandra observation detected no X rays from GD 356, setting the most restrictive upper limit to the X-ray luminosity from any cool white dwarf -- L_{X} < 6.0 x 10^{25} ergs/s, at 99.7% confidence, for a 1-keV thermal-bremsstrahlung spectrum. The corresponding limit to the electron density is n_{0} < 4.4 x 10^{11} cm^{-3}. Our re-analysis of the archival data confirmed the non-detections reported by the original investigators. We discuss the implications of our and prior observations on models for coronal emission from white dwarfs. For magnetic white dwarfs, we emphasize the more stringent constraints imposed by cyclotron radiation. In addition, we describe (in an appendix) a statistical methodology for detecting a source and for constraining the strength of a source, which applies even when the number of source or background events is small.

ROSAT observations found no convincing evidence for X-ray sources located in local dwarf spheroidal galaxies (dSph). Now with more sensitive instruments on board Chandra and XMM-Newton we can reach fainter luminosity levels. We report on an observation of the Sagittarius (Sgr) dSph made using Chandra and an observation of the Carina (Car) dSph made using XMM-Newton. Our observations are sensitive to sources with X-ray luminosities in the 0.1-10 keV band of \~1x10^{32 ergs/s and 3x10^34 ergs/s for the Sgr and Car fields respectively. We have identified a total of 80 sources in the Sgr field and 53 sources in the Car field. Although the source numbers are roughly consistent with the expected number of background AGN, we found a small fraction of X-ray sources which were soft and could be located in the host dSph. Follow-up optical/IR observations may help to identify their optical counterparts and hence determine their nature.

The ratios of hydrogen Balmer emission line intensities in cataclysmic variables are signatures of the physical processes that produce them. To quantify those signatures relative to classifications of cataclysmic variable types, we applied the multivariate statistical analysis methods of principal components analysis and discriminant function analysis to the spectroscopic emission data set of Williams (1983). The two analysis methods reveal two different sources of variation in the ratios of the emission lines. The source of variation seen in the principal components analysis was shown to be correlated with the binary orbital period. The source of variation seen in the discriminant function analysis was shown to be correlated with the equivalent width of the H$\beta$ line. Comparison of the data scatterplot with scatterplots of theoretical models shows that Balmer line emission from T CrB systems is consistent with the photoionization of a surrounding nebula. Otherwise, models that we considered do not reproduce the wide range of Balmer decrements, including "inverted" decrements, seen in the data.

We have carried out a Chandra observation of the globular cluster M54. We detected 7 sources located within the half-mass radius of M54, at a flux limit of 1.5x10^-15 erg/s/cm in the 0.3-8 keV energy band. The spatial distribution and the colour/spectral properties of the 7 sources suggest that they are likely to be cataclysmic variables or LMXBs in the globular cluster. M54 shows the largest number of X-ray sources with luminosities greater than 10^32 ergs/s compared to other globular clusters observed using Chandra and XMM-Newton. We searched for a correlation between the number of sources above this luminosity level with globular cluster parameters. We found evidence that the number of sources peaks at a King concentration parameter c~1.7-1.9, with globular clusters which are core-collapsed or have low-c values having a smaller number of sources. We speculate on possible reasons for this.

We combine 9 ROSAT, 9 Chandra, and 2 XMM-Newton observations of the Cen~A galaxy to obtain the X-ray light curve of 1RXH J132519.8-430312 (=CXOU J132519.9$-$430317) spanning 1990 to 2003. The source reached a peak 0.1-2.4 keV flux F_X>10^{-12} ergs cm^{-2} s^{-1} during a 10~day span in 1995 July. The inferred peak isotropic luminosity of the source therefore exceeded 3 10^{39} ergs s^{-1}, which places the source in the class of ultra-luminous X-ray sources. Coherent pulsations at 13.264 Hz are detected during a second bright episode (F_X >3 times 10^{-13} ergs cm^{-2} s^{-1}) in 1999 December. The source is detected and varies significantly within three additional observations but is below the detection threshold in 7 observations. The X-ray spectrum in 1999 December is best described as a cut-off power law or a disk-blackbody (multi-colored disk). We also detect an optical source, m_F555W ~ 24.1 mag, within the Chandra error circle of 1RXH J132519.8-430312 in HST images taken 195~days before the nearest X-ray observation. The optical brightness of this source is consistent with a late O or early B star at the distance of Cen A. If the optical source is the counterpart, then the X-ray and optical behavior of 1RXH J132519.8-430312 are similar to the transient Be/X-ray pulsar A 0538-66.

We have performed a fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity nearly 0.3c) is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation in the region around r = 3r_S. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving torsional Alfven waves (TAWs). This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk due to the lack of streaming material from an accompanying star.

Before the official first-light images, the Chandra X-Ray Observatory obtained an X-ray image of the field to which its focal plane was first exposed. We describe this historic observation and report our study of the first Chandra field. Chandra's Advanced CCD Imaging Spectrometer (ACIS) detected 15 X-ray sources, the brightest being dubbed ``Leon X-1'' to honor the Chandra Telescope Scientist, Leon Van Speybroeck. Based upon our analysis of the X-ray data and spectroscopy at the European Southern Observatory (ESO; La Silla, Chile), we find that Leon X-1 is a Type-1 (unobscured) active galactic nucleus (AGN) at a redshift z=0.3207. Leon X-1 exhibits strong Fe II emission and a broad-line Balmer decrement that is unusually flat for an AGN. Within the context of the Eigenvector-1 correlation space, these properties suggest that Leon X-1 may be a massive (> 10**{9} solar mass) black hole, accreting at a rate approaching its Eddington limit.

We have constructed a model to describe the optical emission from ultra-luminous X-ray sources (ULXs). We assume a binary model with a black hole accreting matter from a Roche lobe filling companion star. We consider the effects of radiative transport and radiative equilibrium in the irradiated surfaces of both the star and a thin accretion disk. We have developed this model as a tool with which to positively identify the optical counterparts of ULXs, and subsequently derive parameters such as the black hole mass and the luminosity class and spectral type of the counterpart. We examine the dependence of the optical emission on these and other variables. We extend our model to examine the magnitude variation at infrared wavelengths, and we find that observations at these wavelengths may have more diagnostic power than in the optical. We apply our model to existing HST observations of the candidates for the optical counterpart of ULX X-7 in NGC 4559. All candidates could be consistent with an irradiated star alone, but we find that a number of them are too faint to fit with an irradiated star and disk together. Were one of these the optical counterpart to X-7, it would display a significant temporal variation.

We use a two-temperature hydrodynamical formulation to determine the temperature and density structures of the post-shock accretion flows in magnetic cataclysmic variables (mCVs) and calculate the corresponding X-ray spectra. The effects of two-temperature flows are significant for systems with a massive white dwarf and a strong white-dwarf magnetic field. Our calculations show that two-temperature flows predict harder keV spectra than one-temperature flows for the same white-dwarf mass and magnetic field. This result is insensitive to whether the electrons and ions have equal temperature at the shock but depends on the electron-ion exchange rate, relative to the rate of radiative loss along the flow. White-dwarf masses obtained by fitting the X-ray spectra of mCVs using hydrodynamic models including the two-temperature effects will be lower than those obtained using single-temperature models. The bias is more severe for systems with a massive white dwarf.

Field-channelled accretion flows occur in a variety of astrophysical objects, including T Tauri stars,magnetic cataclysmic variables and X-ray pulsars. We consider a curvilinear coordinate system and derive a general hydrodynamic formulation for accretion onto stellar objects confined by a stellar dipole magnetic field. The hydrodynamic equations are solved to determine the velocity, density and temperature profiles of the flow. We use accreting magnetic white-dwarf stars as an illustrative example of astrophysical applications. Our calculations show that the compressional heating due to the field geometry is as important as radiative cooling and gravity in determining the structure of the post-shock flow in accreting white-dwarf stars. The generalisation of the formulation to accretion flows channelled by higher-order fields and the applications to other astrophysical systems are discussed.

The Deep Extragalactic Exploratory Probe (DEEP) is a multi-phase research program dedicated to the study of the formation and evolution of galaxies and of large scale structure in the distant Universe. This paper describes the first five-year phase, denoted DEEP1. A series of ten DEEP1 papers will discuss a range of scientific topics (e.g., the study of photometric and spectral properties of a general distant galaxy survey, the evolution observed in galaxy populations of varied morphologies). The observational basis for these studies is the Groth Survey Strip field, a 127 square arcminute region which has been observed with the Hubble Space Telescope in both broad I-band and V-band optical filters and with the Low Resolution Imaging Spectrograph on the Keck Telescopes. Catalogs of photometric and structural parameters have been constructed for 11,547 galaxies and stars at magnitudes brighter than 29, and spectroscopy has been conducted for a magnitude-color weighted subsample of 818 objects. We evaluate three independent techniques for constructing an imaging catalog for the field from the HST data, and discuss the depth and sampling of the resultant catalogs. The selection of the spectroscopic subsample is discussed, and we describe the multifaceted approach taken to prioritizing objects of interest for a variety of scientific subprograms. A series of Monte Carlo simulations then demonstrates that the spectroscopic subsample can be adequately modeled as a simple function of magnitude and color cuts in the imaging catalog.

We report on the analysis of 100 ks INTEGRAL observations of the Galactic microquasar GRS 1915+105. We focus on INTEGRAL Revolution number 48 when the source was found to exhibit a new type of variability as preliminarily reported in Hannikainen et al. (2003). The variability pattern, which we name $\xi$, is characterized by a pulsing behaviour, consisting of a main pulse and a shorter, softer, and smaller amplitude precursor pulse, on a timescale of 5 minutes in the JEM-X 3-35 keV lightcurve. We also present simultaneous RXTE data. From a study of the individual RXTE/PCA pulse profiles we find that the rising phase is shorter and harder than the declining phase, which is opposite to what has been observed in other otherwise similar variability classes in this source. The position in the colour-colour diagram throughout the revolution corresponds to State A (Belloni et al. 2000) but not to any previously known variability class. We separated the INTEGRAL data into two subsets covering the maxima and minima of the pulses and fitted the resulting two broadband spectra with a hybrid thermal--non-thermal Comptonization model. The fits show the source to be in a soft state characterized by a strong disc component below ~6 keV and Comptonization by both thermal and non-thermal electrons at higher energies.

The E3 giant elliptical galaxy NGC5018 was observed with the Chandra X-ray Observatory on 14 April 2001. Results of analysis of these X-ray data as well as of complementary optical, infrared, and radio data are reported. Seven X-ray point sources, including the nucleus, were detected.If they are intrinsic to NGC5018, then all six non-nuclear sources have luminosities exceeding E+39 ergs/s in the 0.5-8.0 keV energy band; placing them in the class of Ultra-luminous X-ray sources. Comparison of X-ray source positions to archival HST/WFPC2 images reveal four of the six non-nuclear sources are spatially-coincident with bright, absolute-magnitude (V) around -8.6 mag, objects. These four objects have optical magnitudes and (V-I) colors consistent with globular clusters in NGC5018 but they also have X-ray-to-optical flux ratios consistent with background active galactic nuclei. Strong, unpolarized, radio emission has been detected from one of the optically-bright counterparts. Another optically-bright counterpart was observed to vary by one mag in optical observations taken 28 July 1997 and 04 Feb 1999. Extended X-ray emission is detected within a 15"radius of the galaxy center at a luminosity of E+40 ergs/s in the X-ray band. Its thermal X-ray spectrum (kT around 0.4 keV) and its spatial coincidence with strong H_alpha emission are consistent with a hot gas origin. The nucleus itself may be a weak X-ray source, around E+39 ergs/s, that displays a radio spectrum typical of AGN.

We present XMM-Newton observations of the 569 sec period system RX J1914+24 (V407 Vul). This period is believed to represent the binary orbital period making it an ultra-compact binary system. By comparing the phase of the rise to maximum X-ray flux at various epochs (this includes observations made using ROSAT, ASCA and Chandra) we find that the system is spinning up at a rate of 3.17+/-0.07x10^{-12} s/s. We find that the spectra softens as the X-ray flux declines towards the off-phase of the 569 sec period. Further, the spectra are best fitted by an absorbed blackbody component together with a broad emission feature around 0.59keV. This emission feature is most prominent at the peak of the on-phase. We speculate on its origin.

We examine the XMM X-ray spectrum of the LINER-AGN NGC 7213, which is best fit with a power law, K-alpha emission lines from Fe I, Fe XXV and Fe XXVI and a soft X-ray collisionally ionised thermal plasma with kT=0.18 +0.03/-0.01 keV. We find a luminosity of 7x10^(-4) L_Edd, and a lack of soft X-ray excess emission, suggesting a truncated accretion disc. NGC 7213 has intermediate X-ray spectral properties, between those of the weak AGN found in the LINER M81 and higher luminosity Seyfert galaxies. This supports the notion of a continuous sequence of X-ray properties from the Galactic Centre through LINER galaxies to Seyferts, likely determined by the amount of material available for accretion in the central regions.

We present a sample of over 50 luminous field bulges (including ellipticals) found in the Groth Strip Survey (GSS), with 0.73< z < 1.04 and with bulge magnitudes I <= 23. The exponential disk light is removed via decomposition of HST images using GIM2D. We find that 85% of these bulges are nearly as red as local E/S0's and have a shallow slope and a small color dispersion in the color-luminosity relation, suggesting roughly coeval formation. The surface brightnesses of these bulges are about 1 mag higher than local bulges. These results are explained adopting a "drizzling" scenario where a metal-rich early formation is later polluted by small amounts of additional star formation. Almost all disks have the same or bluer colors than their accompanying bulges, regardless of the bulge-disk ratio and bulge luminosity, as expected from semi-analytic hierarchical galaxy formation models. We present evidence that the few blue bulge candidates are not likely to be genuine blue ellipticals or bulges. Our deeper, more extensive, and less disk-contaminated observations challenge prior claims that 30% to 50% of field bulges or ellipticals are in a blue, star-forming phase at z < 1. We conclude that field bulges and ellipticals at z ~ 1, like luminous early- type cluster galaxies at the same redshift, are already dominated by metal-rich, old stellar populations that have been fading from a formation epoch earlier than z ~ 1.5. (abridged)

Compton scattering in the bulk accretion flow of the accretion column in magnetic cataclysmic variables (mCVs) can significantly shift photon energies in the X-ray emission lines resulting from accretion shocks. In particular, Compton recoil can potentially broaden the 6.7 and 6.97 keV Fe K$\alpha$ emission lines produced in the post-shock region, and contaminate the fluorescent 6.4 keV neutral Fe K$\alpha$ line reflected off the white dwarf surface. We present nonlinear Monte Carlo simulations demonstrating these effects, and we discuss the interpretation of observed Fe K$\alpha$ linewidths in mCVs in light of these new results. The implications for other accreting compact objects are also discussed.

The Deep Extragalactic Evolutionary Probe (DEEP) is a series of spectroscopic surveys of faint galaxies, targeted at the properties and clustering of galaxies at redshifts z ~ 1. We present the redshift catalog of the DEEP 1 GSS pilot phase of this project, a Keck/LRIS survey in the HST/WFPC2 Groth Survey Strip. The redshift catalog and data, including reduced spectra, are publicly available through a Web-accessible database. The catalog contains 658 secure galaxy redshifts with a median z=0.65, and shows large-scale structure walls to z = 1. We find a bimodal distribution in the galaxy color-magnitude diagram which persists to z = 1. A similar color division has been seen locally by the SDSS and to z ~ 1 by COMBO-17. For red galaxies, we find a reddening of only 0.11 mag from z ~ 0.8 to now, about half the color evolution measured by COMBO-17. We measure structural properties of the galaxies from the HST imaging, and find that the color division corresponds generally to a structural division. Most red galaxies, ~ 75%, are centrally concentrated, with a red bulge or spheroid, while blue galaxies usually have exponential profiles. However, there are two subclasses of red galaxies that are not bulge-dominated: edge-on disks and a second category which we term diffuse red galaxies (DIFRGs). The distant edge-on disks are similar in appearance and frequency to those at low redshift, but analogs of DIFRGs are rare among local red galaxies. DIFRGs have significant emission lines, indicating that they are reddened mainly by dust rather than age. The DIFRGs in our sample are all at z>0.64, suggesting that DIFRGs are more prevalent at high redshifts; they may be related to the dusty or irregular extremely red objects (EROs) beyond z>1.2 that have been found in deep K-selected surveys. (abridged)

We present an XMM-Newton observation of the Seyfert-LINER galaxy NGC 7213. The RGS soft X-ray spectrum is well fitted with a power law plus soft X-ray collisionally ionised thermal plasma (kT = 0.18 +0.03/-0.01 keV). We confirm the presence of Fe I, XXV and XXVI K-alpha emission in the EPIC spectrum and set tighter constraints on their equivalent widths of 82 +10/-13, 24 +9/-11 and 24 +10/-13 eV respectively. We compare the observed properties together with the inferred mass accretion rate of NGC 7213, to those of other Seyfert and LINER galaxies. We find that NGC 7213 has intermediate X-ray spectral properties lying between those of the weak AGN found in the LINER M81 and higher luminosity Seyfert galaxies. There appears to be a continuous sequence of X-ray properties from the Galactic Centre through LINER galaxies to Seyferts, likely determined by the amount of material available for accretion in the central regions.

Extensive radio data of two Algol systems and two RS CVn binaries were re-analyzed. We found evidence for a new periodicity that we interpret as a superhump in Algol, in which it may have been expected according to its semi-detached nature and low binary mass ratio. This is the first detection of the superhump phenomenon in the radio and the first observation of superhumps in Algol systems. According to our result, the accretion disk in Algol precesses in spite of its non-Keplerian nature and therefore this phenomenon is not restricted to the classical Keplerian disks in compact binaries.We propose that in Algol systems with short orbital periods, the disk is magnetically active as well as the secondary star. The magnetic field in the disk originates from amplification of the seed field in the magnetized material transferred from the secondary. The disk and stellar fields interact with each other, with reconnection of the field lines causing flares and particle acceleration. Relativistic particles are trapped in the field and directed toward the polar regions of the secondary star because of the dipole structure of its magnetic field. Our proposed model for the magnetic activity in Algol systems provides a simple explanation to the observed properties of Algol in the radio wavelengths, and to the presence of quiescent gyrosynchrotron emission near the polar region of the secondary star, where electrons are difficult to be confined if the field lines are open as in normal single magnetic stars. We propose that the superhump variation in the radio is generated by enhanced reconnection when the elongated side of the elliptic disk is the closest to the cool star. This leads to flares and enhancement in particle acceleration and is manifested as stronger gyrosynchrotron radiation.

We have studied the candidate optical counterparts and the stellar population in the star-forming complex around a bright ULX in NGC4559, using HST/WFPC2, XMM-Newton/OM, and ground-based data. We find that the ULX is located near a small group of OB stars. The brightest point source in the Chandra error circle is consistent with a single blue supergiant of mass ~ 20 M_sun and age ~ 10 Myr. A few other stars are resolved inside the error circle: mostly blue and red supergiants with masses ~ 10-15 M_sun and ages ~ 20 Myr. This is consistent with the interpretation of this ULX as a black hole (BH) accreting from a high-mass donor star in its supergiant phase, via Roche-lobe overflow. The observed optical colors and the blue-to-red supergiant ratio suggest a low metal abundance: 0.2 <~ Z/Z_sun <~ 0.4 (Padua tracks), or 0.05 <~ Z/Z_sun <~ 0.2 (Geneva tracks). The age of the star-forming complex is <~ 30 Myr. H-alpha images show that this region has a ring-like appearance. We propose that it is an expanding wave of star formation, triggered by an initial density perturbation, in a region where the gas was only marginally stable to gravitational collapse. A possible trigger was the collision with a satellite dwarf galaxy, visible a few arcsec north-west of the complex, going through the gas-rich outer disk of NGC4559. The X-ray data favour a BH more massive (M > 50 M_sun) than typical Milky Way BH candidates. The optical data favour a ``young'' BH originating in the recent episode of massive star formation; however, they also rule out an association with young massive star clusters. We speculate that other mechanisms may lead to the formation of relatively massive BHs (~ 50-100 M_sun) from stellar evolution processes in low-metallicity environments, or when star formation is triggered by galactic collisions.

We present a general formulation for ray-tracing calculations in curved space-time. The formulation takes full account of relativistic effects in the photon transport and the relative motions of the emitters and the light-of-sight absorbing material. We apply the formulation to calculate the emission from accretion disks and tori around rotating black holes.

One hundred fifty-four discrete non-nuclear Ultra-Luminous X-ray (ULX) sources, with intrinsic X-ray luminosities >1e39 ergs/s, are identified in 82 galaxies observed with Chandra's Advanced CCD Imaging Spectrometer. Source positions, X-ray luminosities, and spectral and timing characteristics are tabulated. Statistical comparisons between these X-ray properties and those of the weaker discrete sources in the same fields are made. Sources above ~1e38 erg/s display similar spatial, spectral, color, and variability distributions. In particular, there is no compelling evidence in the sample for a new and distinct class of X-ray object such as the intermediate-mass black holes. Comparison of the cumulative X-ray luminosity functions of the ULXs to Chandra Deep Field results suggests ~25% of the sources may be background objects including 14% of the ULX candidates in the sample of spiral galaxies and 44% of those in elliptical galaxies implying the elliptical galaxy ULX population is severely compromised by background active galactic nuclei. Correlations with host galaxy properties confirm the number and total X-ray luminosity of the ULXs are associated with recent star formation and with galaxy merging and interactions. The preponderance of ULXs in star-forming galaxies as well as their similarities to less-luminous sources suggest they originate in a young but short-lived population such as the high-mass X-ray binaries with a smaller contribution (based on spectral slope) from recent supernovae. The number of ULXs in elliptical galaxies scales with host galaxy mass and can be explained most simply as the high-luminosity end of the low-mass X-ray binary population.

GRS 1915+105 is being observed as part of an Open Time monitoring program with INTEGRAL. Three out of six observations from the monitoring program are presented here, in addition to data obtained through an exchange with other observers. We also present simultaneous RXTE observations of GRS 1915+105. During INTEGRAL Revolution 48 (2003 March 6) the source was observed to be in a highly variable state, characterized by 5-minute quasi-periodic oscillations. During these oscillations, the rise is faster than the decline, and is harder. This particular type of variability has never been observed before. During subsequent INTEGRAL revolutions (2003 March-May), the source was in a steady or ''plateau'' state (also known as class $\chi$ according to Belloni et al. 2000). Here we discuss both the temporal and spectral characteristics of the source during the first three months of observations. The source was clearly detected with all three gamma-ray and X-ray instruments onboard INTEGRAL.

We present spectral and timing results from a long (130 ks) XMM-NEWTON EPIC observation of the nucleus of the Seyfert/LINER galaxy M 81. During the observation the X-ray flux varied by 20%, but there was no significant change in spectral shape. The 2-10 keV spectrum is well described by a power law continuum and three narrow Fe K emission lines at 6.4, 6.7 and 6.96 keV. The three emission lines have equivalent widths of 39, 47, and 37 eV respectively. The ratios of the three lines are thus more similar to those observed from the Galactic Centre region than to those typically observed from Seyfert galaxies. The high ionization lines most likely originate either from photoionized gas within 0.1 pc of the nucleus of M 81, or from a non-thermal distribution of cosmic-ray electrons interacting with the 0.2-0.6 keV thermal plasma which is found in the bulge of M 81.

We have made a series of snap-shot observations of 37 polars using XMM-Newton. We found that 16 of these systems were in a low, or much reduced, accretion state. Of those, 6 were not detected in X-rays. This suggests that in any survey of polars, around half will be in a low accretion state. We tested if there was a bias towards certain orbital periods: this is not the case. Of the 10 systems which were detected at low, but significant rates in X-rays, 8 showed significant variability in their X-ray light curves. This implies that non-uniform accretion still takes place during low accretion epochs. The bolometric luminosity of these systems is ~10^30 ergs s, two orders of magnitude less than for systems in a high accretion state. The X-ray spectra show no evidence of a distinct soft X-ray component. However, the X-ray and UV data imply that such a low temperature component exists: its temperature is low enough for its flux distribution to move outside the bandpass of the X-ray instruments.

Collimated relativistic outflows, or jets, are amongst the most energetic and relativistic phenomena in the Universe. They are associated with supermassive black holes in distant active galactic nuclei (AGN), accreting black holes and neutron stars in binary systems and are believed to be responsible for gamma-ray bursts. The physics of these jets,however, remains something of a mystery in that their bulk velocities, compositions and energetics remain poorly determined. Here we report the discovery of an ultra-relativistic outflow from a binary accreting neutron star accreting gas within a binary system. The velocity of the outflow is comparable to the fastest-moving flows observed from active galactic nuclei, and its strength is modulated by the rate of accretion of material onto the neutron star. Shocks are energized further downstream in the flow, which are themselves moving with mildly relativistic bulk velocities and are the sites of the observed synchrotron emission from the jet. We conclude that the generation of such highly relativistic outflows does not require properties unique to black holes, such as an event horizon.

We have examined X-ray and optical observations of two ultra-luminous X-ray sources, X7 and X10 in NGC4559, using XMM-Newton, Chandra and HST. The UV/X-ray luminosity of X7 exceeds 2.1e40 erg/s in the XMM-Newton observation, and that of X10 is >1.3e40 erg/s. X7 has both thermal and power-law spectral components. The characteristic temperature of the thermal component is 0.12 keV. The power-law components in the two sources both have slopes with photon index \~2.1. A timing analysis of X7 indicates a break frequency at 28 mHz in the power spectrum, while that for X10 is consistent with an unbroken power law. The luminosity of the blackbody component in the X-ray spectrum of X7 and the nature of its time variability provides evidence that this object is an intermediate mass black hole accreting at sub-Eddington rates, but other scenarios which require high advection efficiencies from a hollowed-out disk might be possible. The emission from X10 can be characterised by a single power-law. This source can be interpreted either as an intermediate mass black hole, or as a stellar-mass black hole with relativistically-beamed Comptonised emission. There are four optical counterparts in the error circle of X7. No counterparts are evident in the error circle for X10.

The spectrum and light curve of the bright X-ray source CG X-1 in the field of the Circinus galaxy are re-examined. Previous analyses have concluded that the source is an accreting black hole of about 50 solar masses although it was noted that the light curve resembles that of an AM Her system. Here we show that the short period and an assumed main sequence companion constrain the mass of the companion to less than one solar mass. Further a possible eclipse seen during one of the Chandra observations and a subsequent XMM-Newton observation constrains the mass of the compact object to less than about 60 solar masses. If such a system lies in the Circinus galaxy, then the accreting object must either radiate anisotropically or strongly violate the Eddington limit. Even if the emission is beamed, then the companion star which intercepts this flux during eclipse will be driven out of thermal equilibrium and evaporate within about 1000 years. We find that the observations cannot rule out an AM Her system in the Milky Way and that such a system can account for the variations seen in the light curve.

We present data from the first of six monitoring Open Time observations of GRS 1915+105 undertaken with the orbiting INTEGRAL satellite. The source was clearly detected with all three X-ray and gamma-ray instruments on board. GRS 1915+105 was in a highly variable state, as demonstrated by the JEM X-2 and ISGRI lightcurves. These and simultaneous RXTE/PCA lightcurves point to a novel type of variability pattern in the source. In addition, we fit the combined JEM X-2 and ISGRI spectrum between 3-300 keV with a disk blackbody + powerlaw model leading to typical parameter values found earlier at similar luminosity levels. A new transient, IGR J19140+098, was discovered during the present observation.

We present soft X-ray spectroscopy of the black hole binary LMC X-3 from the XMM-Newton Reflection Grating Spectrometer. The observations span the full range of spectral states seen in LMC X-3. The spectra are completely dominated by continuum emission, and the neutral absorbing column measured from the O I edge (NH = 3.8 +0.8 -0.7 10^20 cm^-2) is consistent with the Galactic interstellar column density towards LMC X-3. We find no evidence for variability of the neutral absorbing column. We also constrain the ionized column density using the upper limits to the equivalent width of the O II - O VIII K-shell resonance lines: we find that the equivalent hydrogen column density of gas in which O is partially ionized is < 8 10^20 cm^-2. From this upper limit we can rule out a line driven stellar wind as the power source for the X-ray emission of LMC X-3 except when it is faint. At wavelengths longward of the peak emission the spectral shape is well described by a multi-temperature disk-blackbody spectrum; the powerlaw component which dominates at shorter wavelengths does not continue longward of the disk-blackbody peak. This implies that the multi-temperature disk-blackbody component supplies the seed photons which are Compton upscattered in the hot corona, consistent with the standard paradigm for black hole X-ray binary spectral states.

We have identified 127 discrete sources in a Chandra observation of M83, with a detection limit of ~ 3 x 10^{36} erg/s in the 0.3-8.0 keV band. We discuss the individual properties of about 20 bright sources with L_x >~ 10^{38} erg/s, and the statistical properties of the whole sample. About 1/3 of the bright sources show X-ray spectra with a bb component at kT <~ 1 keV, plus a p-law component with Gamma ~ 2.5, typical of X-ray binaries (XRBs) in a soft state; another third have p-law spectra with Gamma ~ 1.5, consistent with XRBs in a hard state. Two bright sources show emission lines on a hard p-law continuum, and are probably XRBs surrounded by a photo-ionized nebula or stellar wind. Among the other bright sources, we also identified two SNR candidates, with optically-thin thermal plasma spectra at kT ~ 0.5 keV. The two brightest supersoft sources have T_bb ~ 70 eV and L_x ~ 10^{38} erg/s. Two candidate X-ray pulsars are detected with periods ~ 200 s. Another source corresponds to the core of a background FRII radio galaxy. The discrete sources can be divided into three groups, based on their spatial, color and luminosity distributions. The first group comprises supersoft sources with no detected emission above 1 keV. The second group consists of soft sources with little or no detected emission above 2 keV. They are strongly correlated with H-alpha emission in the spiral arms and starburst nucleus, tracing a young population. Their relative abundance depends on the current level of star-formation. Most of them are likely to be SNRs. The sources in the third group are mostly XRBs (in a soft or hard state), reaching higher X-ray luminosities than sources in the other two groups. Being a mixture of old low-mass and young high-mass systems, the whole group appears to be of intermediate age when correlated with H-alpha.

RX J0806+15 has recently been identified as the binary system with the shortest known orbital period. We present a series of observations of RX J0806+15 including new optical observations taken one month apart. Using these observations and archival data we find that the period of this system is decreasing over time. Our measurements imply f_dot = 6.11x10^-16 Hz/s, which is in agreement with a rate expected from the gravitational radiation for two white dwarfs orbiting at a given period. However, a smaller value of f_dot = 3.14x10^-16 Hz/s cannot be ruled out. Our result supports the idea that the 321.5 s period is the orbital period and that the system is the shortest period binary known so far and that it is one of the strongest sources of constant gravitational radiation in the sky. Furthermore, the decrease of the period strongly favours the unipolar inductor (or electric star) model rather than the accretion models.

By analogy to Jovian radio emissions powered by the electromagnetic interaction between Jupiter and its moons, we propose that close magnetic-nonmagnetic white-dwarf pairs and white-dwarf planetary systems are strong radio sources. A simple model is developed to predict the flux densities of radio emission generated by a loss-cone-driven electron-cyclotron maser. The radio emission from these systems has high brightness temperatures, is highly polarized, and varies on a periodic cycle following the orbital rotation. Masers from magnetic-nonmagnetic white-dwarf pairs, with orbital periods <10 min, are expected to be detectable over a wide range of radio frequencies. Terrestrial planets in close orbits about magnetic white dwarfs, with orbital periods $\la 30$ hr, can also produce detectable radio emission, thus providing a means to identify Earth-sized extrasolar planets.

At radio frequencies, the current evidence for the microquasar-quasar connection is based on imaging observations showing that relativistic outflows/jets are found in both classes of objects. Some microquasars also display superluminal motion, further strengthening the view that microquasars are in fact Galactic miniatures of quasars. Here we demonstrate that this connection can be extended to incorporate timing and spectral observations. Our argument is based on the striking similarity found in the radio flaring behaviour of the Galactic superluminal source GRO J1655-40 and of the extragalactic sources, such as the blazar 3C 273. We find that the variability of GRO J1655-40 can be explained within the framework of the successful generalised shock model for compact radio sources in which the radio emission arises from shocked plasma in relativistic jets. Specifically, the multifrequency flare amplitudes, time delays and radio polarization position angle measurements are consistent with the predictions of the growth stage of this model.

The X-ray sources observed in nearby galaxies show brightness distributions (log(N(>S))--log(S) curves) which can be described as single or broken power laws. Single power-law distributions are often found for sources in galaxies with vigorous ongoing star-formation activity, while broken power laws are found in elliptical galaxies and in bulges of spiral galaxies. The luminosity break can be caused by a population of X-ray binaries which contain a neutron star accreting at the Eddington limit or by aging of the X-ray binary population. We show that a simple birth-death model can reproduce single and broken power-law log(N(>S))--log(S) curves. We have found that power-law log(N(>S))--log(S) curves are a consequence of smooth continuous formation of X-ray binaries and a luminosity break is the signature of a starburst episode in the recent past. The luminosity break is robust and is determined by the lifespans of the X-ray binaries relative to the look-back time to the starburst epoch. The model successfully explains the different forms of log(N(>S))--log(S) curves for the disk and bulge sources in the spiral galaxy M81.

We review the current observational status of the ROSAT sources RX J1914.4+2456 and RX J0806.3+1527, and the evidence that these are ultra-short period (<10min) binary systems. We argue that an Intermediate Polar interpretation can be ruled out, that they are indeed compact binaries with a degenerate secondary, and that the period seen in the X-ray and optical is the orbital period. A white dwarf primary is preferred, but a neutron star cannot be excluded. We examine the capability of the three current double-degenerate models (Polar, Direct Accretor and Electric Star) to account for the observational characteristics of these systems. All models have difficulties with some aspects of the observations, but none can be excluded with confidence at present. The Electric Star model provides the best description, but the lifetime of this phase requires further investigation. These ultra-short period binaries will be strong gravitational wave emitters in the LISA bandpass, and because of their known source properties will be important early targets for gravitational wave studies.

A hydrodynamic formulation for accretion flow channeled by a dipolar magnetic field is constructed using a curvi-linear coordinate system natural to the field structure. We solve the hydrodynamic equations and determine the velocity, density and temperature profiles of the post-shock accretion flow. The results are applied to accretion flows in intermediate polars. We have found that for systems with massive white dwarfs (~1Msolar) the temperature profiles in the flow can differ significantly to those obtained from models in which the accretion column is assumed to be cylindrical.

We present the soft X-ray spectrum of the LINER galaxy M81 derived from a long observation with the XMM RGS. The spectrum is dominated by continuum emission from the active nucleus, but also contains emission lines from Fe L, and H-like and He-like N, O, and Ne. The emission lines are significantly broader than the RGS point-source spectral resolution; in the cross dispersion direction the emission lines are detected adjacent to, as well as coincident with, the active nucleus. This implies that they originate in a region of a few arc-minutes spatial extent (1 arc-minute ~ 1 kpc in M81). The flux ratios of the OVII triplet suggest that collisional processes are responsible for the line emission. A good fit to the whole RGS spectrum is obtained using a model consisting of an absorbed power law from the active nucleus and a 3 temperature optically thin thermal plasma. Two of the thermal plasma components have temperatures of 0.18+-0.04 keV and 0.64+-0.04 keV, characteristic of the hot interstellar medium produced by supernovae; the combined luminosity of the plasma at these two temperatures accounts for all the unresolved bulge X-ray emission seen in the Chandra observation by Tennant et al. (2001). The third component has a higher temperature (~1.7 keV), and we argue that this, along with some of the 0.64 keV emission, comes from X-ray binaries in the bulge of M81.

We report the detection of circular polarization during the 1994 outburst of the Galactic microquasar GRO J1655-40. The circular polarization is clearly detected at 1.4 and 2.4GHz, but not at 4.8 and 8.4GHz, where its magnitude never exceeds 5 mJy. Both the sign and magnitude of the circular polarization evolve during the outburst. The time dependence and magnitude of the polarized emission can be qualitatively explained by a model based on synchrotron emission from the outbursts, but is most consistent with circular polarization arising from propagation effects through the relativistic plasma surrounding the object.

A Chandra X-Ray Observatory ACIS-S imaging observation is used to study the population of X-ray sources in the nearby Sab galaxy M81 (NGC 3031). A total of 177 sources are detected with 124 located within the D25 isophote to a limiting X-ray luminosity of 3e36 ergs/cm2/s. Source positions, count rates, luminosities in the 0.3-8.0 keV band, limiting optical magnitudes, and potential counterpart identifications are tabulated. Spectral and timing analysis of the 36 brightest sources are reported including the low-luminosity active galactic nucleus, SN 1993J, and the Einstein-discovered ultra-luminous X-ray source X6.

A Chandra ACIS-S imaging observation of the nearby galaxy M81 (NGC 3031) reveals 9 luminous soft X-ray sources. The local environments, X-ray spectral properties, and X-ray light curves of the sources are presented and discussed in the context of prevailing physical models for supersoft sources. It is shown that the sample falls within expectations based on population synthesis models taken from the literature though the high observed luminosities (L~2.0e36 to \~3.0e38 erg/s in the 0.2-2.0 keV band) and equivalent blackbody temperatures T~40 to 80 eV) place the brightest detected M81 objects at the high luminosity end of the class of supersoft sources defined by previous ROSAT and Einstein studies of nearby galaxies. This is interpreted as a natural consequence of the higher sensitivity of Chandra to hotter and more luminous systems. Most of the sources can be explained as canonical \ssss, accreting white dwarfs powered by steady surface nuclear burning, with X-ray spectra well-fit by hot white dwarf local thermodynamic equilibrium atmosphere models. An exceptionally bright source is scrutinized in greater detail as its estimated bolometric luminosity, L~1.5e39 erg/s, greatly exceeds theoretical estimates for supersoft sources. This source may be beyond the stability limit and undergoing a phase of mass outflow under extreme conditions. Alternatively, a model in which the observed X-ray spectrum arises from an accretion disk around a blackhole of mass \~1200/sqrt(cos(i)) solar masses (viewed at an inclination angle i) cannot be excluded.

We present the soft X-ray spectrum of the LINER galaxy M81 from a long observation with the XMM-Newton RGS. As well as continuum emission from the active nucleus, the spectrum shows emission lines from multi temperature collisionally ionized gas. The emission lines are significantly broader than the RGS point-source spectral resolution, and in the cross dispersion direction the emission lines are detected adjacent to, as well as coincident with, the active nucleus. This implies that they originate in a region of more than an arcminute (~ 1 kpc) spatial extent. A good fit to the whole RGS spectrum can be found using a model consisting of an absorbed power law from the active nucleus and a 3 temperature thermal plasma. Two of the thermal plasma components have temperatures of 0.18+-0.04 keV and 0.64+-0.04 keV, characteristic of the hot interstellar medium produced by supernovae; the combined luminosity of these two plasma components accounts for all the unresolved bulge X-ray emission seen in the Chandra observation by Tennant et al (2001). The third component has a higher temperature (~1.7 keV) and we propose X-ray binaries in the bulge of M81 as the most likely source of this emission.

We present our first results on field X-ray sources detected in a deep, 184.7 ks observation with the ACIS-I camera on Chandra. The observations target the Lynx field of SPICES, and contains three known X-ray-emitting clusters out to z=1.27. Not including the known clusters, in the 17'x17' ACIS-I field we detect 132 sources in the 0.5-2 keV (soft) X-ray band down to a limiting flux of \~1.7e-16 erg/cm2/s and 111 sources in the 2-10 keV (hard) X-ray band down to a limiting flux of ~1.3e-15 erg/cm2/s. The combined catalog contains a total of 153 sources, of which 42 are detected only in the soft band and 21 are detected only in the hard band. Confirming previous Chandra results, we find that the fainter sources have harder X-ray spectra, providing a consistent solution to the long-standing `spectral paradox'. From deep optical and near-infrared follow-up data, 77% of the X-ray sources have optical counterparts to I=24 and 71% of the X-ray sources have near-infrared counterparts to K=20. Four of the 24 sources in the near-IR field are associated with extremely red objects (EROs; I-K>4). We have obtained spectroscopic redshifts with the Keck telescopes of 18 of the Lynx Chandra sources. These sources comprise a mix of broad-lined active galaxies, apparently normal galaxies, and two late-type Galactic dwarfs. Intriguingly, one Galactic source is identified with an M7 dwarf exhibiting non-transient, hard X-ray emission. We review non-AGN mechanisms to produce X-ray emission and discuss properties of the Lynx Chandra sample in relation to other samples of X-ray and non-X-ray sources.

We present observations of the proposed double degenerate polar RX J1914+24. Our optical and infrared spectra show no emission lines. This, coupled with the lack of significant levels of polarisation provide difficulties for a double degenerate polar interpretation. Although we still regard the double degenerate polar model as feasible, we have explored alternative scenarios for RX~J1914+24. These include a double degenerate algol system, a neutron star-white dwarf pair and an electrically powered system. The latter model is particularly attractive since it naturally accounts for the lack of both emission lines and detectable polarisation in RX J1914+24. The observed X-ray luminosity is consistent with the predicted power output. If true, then RX J1914+24 would be the first known stellar binary system radiating largely by electrical energy.

We carried out a multiwavelenght study of the black-hole candidate LMC X-3 with XMM-Newton. The system showed a transition to a low-hard state, in which the X-ray spectrum was well fitted by a simple power law. It then returned to a high-soft state, characterised by a strong disk-blackbody component. The line-of-sight absorption column density is <~ 4 x 10^{20} cm^{-2} consistent with the foreground Galactic absorption. This rules out wind accretion. We argue that, despite LMC X-3 being a high-mass X-ray binary, Roche-lobe overflow is the main mechanism of mass transfer. From UV/optical observations in the low-hard state, we determine that the companion is a slightly evolved B5 star with a mass M ~ 4.5 M_sun. This is indeed consistent with the secondary star being close to filling its Roche lobe.

We have studied the face-on, barred spiral M83 (NGC 5236) with Chandra. Eighty-one point sources are detected (above 3.5-sigma) in the ACIS S3 image: 15 of them are within the inner 16 arcsec region (starburst nucleus, resolved for the first time with Chandra), and 23 within the inner 60 arcsec (including the bar). The luminosity distributions of the sources in the inner 60 arcsec region (nucleus and stellar bar) is a single power law, which we interpret as due to continuous, ongoing star formation. Outside this inner region, there is a smaller fraction of bright sources, which we interpret as evidence of an aging population from a past episode of star formation. About 50% of the total emission in the nuclear region is unresolved; of this, about 70% can be attributed to hot thermal plasma, and we interpret the rest as due to unresolved point sources (eg, faint X-ray binaries). The unresolved X-ray emission also shows differences between the nuclear region and the spiral arms. In the nuclear region, the electron temperature of the thermal plasma is = 0.58 keV. In the spiral arms, the thermal component is at kT = 0.31 keV and a power-law component dominates at energies >~1 keV. The high abundance of C, Ne, Mg, Si and S with respect to Fe suggests that the interstellar medium is enriched and heated by core-collapse supernova explosions and winds from massive stars.

IIIZw2 was observed with XMM-Newton in July 2000. Its X-ray spectrum can be described by a powerlaw of photon index Gamma = 1.7 with a Gaussian line at 6.7 keV. There is no significant evidence of intrinsic absorption within the source or of a soft X-ray excess. Multi-wavelength light curves over a period of 25 years show related variations from the radio to the X-rays. We interpret the radio to optical emission as synchrotron radiation, self-absorbed in the radio/millimeter region, and the X-rays as mainly due to Compton up-scattering of low-energy photons by the population of high-energy electrons that give rise to the synchrotron radiation.

Chandra has resolved the starburst nuclear region of the face-on grand-design spiral M83. Eighty-one point sources are detected above 3.5 sigma in the ACIS S3 image, and 15 of them are within the inner 16-arcsec region of the galaxy. A point source with L_x ~ 3 x 10^(38) erg/s in the 0.3--8.0 keV band is found to coincide with the infra-red nuclear photometric peak, one of the two dynamical nuclei of the galaxy. No point-like sources are resolved (at a 2.5-sigma level) at the centre of symmetry of the outer optical isophote ellipses, suspected to be another dynamical nucleus. About 50% of the total emission in the nuclear region is unresolved; of this, about 70% can be attributed to hot thermal plasma, and the rest is probably due to unresolved point sources (eg, faint X-ray binaries). The azimuthally-averaged radial distribution of the unresolved emission has a King-like profile, with no central cusp. Strong emission lines are seen in the spectrum of the optically thin plasma component. The high abundances of C, Ne, Mg, Si and S with respect to Fe suggest that the interstellar medium in the nucleus is enriched and heated by type-II supernova explosions and winds from massive stars. The cumulative luminosity distribution of the discrete X-ray sources is neither a single nor a broken power law. Separating the sources in the nuclear region (within a distance of 60 arcsec from the X-ray centre) from the rest reveals that the two groups have different luminosity distributions. The log[N(>S)]--log S curve of the sources in the inner region (nucleus and stellar bar) is a single power law, which we interpret as due to continuous, ongoing star formation. Outside the central region, there is a smaller fraction of sources brighter than the Eddington limit for an accreting neutron star.

We report on observations of a TypeII quasar at redshift z=3.288, identified as a hard X-ray source in a 185 ks observation with the Chandra X-ray Observatory and as a high-redshift photometric candidate from deep, multiband optical imaging. CXOJ084837.9+445352 (hereinafter CXO52) shows an unusually hard X-ray spectrum from which we infer an absorbing column density N(H) = (4.8+/-2.1)e23 / cm2 (90% confidence) and an implied unabsorbed 2-10 keV rest-frame luminosity of L(2-10) = 3.3e44 ergs/s, well within the quasar regime. Hubble Space Telescope imaging shows CXO52 to be elongated with slight morphological differences between the WFPC2 F814W and NICMOS F160W bands. Optical and near-infrared spectroscopy of CXO52 show high-ionization emission lines with velocity widths ~1000 km/s and flux ratios similar to a Seyfert2 galaxy or radio galaxy. The latter are the only class of high-redshift TypeII luminous AGN which have been extensively studied to date. Unlike radio galaxies, however, CXO52 is radio quiet, remaining undetected at radio wavelengths to fairly deep limits, f(4.8GHz) < 40 microJy. High-redshift TypeII quasars, expected from unification models of active galaxies and long-thought necessary to explain the X-ray background, are poorly constrained observationally with few such systems known. We discuss recent observations of similar TypeII quasars and detail search techniques for such systems: namely (1) X-ray selection, (2) radio selection, (3) multi-color imaging selection, and (4) narrow-band imaging selection. Such studies are likely to begin identifying luminous, high-redshift TypeII systems in large numbers. We discuss the prospects for these studies and their implications to our understanding of the X-ray background.

We propose a model for stellar binary systems consisting of a magnetic and a non-magnetic white-dwarf pair which is powered principally by electrical energy. In our model the luminosity is caused by resistive heating of the stellar atmospheres due to induced currents driven within the binary. This process is reminiscent of the Jupiter-Io system, but greatly increased in power because of the larger companion and stronger magnetic field of the primary. Electrical power is an alternative stellar luminosity source, following on from nuclear fusion and accretion. We find that this source of heating is sufficient to account for the observed X-ray luminosity of the 9.5-min binary RX J1914+24, and provides an explanation for its puzzling characteristics.

We present a 10-year X-ray light curve and the spectra of a peculiar X-ray transient in the spiral galaxy M81. The source was below the detection limit of ROSAT PSPC before 1993, but it brightened substantially in 1993, with luminosities exceeding the Eddington limit of a 1.5 solar-mass compact accretor. It then faded and was not firmly detected in the ROSAT HRI and PSPC observations after 1994. The Chandra image obtained in 2000 May, however, shows an X-ray source at its position within the instrumental uncertainties. The Chandra source is coincident with a star-like object in the Digitized-Sky-Survey. A Hubble image suggests that the optical object may be extended. While these three observations could be of the same object, which may be an X-ray binary containing a black-hole candidate, the possibility that the ROSAT and Chandra sources are two different objects in a dense stellar environment cannot be ruled out. The Hubble data suggests that the optical object may be a globular cluster yet to be identified.

The 1998 September outburst of the black-hole X-ray binary XTE J1550-564 was monitored at X-ray, optical and radio wavelengths. We divide the outburst sequence into five phases and discuss their multiwavelength properties. The outburst starts with a hard X-ray spike, while the soft X-ray flux rises with a longer timescale. We suggest that the onset of the outburst is determined by an increased mass transfer rate from the companion star, but the outburst morphology is determined by the distribution of specific angular momentum in the accreting matter. The companion in XTE J1550-564 is likely to be an active magnetic star, with a surface field strong enough to influence the dynamics of mass transfer. We suggest that its magnetic field can create a magnetic bag capable of confining gas inside the Roche lobe of the primary. The impulsive rise in the hard X-rays is explained by the inflow of material with low angular momentum onto the black hole, on a free-fall timescale, when the magnetic confinement breaks down. At the same time, high angular momentum matter, transferred via ordinary Roche-lobe overflow, is responsible for the formation of a disk.

The observed luminosity distributions of X-ray sources indicate the presence of several populations of X-ray binaries in the nearby galaxies. Each population has its formation and evolutionary history, depending on the host environment. The features seen in the log N($>$S) -- log S curves for different types of galaxies and for different galactic components can be reproduced by a birth-death model, in which the lifespans of the binaries are inversely proportional to their X-ray brightness. Conversely, the dynamical history of a galaxy can be inferred from the luminosity distributions of its X-ray binary populations.

We have observed variations in the optical emission lines from the X-ray binary Circinus X-1. These variations may be attributed both to orbital variations and to long term secular changes in line strength. We have detected double-peaked H-alpha emission lines on two occasions, providing the first direct evidence for an accretion disk in the system. The separation of the peaks was different on the two occasions, suggesting that the disk might have a different size. The equivalent width of the emission lines dropped by more than a factor of three between 1999 and 2000; this continues the trend seen in earlier data, so that the H-alpha equivalent width has now declined by a factor of twenty since 1976. The emission lines do not appear to show signature of orbital motion, except for the data taken near phase 0, which show a significant velocity shift. We have observed an absorption component to the He I lines on one occasion. We suggest that, unlike the P Cygni profiles seen in X-ray spectra, this absorption does not arise in the accelerating zone of a radiatively driven wind. Instead, the absorption arises in material previously ejected from the system. It was only seen on this one occasion because the strength of the emission line had dropped dramatically.

The 1998 outburst of XTE J1550-564 started with a hard X-ray spike, rising in less than a day and declining after 3-4 days; at the same time, the soft X-ray flux was rising with a longer timescale (~ 10 days). We suggest that the soft and the initial hard X-ray emission are produced by two different components of the accretion flow: a higher angular momentum flow, which forms the disk, and a lower angular momentum flow feeding the hot inner region. Thus, we argue that the onset of the outburst is determined by an increased mass transfer rate from the companion star, but the outburst morphology is also determined by the distribution of specific angular momentum in the accreting matter. In XTE J1550-564, we attribute the initial, impulsive low angular momentum accretion to the breaking down of magnetic confinement by the magnetically active secondary star. We show that a hard X-ray spike is seen at the onset of an outburst in other BHCs.

We have investigated the ionisation structure of the post-shock regions of magnetic cataclysmic variables using an analytic density and temperature structure model in which effects due to bremsstrahlung and cyclotron cooling are considered. We find that in the majority of the shock-heated region where H- and He-like lines of the heavy elements are emitted, the collisional-ionisation and corona-condition approximations are justified. We have calculated the line emissivity and ionisation profiles for Iron as a function of height within the post-shock flow. For low-mass white dwarfs, line emission takes place near the shock. For high-mass white dwarfs, most of the line emission takes place in regions well below the shock and hence it is less sensitive to the shock temperature. Thus, the line ratios are useful to determine the white-dwarf masses for the low-mass white dwarfs, but the method is less reliable when the white dwarfs are massive. Line spectra can, however, be used to map the hydrodynamic structure of the post-shock accretion flow.

Using a parameterised function for the mass loss at the base of the post-shock region, we have constructed a formulation for magnetically confined accretion flows which avoids singularities, such as the infinity in density, at the base associated with all previous formulations. With the further inclusion of a term allowing for the heat input into the base from the accreting white dwarf we are able also to obtain the hydrodynamic variables to match the conditions in the stellar atmosphere. (We do not, however, carry out a mutually consistent analysis for the match). Changes to the emitted X-ray spectra are negligible unless the thickness of mass leakage region at the base approaches or exceeds one percent of the height of the post-shock region. In this case the predicted spectra from higher-mass white dwarfs will be harder, and fits to X-ray data will predict lower white-dwarf masses than previous formulations.

Sharp dips observed in the pulse profiles of three X-ray pulsars (GX 1+4, RX J0812.4-3114 and A 0535+26) have previously been suggested to arise from partial eclipses of the emission region by the accretion column occurring once each rotation period. We present pulse-phase spectroscopy from Rossi X-ray Timing Explorer satellite observations of GX 1+4 and RX J0812.4-3114 which for the first time confirms this interpretation. The dip phase corresponds to the closest approach of the column axis to the line of sight, and the additional optical depth for photons escaping from the column in this direction gives rise to both the decrease in flux and increase in the fitted optical depth measured at this phase. Analysis of the arrival time of individual dips in GX~1+4 provides the first measurement of azimuthal wandering of a neutron star accretion column. The column longitude varies stochastically with standard deviation 2-6 degrees depending on the source luminosity. Measurements of the phase width of the dip both from mean pulse profiles and individual eclipses demonstrates that the dip width is proportional to the flux. The variation is consistent with that expected if the azimuthal extent of the accretion column depends only upon the Keplerian velocity at the inner disc radius, which varies as a consequence of the accretion rate Mdot.

We present a general formulation for stability analyses of radiative shocks with multiple cooling processes, longitudinal and transverse perturbations, and unequal electron and ion temperatures. Using the accretion shocks of magnetic cataclysmic variables as an illustrative application, we investigate the shock instabilities by examining the eigenfunctions of the perturbed hydrodynamic variables. We also investigate the effects of varying the condition at the lower boundary of the post-shock flow from a zero-velocity fixed wall to several alternative type of boundaries involving the perturbed hydrodynamic variables, and the variations of the emission from the post-shock flow under different modes of oscillations. We found that the stability properties for flow with a stationary-wall lower boundary are not significantly affected by perturbing the lower boundary condition, and they are determined mainly by the energy-transport processes. Moreover, there is no obvious correlation between the amplitude or phase of the luminosity response and the stability properties of the system. Stability of the system can, however, be modified in the presence of transverse perturbation. The luminosity responses are also altered by transverse perturbation.

We report multifrequency radio observations of XTE J1550-564 obtained with the Molonglo Observatory Synthesis Telescope and the Australia Telescope Compact Array at the time of its discovery and subsequent hard and soft X-ray outburst in 1998 September. A large radio flare was observed, peaking about 1.8 days after the X-ray flare. In addition, we present Australian Long Baseline Array images obtained shortly after the maximum of the radio flare which show evolving structure. The apparent separation velocity of the two outermost components is v>2c.

The HI Balmer emission lines of black-hole binaries show double-peaked profiles during the high-soft state and the quiescent state. In the high-soft state the profiles are asymmetric with a stronger red peak, but the profiles are symmetric in the quiescent state. We suggest that in the high-soft state the emission lines originate from the temperature-inversion layer caused by irradiative heating of an optically thick accretion disk. Irradiative heating also causes the formation of a disk wind, which mildly absorbs the blue peak of the lines. The double-peaked lines seen in the quiescent state arise from an optically thin disk. In the absence of a disk wind, the lines are unabsorbed and so the symmetry of the line profiles is preserved.

The Chandra X-ray Observatory obtained a 50-ks observation of the central region of M81 using the ACIS-S in imaging mode. The global properties of the 97 x-ray sources detected in the inner 8.3x8.3 arcmin field of M81 are examined. Roughly half the sources are concentrated within the central bulge. The remainder are distributed throughout the disk with the brightest disk sources lying preferentially along spiral arms. The average hardness ratios of both bulge and disk sources are consistent with power law spectra of index Gamma~1.6 indicative of a population of x-ray binaries. A group of much softer sources are also present. The background source-subtracted logN-logS distribution of the disk follows a power law of index ~ -0.5 with no change in slope over three decades in flux. The logN-logS distribution of the bulge follows a similar shape but with a steeper slope above ~4.0e+37 ergs/s. There is unresolved x-ray flux from the bulge with a radial profile similar to that of the bulge sources. This unresolved flux is softer than the average of the bulge sources and extrapolating the bulge logN-logS distribution towards weaker sources can only account for 20% of the unresolved flux. No strong time variability was observed for any source with the exception of one bright, soft source.

We present first results from SPICES, the Spectroscopic, Photometric, Infrared-Chosen Extragalactic Survey. SPICES is comprised of four ~30 square arcminute high Galactic latitude fields with deep BRIzJK imaging reaching depths of ~25th magnitude (AB) in the optical and ~23rd magnitude (AB) in the near-infrared. To date we have 626 spectroscopic redshifts for infrared-selected SPICES sources with K<20 (Vega). The project is poised to address galaxy formation and evolution to redshift z~2. We discuss initial results from the survey, including the surface density of extremely red objects and the fraction of infrared sources at z>1. One of the SPICES fields has been the target of a deep 190 ksec Chandra exposure; we discuss initial results from analysis of that data set. Finally, we briefly discuss a successful campaign to identify high-redshift sources in the SPICES fields.

We report the results of preliminary analysis of the XMM_Newton EPIC and RGS observations of the candidate black-hole binary LMC X-3 between February and June 2000. The observations covered both the soft and the hard X-ray spectral states. The hard-state spectra were dominated by a power-law component with a photon index Gamma = 1.9 +/- 0.1. The soft-state spectra consisted of a thermal component with a multi-colour disk temperature T_in = 0.9 keV and a power-law tail with Gamma ~ 2.5--2.7. The model in which the X-rays from LMC X-3 in the high-soft state are powered by a strong stellar wind from a massive companion is not supported by the small line-of-sight absorption (n_H <~ 10^{21} cm^{-2}) deduced from the RGS data. The transition from the soft to the hard state appears to be a continuous process associated with the changes in the mass-transfer rate.

In 1998 September, the X-ray transient XTE J1550-564 underwent a major outburst in soft and hard X-rays, followed by a radio flare. Australian Long Baseline Array images obtained shortly after the peak in the radio flare showed evolving structure. The components observed have an apparent separation velocity of >2c.

We study the optical counterpart of the black-hole X-ray binary LMC X-3, by using XMM-Newton/OM observations carried out during a low-hard X-ray state. We derive a better constraint for the temperature, mass and radius of the companion star, and we show that the star is likely to be a ~ B5 subgiant filling its Roche lobe. Taking into account the effect of X-ray irradiation, we suggest a value f_M = (1.5 +/- 0.3) M_sun for the mass function in this system, lower than previously thought; we provide a more accurate lower limit to the mass of the compact object.

We have carried out observations of the X-ray transient GX339-4 during its high-soft and low-hard X-ray spectral states. Our high-resolution spectroscopic observation in 1999 April suggests that the H-alpha line has a single-peaked profile in the low-hard state as speculated in our previous paper. The HeII 4686 line, however, has a double-peaked profile in both the high-soft and low-hard states. This suggests that the line-emission mechanism is different in the two states. Our interpretation is that double-peaked lines are emitted from a temperature-inversion layer on the accretion-disk surface when it is irradiatively heated by soft X-rays. Single-peaked lines may be emitted from outflow/wind matter driven by hard X-ray heating. We have constructed a simple plane-parallel model and we use it to illustrate that a temperature-inversion layer can be formed at the disk surface under X-ray illumination. We also discuss the conditions required for the formation of temperature inversion and line emission. Based on the velocity separations measured for the double-peaked lines in the high-soft state, we propose that GX339-4 is a low-inclination binary system. The orbital inclination is about 15 deg if the orbital period is 14.8 hours.

The strong deviation in the properties of X-ray clusters from simple scaling laws highlights the importance of non-gravitational heating and cooling processes in the evolution of proto-cluster gas. We investigate this from two directions: by finding the amount of `excess energy' required in intra-cluster gas in order to reproduce the observed X-ray cluster properties, and by studying the excess energies obtained from supernovae in a semi-analytic model of galaxy formation. Using the insights obtained from the model, we then critically discuss possible ways of achieving the high excess specific energies required in clusters. These include heating by supernovae and active galactic nuclei, the role of entropy, and the effect of removing gas through radiative cooling. Our model self-consistently follows the production of excess energy as well as its effect on star formation in galaxies. Excess energy is retained in gas as gravitational, kinetic and/or thermal energy. The gas distribution of virialized halos are modelled with a 2-parameter family of gas density profiles; gas profiles are then chosen according to the total energy of the gas. We investigate four contrasting ways of modifying gas density profiles in the presence of excess energy, in order to study the sensitivity of results to the model. In addition, we estimate the minimum excess specific energy required when all available gas profiles are considered, using a fiducial cluster with a temperature of around 2 keV. We find that the excess specific energies required in clusters lie roughly in the range 1--3 keV/particle. (abridged)

We report multifrequency radio observations of GRO J1655-40 obtained with the Australia Telescope Compact Array, the Molonglo Observatory Synthesis Telescope and the Hartebeesthoek Radio Astronomy Observatory at the time of the major hard X-ray and radio outbursts in 1994 August-September. The radio emission reached levels of the order of a few Jy and was found to be linearly polarized by up to 10%, indicating a synchrotron origin. The light curves are in good agreement with those measured with the VLA, but our closer time sampling has revealed two new short-lived events and significant deviations from a simple exponential decay. The polarization data show that the magnetic field is well ordered and aligned at right angles to the radio jets for most of the monitoring period. The time evolution of the polarization cannot be explained solely in terms of a simple synchrotron bubble model, and we invoke a hybrid `core-lobe' model with a core which contributes both synchrotron and free-free emission and `lobes' which are classical synchrotron emitters.

We have obtained optical spectra of the soft X-ray transient GRO J1655-40 during different X-ray spectral states (quiescence, high-soft and hard outburst) between 1994 Aug and 1997 Jun. Characteristic features observed during the 1996-97 high-soft state were: a) broad absorption lines at Halpha and Hbeta, probably formed in the inner disk; b) double-peaked HeII 4686 emission lines, formed in a temperature-inversion layer on the disk surface, created by the soft X-ray irradiation; c) double-peaked Halpha emission, with a strength associated with the hard X-ray flux, suggesting that it was probably emitted from deeper layers than He II 4686. The observed rotational velocities of all the double-peaked lines suggest that the disk was extended slightly beyond its tidal radius. Three classes of lines were identified in the spectra taken in 1994 Aug-Sep, during a period of low X-ray activity between two strong X-ray flares: broad absorption, broad (flat-topped) emission and narrow emission. We have found that the narrow emission lines (single-peaked or double-peaked) cannot be explained by a thin-disk model. We propose that the system was in a transient state, in which the accretion disk might have had an extended optically thin cocoon and significant matter outflow. After the onset of a hard X-ray flare the disk signatures disappeared, and strong single-peaked Halpha and Paschen emission was detected, suggesting that the cocoon became opaque to optical radiation. High-ionisation lines disappeared or weakened. Two weeks after the end of the flare, the cocoon appeared to be once again optically thin.

Pulse profiles from X-ray pulsars often exhibit strong energy dependence and both periodic and aperiodic variations with time. The great variety of profiles observed in various sources, and even from individual sources, makes it difficult to separate the numerous factors influencing the phase-dependence of the X-ray emission. These factors include the system geometry and particularly the photon energy and angle dependence of emission about the neutron star poles. Comptonisation may play an important role in determining beam patterns and hence pulse profiles in X-ray pulsars. A Monte Carlo simulation is used to investigate the beaming due to Comptonisation in a simple accretion column geometry. We apply the model to the extremely variable pulse profiles of the high-magnetic field pulsar GX 1+4.

Almost all of the extragalactic X-ray background (XRB) at 0.25 keV can be accounted for by radio-quiet quasars, allowing us to derive an upper limit of 4 \bgunit\ for the remaining background at 0.25 keV. However, the XRB from the gas halos of groups of galaxies, with gas removal due to cooling accounted for, exceeds this upper limit by an order of magnitude if non-gravitational heating is not included. We calculate this using simulations of halo merger trees and realistic gas density profiles, which we require to reproduce the observed gas fractions and abundances of X-ray clusters. In addition, we find that the entire mass range of groups, from $\sim 5\times 10^{12}$ to $\sim 10^{14}$\Ms, contributes to the 0.25 keV background in this case. In a further study, we reduce the luminosities of groups by maximally heating their gas halos while maintaining the same gas fractions. This only reduces the XRB by a factor of 2 or less. We thus argue that most of the gas associated with groups must be outside their virial radii. This conclusion is supported by X-ray studies of individual groups. The properties of both groups and X-ray clusters can be naturally explained by a model in which the gas is given excess specific energies of $\sim 1$ keV/particle by non-gravitational heating. With this excess energy, the gas is gravitationally unbound from groups, but recollapses with the formation of a cluster of temperature $\ga 1$ keV. This is similar to a model proposed by Pen, but is contrary to the evolution of baryons described by Cen \& Ostriker. (abridged)

The structure of the hot downstream region below a radiative accretion shock, such as that of an accreting compact object, may oscillate due to a global thermal instability. The oscillatory behaviour depends on the functional forms of the cooling processes, the energy exchanges of electrons and ions in the shock-heated matter, and the boundary conditions. We analyse the stability of a shock with unequal electron and ion temperatures, where the cooling consists of thermal bremsstrahlung radiation which promotes instability, plus a competing process which tends to stabilize the shock. The effect of transverse perturbations is considered also. As an illustration, we study the special case in which the stabilizing cooling process is of order 3/20 in density and 5/2 in temperature, which is an approximation for the effects of cyclotron cooling in magnetic cataclysmic variables. We vary the efficiency of the second cooling process, the strength of the electron-ion exchange and the ratio of electron and ion pressures at the shock, to examine particular effects on the stability properties and frequencies of oscillation modes.

We have analyzed the kinematical parameters of Cir X-1 to constrain the nature of its companion star, the eccentricity of the binary and the pre-supernova parameter space. We argue that the companion is most likely to be a low-mass (< 2.0 M_sun) unevolved star and that the eccentricity of the orbit is 0.94 +/- 0.04. We have evaluated the dynamical effects of the supernova explosion and we find it must have been asymmetric. On average, we find that a kick of 740 km/s is needed to account for the recently measured radial velocity of +430 km/s (Johnston, Fender & Wu) for this extreme system. The corresponding minimum kick velocity is 500 km/s. This is the largest kick needed to explain the motion of any observed binary system. If Cir X-1 is associated with the supernova remnant G321.9-0.3 then we find a limiting minimum age of this remnant of 60000 yr. Furthermore, we predict that the companion star has lost 10% of its mass as a result of stripping and ablation from the impact of the supernova shell shortly after the explosion.

We have detected the optical counterpart of the proposed double degenerate polar RX J1914+24. The I band light curve is modulated on the 9.5 min period seen in X-rays. There is no evidence for any other periods. No significant modulation is seen in J. The infrared colours of RX J1914+24 are not consistent with a main sequence dwarf secondary star. Our ASCA spectrum of RX J1914+24 is typical of a heavily absorbed polar and our ASCA light curve also shows only the 9.5 min period. We find that the folded I band and X-ray light curves are out of phase. We attribute the I band flux to the irradiated face of the donor star. The long term X-ray light curve shows a variation in the observed flux of up to an order of magnitude. These observations strengthen the view that RX J1914+24 is indeed the first double degenerate polar to be detected. In this light, we discuss the synchronising mechanisms in such a close binary and other system parameters.

We re-examine the vertical structure of the post-shock flow in the accretion region of mCVs, and the X-ray emission as a function of height. We then predict X-ray light curves and phase-resolved spectra, taking into account the vertical structure, examine the implications and check whether the predicted heights are compatible with observation.

We carried out spectroscopic observations of the candidate black-hole binary GX339-4 during its low-hard and high-soft X-ray states, in 1997 and 1998 respectively. We have found that the spectrum is dominated by emission lines of neutral elements with asymmetric, round-topped profiles in the low-hard state (probably due to a dense matter outflow from an inflated outer disk). In the high-soft state, however, the emission lines from both neutral and ionized elements have unambiguously resolved double-peaked profiles; this indicates the presence of an irradiatively-heated thin accretion disk. Higher-ionization lines have larger peak-to-peak velocity separations. Our data do not show velocity modulations of the line centres due to the orbital motion of the compact object; there are no detectable lines from the companion star. These features are consistent with those of a system with a low-mass companion star and low orbital inclination. More results (from the 1999 low-hard state) coming soon in Paper II, stay tuned...

We carried out I, R, V and B photometric observations of the neutron star X-ray binary RXTE J2123-058 shortly after the end of the X-ray outburst in mid-1998. We adopt the low mass binary model to interpret our observations. After folding our data on the 0.24821-d orbital period, and correcting for the steady brightness decline following the outburst, we observe sinusoidal oscillations with hints of ellipsoidal modulations which became progressively more evident. Our data also show that the decline in brightness was faster in the V band than in the R and I bands. This suggests both the cooling of an irradiation-heated secondary star and the fading of an accretion disc over the nights of our observations.

We present new optical and infrared (IR) observations of Cir X-1 taken near apastron. Both sets of spectra show asymmetric emission lines. Archival optical observations show that an asymmetric H-alpha emission line has been in evidence for the past twenty years, although the shape of the line has changed significantly. We present an eccentric (e ~ 0.7-0.9) low mass binary model, where the system consists of a neutron star orbiting around a (sub-)giant companion star of 3-5 solar masses. We suggest that the broad components of the emission lines arise in a high-velocity, optically thick flow near the neutron star; while the narrow components of the optical and the IR lines arise near the companion star and a heated ejecta-shell surrounding the binary respectively. In this model, the velocity of the narrow component reflects the space velocity of the binary; the implied radial velocity (+430 km/s after correcting for Galactic rotation) is the highest velocity known for an X-ray binary.

We have calculated the temperature and density structure of the hot postshock plasma in magnetically confined accretion flows, including the gravitational potential. This avoids the inconsistency of previous calculations which assume that the height of the shock is negligible. We assume a stratified accretion column with 1-d flow along the symmetry axis. We find that the calculations predict a lower shock temperature than previous calculations, with a flatter temperature profile with height. We have revised previous determinations of the masses of the white dwarf primary stars and find that for higher mass white dwarfs there is a general reduction in derived masses when the gravitational potential is included. This is because the spectrum from such flows is harder than that of previous prescriptions at intermediate energies.

As part of the Deep Extragalactic Evolutionary Probe (DEEP) survey, a sample of 190 field galaxies (I_{814} <= 23.5) in the ``Groth Survey Strip'' has been used to analyze the magnitude-size relation over the range 0.1 < z < 1.1. The survey is statistically complete to this magnitude limit. All galaxies have photometric structural parameters, including bulge fractions (B/T), from Hubble Space Telescope images, and spectroscopic redshifts from the Keck Telescope. The analysis includes a determination of the survey selection function in the magnitude-size plane as a function of redshift, which mainly drops faint galaxies at large distances. Our results suggest that selection effects play a very important role. A first analysis treats disk-dominated galaxies with B/T < 0.5. If selection effects are ignored, the mean disk surface brightness (averaged over all galaxies) increases by ~1.3 mag from z = 0.1 to 0.9. However, most of this change is plausibly due to comparing low luminosity galaxies in nearby redshift bins to high luminosity galaxies in distant bins. If this effect is allowed for, no discernible evolution remains in the disk surface brightness of bright (M_B < -19) disk-dominated galaxies. A second analysis treats all galaxies by substituting half-light radius for disk scale length, with similar conclusions. Indeed, at all redshifts, the bulk of galaxies is consistent with the magnitude-size envelope of local galaxies, i.e., with little or no evolution in surface brightness. In the two highest redshift bins (z > 0.7), a handful of luminous, high surface brightness galaxies appears that occupies a region of the magnitude-size plane rarely populated by local galaxies. Their wide range of colors and bulge fractions points to a variety of possible origins.

Models of galaxy formation should be able to predict the properties of clusters of galaxies, in particular their gas fractions, metallicities, X-ray luminosity-temperature relation, temperature function and mass-deposition-rate function. Fitting these properties places important constaints on galaxy formation on all scales. By following gas processes in detail, our semi-analytic model (based on that of Nulsen & Fabian 1997) is the only such model able to predict all of the above cluster properties. We use realistic gas fractions and gas density profiles, and as required by observations we break the self-similarity of cluster structure by including supernova heating of intracluster gas, the amount of which is indicated by the observed metallicities. We also highlight the importance of the mass-deposition-rate function as an independent and very sensitive probe of cluster structure.

The post-shock flow in the accretion region of MCVs is considered to be a magnetically confined cooling flow. The X-ray spectrum from this hot plasma is emitted from a range of temperatures and densities. Extracting information from such spectra obtained with the next generation of X-ray satellites will rely on models for the structure of this flow which include appropriate treatment of the cooling mechanisms, boundary conditions and of absorptions and reflection components. This paper summarises the current situation and explores some of the difficulties arising in the exploitation of the line emission from these systems. It concludes by looking at improvements to the models now in the pipeline.

In observations with the Australia Telescope Compact Array we have resolved the radio counterpart of the unusual X-ray binary Cir X-1 into an asymmetric extended structure on arcsecond scales. In order to quantify the asymmetry we have redetermined as accurately as possible both the optical and radio coordinates of the source. The extended emission can be understood as a compact, absorbed core at the location of the X-ray binary, and extended emission up to 2 arcsec to the southeast of the core. The arcsec-scale extended emission aligns with the larger, more symmetric arcmin-scale collimated structures in the surrounding synchrotron nebula. This suggests that the transport of mass and/or energy from the X-ray binary to the synchroton nebula is occurring via the arcsec-scale structures. The ratio of extended flux from the southeast to that from the northwest of the core is at least 2:1. Interpreted as relativistic aberration of an intrinsically symmetric jet from the source, this implies a minimum outflow velocity of 0.1 c. Alternatively, the emission may be intrinsically asymmetric, perhaps as a result of the high space velocity of the system.

New measurements of galaxy clustering and background radiations provide improved constraints on the isotropy and homogeneity of the Universe on large scales. In particular, the angular distribution of radio sources and the X-Ray Background probe density fluctuations on scales intermediate between those explored by galaxy surveys and the Cosmic Microwave Background experiments. On a scale of $\sim 100 \Mpc$ the rms density fluctuations are at the level of $\sim 10%$ and on scales larger than $300 \Mpc$ the distribution of both mass and luminous sources safely satisfies the `Cosmological Principle' of isotropy and homogeneity. The transition with scale from clumpiness to homogeneity can be phrased in terms of the fractal dimension of the galaxy and mass distributions.

We report optical spectroscopic observations of the Low Mass X-ray Binary GRO J1655-40 during outbursts in August-September 1994 and in June 1996. The spectra show emission and absorption lines from the accretion disk surrounding the X-ray source, which have enabled us to measure the motion of the black-hole candidate in this system. It is the first time among black-hole candidates that an emission-line radial velocity curve is seen in phase with the expected motion of the primary. The projected radial velocity semi-amplitude determined for the primary [K(X) = 76.2 +/- 7.5 km/s], combined with estimates of the projected radial velocity semi-amplitude of the companion star from previous observations [K(C) = 228.2 +/- 3.0 km/s], yields a 95% confidence limit, derived purely from kinematics, of M(X) > 5.1 M(sun) for the mass of the compact object in GRO J1655-40, confirming its black hole nature.

We report results based on 35 new spectroscopic redshifts obtained with the Keck Telescope for field galaxies that also have photometry and morphology from survey images taken by the refurbished HST. A sample of 24 redshifts for galaxies fainter than I = 22 has a median redshift of z ~ 0.81. This result is inconsistent with the lower median redshift of z ~ 0.6 predicted by the ``maximal merger models'' of Carlberg (1996), which otherwise fit existing data. The data match an extrapolation of the CFRS, as well as predictions of certain mild luminosity-evolution models. Nearly half of the redshifts lie in two structures at z ~ 0.81 and z ~ 1.0, showing the presence of high density concentrations spanning scales of ~ 1/h Mpc, i.e., the size of groups. We find emission lines or the presence of possible neighbors in 7 of 9 otherwise luminous galaxies with red central regions at redshifts beyond z ~ 0.7. We also note a diversity of morphological types among blue galaxies at z ~ 1, including small compact galaxies, ``chains,'' and ``blue nucleated galaxies.'' These morphologies are found among local, but generally less luminous, galaxies. Distant blue galaxies also include apparently normal late-type spirals. These findings could imply modest bursts of star formation caused by mergers or interactions of small, gas-rich galaxies with each other or with larger, well-formed galaxies. This first glimpse of very faint z ~ 1 field galaxies of diverse colors and morphologies suggests that a mixture of physical processes is at work in the formation and evolution of faint field galaxies.