Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

Pulsar timing experiments typically generate a phase-connected timing solution from a sequence of times-of-arrival (TOAs) by absolute pulse numbering, i.e. by fitting an integer number of pulses between TOAs in order to minimize the residuals with respect to a parametrized phase model. In this observing mode, rotational glitches are discovered, when the residuals of the no-glitch phase model diverge after some epoch, and glitch parameters are refined by Bayesian follow-up. Here an alternative, complementary approach is presented which tracks the pulse frequency $f$ and its time derivative $df/dt$ with a hidden Markov model (HMM), whose dynamics include stochastic spin wandering (timing noise) and impulsive jumps in $f$ and $df/dt$ (glitches). The HMM tracks spin wandering explicitly, as a specific realization of a discrete-time Markov chain. It discovers glitches by comparing the Bayes factor for glitch and no-glitch models. It ingests standard TOAs for convenience and, being fully automated, allows performance bounds to be calculated quickly via Monte Carlo simulations. Practical, user-oriented plots are presented of the false alarm probability and detection threshold (e.g. minimum resolvable glitch size) versus observational scheduling parameters (e.g. TOA uncertainty, mean delay between TOAs) and glitch parameters (e.g. transient and permanent jump sizes, exponential recovery time-scale). The HMM is also applied to $\sim 1$ yr of real data bracketing the 2016 December 12 glitch in PSR J0835-4510 as a proof of principle. It detects the known glitch and confirms that no other glitch exists in the same data with size $> 10^{-7} f$.

The past four years have seen a scientific revolution through the birth of a new field: gravitational-wave astronomy. The first detection of gravitational waves---recognised by the 2017 Nobel Prize in Physics---provided unprecedented tests of general relativity while unveiling a previously unknown class of massive black holes, thirty times more massive than the Sun. The subsequent detection of gravitational waves from a merging binary neutron star confirmed the hypothesised connection between binary neutron stars and short gamma-ray bursts while providing an independent measurement of the expansion of the Universe. The discovery enabled precision measurement of the speed of gravity while shedding light on the origin of heavy elements. At the time of writing, the Laser Interferometer Gravitational-wave Observatory (LIGO) and its European partner, Virgo, have published the detection of eleven gravitational-wave events. New, not-yet-published detections are announced on a nearly weekly basis. This fast-growing catalogue of gravitational-wave transients is expected to yield insights into a number of topics, from the equation of state of matter at supra-nuclear densities to the fate of massive stars. The science potential of 3G observatories is enormous, enabling measurements of gravitational waves from the edge of the Universe and precise determination of the neutron star equation of state. Australia is well-positioned to help develop the required technology. The Mid-term Review for the Decadal plan for Australian astronomy 2016-2025 should consider investment in a scoping study for an Australian Gravitational-Wave Pathfinder that develops and validates core technologies required for the global 3G detector network.

The dust mass absorption coefficient, $\kappa_{d}$, is the conversion function used to infer physical dust masses from observations of dust emission. However, it is notoriously poorly constrained, and it is highly uncertain how it varies, either between or within galaxies. Here we present the results of a proof-of concept study, using the DustPedia data for two nearby face-on spiral galaxies M74 (NGC 628) and M83 (NGC 5236), to create the first ever maps of $\kappa_{d}$ in galaxies. We determine $\kappa_{d}$ using an empirical method that exploits the fact that the dust-to-metals ratio of the interstellar medium is constrained by direct measurements of the depletion of gas-phase metals. We apply this method pixel-by-pixel within M74 and M83, to create maps of $\kappa_{d}$. We also demonstrate a novel method of producing metallicity maps for galaxies with irregularly-sampled measurements, using the machine learning technique of Gaussian process regression. We find strong evidence for significant variation in $\kappa_{d}$. We find values of $\kappa_{d}$ at 500 $\mu$m spanning the range 0.11-0.25 ${\rm m^{2}\,kg^{-1}}$ in M74, and 0.15-0.80 ${\rm m^{2}\,kg^{-1}}$ in M83. Surprisingly, we find that $\kappa_{d}$ shows a distinct inverse correlation with the local density of the interstellar medium. This inverse correlation is the opposite of what is predicted by standard dust models. However, we find this relationship to be robust against a large range of changes to our method - only the adoption of unphysical or highly unusual assumptions would be able to suppress it.

We use a sub-set of the DustPedia galaxy sample (461 galaxies) to investigate the effect the environment has had on galaxies. We consider Virgo cluster and field samples and also assign a density contrast parameter to each galaxy, as defined by the local density of SDSS galaxies. We consider their chemical evolution (using M_{Dust}/M_{Baryon} and M_{Gas}/M_{Baryon}), their specific star formation rate (SFR/M_{Stars}), star formation efficiency (SFR/M_{Gas}), stars-to-dust mass ratio (M_{Stars}/M_{Dust}), gas-to-dust mass ratio (M_{Gas}/M_{Dust}) and the relationship between star formation rate per unit mass of dust and dust temperature (SFR/M_{Dust} and T_{Dust}). Late type galaxies (later than Sc) in all of the environments can be modelled using simple closed box chemical evolution and a simple star formation history (SFR(t) \propto t\exp{-t/\tau}). For earlier type galaxies the physical mechanisms that give rise to their properties are clearly much more varied and require a more complicated model (mergers, gas in or outflow). However, we find little or no difference in the properties of galaxies of the same morphological type within the cluster, field or with different density contrasts. It appears that it is morphology, how and whenever this is laid down, and consistent internal physical processes that primarily determine the derived properties of galaxies in the DustPedia sample and not processes related to differences in the local environment.

Astronomy plays a major role in the scientific landscape of Namibia and Southern Africa. Considerable progress has been achieved scientifically as well as in terms of human capacity development in the field. In all wavelength regimes accessible with ground-based instruments, the largest of those instruments are situated in Southern Africa: MeerKAT, the Southern African Large Telescope, and the High Energy Stereoscopic System. Because of the excellent observing conditions from Namibian soil, further large-scale projects such as the Cherenkov Telescope Array considered sites in Namibia and the Africa Millimetre Telescope will eventually be built there. Against this background, the current situation of astronomical research and education in Namibia is reviewed, focusing on optical, radio and gamma-ray astronomy and also including smaller scale projects. Further, the role of astronomy, with particular focus on developmental aspects in the African context is outlined and the progress in human capacity development is summarized.

We present the first data release of high-resolution ($\leq0.2$ arcsec) 1.5-GHz radio images of 103 nearby galaxies from the Palomar sample, observed with the eMERLIN array, as part of the LeMMINGs survey. This sample includes galaxies which are active (LINER and Seyfert) and quiescent (HII galaxies and Absorption line galaxies, ALG), which are reclassified based upon revised emission-line diagrams. We detect radio emission $\gtrsim$ 0.2 mJy for 47/103 galaxies (22/34 for LINERS, 4/4 for Seyferts, 16/51 for HII galaxies and 5/14 for ALGs) with radio sizes typically of $\lesssim$100 pc. We identify the radio core position within the radio structures for 41 sources. Half of the sample shows jetted morphologies. The remaining half shows single radio cores or complex morphologies. LINERs show radio structures more core-brightened than Seyferts. Radio luminosities of the sample range from 10$^{32}$ to 10$^{40}$ erg s$^{-1}$: LINERs and HII galaxies show the highest and the lowest radio powers respectively, while ALGs and Seyferts have intermediate luminosities. We find that radio core luminosities correlate with black hole (BH) mass down to $\sim$10$^{7}$ M$_{\odot}$, but a break emerges at lower masses. Using [O III] line luminosity as a proxy for the accretion luminosity, active nuclei and jetted HII galaxies follow an optical fundamental plane of BH activity, suggesting a common disc-jet relationship. In conclusion, LINER nuclei are the scaled-down version of FR I radio galaxies; Seyferts show less collimated jets; HII galaxies may host weak active BHs and/or nuclear star-forming cores; and recurrent BH activity may account for ALG properties.

The New IRAM KID Arrays 2 (NIKA2) consortium has just finished installing and commissioning a millimetre camera on the IRAM 30 m telescope. It is a dual-band camera operating with three frequency multiplexed kilo-pixels arrays of Lumped Element Kinetic Inductance Detectors (LEKID) cooled at 150 mK, designed to observe the intensity and polarisation of the sky at 260 and 150 GHz (1.15 and 2 mm). NIKA2 is today an IRAM resident instrument for millimetre astronomy, such as Intra Cluster Medium from intermediate to distant clusters and so for the follow-up of Planck satellite detected clusters, high redshift sources and quasars, early stages of star formation and nearby galaxies emission. We present an overview of the instrument performance as it has been evaluated at the end of the commissioning phase.

Searches for persistent gravitational radiation from nonpulsating neutron stars in young supernova remnants (SNRs) are computationally challenging because of rapid stellar braking. We describe a practical, efficient, semi-coherent search based on a hidden Markov model (HMM) tracking scheme, solved by the Viterbi algorithm, combined with a maximum likelihood matched filter, the $\mathcal{F}$-statistic. The scheme is well suited to analyzing data from advanced detectors like the Advanced Laser Interferometer Gravitational Wave Observatory (Advanced LIGO). It can track rapid phase evolution from secular stellar braking and stochastic timing noise torques simultaneously without searching second- and higher-order derivatives of the signal frequency, providing an economical alternative to stack-slide-based semi-coherent algorithms. One implementation tracks the signal frequency alone. A second implementation tracks the signal frequency and its first time derivative. It improves the sensitivity by a factor of a few upon the first implementation, but the cost increases by two to three orders of magnitude.

Millimeter-wave continuum astronomy is today an indispensable tool for both general Astrophysics studies and Cosmology. General purpose, large field-of-view instruments are needed to map the sky at intermediate angular scales not accessible by the high-resolution interferometers and by the coarse angular resolution space-borne or ground-based surveys. These instruments have to be installed at the focal plane of the largest single-dish telescopes. In this context, we have constructed and deployed a multi-thousands pixels dual-band (150 and 260 GHz, respectively 2mm and 1.15mm wavelengths) camera to image an instantaneous field-of-view of 6.5arc-min and configurable to map the linear polarization at 260GHz. We are providing a detailed description of this instrument, named NIKA2 (New IRAM KID Arrays 2), in particular focusing on the cryogenics, the optics, the focal plane arrays based on Kinetic Inductance Detectors (KID) and the readout electronics. We are presenting the performance measured on the sky during the commissioning runs that took place between October 2015 and April 2017 at the 30-meter IRAM (Institut of Millimetric Radio Astronomy) telescope at Pico Veleta. NIKA2 has been successfully deployed and commissioned, performing in-line with the ambitious expectations. In particular, NIKA2 exhibits FWHM angular resolutions of around 11 and 17.5 arc-seconds at respectively 260 and 150GHz. The NEFD (Noise Equivalent Flux Densities) demonstrated on the maps are, at these two respective frequencies, 33 and 8 mJy*sqrt(s). A first successful science verification run has been achieved in April 2017. The instrument is currently offered to the astronomical community during the coming winter and will remain available for at least the next ten years.

A hidden Markov model (HMM) scheme for tracking continuous-wave gravitational radiation from neutron stars in low-mass X-ray binaries (LMXBs) with wandering spin is extended by introducing a frequency-domain matched filter, called the J-statistic, which sums the signal power in orbital sidebands coherently. The J-statistic is similar but not identical to the binary-modulated F-statistic computed by demodulation or resampling. By injecting synthetic LMXB signals into Gaussian noise characteristic of the Advanced Laser Interferometer Gravitational-wave Observatory (Advanced LIGO), it is shown that the J-statistic HMM tracker detects signals with characteristic wave strain $h_0 \geq 2 \times 10^{-26}$ in 370 d of data from two interferometers, divided into 37 coherent blocks of equal length. When applied to data from Stage I of the Scorpius X-1 Mock Data Challenge organised by the LIGO Scientific Collaboration, the tracker detects all 50 closed injections ($h_0 \geq 6.84 \times 10^{-26}$), recovering the frequency with a root-mean-square accuracy of $\leq 1.95\times10^{-5}$ Hz. Of the 50 injections, 43 (with $h_0 \geq 1.09 \times 10^{-25}$) are detected in a single, coherent 10-d block of data. The tracker employs an efficient, recursive HMM solver based on the Viterbi algorithm, which requires $\sim 10^5$ CPU-hours for a typical, broadband (0.5-kHz), LMXB search.

The DustPedia project is capitalising on the legacy of the Herschel Space Observatory, using cutting-edge modelling techniques to study dust in the 875 DustPedia galaxies - representing the vast majority of extended galaxies within 3000 km s$^{-1}$ that were observed by Herschel. This work requires a database of multiwavelength imagery and photometry that greatly exceeds the scope (in terms of wavelength coverage and number of galaxies) of any previous local-Universe survey. We constructed a database containing our own custom Herschel reductions, along with standardised archival observations from GALEX, SDSS, DSS, 2MASS, WISE, Spitzer, and Planck. Using these data, we performed consistent aperture-matched photometry, which we combined with external supplementary photometry from IRAS and Planck. We present our multiwavelength imagery and photometry across 42 UV-microwave bands for the 875 DustPedia galaxies. Our aperture-matched photometry, combined with the external supplementary photometry, represents a total of 21,857 photometric measurements. A typical DustPedia galaxy has multiwavelength photometry spanning 25 bands. We also present the Comprehensive & Adaptable Aperture Photometry Routine (CAAPR), the pipeline we developed to carry out our aperture-matched photometry. CAAPR is designed to produce consistent photometry for the enormous range of galaxy and observation types in our data. In particular, CAAPR is able to determine robust cross-compatible uncertainties, thanks to a novel method for reliably extrapolating the aperture noise for observations that cover a very limited amount of background. Our rich database of imagery and photometry is being made available to the community

The purpose of this work is the characterization of the radial distribution of dust, stars, gas, and star-formation rate (SFR) in a sub-sample of 18 face-on spiral galaxies extracted from the DustPedia sample. This study is performed by exploiting the multi-wavelength, from UV to sub-mm bands, DustPedia database, in addition to molecular (12CO) and atomic (HI) gas maps and metallicity abundance information available in the literature. We fitted the surface brightness profiles of the tracers of dust and stars, the mass surface density profiles of dust, stars, molecular gas, and total gas, and the SFR surface density profiles with an exponential curve and derived their scale-lengths. We also developed a method to solve for the CO-to-H2 conversion factor (alpha_CO) per galaxy by using dust and gas mass profiles. Although each galaxy has its own peculiar behaviour, we identified a common trend of the exponential scale-lengths vs. wavelength. On average, the scale-lengths normalized to the B-band 25 mag/arcsec^2 radius decrease from UV to 70 micron, from 0.4 to 0.2, and then increase back up to 0.3 at 500 microns. The main result is that, on average, the dust mass surface density scale-length is about 1.8 times the stellar one derived from IRAC data and the 3.6 micron surface brightness, and close to that in the UV. We found a mild dependence of the scale-lengths on the Hubble stage T: the scale-lengths of the Herschel bands and the 3.6 micron scale-length tend to increase from earlier to later types, the scale-length at 70 micron tends to be smaller than that at longer sub-mm wavelength with ratios between longer sub-mm wavelengths and 70 micron that decrease with increasing T. The scale-length ratio of SFR and stars shows a weak increasing trend towards later types.

The European Space Agency has invested heavily in two cornerstones missions; Herschel and Planck. The legacy data from these missions provides us with an unprecedented opportunity to study cosmic dust in galaxies so that we can answer fundamental questions about, for example: the origin of the chemical elements, physical processes in the interstellar medium (ISM), its effect on stellar radiation, its relation to star formation and how this relates to the cosmic far infrared background. In this paper we describe the DustPedia project, which is enabling us to develop tools and computer models that will help us relate observed cosmic dust emission to its physical properties (chemical composition, size distribution, temperature), to its origins (evolved stars, super novae, growth in the ISM) and the processes that destroy it (high energy collisions and shock heated gas). To carry out this research we will combine the Herschel/Planck data with that from other sources of data, providing observations at numerous wavelengths (< 41) across the spectral energy distribution, thus creating the DustPedia database. To maximise our spatial resolution and sensitivity to cosmic dust we limit our analysis to 4231 local galaxies (v < 3000 km/s) selected via their near infrared luminosity (stellar mass). To help us interpret the data we have developed a new physical model for dust (THEMIS), a new Bayesian method of fitting and interpreting spectral energy distributions (HerBIE) and a state-of-the-art Monte Carlo photon tracing radiative transfer model (SKIRT). In this the first of the DustPedia papers we describe the project objectives, data sets used and provide an insight into the new scientific methods we plan to implement.

Gravitational wave searches for continuous-wave signals from neutron stars are especially challenging when the star's spin frequency is unknown a priori from electromagnetic observations and wanders stochastically under the action of internal (e.g. superfluid or magnetospheric) or external (e.g. accretion) torques. It is shown that frequency tracking by hidden Markov model (HMM) methods can be combined with existing maximum likelihood coherent matched filters like the F-statistic to surmount some of the challenges raised by spin wandering. Specifically it is found that, for an isolated, biaxial rotor whose spin frequency walks randomly, HMM tracking of the F-statistic output from coherent segments with duration T_drift = 10d over a total observation time of T_obs = 1yr can detect signals with wave strains h0 > 2e-26 at a noise level characteristic of the Advanced Laser Interferometer Gravitational Wave Observatory (Advanced LIGO). For a biaxial rotor with randomly walking spin in a binary orbit, whose orbital period and semi-major axis are known approximately from electromagnetic observations, HMM tracking of the Bessel-weighted F-statistic output can detect signals with h0 > 8e-26. An efficient, recursive, HMM solver based on the Viterbi algorithm is demonstrated, which requires ~10^3 CPU-hours for a typical, broadband (0.5-kHz) search for the low-mass X-ray binary Scorpius X-1, including generation of the relevant F-statistic input. In a "realistic" observational scenario, Viterbi tracking successfully detects 41 out of 50 synthetic signals without spin wandering in Stage I of the Scorpius X-1 Mock Data Challenge convened by the LIGO Scientific Collaboration down to a wave strain of h0 = 1.1e-25, recovering the frequency with a root-mean-square accuracy of <= 4.3e-3 Hz.

The process by which the Sun affects the terrestrial environment on short timescales is predominately driven by the amount of magnetic reconnection between the solar wind and Earth's magnetosphere. Reconnection occurs most efficiently when the solar wind magnetic field has a southward component. The most severe impacts are during the arrival of a coronal mass ejection (CME) when the magnetosphere is both compressed and magnetically connected to the heliospheric environment. Unfortunately, forecasting magnetic vectors within coronal mass ejections remains elusive. Here we report how, by combining a statistically robust helicity rule for a CME's solar origin with a simplified flux rope topology the magnetic vectors within the Earth-directed segment of a CME can be predicted. In order to test the validity of this proof-of-concept architecture for estimating the magnetic vectors within CMEs, a total of eight CME events (between 2010 and 2014) have been investigated. With a focus on the large false alarm of January 2014, this work highlights the importance of including the early evolutionary effects of a CME for forecasting purposes. The angular rotation in the predicted magnetic field closely follows the broad rotational structure seen within the in situ data. This time-varying field estimate is implemented into a process to quantitatively predict a time-varying Kp index that is described in detail in paper II. Future statistical work, quantifying the uncertainties in this process, may improve the more heuristic approach used by early forecasting systems.

Accurate space weather forecasting requires knowledge of the trajectory of coronal mass ejections (CMEs), including any deflections close to the Sun or through interplanetary space. Kay et al. 2013 introduced ForeCAT, a model of CME deflection resulting from the background solar magnetic field. For a magnetic field solution corresponding to Carrington Rotation (CR) 2029 (declining phase, April-May 2005), the majority of the CMEs deflected to the Heliospheric Current Sheet (HCS), the minimum in magnetic pressure on global scales. Most of the deflection occurred below 4 Rs. Here we extend ForeCAT to include a three dimensional description of the deflecting CME. We attempt to answer the following questions: a) Do all CMEs deflect to the magnetic minimum? and b) Does most deflection occur within the first few solar radii (~4 Rs)? Results for solar minimum and declining phase CMEs show that not every CME deflects to the magnetic minimum and that the deflection is typically determined below 2 Rs. Slow, wide, low mass CMEs in declining phase solar backgrounds with strong magnetic field and magnetic gradients exhibit the largest deflections. Local gradients related to active regions tend to cause the largest deviations from the deflection predicted by global magnetic gradients, but variations can also be seen for CMEs in the quiet sun regions of the declining phase CR. We show the torques due to differential forces along the CME can cause rotation about the CME's toroidal axis.

The acceleration of protons and electrons to high (sometimes GeV/nucleon) energies by solar phenomena is a key component of space weather. These solar energetic particle (SEP) events can damage spacecraft and communications, as well as present radiation hazards to humans. In-depth particle acceleration simulations have been performed for idealized magnetic fields for diffusive acceleration and particle propagation, and at the same time the quality of MHD simulations of coronal mass ejections (CMEs) has improved significantly. However, to date these two pieces of the same puzzle have remained largely decoupled. Such structures may contain not just a shock but also sizable sheath and pileup compression regions behind it, and may vary considerably with longitude and latitude based on the underlying coronal conditions. In this work, we have coupled results from a detailed global three-dimensional MHD time-dependent CME simulation to a global proton acceleration and transport model, in order to study time-dependent effects of SEP acceleration between 1.8 and 8 solar radii in the 2005 May 13 CME. We find that the source population is accelerated to at least 100 MeV, with distributions enhanced up to six orders of magnitude. Acceleration efficiency varies strongly along field lines probing different regions of the dynamically evolving CME, whose dynamics is influenced by the large-scale coronal magnetic field structure. We observe strong acceleration in sheath regions immediately behind the shock.

To predict whether a coronal mass ejection (CME) will impact Earth, the effects of the background on the CME's trajectory must be taken into account. We develop a model, ForeCAT (Forecasting a CME's Altered Trajectory), of CME deflection due to magnetic forces. ForeCAT includes CME expansion, a three-part propagation model, and the effects of drag on the CME's deflection. Given the background solar wind conditions, the launch site of the CME, and the properties of the CME (mass, final propagation speed, initial radius, and initial magnetic strength), ForeCAT predicts the deflection of the CME. Two different magnetic backgrounds are considered: a scaled background based on type II radio burst profiles and a Potential Field Source Surface (PFSS) background. For a scaled background where the CME is launched from an active region located between a CH and streamer region the strong magnetic gradients cause a deflection of 8.1 degrees in latitude and 26.4 degrees in longitude for a 1e15 g CME propagating out to 1 AU. Using the PFSS background, which captures the variation of the streamer belt position with height, leads to a deflection of 1.6 degrees in latitude and 4.1 degrees in longitude for the control case. Varying the CME's input parameters within observed ranges leads to the majority of CMEs reaching the streamer belt within the first few solar radii. For these specific backgrounds, the streamer belt acts like a potential well that forces the CME into an equilibrium angular position.

We present a new global model of the solar corona, including the low corona, the transition region and the top of chromosphere. The realistic 3D magnetic field is simulated using the data from the photospheric magnetic field measurements. The distinctive feature of the new model is incorporating the MHD Alfven wave turbulence. We assume this turbulence and its non-linear dissipation to be the only momentum and energy source for heating the coronal plasma and driving the solar wind. The difference between the turbulence dissipation efficiency in coronal holes and that in closed field regions is because the non-linear cascade rate degrades in strongly anisotropic turbulence in coronal holes (no inward propagating wave), thus resulting in colder coronal holes with the bi-modal solar wind originating from them. The detailed presentation of the theoretical model is illustrated with the synthetic images for multi-wavelength EUV emission compared with the observations from SDO AIA and Stereo EUVI instruments for the Carrington rotation 2107.

Here we investigate the contribution of surface Alfven wave damping to the heating of the solar wind in minima conditions. These waves are present in regions of strong inhomogeneities in density or magnetic field (e. g., the border between open and closed magnetic field lines). Using a 3-dimensional Magnetohydrodynamics (MHD) model, we calculate the surface Alfven wave damping contribution between 1-4 solar radii, the region of interest for both acceleration and coronal heating. We consider waves with frequencies lower than those that are damped in the chromosphere and on the order of those dominating the heliosphere. In the region between open and closed field lines, within a few solar radii of the surface, no other major source of damping has been suggested for the low frequency waves we consider here. This work is the first to study surface Alfven waves in a 3D environment without assuming a priori a geometry of field lines or magnetic and density profiles. We determine that waves with frequencies >2.8x10^-4 Hz are damped between 1-4 solar radii. In quiet sun regions, surface Alfven waves are damped at further distances compared to active regions, thus carrying additional wave energy into the corona. We compare the surface Alfven wave contribution to the heating by a variable polytropic index and find that it an order of magnitude larger than needed for quiet sun regions. For active regions the contribution to the heating is twenty percent. As it has been argued that a variable gamma acts as turbulence, our results indicate that surface Alfven wave damping is comparable to turbulence in the lower corona. This damping mechanism should be included self consistently as an energy driver for the wind in global MHD models.

We present observations that were carried out with the Two HUndred Micron PhotometER (THUMPER) mounted on the James Clerk Maxwell Telescope (JCMT) in Hawaii, at a wavelength of 200 um (frequency 1.5 THz). The observations utilise a small atmospheric window that opens up at this wavelength under very dry conditions at high-altitude observing sites. The atmosphere was calibrated using the sky-dipping method and a relation was established between the optical depth, tau, at 1.5 THz and that at 225 GHz: tau_1.5THz = (95 +/- 10)*tau_225GHz. Mars and Jupiter were mapped from the ground at this wavelength for the first time, and the system characteristics measured. A noise equivalent flux density (NEFD) of ~65 +/- 10 Jy (1 sigma 1 second) was measured for the THUMPER-JCMT combination, consistent with predictions based upon our laboratory measurements. The main-beam resolution of 14 arcsec was confirmed and an extended error-beam detected at roughly two-thirds of the magnitude of the main beam. Measurements of the Sun allow us to estimate that the fraction of the power in the main beam is ~15%, consistent with predictions based on modelling the dish surface accuracy. It is therefore shown that the sky over Mauna Kea is suitable for astronomy at this wavelength under the best conditions. However, higher or drier sites should have a larger number of useable nights per year.

The Monte Carlo method is a powerful tool for performing radiative equilibrium calculations, even in complex geometries. The main drawback of the standard Monte Carlo radiative equilibrium methods is that they require iteration, which makes them numerically very demanding. Bjorkman & Wood recently proposed a frequency distribution adjustment scheme, which allows radiative equilibrium Monte Carlo calculations to be performed without iteration, by choosing the frequency of each re-emitted photon such that it corrects for the incorrect spectrum of the previously re-emitted photons. Although the method appears to yield correct results, we argue that its theoretical basis is not completely transparent, and that it is not completely clear whether this technique is an exact rigorous method, or whether it is just a good and convenient approximation. We critically study the general problem of how an already sampled distribution can be adjusted to a new distribution by adding data points sampled from an adjustment distribution. We show that this adjustment is not always possible, and that it depends on the shape of the original and desired distributions, as well as on the relative number of data points that can be added. Applying this theorem to radiative equilibrium Monte Carlo calculations, we provide a firm theoretical basis for the frequency distribution adjustment method of Bjorkman & Wood, and we demonstrate that this method provides the correct frequency distribution through the additional requirement of radiative equilibrium. We discuss the advantages and limitations of this approach, and show that it can easily be combined with the presence of additional heating sources and the concept of photon weighting. However, the method may fail if small dust grains are included... (abridged)

The HIDEEP survey (Minchin et al. 2003) was done in an attempt to find objects having low inferred neutral hydrogen column densities, yet they found a distribution which was strongly peaked at 10^20.65 cm^-2. In an attempt to understand this distribution and similar survey results, we model HI profiles of gas discs and use simple simulations of objects having a wide range of HI properties in the presence of an ionizing background. We find that inferred column density (N_HI^o) values, which are found by averaging total HI masses over some disc area, do not vary strongly with central column density (N_max) for detectable objects, so that even a population having a wide range of N_max values will give rise to a strongly peaked distribution of N_HI^o values. We find that populations of objects, having a wide range of model parameters, give rise to inferred column density distributions around 10^20.6+/-0.3 cm^-2. However, populations of fairly massive objects having a wide range of central column densities work best in reproducing the HIDEEP data, and these populations are also consistent with observed Lyman limit absorber counts. It may be necessary to look two orders of magnitude fainter than HIDEEP limits to detect ionized objects having central column densities <10^20 cm^-2, but the inferred column densities of already detected objects might be lower if their radii could be estimated more accurately.

According to the Cold Dark Matter (CDM) hierarchical clustering theory of galaxy and large scale structure formation, there should be numerous low mass dark matter haloes present in the Universe today. If these haloes contain sufficient stars they should be detectable as low luminosity stellar systems or dwarf galaxies. We have previously described a new detection method for faint low surface brightness objects and shown that there are relatively large numbers of very faint dwarf galaxies in the nearby Virgo cluster. In this paper we present results from a similar survey carried out on the Millennium Galaxy strip which runs along the celestial equator and samples a very different galaxy environment. We show that the dwarf-to-giant galaxy number ratio along this strip ranges from 0.7:1 to, at most, 6:1, corresponding to a flat luminosity function (alpha approx -0.8 to -1.0). This is very different to our value of 20:1 for the Virgo cluster. There is no population of low surface brightness dwarf galaxies in the field that have gone undetected by the redshift surveys. This result is exactly opposite to what CDM models predict for the environmental dependence of the dark matter mass function which is that there are proportionally more small dark matter haloes in lower density environments.

We have carried out a fully sampled large area ($4^{\circ} \times 8^{\circ}$) 21cm \HI line survey of part of the Virgo cluster using the Jodrell Bank multi-beam instrument. The survey has a sensitivity some 3 times better than the standard HIJASS and HIPASS surveys. We detect 31 galaxies, 27 of which are well known cluster members. The four new detections have been confirmed in the HIPASS data and by follow up Jodrell Bank pointed observations. One object lies behind M86, but the other 3 have no obvious optical counter parts upon inspection of the digital sky survey fields. These 3 objects were mapped at Arecibo with a smaller \am{3}{6} HPBW and a 4 times better sensitivity than the Jodrell Bank data, which allow an improved determination of the dimensions and location of two of the objects, but surprisingly failed to detect the third. The two objects are resolved by the Arecibo beam giving them a size far larger than any optical images in the nearby field. To our mass limit of $5 \times 10^{7}$ $\frac{\Delta v}{50 km s^{-1}}$ $M_{\odot}$ and column density limit of $3 \times 10^{18}$ $\frac{\Delta v}{50 km s^{-1}}$ atoms cm$^{-2}$ these new detections represent only about 2% of the cluster atomic hydrogen mass. Our observations indicate that the \HI mass function of the cluster turns down at the low mass end making it very different to the field galaxy \HI mass function. This is quite different to the Virgo cluster optical luminosity function which is much steeper than that in the general field. Many of the sample galaxies are relatively gas poor compared to \HI selected samples of field galaxies, confirming the 'anaemic spirals' view of Virgo cluster late type galaxies.

We present SKIRT (Stellar Kinematics Including Radiative Transfer), a new Monte Carlo radiative transfer code that allows the calculation of the observed stellar kinematics of a dusty galaxy. The code incorporates the effects of both absorption and scattering by interstellar dust grains, and calculates the Doppler shift of the emerging radiation exactly by taking into account the velocities of the emitting stars and the individual scattering dust grains. The code supports arbitrary distributions of dust through a cellular approach, whereby the integration through the dust is optimized by means of a novel efficient trilinear interpolation technique. We apply our modelling technique to calculate the observed kinematics of realistic models for dusty disc galaxies. We find that the effects of dust on the mean projected velocity and projected velocity dispersion are severe for edge-on galaxies. For galaxies which deviate more than a few degrees from exactly edge-on, the effects are already strongly reduced. As a consequence, dust attenuation cannot serve as a possible way to reconcile the discrepancy between the observed shallow slopes of the inner rotation curves of LSB galaxies and the predictions of CDM cosmological models. For face-on galaxies, the velocity dispersion increases with increasing dust mass due to scattering, but the effects are limited, even for extended dust distributions. Finally, we show that serious errors can be made when the individual velocities of the dust grains are neglected in the calculations.

Some Lyman continuum photons are likely to escape from most galaxies, and these can play an important role in ionizing gas around and between galaxies, including gas that gives rise to Lyman alpha absorption. Thus the gas surrounding galaxies and in the intergalactic medium will be exposed to varying amounts of ionizing radiation depending upon the distances, orientations, and luminosities of any nearby galaxies. The ionizing background can be recalculated at any point within a simulation by adding the flux from the galaxies to a uniform quasar contribution. Normal galaxies are found to almost always make some contribution to the ionizing background radiation at redshift zero, as seen by absorbers and at random points in space. Assuming that about 2 percent of ionizing photons escape from a galaxy like the Milky Way, we find that normal galaxies make a contribution of at least 30 to 40 percent of the assumed quasar background. Lyman alpha absorbers with a wide range of neutral column densities are found to be exposed to a wide range of ionization rates, although the distribution of photoionization rates for absorbers is found to be strongly peaked. On average, less highly ionized absorbers are found to arise farther from luminous galaxies, while local fluctuations in the ionization rate are seen around galaxies having a wide range of properties.

We have computed new estimates of the local rates of supernovae (SNe) adding the updated log of Evans' visual search to our SN search database. In this way, we have accumulated the largest SN statistics ever assembled for this purpose. The new SN rates are corrected on an empirical basis for the bias in the inner regions of galaxies and that in inclined spirals. We also tested an alternative approach based on the simple model proposed by Hatano et al. (1998) for the SN and dust distribution in spirals. It turns out that, although the two approaches give similar average rates, the Hatano et al. model appears to overcorrect the SN rate of distant galaxies. We used these updated statistics to probe the SN rates with different tracers of the star formation activity in galaxies, namely integrated colors, infrared luminosities and nuclear activities. We found a clear relation between the core-collapse SN rate and the integrated galaxy color, which appears consistent with the prediction of galaxy evolutionary models. We also compared SN rates in galaxies with different L_FIR with unfavorable outcome, and we argue that L_FIR is not a universal measurement of SFR. Finally, we confirm that the SN rate is not enhanced in AGN host galaxies which indicates that the nuclear engine does not significantly stimulate the extranuclear SF.

We have used the Wide Field Planetary Camera 2 aboard the Hubble Space Telescope to image the compact star cluster nuclei of the nearby, late-type, low-luminosity spiral galaxies NGC 4395, NGC 4242, and ESO 359-029. We also analyze archival WFPC2 observations of the compact star cluster nucleus of M33. A comparative analysis of the structural and photometric properties of these four nuclei is presented. All of the nuclei are very compact, with luminosity densities increasing at small radii to the resolution limit of our data. NGC 4395 contains a Seyfert 1 nucleus with a distinct bipolar structure and bright associated filaments which are likely due to [OIII] emission. The M33 nucleus has a complex structure, with elongated isophotes and possible signatures of weak activity, including a jet-like component. The other two nuclei are not known to be active, but share similar physical size scales and luminosities to the M33 and NGC 4395 nuclei. The circumnuclear environments of all four of our program galaxies are extremely diffuse, have only low-to-moderate star formation, and appear to be devoid of large quantities of dust. The central gravitational potentials of the galaxies are also quite shallow, making the origin of these types of `naked' nuclei problematic.

We present HST photometry for three fields in the outer disk of the LMC extending approximately four magnitudes below the faintest main sequence turnoff. We cannot detect any strongly significant differences in the stellar populations of the three fields based on the morphologies of the color-magnitude diagrams, the luminosity functions, and the relative numbers of stars in different evolutionary stages. Our observations therefore suggest similar star formation histories in these regions, although some variations are certainly allowed. The fields are located in two regions of the LMC: one is in the north-east field and two are located in the north-west. Under the assumption of a common star formation history, we combine the three fields with ground-based data at the same location as one of the fields to improve statistics for the brightest stars. We compare this stellar population with those predicted from several simple star formation histories suggested in the literature, using a combination of the R-method of Bertelli et al (1992) and comparisons with the observed luminosity function. The only model which we consider that is not rejected by the observations is one in which the star formation rate is roughly constant for most of the LMC's history and then increases by a factor of three about 2 Gyr ago. Such a model has roughly equal numbers of stars older and younger than 4 Gyr, and thus is not dominated by young stars. This star formation history, combined with a closed box chemical evolution model, is consistent with observations that the metallicity of the LMC has doubled in the past 2 Gyr.

We present an F606W-F814W color-magnitude diagram for the Draco dwarf spheroidal galaxy based on Hubble Space Telescope WFPC2 images. The luminosity function is well-sampled to 3 magnitudes below the turn-off. We see no evidence for multiple turnoffs and conclude that, at least over the field of the view of the WFPC2, star formation was primarily single-epoch. If the observed number of blue stragglers is due to extended star formation, then roughly 6% (upper limit) of the stars could be half as old as the bulk of the galaxy. The color difference between the red giant branch and the turnoff is consistent with an old population and is very similar to that observed in the old, metal-poor Galactic globular clusters M68 and M92. Despite its red horizontal branch, Draco appears to be older than M68 and M92 by 1.6 +/- 2.5 Gyrs, lending support to the argument that the ``second parameter'' which governs horizontal branch morphology must be something other than age. Draco's observed luminosity function is very similar to that of M68, and the derived initial mass function is consistent with that of the solar neighborhood.

With the purpose to obtain new estimates of the rate of supernovae we joined the logs of five SN searches, namely the Asiago, Crimea, Cal{\'a}n-Tololo and OCA photographic surveys and the visual search by Evans (the sample counts 110 SNe). We found that the most prolific galaxies are late spirals in which most SNe are of type II (0.88 SNu). SN Ib/c are rarer than SN Ia (0.16 and 0.24 SNu, respectively), ruling out previous claims of a very high rate of SNIb/c. We also found that the rate of SN Ia in ellipticals (0.13 SNu) is smaller than in spirals, supporting the hypothesis of different ages of the progenitor systems in early and late type galaxies. Finally, we estimated that even assuming that separate classes of faint SN Ia and SN II do exist (SNe 1991bg and 1987A could be the respective prototypes) the overall SN rate is raised only by 20-30%, therefore excluding that faint SNe represent the majority of SN explosions. Also, the bright SNIIn are intrinsically very rare (2 to 5% of all SNII in spirals).