Chemically tagging groups of stars born in the same birth cluster is a major goal of spectroscopic surveys. To investigate the feasibility of such strong chemical tagging, we perform a blind chemical tagging experiment on abundances measured from APOGEE survey spectra. We apply a density-based clustering algorithm to the eight dimensional chemical space defined by [Mg/Fe], [Al/Fe], [Si/Fe], [K/Fe], [Ti/Fe], [Mn/Fe], [Fe/H], and [Ni/Fe], abundances ratios which together span multiple nucleosynthetic channels. In a high quality sample of 182,538 giant stars, we detect twenty-one candidate clusters with more than fifteen members. Our candidate clusters are more chemically homogeneous than a population of non-member stars with similar [Mg/Fe] and [Fe/H], even in abundances not used for tagging. Group members are consistent with having the same age and fall along a single stellar-population track in logg vs. Teff space. Each group's members are distributed over multiple kpc, and the spread in their radial and azimuthal actions increases with age. We qualitatively reproduce this increase using N-body simulations of cluster dissolution in Galactic potentials that include transient winding spiral arms. Observing our candidate birth clusters with high-resolution spectroscopy in other wavebands to investigate their chemical homogeneity in other nucleosynthetic groups will be essential to confirming the efficacy of strong chemical tagging. Our initially spatially-compact but now widely dispersed candidate clusters will provide novel limits on chemical evolution and orbital diffusion in the Galactic disc, and constraints on star formation in loosely-bound groups.

We compare the results of high-resolution simulations of individual dark matter subhalos evolving in external tidal fields with and without baryonic bulge and disk components, where the average dark matter particle mass is three orders of magnitude smaller than cosmological zoom-in simulations of galaxy formation. The Via Lactea II simulation is used to setup our initial conditions and provides a basis for our simulations of subhalos in a dark matter-only tidal field, while an observationally motivated model for the Milky Way is used for the tidal field that is comprised of a dark matter halo, a stellar disk, and a stellar bulge. Our simulations indicate that including stellar components in the tidal field results in the number of subhalos in Milky Way-like galaxies being only $65\%$ of what is predicted by $\Lambda$ Cold Dark Matter ($\Lambda$CDM). For subhalos with small pericentres $(r_p \lesssim 25$ kpc), the subhalo abundance is reduced further to $40\%$, with the surviving subhalos being less dense and having a tangentially-anisotropic orbital distribution. Conversely, subhalos with larger pericentres are minimally affected by the inclusion of a stellar component in the tidal field, with the total number of outer subhalos $\approx 75\%$ of the $\Lambda$CDM prediction. The densities of outer subhalos are comparable to predictions from $\Lambda$CDM, with the subhalos having an isotropic distribution of orbits. These ratios are higher than those found in previous studies that include the effects baryonic matter, which are affected by spurious disruption caused by low resolution.

Massive globular clusters lose stars via internal and external processes. Internal processes include mainly two-body relaxation, while external processes include interactions with the Galactic tidal field. We perform a suite of N-body simulations of such massive clusters using three different direct-summation N-body codes, exploring different Galactic orbits and particle numbers. By inspecting the rate at which a star's energy changes as it becomes energetically unbound from the cluster, we can neatly identify two populations we call kicks and sweeps, that escape through two-body encounters internal to the cluster and the external tidal field, respectively. We find that for a typical halo globular cluster on a moderately eccentric orbit, sweeps are far more common than kicks but the total mass loss rate is so low that these clusters can survive for hundreds of Hubble times. The different N-body codes give largely consistent results, but we find that numerical artifacts may arise in relation to the time-step parameter of the Hermite integration scheme, namely that the value required for convergent results is sensitive to the number of particles.

We make use of APOGEE and $Gaia$ data to identify stars that are consistent with being born in the same association or star cluster as the Sun. We limit our analysis to stars that match solar abundances within their uncertainties, as they could have formed from the same Giant Molecular Cloud (GMC) as the Sun. We constrain the range of orbital actions that solar siblings can have with a suite of simulations of solar birth clusters evolved in static and time-dependent tidal fields. The static components of each galaxy model are the bulge, disk, and halo, while the various time-dependent components include a bar, spiral arms, and GMCs. In galaxy models without GMCs, simulated solar siblings all have $J_R < 122$ km $\rm s^{-1}$ kpc, $990 < L_z < 1986$ km $\rm s^{-1}$ kpc, and $0.15 < J_z < 0.58$ km $\rm s^{-1}$ kpc. Given the actions of stars in APOGEE and $Gaia$, we find 104 stars that fall within this range. One candidate in particular, Solar Sibling 1, has both chemistry and actions similar enough to the solar values that strong interactions with the bar or spiral arms are not required for it to be dynamically associated with the Sun. Adding GMCs to the potential can eject solar siblings out of the plane of the disk and increase their $J_z$, resulting in a final candidate list of 296 stars. The entire suite of simulations indicate that solar siblings should have $J_R < 122$ km $\rm s^{-1}$ kpc, $353 < L_z < 2110$ km $\rm s^{-1}$ kpc, and $J_z < 0.8$ km $\rm s^{-1}$ kpc. Given these criteria, it is most likely that the association or cluster that the Sun was born in has reached dissolution and is not the commonly cited open cluster M67.

Observations of the redshift z=7.085 quasar J1120+0641 have been used to search for variations of the fine structure constant, alpha, over the redshift range 5.5 to 7.1. Observations at z=7.1 probe the physics of the universe when it was only 0.8 billion years old. These are the most distant direct measurements of alpha to date and the first measurements made with a near-IR spectrograph. A new AI analysis method has been employed. Four measurements from the X-SHOOTER spectrograph on the European Southern Observatory's Very Large Telescope (VLT) directly constrain any changes in alpha relative to the value measured on Earth (alpha_0). The weighted mean strength of the electromagnetic force over this redshift range in this location in the universe is da/a = (alpha_z - alpha_0)/alpha_0 = (-2.18 +/- 7.27) X 10^{-5}, i.e. we find no evidence for a temporal change from the 4 new very high redshift measurements. When the 4 new measurements are combined with a large existing sample of lower redshift measurements, a new limit on possible spatial variation of da/a is marginally preferred over a no-variation model at the 3.7 sigma level.

Astrocombs are ideal spectrograph calibrators whose limiting precision can be derived using a second, independent, astrocomb system. We therefore analyse data from two astrocombs (one 18 GHz and one 25 GHz) used simultaneously on the HARPS spectrograph at the European Southern Observatory. The first aim of this paper is to quantify the wavelength repeatability achieved by a particular astrocomb. The second aim is to measure wavelength calibration consistency between independent astrocombs, that is to place limits or measure any possible zero-point offsets. We present three main findings, each with important implications for exoplanet detection, varying fundamental constant and redshift drift measurements. Firstly, wavelength calibration procedures are important: using multiple segmented polynomials within one echelle order results in significantly better wavelength calibration compared to using a single higher-order polynomial. Segmented polynomials should be used in all applications aimed at precise spectral line position measurements. Secondly, we found that changing astrocombs causes significant zero-point offsets ($\approx 60{\rm cms}^{-1}$ in our raw data) which were removed. Thirdly, astrocombs achieve a precision of $\lesssim 4{\rm cms}^{-1}$ in a single exposure ($\approx 10\% $ above the measured photon-limited precision) and $1 {\rm cms}^{-1}$ when time-averaged over a few hours, confirming previous results. Astrocombs therefore provide the technological requirements necessary for detecting Earth-Sun analogues, measuring variations of fundamental constants and the redshift drift.

For around 100 years, hydrogen spectral modelling has been based on Voigt profile fitting. The semi-classical Voigt profile is based on a 2-level atom approximation. Whilst the Voigt profile is excellent for many circumstances, the accuracy is insufficient for very high column density damped Lyman-$\alpha$ absorption systems. We have adapted the quantum-mechanical Kramers-Heisenberg model to include thermal broadening, producing a new profile, the KHT profile. Interactions involving multiple discrete atomic levels and continuum terms, not accounted for in the Voigt model, generate asymmetries in the Lyman line wings. If not modelled, this can lead to significant systematics in parameter estimation when modelling real data. There are important ramifications in particular for measurements of the primordial deuterium abundance. However, the KHT model is complicated. We therefore present a simplified formulation based on Taylor series expansions and look-up tables, quantifying the impact of the approximations made. The KHT profile has been implemented within the widely-used VPFIT code.

We present the results of a detailed search for members of the Pal 5 tidal tail system in Gaia Data Release 2 (DR2). Tidal tails provide a sensitive method for measuring the current and past gravitational potential of their host galaxy as well as for testing predictions for the abundance of dark matter subhalos. The Pal 5 globular cluster and its associated tails are an excellent candidate for such analysis; however, only ~23 degrees of arc are currently known, with in particular the leading tail much shorter than the trailing. Using Gaia DR2 and its precise astrometry, we extend the known extent of the Pal 5 tail to ~30 degrees, 7 degrees of which are newly detected along the leading arm. The detected leading and trailing arms are symmetric in length and remain near constant width. This detection constrains proposed models in which the Galactic bar truncates Pal 5's leading arm. Follow-up spectroscopic observations are necessary to verify the candidate stream stars are consistent with the known tidal tails. If confirmed, this Pal 5 stream extension opens up new possibilities to constrain the Galactic potential.

The typical blazar S5 0716$+$714 is very interesting due to its rapid and large amplitude variability and high duty cycle of micro-variability in optical band. We analyze the observations in I, R and V bands obtained with the $1.0m$ telescope at Weihai observatory of Shandong University from 2011 to 2018. The model of synchrotron radiation from turbulent cells in a jet has been proposed as a mechanism for explaining micro-variability seen in blazar light curves. Parameters such as the sizes of turbulent cells, the enhanced particle densities, and the location of the turbulent cells in the jet can be studied using this model. The model predicts a time lag between variations as observed in different frequency bands. Automatic model fitting method for micro-variability is developed, and the fitting results of our multi-frequency micro-variability observations support the model. The results show that both the amplitude and duration of flares decomposed from the micro-variability light curves confirm to the log-normal distribution. The turbulent cell size is within the range of about 5 to 55 AU, and the time lags of the micro-variability flares between the I-R and R-V bands should be several minutes. The time lags obtained from the turbulence model are consistent with the fitting statistical results, and the time lags of flares are correlated with the time lags of the whole light curve.

We report on the extent of the effects of the Milky Way's gravitational field in shaping the structural parameters and internal dynamics of its globular cluster population. We make use of a homogeneous, up-to-date data set with kinematics, structural properties, current and initial masses of 156 globular clusters. In general, cluster radii increase as the Milky Way potential weakens; with the core and Jacobi radii being those which increase at the slowest and fastest rate respectively. We interpret this result as the innermost regions of globular clusters being less sensitive to changes in the tidal forces with the Galactocentric distance. The Milky Way's gravitational field also seems to have differentially accelerated the internal dynamical evolution of individual clusters, with those toward the bulge appearing dynamically older. Finally we find a sub-population consisting of both compact and extended globular clusters (as defined by their rh/rJ ratio) beyond 8 kpc that appear to have lost a large fraction of their initial mass lost via disruption. Moreover, we identify a third group with rh/rJ > 0.4, which have lost an even larger fraction of their initial mass by disruption. In both cases the high fraction of mass lost is likely due to their large orbital eccentricities and inclination angles, which lead to them experiencing more tidal shocks at perigalacticon and during disc crossings. Comparing the structural and orbital parameters of individual clusters allows for constraints to be placed on whether or not their evolution was relaxation or tidally dominated.

The Cosmological Principle, the combined assumptions of cosmological isotropy and homogeneity, underpins the standard model of Big Bang cosmology with which we interpret astronomical observations. A new test of isotropy over the redshift range $2<z<4$ and across large angular scales on the sky is presented. We use the cosmological distribution of neutral hydrogen, as probed by the Ly$\alpha$ forest seen towards distant quasars. The Sloan Digital Sky Survey provides the largest dataset of quasar spectra available to date. We use combined information from Data Releases 12 and 14 to select a sample of 142,661 quasars most suitable for this purpose. The scales covered by the data extend beyond post-inflation causality scales, thus probing initial conditions in the early universe. We identify significant spatially correlated systematic effects that can emulate cosmological anisotropy. Once these systematics have been accounted for, the data are found to be consistent with isotropy, providing an important independent check on the standard model, consistent with results from cosmic microwave background data.

We present direct $N$-body simulations of tidally filling 30,000 ${\rm M}_\odot$ star clusters orbiting between 10 kpc and 100 kpc in galaxies with a range of dark matter substructure properties. The time-dependent tidal force is determined based on the combined tidal tensor of the galaxy's smooth and clumpy dark matter components, the latter of which causes fluctuations in the tidal field that can heat clusters. The strength and duration of these fluctuations are sensitive to the local dark matter density, substructure fraction, sub-halo mass function, and the sub-halo mass-size relation. Based on the cold dark matter framework, we initially assume sub-halos are Hernquist spheres following a power-law mass function between $10^5$ and $10^{11} {\rm M}_\odot$ and find that tidal fluctuations are too weak and too short to affect star cluster evolution. Treating sub-halos as point masses, to explore how denser sub-halos affect clusters, we find that only sub-halos with masses greater than $10^{6} {\rm M}_\odot$ will cause cluster dissolution times to decrease. These interactions can also decrease the size of a cluster while increasing the velocity dispersion and tangential anisotropy in the outer regions via tidal heating. Hence increased fluctuations in the tidal tensor, especially fluctuations that are due to low-mass halos, do not necessarily translate into mass loss. We further conclude that the tidal approximation can be used to model cluster evolution in the tidal fields of cosmological simulations with a minimum cold dark matter sub-halo mass of $10^{6} {\rm M}_\odot$, as the effect of lower-mass sub-halos on star clusters is negligible.

Open-access telescopes of all apertures are needed to operate a competitive and efficient national science program. While larger facilities contribute light-gathering power and angular resolution, smaller ones dominate for field of view, time-resolution, and especially, total available observing time, thereby enabling our entire, diversely-expert community. Smaller aperture telescopes therefore play a critical and indispensable role in advancing science. Thus, the divestment of NSF support for modest-aperture (1 - 4 m) public telescopes poses a serious threat to U.S. scientific leadership, which is compounded by the unknown consequences of the shift from observations driven by individual investigators to survey-driven science. Given the much higher cost efficiency and dramatic science returns for investments in modest aperture telescopes, it is hard to justify funding only the most expensive facilities. We therefore urge the Astro2020 panel to explicitly make the case for modest aperture facilities, and to recommend enhancing this funding stream to support and grow this critical component of the OIR System. Further study is urgently needed to prioritize the numerous exciting potential capabilities of smaller facilities,and to establish sustainable, long-term planning for the System.

In this State of the Profession Consideration, we will discuss the state of hands-on observing within the profession, including: information about professional observing trends; student telescope training, beginning at the undergraduate and graduate levels, as a key to ensuring a base level of technical understanding among astronomers; the role that amateurs can take moving forward; the impact of telescope training on using survey data effectively; and the need for modest investments in new, standard instrumentation at mid-size aperture telescope facilities to ensure their usefulness for the next decade.

In this paper, we continue our study on the evolution of black holes (BHs) that receive velocity kicks at the origin of their host star cluster potential. We now focus on BHs in rotating clusters that receive a range of kick velocities in different directions with respect to the rotation axis. We perform N-body simulations to calculate the trajectories of the kicked BHs and develop an analytic framework to study their motion as a function of the host cluster and the kick itself. Our simulations indicate that for a BH that is kicked outside of the cluster's core, as its orbit decays in a rotating cluster the BH will quickly gain angular momentum as it interacts with stars with high rotational frequencies. Once the BH decays to the point where its orbital frequency equals that of local stars, its orbit will be circular and dynamical friction becomes ineffective since local stars will have low relative velocities. After circularization, the BH's orbit decays on a longer timescale than if the host cluster was not rotating. Hence BHs in rotating clusters will have longer orbital decay times. The timescale for orbit circularization depends strongly on the cluster's rotation rate and the initial kick velocity, with kicked BHs in slowly rotating clusters being able to decay into the core before circularization occurs. The implication of the circularization phase is that the probability of a BH undergoing a tidal capture event increases, possibly aiding in the formation of binaries and high-mass BHs.

Using newly obtained high-resolution data ($R\sim{1\times{10}^{5}}$) from the \textit{Hubble Space Telescope}, and archival UV data from the \textit{Far Ultraviolet Spectroscopic Explorer} we have conducted a detailed UV survey of the three hot, metal-polluted white dwarfs WD0455-282, WD0621-376, and WD2211-495. Using bespoke model atmospheres we measured $T_{\mathrm{eff}}$, log $g$, and photospheric abundances for these stars. In conjunction with data from Gaia we measured masses, radii, and gravitational redshift velocities for our sample of objects. We compared the measured photospheric abundances with those predicted by radiative levitation theory, and found that the observed Si abundances in all three white dwarfs, and the observed Fe abundances in WD0621-376 and WD2211-495, were larger than those predicted by an order of magnitude. These findings imply not only an external origin for the metals, but also ongoing accretion, as the metals not supported by radiative levitation would sink on extremely short timescales. We measured the radial velocities of several absorption features along the line of sight to the three objects in our sample, allowing us to determine the velocities of the photospheric and interstellar components along the line of sight for each star. Interestingly, we made detections of circumstellar absorption along the line of sight to WD0455-282 with three velocity components. To our knowledge, this is the first such detection of multi-component circumstellar absorption along the line of sight to a white dwarf.

Understanding the evolution of stellar clusters in an evolving tidal field is critical for studying the disruption of stellar clusters in a cosmological context. We systematically characterise the response of stellar clusters to tidal shocks using controlled N-body simulations of clusters with various properties that are subjected to different types of shocks. We find that the strength of the shock and the density of the cluster within the half-mass radius are the dominant properties that drive the amount of mass lost by the cluster, with the shape of the cluster profile being of minor influence. When the shock is applied as two separate sub-shocks, the amount of mass loss during the second sub-shock is sensitive to the gap time between them. Clusters that experience successive sub-shocks separated by less than their crossing time attain the same masses and sizes at the end of the simulation. However, clusters subjected to sub-shocks separated by more than a crossing time experience different evolutionary histories. The amount of mass lost in the N-body models and its scaling with shock and cluster properties differs from that predicted by classical tidal disruption theory. We demonstrate that the discrepancy is alleviated by including a dependence on the escape time-scale of unbound stars, analogously to mass loss driven by two-body relaxation. With our new theoretical model for shock-driven mass loss, the predicted relative amounts of mass loss agree with the results of the N-body simulations to ~0.3 dex across the full suite of simulations.

Gaia Data Release 2 provides a wealth of data to study the internal structure of nearby globular clusters. We use this data to investigate the internal kinematics of 11 nearby globular clusters, with a particular focus on their poorly-studied outer regions. We apply a strict set of selection criteria to remove contaminating sources and create pure cluster-member samples over a significant fraction of the radial range of each cluster. We confirm previous measurements of rotation (or a lack thereof) in the inner regions of several clusters, while extending the detection of rotation well beyond where it was previously measured and finding a steady decrease in rotation with radius. We also determine the orbital anisotropy profile and determine that clusters have isotropic cores, are radially anisotropic out to $\approx$ 4 half-light radii or $35\%$ of their limiting radii, and are then isotropic out to the limits of our datasets. We detect for the first time the presence of radial anisotropy in M 22, while confirming previous detections of radial anisotropy in 47 Tuc, M 3, M 13, M 15, and $\omega$ Cen's innermost regions. The implications of these measurements are that clusters can be separated into two categories: 1) clusters with observed radial anisotropy that likely formed tidally under-filling or are dynamically young, and 2) clusters that are primarily isotropic that likely formed tidally filling or are dynamically old.

We perform a large suite of direct N-body simulations aimed at revealing the location of the progenitor, or its remnant, of the GD-1 stream. Data from \gaia\ DR2 reveals the GD-1 stream extends over $\approx 100^\circ$, allowing us to determine the stream's leading and trailing ends. Our models suggest the length of the stream is consistent with a dynamical age of between 2-3 Gyr and the exact length, width and location of the GD-1 stream correspond to the stream's progenitor being located between $-30^\circ < \phi_{1,\mathrm{pro}} < -45^{\circ}$ in the standard GD-1 coordinate system. The model stream density profiles reveal that intact progenitors leave a strong over-density, recently-dissolved progenitors appear as gaps in the stream as escaped stars continue to move away from the remnant progenitor's location, and long-dissolved progenitors leave no observational signature on the remaining stream. Comparing our models to the GD-1 stream yields two possible scenarios for its progenitor's history: a) the progenitor reached dissolution approximately 500 Myr ago during the cluster's previous perigalactic pass and is both located at and responsible for the observed gap at $\phi_1=-40^{\circ}$ or b) the progenitor reached dissolution over 2.5 Gyr ago, the fully-dissolved remnant is at $-30^\circ < \phi_1 < -45^{\circ}$, and an observational signature of its location no longer exists. That the dissolved progenitor is in the range $-30^\circ < \phi_1 < -45^{\circ}$ implies that density fluctuations outside of this range, e.g., a deep gap at $\phi_1 \approx -20^\circ$, are likely produced by compact baryonic or dark-matter perturbers.

NGC 288 is a Galactic globular cluster having observed extra tidal structure, without confirmed tidal tails. $Gaia$ DR2 provides photometric and astrometric data for many of the stars in NGC 288 and its extra tidal structure. To compare with $Gaia$ data, we simulate an $N$-body model of a star cluster with the same orbit as NGC 288 in a Milky Way potential. The simulation shows that the cluster forms tidal tails that are compressed along the cluster's orbit when it is at apocentre and are expected to be a diffuse bipolar structure. In this letter, we present a comparison between the simulation and observations from $Gaia$ DR2. We find that both the simulation and the observations share comparable trends in the position on the sky and proper motions of the extra-tidal stars, supporting the presence of tidal tails around NGC 288.

White dwarf atmospheres are subjected to gravitational potentials around $10^5$ times larger than occur on Earth. They provide a unique environment in which to search for any possible variation in fundamental physics in the presence of strong gravitational fields. However, a sufficiently strong magnetic field will alter absorption line profiles and introduce additional uncertainties in measurements of the fine structure constant. Estimating the magnetic field strength is thus essential in this context. Here we model the absorption profiles of a large number of atomic transitions in the white dwarf photosphere, including first-order Zeeman effects in the line profiles, varying the magnetic field as a free parameter. We apply the method to a high signal-to-noise, high-resolution, far-ultraviolet HST/STIS spectrum of the white dwarf G191-B2B. The method yields a sensitive upper limit on its magnetic field of $B < 2300$ Gauss at the $3\sigma$ level. Using this upper limit we find that the potential impact of quadratic Zeeman shifts on measurements of the fine structure constant in G191-B2B is 4 orders of magnitude below laboratory wavelength uncertainties.

We present the results of a study aimed at investigating the effects of dynamical evolution on the spatial distribution and mixing of primordial binary stars in multiple-population globular clusters. Multiple stellar population formation models predict that second-generation (SG) stars form segregated in the inner regions of a more extended first-generation (FG) cluster. Our study, based on the results of a survey of N-body simulations, shows that the spatial mixing process for binary stars is more complex than that of single stars since additional processes such as binary ionization, recoil and ejection following binary-single and binary-binary interactions play a key role in determining the spatial distribution of the population of surviving binaries. The efficiency and relative importance of these additional effects depends on the binary binding energy and determines the timescale of the spatial mixing of FG and SG binaries. Our simulations illustrate the role of ionization, recoil and ejection combined with the effects of mass segregation driven by two-body relaxation and show that the complex interplay of all these processes results in a significant extension of the time needed for the complete spatial mixing of FG and SG binaries compared to that of single stars. Clusters in which FG and SG single stars have already reached complete spatial mixing might be characterized by a significant radial gradient in the ratio of the FG-to-SG binary fraction. The implications of the delayed mixing of FG and SG binaries for the differences between the kinematics of the two populations are discussed.

We incorporate a semi-analytic formula for the main sequence lifetime of helium-rich stars in N-body simulations of multiple population globular clusters to investigate how the enriched helium stars impact the dynamics of globular clusters. We show that a globular cluster with a helium-rich concentrated population will be slightly smaller than a globular cluster with a normal-helium second generation, with the largest difference seen in the extended normal-helium population. This effect is shown both for a cluster in isolation and one in a realistic Milky Way tidal field. We show that this effect is a result of mass segregation, and the earlier loss of more massive stars in a helium-rich concentrated population due to their decreased main sequence lifetime. The two populations will therefore become dynamically mixed at a slightly earlier time than if they have the same helium abundance. Furthermore, we find that it is possible for the helium-enriched population to become more extended than the normal-helium population if it forms with a low initial concentration or the cluster is able to evolve for a large number of relaxation times. We conclude that the dynamical effects of helium on stellar masses are modest, and that the initial concentration of the two populations and the strength of the Milky Way tidal field are more important in determining the relative radial distributions of multiple populations.

We present the results of our power spectral density analysis for the BL Lac object OJ\,287, utilizing the {\it Fermi}-LAT survey at high-energy $\gamma$-rays, {\it Swift}-XRT in X-rays, several ground-based telescopes and the {\it Kepler} satellite in the optical, and radio telescopes at GHz frequencies. The light curves are modeled in terms of continuous-time auto-regressive moving average (CARMA) processes. Owing to the inclusion of the {\it Kepler} data, we were able to construct \emph{for the first time} the optical variability power spectrum of a blazar without any gaps across $\sim6$ dex in temporal frequencies. Our analysis reveals that the radio power spectra are of a colored-noise type on timescales ranging from tens of years down to months, with no evidence for breaks or other spectral features. The overall optical power spectrum is also consistent with a colored noise on the variability timescales ranging from 117 years down to hours, with no hints of any quasi-periodic oscillations. The X-ray power spectrum resembles the radio and optical power spectra on the analogous timescales ranging from tens of years down to months. Finally, the $\gamma$-ray power spectrum is noticeably different from the radio, optical, and X-ray power spectra of the source: we have detected a characteristic relaxation timescale in the {\it Fermi}-LAT data, corresponding to $\sim 150$\,days, such that on timescales longer than this, the power spectrum is consistent with uncorrelated (white) noise, while on shorter variability timescales there is correlated (colored) noise.

Through a suite of direct N-body simulations, we explore how the structural and kinematic evolution of a star cluster located at the center of a dwarf galaxy is affected by the shape of its host's dark matter density profile. The stronger central tidal fields of cuspier halos minimize the cluster's ability to expand in response to mass loss due to stellar evolution during its early evolutionary stages and during its subsequent long-term evolution driven by two-body relaxation. Hence clusters evolving in cuspier dark matter halos are characterized by more compact sizes, higher velocity dispersions and remain approximately isotropic at all clustercentric distances. Conversely, clusters in cored halos can expand more and develop a velocity distribution profile that becomes increasingly radially anisotropic at larger clustercentric distances. Finally, the larger velocity dispersion of clusters evolving in cuspier dark matter profiles results in them having longer relaxation times. Hence clusters in cuspy galaxies relax at a slower rate and, consequently, they are both less mass segregated and farther from complete energy equipartition than cluster's in cored galaxies. Application of this work to observations allows for star clusters to be used as tools to measure the distribution of dark matter in dwarf galaxies and to distinguish isolated star clusters from ultra-faint dwarf galaxies.

We present the first paper of the series Origin of Metals around Galaxies (OMG) aimed to study the origin of the metals observed in the circumgalactic and intergalactic media. In this work we extract and build the catalogs of metal absorbers that will be used in future analyses, and make our results publicly available to the community. We design a fully automatic algorithm to search for absorption metal-line doublets of the species CIV, NV, SiIV and MgII in high-resolution ($R\gtrsim30\,000$) quasar spectra without human intervention, and apply it to the high-resolution and signal-to-noise ratio spectra of 690 quasars, observed with the UVES and HIRES instruments. We obtain $5\,656$ CIV doublets, $7\,919$ doublets of MgII, $2\,258$ of SiIV, and 239 of NV, constituting the largest high-resolution metal-doublet samples to date, and estimate the dependence of their completeness and purity on various doublet parameters such as equivalent width and redshift, using real and artificial quasar spectra. The catalogs include doublets with rest-frame line equivalent widths down to a few ${\rm m\AA}$, all detected at a significance above 3$\sigma$, and covering the redshifts between $1<z \lesssim 5$, properties that make them useful for a wide range of chemical evolution studies.

The spectrum of the $z_{\rm em} = 2.63$ quasar Q1009+2956 has been observed extensively on the Keck telescope. The Lyman limit absorption system $z_{\rm abs} = 2.504$ was previously used to measure D/H by Burles & Tytler using a spectrum with signal to noise approximately 60 per pixel in the continuum near Ly {\alpha} at $z_{\rm abs} = 2.504$. The larger dataset now available combines to form an exceptionally high signal to noise spectrum, around 147 per pixel. Several heavy element absorption lines are detected in this LLS, providing strong constraints on the kinematic structure. We explore a suite of absorption system models and find that the deuterium feature is likely to be contaminated by weak interloping Ly {\alpha} absorption from a low column density H I cloud, reducing the expected D/H precision. We find D/H = $2.48^{+0.41}_{-0.35}\times10^{-5}$ for this system. Combining this new measurement with others from the literature and applying the method of Least Trimmed Squares to a statistical sample of 15 D/H measurements results in a "reliable" sample of 13 values. This sample yields a primordial deuterium abundance of (D/H)$_{\rm p} = (2.545 \pm 0.025)\times10^{-5}$. The corresponding mean baryonic density of the Universe is $\Omega_{\rm b}h^2 = 0.02174\pm0.00025$. The quasar absorption data is of the same precision as, and marginally inconsistent with, the 2015 CMB Planck (TT+lowP+lensing) measurement, $\Omega_{\rm b}h^2 = 0.02226\pm0.00023$. Further quasar and more precise nuclear data are required to establish whether this is a random fluctuation.

We present a detailed analysis of the radial distribution of light-element multiple populations (LE-MPs) in the massive and dense globular cluster M80 based on the combination of UV and optical Hubble Space Telescope data. Surprisingly, we find that first generation stars (FG) are significantly more centrally concentrated than extreme second generation ones (SG) out to $\sim 2.5 r_h$ from the cluster center. To understand the origin of such a peculiar behavior, we used a set of $N$-body simulations following the long-term dynamical evolution of LE-MPs. We find that, given the advanced dynamical state of the cluster, the observed difference does not depend on the primordial relative distributions of FG and SG stars. On the contrary, a difference of $\sim 0.05-0.10 M_{\odot}$ between the average masses of the two sub-populations is needed to account for the observed radial distributions. We argue that such a mass difference might be the result of the higher He abundance of SG stars (of the order of $\Delta Y\sim 0.05-0.06$) with respect to FG. Interestingly, we find that a similar He variation is necessary to reproduce the horizontal branch morphology of M80. These results demonstrate that differences in mass among LE-MPs, due to different He content, should be properly taken into account for a correct interpretation of their radial distribution, at least in dynamically evolved systems.

We present the results of a survey of N-body simulations aimed at studying the effects of the long-term dynamical evolution on the stellar mass function (MF) of multiple stellar populations in globular clusters. Our simulations show that if first-(1G) and second-generation (2G) stars have the same initial MF (IMF), the global MFs of the two populations are affected similarly by dynamical evolution and no significant differences between the 1G and the 2G MFs arise during the cluster's evolution. If the two populations have different IMFs, dynamical effects do not completely erase memory of the initial differences. Should observations find differences between the global 1G and 2G MF, these would reveal the fingerprints of differences in their IMFs. Irrespective of whether the 1G and 2G populations have the same global IMF or not, dynamical effects can produce differences between the local (measured at various distances from the cluster centre) 1G and 2G MFs; these differences are a manifestation of the process of mass segregation in populations with different initial structural properties. In dynamically old and spatially mixed clusters, however, differences between the local 1G and 2G MFs can reveal differences between the 1G and 2G global MFs. In general, for clusters with any dynamical age, large differences between the local 1G and 2G MFs are more likely to be associated with differences in the global MF. Our study also reveals a dependence of the spatial mixing rate on the stellar mass, another dynamical consequence of the multiscale nature of multiple-population clusters.

Dynamical evolution drives globular clusters toward core collapse, which strongly shapes their internal properties. Diagnostics of core collapse have so far been based on photometry only, namely on the study of the concentration of the density profiles. Here we present a new method to robustly identify core-collapsed clusters based on the study of their stellar kinematics. We introduce the \textit{kinematic concentration} parameter, $c_k$, the ratio between the global and local degree of energy equipartition reached by a cluster, and show through extensive direct $N$-body simulations that clusters approaching core collapse and in the post-core collapse phase are strictly characterized by $c_k>1$. The kinematic concentration provides a suitable diagnostic to identify core-collapsed clusters, independent from any other previous methods based on photometry. We also explore the effects of incomplete radial and stellar mass coverage on the calculation of $c_k$ and find that our method can be applied to state-of-art kinematic datasets.

The z_abs = 2.504 Lyman Limit absorption system (LLS) towards Q1009+2956 has previously been used to estimate the primordial deuterium abundance. Since the initial measurement by Burles & Tytler, this quasar has been observed extensively with the Keck telescope, providing a substantial increase in signal-to-noise (from 60 to 147 at continuum level of Ly-alpha at z_abs=2.504). We explore a set of different models for the absorption system and find that the deuterium feature is contaminated by Ly-alpha absorption from a low column density H I cloud. This significantly limits precision to which one can estimate the D/H ratio in this LLS. Our final result for this system D/H =2.48^{+0.41}_{-0.35}*10^{-5} has the same relative uncertainty of 17% as the previous estimate by Burles & Tytler despite the far higher signal-to-noise of our dataset. A weighted mean of 13 D/H measurements available in the literature (including our result) gives the primordial deuterium abundance of (D/H)_p = (2.545 +/- 0.025)*10^{-5} and the corresponding baryon density of the Universe of Omega_b h^2 = 0.02174 +/- 0.00025 marginally inconsistent with the 2015 Planck CMB data.

We consider how dynamical friction acts on black holes that receive a velocity kick while located at the center of a gravitational potential, analogous to a star cluster, due to either a natal kick or the anisotropic emission of gravitational waves during a black hole-black hole merger. Our investigation specifically focuses on how well various Chandrasekhar-based dynamical friction models can predict the orbital decay of kicked black holes with $m_{bh} \lesssim 100 M_\odot$ due to an inhomogeneous background stellar field. In general, the orbital evolution of a kicked black hole follows that of a damped oscillator where two-body encounters and dynamical friction serve as sources of damping. However, we find models for approximating the effects of dynamical friction do not accurately predict the amount of energy lost by the black hole if the initial kick velocity $v_{k}$ is greater than the stellar velocity dispersion $\sigma$. For all kick velocities, we also find that two-body encounters with nearby stars can cause the energy evolution of a kicked BH to stray significantly from standard dynamical friction theory as encounters can sometimes lead to an energy gain. For larger kick velocities, we find the orbital decay of a black hole departs from classical theory completely as the black hole's orbital amplitude decays linearly with time as opposed to exponentially. Therefore, we have developed a linear decay formalism which scales linearly with black hole mass and $\frac{v_{k}}{\sigma}$ in order to account for the variations in the local gravitational potential.

In this work, we present the results of our multi-band microvariability study of the famous blazar BL Lac. We performed microvariablity observations of the source in the optical VRI bands for 4 nights in the year 2016. We studied the intranight flux and spectral variability of the source in detail with an objective to characterize microvariability in the blazars, a frequently observed phenomenon in blazars. The results show that the source often displays a fast flux variability with an amplitude as large as ~ 0.2 magnitude within a few hours, and that the color variability in the similar time scales can be characterized as "bluer-when-brighter" trend. We also observed markedly curved optical spectrum during one of the nights. Furthermore, the correlation between multi-band emission shows that in general the emission in all the bands are highly correlated; and in one of the nights V band emission was found to lead the I band emission by ~ 20 minutes. The search for characteristic timescale using auto-correlation function and the structure function analyses reveals characteristic timescale of ~ 48 minutes in one of the R band observations. We try to explain the observed results in the context of the passage of shock waves through the relativistic outflows in blazars.

We measure how the slope $\alpha$ of the stellar mass function (MF) changes as a function of clustercentric distance $r$ in five Galactic globular clusters and compare $\alpha(r)$ to predictions from direct $N$-body star cluster simulations. Theoretical studies predict that $\alpha(r)$ (which traces the degree of mass segregation in a cluster) should steepen with time as a cluster undergoes two-body relaxation and that the amount by which the global MF can evolve from its initial state due to stellar escape is directly linked to $\alpha(r)$. We find that the amount of mass segregation in M10, NGC 6218, and NGC 6981 is consistent with their dynamical ages, but only the global MF of M10 is consistent with its degree of mass segregation as well. NGC 5466 and NGC 6101 on the other hand appear to be less segregated than their dynamical ages would indicate. Furthermore, despite the fact that the escape rate of stars in non-segregated clusters is independent of stellar mass, both NGC 5466 and NGC 6101 have near-flat MFs. We discuss various mechanisms which could produce non-segregated clusters with near-flat MFs, including higher mass-loss rates and black hole retention, but argue that for some clusters (NGC 5466 and NGC 6101) explaining the present-day properties might require either a non-universal IMF or a much more complex dynamical history.

We monitored BL Lacertae simultaneously in the optical B, V, R and I bands for 13 nights during the period 2012-2016. The variations were well correlated in all bands and the source showed significant intraday variability (IDV). We also studied its optical flux and colour behaviour, and searched for inter-band time lags. A strong bluer-when-brighter chromatism was found on the intra-night time-scale. The spectral changes are not sensitive to the host galaxy contribution. Cross-correlation analysis revealed possible time delay of about 10 min between variations in the V and R bands. We interpreted the observed flares in terms of the model consisting of individual synchrotron pulses.

We carried out a pilot campaign of radio and optical band intra-day variability (IDV) observations of five blazars (3C66A, S5 0716+714, OJ287, B0925+504, and BL Lacertae) on December 18--21, 2015 by using the radio telescope in Effelsberg (Germany) and several optical telescopes in Asia, Europe, and America. After calibration, the light curves from both 5 GHz radio band and the optical R band were obtained, although the data were not smoothly sampled over the sampling period of about four days. We tentatively analyse the amplitudes and time scales of the variabilities, and any possible periodicity. The blazars vary significantly in the radio (except 3C66A and BL Lacertae with only marginal variations) and optical bands on intra- and inter-day time scales, and the source B0925+504 exhibits a strong quasi-periodic radio variability. No significant correlation between the radio- and optical-band variability appears in the five sources, which we attribute to the radio IDV being dominated by interstellar scintillation whereas the optical variability comes from the source itself. However, the radio- and optical-band variations appear to be weakly correlated in some sources and should be investigated based on well-sampled data from future observations.

We recently presented a new "artificial intelligence" method for the analysis of high-resolution absorption spectra (Bainbridge and Webb, Mon. Not. R. Astron. Soc. 2017, 468,1639-1670). This new method unifies three established numerical methods: a genetic algorithm (GVPFIT); non-linear least-squares optimisation with parameter constraints (VPFIT); and Bayesian Model Averaging (BMA). In this work, we investigate the performance of GVPFIT and BMA over a broad range of velocity structures using synthetic spectra. We found that this new method recovers the velocity structures of the absorption systems and accurately estimates variation in the fine structure constant. Studies such as this one are required to evaluate this new method before it can be applied to the analysis of large sets of absorption spectra. This is the first time that a sample of synthetic spectra has been utilised to investigate the analysis of absorption spectra. Probing the variation of nature's fundamental constants (such as the fine structure constant), through the analysis of absorption spectra, is one of the most direct ways of testing the universality of physical laws. This "artificial intelligence" method provides a way to avoid the main limiting factor, i.e., human interaction, in the analysis of absorption spectra.

Spectra observed with the Ultraviolet and Visual Echelle Spectrograph (UVES) on the European Southern Observatory's VLT exhibit long-range wavelength distortions. These distortions impose a systematic error on high-precision measurements of the fine-structure constant, $\alpha$, derived from intervening quasar absorption systems. If the distortion is modelled using a model that is too simplistic, the resulting bias in $\Delta\alpha/\alpha$ away from the true value can be larger than the statistical uncertainty on the $\alpha$ measurement. If the effect is ignored altogether, the same is true. If the effect is modelled properly, accounting for the way in which final spectra are generally formed from the co-addition of exposures made at several different instrumental settings, the effect can be accurately removed and the correct $\Delta\alpha/\alpha$ recovered.

Hot white dwarf stars are the ideal probe for a relationship between the fine-structure constant and strong gravitational fields, providing us with an opportunity for a direct observational test. We study a sample of hot white dwarf stars, combining far-UV spectroscopic observations, atomic physics, atmospheric modelling and fundamental physics, in the search for variation in the fine structure constant. This variation manifests as shifts in the observed wavelengths of absorption lines, such as quadruply ionized iron (FeV) and quadruply ionized nickel (NiV), when compared to laboratory wavelengths. Berengut et al. (Phys. Rev. Lett. 2013, 111, 010801) demonstrated the validity of such an analysis using high-resolution Hubble Space Telescope (HST)/Space Telescope Imaging Spectrograph (STIS) spectra of G191-B2B. We have made three important improvements by: (a) using three new independent sets of laboratory wavelengths, (b) analysing a sample of objects, and (c) improving the methodology by incorporating robust techniques from previous studies towards quasars (the Many Multiplet method). A successful detection would be the first direct measurement of a gravitational field effect on a bare constant of nature. Here we describe our approach and present preliminary results from nine objects using both FeV and NiV.

The primordial deuterium abundance probes fundamental physics during the Big Bang Nucleosynthesis and can be used to infer cosmological parameters. Observationally, the abundance can be measured using absorbing clouds along the lines of sight to distant quasars. Observations of the quasar PKS1937--101 contain two absorbers for which the deuterium abundance has previously been determined. Here we focus on the higher redshift one at $z_{abs} = 3.572$. We present new observations with significantly increased signal-to-noise ratio which enable a far more precise and robust measurement of the deuterium to hydrogen column density ratio, resulting in D/H = $2.62\pm0.05\times10^{-5}$. This particular measurement is of interest because it is among the most precise assessments to date and it has been derived from the second lowest column-density absorber [N(HI) $=17.9\mathrm{cm}^{-2}$] that has so-far been utilised for deuterium abundance measurements. The majority of existing high-precision measurements were obtained from considerably higher column density systems [i.e. N(HI) $>19.4\mathrm{cm}^{-2}$]. This bodes well for future observations as low column density systems are more common.

Spectral lines belonging to the short-lifetime heavy radioactive elements up to Es ($Z$=99) have been found in the spectra of the Przybylski's star. We suggest that these unstable elements may be decay products of a "magic" metastable nucleus belonging to the the island of stability where the nuclei have a magic number of neutrons $N=184$. The laboratory-produced nuclei have a significantly smaller number of neutrons. To identify spectra of the $N=184$ isotopes of these nuclei and their neutron-reach superheavy decay products in astrophysical data we calculate the isotope shift which should be added to the laboratory - measured wavelenghs. The results for the isotopic shifts in the strongest optical electromagnetic transitions in No, Lr, Nh, Fl,and $Z$=120 elements are presented.

We present N-body simulations of star clusters that initially evolve within a strong compressive tidal field and then transition into extensive tidal fields of varying strengths. While subject to compressive tides, clusters can undergo significant heating due to two-body interactions and mass loss due to stellar evolution. When the cluster transitions into an extensive tidal field it is super-virialized, which leads to a rapid expansion and significant mass loss before the cluster reaches virial equilibrium. After the transition, clusters are significantly less massive, more extended and therefore more tidally filling than clusters which have spent their entire lifetime in a similar extensive tidal field.

A new and automated method is presented for the analysis of high-resolution absorption spectra. Three established numerical methods are unified into one "artificial intelligence" process: a genetic algorithm (GVPFIT); non-linear least-squares with parameter constraints (VPFIT); and Bayesian Model Averaging (BMA). The method has broad application but here we apply it specifically to the problem of measuring the fine structure constant at high redshift. For this we need objectivity and reproducibility. GVPFIT is also motivated by the importance of obtaining a large statistical sample of measurements of $\Delta\alpha/\alpha$. Interactive analyses are both time consuming and complex and automation makes obtaining a large sample feasible. In contrast to previous methodologies, we use BMA to derive results using a large set of models and show that this procedure is more robust than a human picking a single preferred model since BMA avoids the systematic uncertainties associated with model choice. Numerical simulations provide stringent tests of the whole process and we show using both real and simulated spectra that the unified automated fitting procedure out-performs a human interactive analysis. The method should be invaluable in the context of future instrumentation like ESPRESSO on the VLT and indeed future ELTs. We apply the method to the $z_{abs} = 1.8389$ absorber towards the $z_{em} = 2.145$ quasar J110325-264515. The derived constraint of $\Delta\alpha/\alpha = 3.3 \pm 2.9 \times 10^{-6}$ is consistent with no variation and also consistent with the tentative spatial variation reported in Webb et al (2011) and King et al (2012).

We make use of $N$-body simulations to determine the relationship between two observable parameters that are used to quantify mass segregation and energy equipartition in star clusters. Mass segregation can be quantified by measuring how the slope of a cluster's stellar mass function $\alpha$ changes with clustercentric distance r, and then calculating $\delta_\alpha = \frac{d \alpha(r)}{d ln(r/r_m)}$ where $r_m$ is the cluster's half-mass radius. The degree of energy equipartition in a cluster is quantified by $\eta$, which is a measure of how stellar velocity dispersion $\sigma$ depends on stellar mass m via $\sigma(m) \propto m^{-\eta}$. Through a suite of $N$-body star cluster simulations with a range of initial sizes, binary fractions, orbits, black hole retention fractions, and initial mass functions, we present the co-evolution of $\delta_\alpha$ and $\eta$. We find that measurements of the global $\eta$ are strongly affected by the radial dependence of $\sigma$ and mean stellar mass and the relationship between $\eta$ and $\delta_\alpha$ depends mainly on the cluster's initial conditions and the tidal field. Within $r_m$, where these effects are minimized, we find that $\eta$ and $\delta_\alpha$ initially share a linear relationship. However, once the degree of mass segregation increases such that the radial dependence of $\sigma$ and mean stellar mass become a factor within $r_m$, or the cluster undergoes core collapse, the relationship breaks down. We propose a method for determining $\eta$ within $r_m$ from an observational measurement of $\delta_\alpha$. In cases where $\eta$ and $\delta_\alpha$ can be measured independently, this new method offers a way of measuring the cluster's dynamical state.

A number of recent observational studies of Galactic globular clusters have measured the variation in the slope of a cluster's stellar mass function $\alpha$ with clustercentric distance $r$. In order to gather a deeper understanding of the information contained in such observations, we have explored the evolution of $\alpha(r)$ for star clusters with a variety of initial conditions using a large suite of $N$-body simulations. We have specifically studied how the time evolution of $\alpha(r)$ is affected by initial size, mass, binary fraction, primordial mass segregation, black hole retention, an external tidal field, and the initial mass function itself. Previous studies have shown that the evolution of $\alpha_G$ is closely related to the amount of mass loss suffered by a cluster. Hence for each simulation we have also followed the evolution of the slope of the cluster's global stellar mass function, $\alpha_G$, and have shown that clusters follow a well-defined track in the $\alpha_G$-$d\alpha(r)/d(ln(r/r_m))$ plane. The location of a cluster on the $\alpha_G-d\alpha(r)/d(ln(r/r_m))$ plane can therefore constrain its dynamical history and, in particular, constrain possible variations in the stellar initial mass function. The $\alpha_G$-$d\alpha(r)/d(ln(r/r_m))$ plane thus serves as a key tool for fully exploiting the information contained in wide field studies of cluster stellar mass functions.

We describe the remote facilities operated by the Southeastern Association for Research in Astronomy (SARA), a consortium of colleges and universities in the US partnered with Lowell Observatory, the Chilean National Telescope Allocation Committee, and the Instituto de Astrofisica de Canarias. SARA observatories comprise a 0.96m telescope at Kitt Peak, Arizona; a 0.6m instrument on Cerro Tololo, Chile; and the 1m Jacobus Kapteyn Telescope at the Roque de los Muchachos, La Palma, Spain. All are operated using standard VNC or Radmin protocols communicating with on-site PCs. Remote operation offers considerable flexibility in scheduling, allowing long-term observational cadences difficult to achieve with classical observing at remote facilities, as well as obvious travel savings. Multiple observers at different locations can share a telescope for training, educational use, or collaborative research programs. Each telescope has a CCD system for optical imaging, using thermoelectric cooling to avoid the need for frequent local service, and a second CCD for offset guiding. The Arizona and Chile instruments also have fiber-fed echelle spectrographs. Switching between imaging and spectroscopy is very rapid, so a night can easily accommodate mixed observing modes. We present some sample observational programs. For the benefit of other groups organizing similar consortia, we describe the operating structure and principles of SARA, as well as some lessons learned from almost 20 years of remote operations.

Here we report on the results of the WEBT photo-polarimetric campaign targeting the blazar S5~0716+71, organized in March 2014 to monitor the source simultaneously in BVRI and near IR filters. The campaign resulted in an unprecedented dataset spanning $\sim 110$\,h of nearly continuous, multi-band observations, including two sets of densely sampled polarimetric data mainly in R filter. During the campaign, the source displayed pronounced variability with peak-to-peak variations of about $30\%$ and "bluer-when-brighter" spectral evolution, consisting of a day-timescale modulation with superimposed hourlong microflares characterized by $\sim 0.1$\,mag flux changes. We performed an in-depth search for quasi-periodicities in the source light curve; hints for the presence of oscillations on timescales of $\sim 3$\,h and $\sim 5$\,h do not represent highly significant departures from a pure red-noise power spectrum. We observed that, at a certain configuration of the optical polarization angle relative to the positional angle of the innermost radio jet in the source, changes in the polarization degree led the total flux variability by about 2\,h; meanwhile, when the relative configuration of the polarization and jet angles altered, no such lag could be noted. The microflaring events, when analyzed as separate pulse emission components, were found to be characterized by a very high polarization degree ($> 30\%$) and polarization angles which differed substantially from the polarization angle of the underlying background component, or from the radio jet positional angle. We discuss the results in the general context of blazar emission and energy dissipation models.

We investigate the shallow increase in globular cluster half-light radii with projected galactocentric distance $R_{gc}$ observed in the giant galaxies M87, NGC 1399, and NGC 5128. To model the trend in each galaxy, we explore the effects of orbital anisotropy and tidally under-filling clusters. While a strong degeneracy exists between the two parameters, we use kinematic studies to help constrain the distance $R_\beta$ beyond which cluster orbits become anisotropic, as well as the distance $R_{f\alpha}$ beyond which clusters are tidally under-filling. For M87 we find $R_\beta > 27$ kpc and $20 < R_{f\alpha} < 40$ kpc and for NGC 1399 $R_\beta > 13$ kpc and $10 < R_{f\alpha} < 30$ kpc. The connection of $R_{f\alpha}$ with each galaxy's mass profile indicates the relationship between size and $R_{gc}$ may be imposed at formation, with only inner clusters being tidally affected. The best fitted models suggest the dynamical histories of brightest cluster galaxies yield similar present-day distributions of cluster properties. For NGC 5128, the central giant in a small galaxy group, we find $R_\beta > 5$ kpc and $R_{f\alpha} > 30$ kpc. While we cannot rule out a dependence on $R_{gc}$, NGC 5128 is well fitted by a tidally filling cluster population with an isotropic distribution of orbits, suggesting it may have formed via an initial fast accretion phase. Perturbations from the surrounding environment may also affect a galaxy's orbital anisotropy profile, as outer clusters in M87 and NGC 1399 have primarily radial orbits while outer NGC 5128 clusters remain isotropic.

OJ287 is a quasi-periodic quasar with roughly 12 year optical cycles. It displays prominent outbursts which are predictable in a binary black hole model. The model predicted a major optical outburst in December 2015. We found that the outburst did occur within the expected time range, peaking on 2015 December 5 at magnitude 12.9 in the optical R-band. Based on Swift/XRT satellite measurements and optical polarization data, we find that it included a major thermal component. Its timing provides an accurate estimate for the spin of the primary black hole, chi = 0.313 +- 0.01. The present outburst also confirms the established general relativistic properties of the system such as the loss of orbital energy to gravitational radiation at the 2 % accuracy level and it opens up the possibility of testing the black hole no-hair theorem with a 10 % accuracy during the present decade.

We present an analysis of 23 absorption systems along the lines of sight towards 18 quasars in the redshift range of $0.4 \leq z_{abs} \leq 2.3$ observed on the Very Large Telescope (VLT) using the Ultraviolet and Visual Echelle Spectrograph (UVES). Considering both statistical and systematic error contributions we find a robust estimate of the weighted mean deviation of the fine-structure constant from its current, laboratory value of $\Delta\alpha/\alpha=\left(0.22\pm0.23\right)\times10^{-5}$, consistent with the dipole variation reported in Webb et al. and King et al. This paper also examines modelling methodologies and systematic effects. In particular we focus on the consequences of fitting quasar absorption systems with too few absorbing components and of selectively fitting only the stronger components in an absorption complex. We show that using insufficient continuum regions around an absorption complex causes a significant increase in the scatter of a sample of $\Delta\alpha/\alpha$ measurements, thus unnecessarily reducing the overall precision. We further show that fitting absorption systems with too few velocity components also results in a significant increase in the scatter of $\Delta\alpha/\alpha$ measurements, and in addition causes $\Delta\alpha/\alpha$ error estimates to be systematically underestimated. These results thus identify some of the potential pitfalls in analysis techniques and provide a guide for future analyses.

We perform $N$-body simulations of star clusters in time-dependant galactic potentials. Since the Milky Way was built-up through mergers with dwarf galaxies, its globular cluster population is made up of clusters formed both during the initial collapse of the Galaxy and in dwarf galaxies that were later accreted. Throughout a dwarf-Milky Way merger, dwarf galaxy clusters are subject to a changing galactic potential. Building on our previous work, we investigate how this changing galactic potential affects the evolution of a cluster's half mass radius. In particular, we simulate clusters on circular orbits around a dwarf galaxy that either falls into the Milky Way or evaporates as it orbits the Milky Way. We find that the dynamical evolution of a star cluster is determined by whichever galaxy has the strongest tidal field at the position of the cluster. Thus, clusters entering the Milky Way undergo changes in size as the Milky Way tidal field becomes stronger and that of the dwarf diminishes. We find that ultimately accreted clusters quickly become the same size as a cluster born in the Milky Way on the same orbit. Assuming their initial sizes are similar, clusters born in the Galaxy and those that are accreted cannot be separated based on their current size alone.

We investigate how the Milky Way tidal field can affect the spatial mixing of multiple stellar populations in the globular cluster NGC 6362. We use $N$-body simulations of multiple population clusters on the orbit of this cluster around the Milky Way. Models of the formation of multiple populations in globular clusters predict that the second population should initially be more centrally concentrated than the first. However, NGC 6362 is comprised of two chemically distinct stellar populations having the same radial distribution. We show that the high mass loss rate experienced on this cluster's orbit significantly accelerates the spatial mixing of the two populations expected from two body relaxation. We also find that for a range of initial second population concentrations, cluster masses, tidal filling factors and fraction of first population stars, a cluster with two populations should be mixed when it has lost 70-80 per cent of its initial mass. These results fully account for the complete spatial mixing of NGC 6362, since, based on its shallow present day mass function, independent studies estimate that the cluster has lost 85 per cent of its initial mass.

We use N-body simulations to model the 12 Gyr evolution of a suite of star clusters with identical initial stellar mass functions over a range of initial cluster masses, sizes, and orbits. Our models reproduce the distribution of present-day global stellar mass functions that is observed in the Milky Way globular cluster population. We find that the slope of a star cluster's stellar mass function is strongly correlated with the fraction of mass that the cluster has lost, independent of the cluster's initial mass, and nearly independent of its orbit and initial size. Thus, the mass function - initial mass relation can be used to determine a Galactic cluster's initial total stellar mass, if the initial stellar mass function is known. We apply the mass function - initial mass relation presented here to determine the initial stellar masses of 33 Galactic globular clusters, assuming an universal Kroupa initial mass function. Our study suggests that globular clusters had initial masses that were on average a factor of 4.5 times larger than their present day mass, with three clusters showing evidence for being 10 times more massive at birth.

The occurrence of low-amplitude flux variations in blazars on hourly timescales, commonly known as microvariability, is still a widely debated subject in high-energy astrophysics. Several competing scenarios have been proposed to explain such occurrences, including various jet plasma instabilities leading to the formation of shocks, magnetic reconnection sites, and turbulence. In this letter we present the results of our detailed investigation of a prominent, five-hour-long optical microflare detected during recent WEBT campaign in 2014, March 2-6 targeting the blazar 0716+714. After separating the flaring component from the underlying base emission continuum of the blazar, we find that the microflare is highly polarized, with the polarization degree $\sim (40-60)\%$$\pm (2-10)\%$, and the electric vector position angle $\sim (10 - 20)$deg$\pm (1-8)$deg slightly misaligned with respect to the position angle of the radio jet. The microflare evolution in the $(Q,\,U)$ Stokes parameter space exhibits a looping behavior with a counter-clockwise rotation, meaning polarization degree decreasing with the flux (but higher in the flux decaying phase), and approximately stable polarization angle. The overall very high polarization degree of the flare, its symmetric flux rise and decay profiles, and also its structured evolution in the $Q-U$ plane, all imply that the observed flux variation corresponds to a single emission region characterized by a highly ordered magnetic field. As discussed in the paper, a small-scale but strong shock propagating within the outflow, and compressing a disordered magnetic field component, provides a natural, though not unique, interpretation of our findings.

The primordial deuterium abundance is an important tracer of the fundamental physics taking place during Big Bang Nucleosynthesis. It can be determined from absorption features along the line of sight to distant quasars. The quasar PKS1937-1009 contains two absorptions systems that have been used to measure the primordial deuterium abundance, the lower redshift one being at z_abs = 3.256. New observations of this absorber are of a substantially higher signal-to-noise and thus permit a significantly more robust estimate of the primordial deuterium abundance, leading to a D/H ratio of 2.45+/-0.28 x10^-5. Whilst the precision of the new measurement presented here is below that obtained from the recent cosmological parameter measurements by Planck, our analysis illustrates how a statistical sample obtained using similarly high spectral signal-to-noise can make deuterium a competitive and complementary cosmological parameter estimator and provide an explanation for the scatter seen between some existing deuterium measurements.

In this paper, we discuss the origin of the observed correlation between cluster concentration c and present-day mass function (PDMF) slope {\alpha} reported by De Marchi, Paresce & Pulone. This relation can either be reproduced from universal initial conditions combined with some dynamical mechanism(s) that alter(s) the cluster structure and mass function over time, or it must arise early on in the cluster lifetime, such as during the gas-embedded phase of cluster formation. Using a combination of Monte Carlo and N-body models for globular cluster evolution performed with the MOCCA and NBODY6 codes, respectively, we explore a number of dynamical mechanisms that could affect the observed relation. For the range of initial conditions considered here, our results are consistent with an universal initial binary fraction \approx 10% (which does not, however, preclude 100%) and an universal initial stellar mass function resembling the standard Kroupa distribution. Most of the dispersion observed in the c-{\alpha} relation can be attributed to two-body relaxation and Galactic tides. However, dynamical processes alone could not have reproduced the dispersion in concentration, and we require at least some correlation between the initial concentration and the total cluster mass. We argue that the origin of this trend could be connected to the gas-embedded phase of cluster evolution.

Recently, white dwarf stars have found a new use in the fundamental physics community. Many prospective theories of the fundamental interactions of Nature allow traditional constants, like the fine structure constant $\alpha$, to vary in some way. A study by Berengut et al. (2013) used the Fe/Ni V line measurements made by Preval et al. (2013) from the hot DA white dwarf G191-B2B, in an attempt to detect any variation in $\alpha$. It was found that the Fe V lines indicated an increasing alpha, whereas the Ni V lines indicated a decreasing alpha. Possible explanations for this could be misidentification of the lines, inaccurate atomic data, or wavelength dependent distortion in the spectrum. We examine the first two cases by using a high S/N reference spectrum from the hot sdO BD+28$^{\circ}$4211 to calibrate the Fe/Ni V atomic data. With this new data, we re-evaluate the work of Berengut et al. (2013) to derive a new constraint on the variation of alpha in a gravitational field.

(abridged) In this paper, we constrain the properties of primordial binary populations in Galactic globular clusters using the MOCCA Monte Carlo code for cluster evolution. Our results are compared to the observations of Milone et al. (2012) using the photometric binary populations as proxies for the true underlying distributions, in order to test the hypothesis that the data are consistent with an universal initial binary fraction near unity and the binary orbital parameter distributions of Kroupa (1995). With the exception of a few possible outliers, we find that the data are to first-order consistent with the universality hypothesis. Specifically, the present-day binary fractions inside the half-mass radius r$_{\rm h}$ can be reproduced assuming either high initial binary fractions near unity with a dominant soft binary component as in the Kroupa distribution combined with high initial densities (10$^4$-10$^6$ M$_{\odot}$ pc$^{-3}$), or low initial binary fractions ($\sim$ 5-10%) with a dominant hard binary component combined with moderate initial densities near their present-day values (10$^2$-10$^3$ M$_{\odot}$ pc$^{-3}$). This apparent degeneracy can be broken using the binary fractions outside r$_{\rm h}$- only high initial binary fractions with a significant soft component combined with high initial densities can contribute to reproducing the observed anti-correlation between the binary fractions outside r$_{\rm h}$ and the total cluster mass. We further illustrate using the simulated present-day binary orbital parameter distributions and the technique introduced in Leigh et al. (2012) that the relative fractions of hard and soft binaries can be used to further constrain the initial cluster density and mass-density relation. Our results favour an initial mass-density relation of the form r$_{\rm h} \propto$ M$_{\rm clus}^{\alpha}$ with $\alpha <$ 1/3.

We perform N-body simulations of a cluster that forms in a dwarf galaxy and is then accreted by the Milky Way to investigate how a cluster's structure is affected by a galaxy merger. We find that the cluster's half mass radius will respond quickly to this change in potential. When the cluster is placed on an orbit in the Milky Way with a stronger tidal field the cluster experiences a sharp decrease in size in response to increased tidal forces. Conversely, when placed on an orbit with a weaker tidal field the cluster expands since tidal forces decrease and stars moving outwards due to internal effects remain bound at further distances than before. In all cases, we find that the cluster's half mass radius will eventually be indistinguishable from a cluster that has always lived in the Milky Way on that orbit. These adjustments occur within 1-2 half mass relaxation times of the cluster in the dwarf galaxy. We also find this effect to be qualitatively independent of the time that the cluster is taken from the dwarf galaxy. In contrast to the half mass radius, we show the core radius of the cluster is not affected by the potential the cluster lives in. Our work suggests that structural properties of accreted clusters are not distinct from clusters born in the Milky Way. Other cluster properties, such as metallicity and horizontal branch morphology, may be the only way to identify accreted star clusters in the Milky Way.

We have performed N-body simulations of tidally filling star clusters with a range of orbits in a Milky Way-like potential to study the effects of orbital inclination and eccentricity on their structure and evolution. At small galactocentric distances Rgc, a non-zero inclination results in increased mass loss rates. Tidal heating and disk shocking, the latter sometimes consisting of two shocking events as the cluster moves towards and away from the disk, help remove stars from the cluster. Clusters with inclined orbits at large Rgc have decreased mass loss rates than the non-inclined case, since the strength the disk potential decreases with Rgc. Clusters with inclined and eccentric orbits experience increased tidal heating due to a constantly changing potential, weaker disk shocks since passages occur at higher Rgc, and an additional tidal shock at perigalacticon. The effects of orbital inclination decrease with orbital eccentricity, as a highly eccentric cluster spends the majority of its lifetime at a large Rgc. The limiting radii of clusters with inclined orbits are best represented by the rt of the cluster when at its maximum height above the disk, where the cluster spends the majority of its lifetime and the rate of change in rt is a minimum. Conversely, the effective radius is independent of inclination in all cases.

We use N-body simulations to explore the influence of orbital eccentricity on the dynamical evolution of star clusters. Specifically we compare the mass loss rate, velocity dispersion, relaxation time, and the mass function of star clusters on circular and eccentric orbits. For a given perigalactic distance, increasing orbital eccentricity slows the dynamical evolution of a cluster due to a weaker mean tidal field. However, we find that perigalactic passes and tidal heating due to an eccentric orbit can partially compensate for the decreased mean tidal field by energizing stars to higher velocities and stripping additional stars from the cluster, accelerating the relaxation process. We find that the corresponding circular orbit which best describes the evolution of a cluster on an eccentric orbit is much less than its semi-major axis or time averaged galactocentric distance. Since clusters spend the majority of their lifetimes near apogalacticon, the properties of clusters which appear very dynamically evolved for a given galactocentric distance can be explained by an eccentric orbit. Additionally we find that the evolution of the slope of the mass function within the core radius is roughly orbit-independent, so it could place additional constraints on the initial mass and initial size of globular clusters with solved orbits. We use our results to demonstrate how the orbit of Milky Way globular clusters can be constrained given standard observable parameters like galactocentric distance and the slope of the mass function. We then place constraints on the unsolved orbits of NGC 1261,NGC 6352, NGC 6496, and NGC 6304 based on their positions and mass functions.

We present new Hubble Space Telescope imaging of the outer regions of M87 in order to study its globular cluster (GC) population out to large galactocentric distances. We discuss particularly the relationship between GC effective radii $r_h$ and projected galactocentric distance $R_{gc}$. The observations suggest a shallow trend $r_h \propto R_{gc}^{0.14}$ out to $R_{gc} \sim 100$ kpc, in agreement with studies of other giant elliptical galaxies. To theoretically reproduce this relationship we simulate GC populations with various distributions of orbits. For an isotropic distribution of cluster orbits we find a steeper trend of $r_h \propto R_{gc}^{0.4}$. Instead we suggest that (a) if the cluster system has an orbital anisotropy profile, where orbits become preferentially radial with increasing galactocentric distance, AND (b) if clusters become more tidally under-filling with galactocentric distance, the observed relationship can be recovered. We also apply this approach to the red and blue GC populations separately and predict that red clusters are preferentially under-filling at large $R_{gc}$ and have a more isotropic distribution of orbits than blue clusters.

Context. The international whole earth blazar telescope (WEBT) consortium planned and carried out three days of intensive micro-variability observations of S5 0716+714 from February 22, 2009 to February 25, 2009. This object was chosen due to its bright apparent magnitude range, its high declination, and its very large duty cycle for micro-variations. Aims. We report here on the long continuous optical micro-variability light curve of 0716+714 obtained during the multi-site observing campaign during which the Blazar showed almost constant variability over a 0.5 magnitude range. The resulting light curve is presented here for the first time. Observations from participating observatories were corrected for instrumental differences and combined to construct the overall smoothed light curve. Methods. Thirty-six observatories in sixteen countries participated in this continuous monitoring program and twenty of them submitted data for compilation into a continuous light curve. The light curve was analyzed using several techniques including Fourier transform, Wavelet and noise analysis techniques. Those results led us to model the light curve by attributing the variations to a series of synchrotron pulses. Results. We have interpreted the observed microvariations in this extended light curve in terms of a new model consisting of individual stochastic pulses due to cells in a turbulent jet which are energized by a passing shock and cool by means of synchrotron emission. We obtained an excellent fit to the 72-hour light curve with the synchrotron pulse model.

We propose a new probe of the dependence of the fine structure constant, alpha, on a strong gravitational field using metal lines in the spectra of white dwarf stars. Comparison of laboratory spectra with far-UV astronomical spectra from the white dwarf star G191-B2B recorded by the Hubble Space Telescope Imaging Spectrograph gives limits on the fractional variation of alpha of (Delta alpha/alpha)=(4.2 +- 1.6)x10^(-5) and (-6.1 +- 5.8)x10^(-5) from Fe V and Ni V spectra, respectively, at a dimensionless gravitational potential relative to Earth of (Delta phi) ~ 5x10^(-5). With better determinations of the laboratory wavelengths of the lines employed these results could be improved by up to two orders of magnitude.

We report the detection of a strong HI 21-cm absorption system at z=0.5344, as well as a candidate system at z=0.3389, in the sight-line towards the z=2.64 quasar MG J0414+0534. This, in addition to the absorption at the host redshift and the other two intervening absorbers, takes the total to four (possibly five). The previous maximum number of 21-cm absorbers detected along a single sight-line is two and so we suspect that this number of gas-rich absorbers is in some way related to the very red colour of the background source. Despite this, no molecular gas (through OH absorption) has yet been detected at any of the 21-cm redshifts, although, from the population of 21-cm absorbers as a whole, there is evidence for a weak correlation between the atomic line strength and the optical--near-infrared colour. In either case, the fact that so many gas-rich galaxies (likely to be damped Lyman-alpha absorption systems) have been found along a single sight-line towards a highly obscured source may have far reaching implications for the population of faint galaxies not detected in optical surveys, a possibility which could be addressed through future wide-field absorption line surveys with the Square Kilometre Array.

We have performed N-body simulations of star clusters orbiting in a spherically symmetric smooth galactic potential. The model clusters cover a range of initial half-mass radii and orbital eccentricities in order to test the historical assumption that the tidal radius of a cluster is imposed at perigalacticon. The traditional assumption for globular clusters is that since the internal relaxation time is larger than its orbital period, the cluster is tidally stripped at perigalacticon. Instead, our simulations show that a cluster with an eccentric orbit does not need to fully relax in order to expand. After a perigalactic pass, a cluster re-captures previously unbound stars, and the tidal shock at perigalacticon has the effect of energizing inner region stars to larger orbits. Therefore, instead of the limiting radius being imposed at perigalacticon, it more nearly traces the instantaneous tidal radius of the cluster at any point in the orbit. We present a numerical correction factor to theoretical tidal radii calculated at perigalacticon which takes into consideration both the orbital eccentricity and current orbital phase of the cluster.

Metal-rich (red) globular clusters in massive galaxies are, on average, smaller than metal-poor (blue) globular clusters. One of the possible explanations for this phenomenon is that the two populations of clusters have different spatial distributions. We test this idea by comparing clusters observed in unusually deep, high signal-to-noise images of M87 with a simulated globular cluster population in which the red and blue clusters have different spatial distributions, matching the observations. We compare the overall distribution of cluster effective radii as well as the relationship between effective radius and galactocentric distance for both the observed and simulated red and blue subpopulations. We find that the different spatial distributions does not produce a significant size difference between the red and blue subpopulations as a whole, or at a given galactocentric distance. These results suggest that the size difference between red and blue globular clusters is likely due to differences during formation or later evolution

We report results from a 1 week multi-wavelength campaign to monitor the BL Lac object S5 0716+714 (on December 9-16, 2009). In the radio bands the source shows rapid (~ (0.5-1.5) day) intra-day variability with peak amplitudes of up to ~ 10 %. The variability at 2.8 cm leads by about 1 day the variability at 6 cm and 11 cm. This time lag and more rapid variations suggests an intrinsic contribution to the source's intraday variability at 2.8 cm, while at 6 cm and 11 cm interstellar scintillation (ISS) seems to predominate. Large and quasi-sinusoidal variations of ~ 0.8 mag were detected in the V, R and I-bands. The X-ray data (0.2-10 keV) do not reveal significant variability on a 4 day time scale, favoring reprocessed inverse-Compton over synchrotron radiation in this band. The characteristic variability time scales in radio and optical bands are similar. A quasi-periodic variation (QPO) of 0.9 - 1.1 days in the optical data may be present, but if so it is marginal and limited to 2.2 cycles. Cross-correlations between radio and optical are discussed. The lack of a strong radio-optical correlation indicates different physical causes of variability (ISS at long radio wavelengths, source intrinsic origin in the optical), and is consistent with a high jet opacity and a compact synchrotron component peaking at ~= 100 GHz in an ongoing very prominent flux density outburst. For the campaign period, we construct a quasi-simultaneous spectral energy distribution (SED), including gamma-ray data from the FERMI satellite. We obtain lower limits for the relativistic Doppler-boosting of delta >= 12-26, which for a BL\,Lac type object, is remarkably high.

(abridged) We present a new analysis of a large sample of quasar absorption-line spectra obtained using UVES (the Ultraviolet and Visual Echelle Spectrograph) on the VLT (Very Large Telescope) in Chile. In the VLT sample (154 absorbers), we find evidence that alpha increases with increasing cosmological distance from Earth. However, as previously shown, the Keck sample (141 absorbers) provided evidence for a smaller alpha in the distant absorption clouds. Upon combining the samples an apparent variation of alpha across the sky emerges which is well represented by an angular dipole model pointing in the direction RA=(17.3 +/- 1.0) hr, dec. = (-61 +/- 10) deg, with amplitude (0.97 +0.22/-0.20) x 10^(-5). The dipole model is required at the 4.1 sigma statistical significance level over a simple monopole model where alpha is the same across the sky (but possibly different to the current laboratory value). The data sets reveal a number of remarkable consistencies: various data cuts are consistent and there is consistency in the overlap region of the Keck and VLT samples. Assuming a dipole-only (i.e. no-monopole) model whose amplitude grows proportionally with `lookback-time distance' (r=ct, where t is the lookback time), the amplitude is (1.1 +/- 0.2) x 10^(-6) GLyr^(-1) and the model is significant at the 4.2 sigma confidence level over the null model [Delta alpha]/alpha = 0). We apply robustness checks and demonstrate that the dipole effect does not originate from a small subset of the absorbers or spectra. We present an analysis of systematic effects, and are unable to identify any single systematic effect which can emulate the observed variation in alpha.

We present the results of a z>2.9 survey for HI 21-cm and molecular absorption in the hosts of radio quasars using the GMRT and the Tidbinbilla 70-m telescope. Previously published searches, which are overwhelmingly at redshifts of z<1, exhibit a 42% detection rate (31 out of 73 sources), and the inclusion of our survey yields a 17% detection rate (2 out of 12 sources) at z>2.5. We therefore believe that our high redshift selection is responsible for our exclusive non-detections, and find that at ultra-violet luminosities of >10e23 W/Hz, 21-cm absorption has never been detected. We also find this to not only apply to our targets, but also those at low redshift exhibiting similar luminosities, giving zero detections out of a total of 16 sources over z=0.24 to 3.8. This is in contrast to the < 10e23 W/Hz sources where there is a near 50% detection rate of 21-cm absorption. The mix of 21-cm detections and non-detections is currently attributed to orientation effects, where according to unified schemes of active galactic nuclei, 21-cm absorption is more likely to occur in sources designated as radio galaxies (type-2 objects, where the nucleus is viewed through dense obscuring circumnuclear gas) than in quasars(type-1 objects, where we have a direct view to the nucleus). However, due to the exclusively high ultra-violet luminosities of our targets it is not clear whether orientation effects alone can wholly account for the distribution, although there exists the possibility that the large luminosities are indicative of a changing demographic of galaxy types. We also find that below luminosities of ~10e23 W/Hz, both type-1 and type-2 objects have a 50% likelihood of exhibiting 21-cm absorption.

Globular clusters have linear sizes (tidal radii) which theory tells us are determined by their masses and by the gravitational potential of their host galaxy. To explore the relationship between observed and expected radii, we utilize the globular cluster population of the Virgo giant M87. Unusually deep, high signal-to-noise images of M87 are used to measure the effective and limiting radii of approximately 2000 globular clusters. To compare with these observations, we simulate a globular cluster population that has the same characteristics as the observed M87 cluster population. Placing these simulated clusters in the well-studied tidal field of M87, the orbit of each cluster is solved and the theoretical tidal radius of each cluster is determined. We compare the predicted relationship between cluster size and projected galactocentric distance to observations. We find that for an isotropic distribution of cluster velocities, theoretical tidal radii are approximately equal to observed limiting radii for Rgc < 10 kpc. However, the isotropic simulation predicts a steep increase in cluster size at larger radii, which is not observed in large galaxies beyond the Milky Way. To minimize the discrepancy between theory and observations, we explore the effects of orbital anisotropy on cluster sizes, and suggest a possible orbital anisotropy profile for M87 which yields a better match between theory and observations. Finally, we suggest future studies which will establish a stronger link between theoretical tidal radii and observed radii.

We previously reported Keck telescope observations suggesting a smaller value of the fine structure constant, alpha, at high redshift. New Very Large Telescope (VLT) data, probing a different direction in the universe, shows an inverse evolution; alpha increases at high redshift. Although the pattern could be due to as yet undetected systematic effects, with the systematics as presently understood the combined dataset fits a spatial dipole, significant at the 4.2-sigma level, in the direction right ascension 17.5 +/- 0.9 hours, declination -58 +/- 9 degrees. The independent VLT and Keck samples give consistent dipole directions and amplitudes, as do high and low redshift samples. A search for systematics, using observations duplicated at both telescopes, reveals none so far which emulate this result.

We compare the redshifts of neutral carbon and carbon monoxide in the redshifted sources in which the fine structure transition of neutral carbon, [CI], has been detected, in order to measure space-time variation of the fundamental constants. Comparison with the CO rotational lines measures gives the same combination of constants obtained from the comparison fine structure line of singly ionised carbon, [CII]. However, neutral carbon has the distinct advantage that it may be spatially coincident with the carbon monoxide, whereas [CII] could be located in the diffuse medium between molecular clouds, and so any comparison with CO could be dominated by intrinsic velocity differences. Using [CI], we obtain a mean variation of dF/F = (-3.6 +/- 8.5) x 10^-5, over z = 2.3 - 4.1, for the eight [CI] systems, which degrades to (-1.5+/- 11) x 10^-5, over z = 2.3 - 6.4 when the two [CII] systems are included. That is, zero variation over look-back times of 10.8-12.8 Gyr. However, the latest optical results indicate a spatial variation in alpha, which describes a dipole and we see the same direction in dF/F. This trend is, however, due to a single source for which the [CI] spectrum is of poor quality. This also applies to one of the two [CII] spectra previously used to find a zero variation in alpha^2/mu. Quantifying this, we find an anti-correlation between |dF/F| and the quality of the carbon detection, as measured by the spectral resolution, indicating that the typical values of >50 km/s, used to obtain a detection, are too coarse to reliably measure changes in the constants. From the fluxes of the known z > 1 CO systems, we predict that current instruments are incapable of the sensitivities required to measure changes in the constants through the comparison of CO and carbon lines. We therefore discuss in detail the use of ALMA for such an undertaking ... ABRIDGED

Molecular hydrogen transitions in quasar spectra can be used to constrain variation in the proton-to-electron mass ratio, $\mu\equiv m_p/m_e$, at high redshifts ($z\gtrsim 2$). We present here an analysis of a new spectrum of the quasar Q0528$-$250, obtained on VLT/UVES (the Ultraviolet and Visual Echelle Spectrograph, on the Very Large Telescope), and analyse the well-known H$_2$ absorber at $z=2.811$ in this spectrum. For the first time we detect HD (deuterated molecular hydrogen) in this system with a column density of $\log_{10}(N/\mathrm{cm^{-2}})=13.27 \pm 0.07$; HD is sensitive to variation in $\mu$, and so we include it in our analysis. Using 76 H$_2$ and 7 HD transitions we constrain variation in $\mu$ from the current laboratory value to be $\Delta\mu/\mu = (0.3\pm 3.2_\mathrm{stat} \pm 1.9_\mathrm{sys})\times 10^{-6}$, which is consistent with no cosmological variation in $\mu$, as well as with previous results from other H$_2$/HD absorbers. The main sources of systematic uncertainty relate to accurate wavelength calibration of the spectra and the re-dispersion of multiple telescope exposures onto the one pixel grid.

We have undertaken a search for millimetre-wave band absorption (through the CO and HCO+ rotational transitions) in the host galaxies of reddened radio sources (z = 0.405-1.802). Despite the colour selection (optical-near infrared colours of V - K > 5 in all but one source), no absorption was found in any of the eight quasars for which the background continuum flux was detected. On the basis of the previous (mostly intervening) H2 and OH detections, the limits reached here and in some previous surveys should be deep enough to detect molecular absorption according to their V - K colours. However, our survey makes the assumption that the reddening is associated with dust close to the emission redshift of the quasar and that the narrow millimetre component of this emission is intercepted by the compact molecular cores. By using the known millimetre absorbers to define the colour depth and comparing this with the ultra-violet luminosities of the sources, we find that, even if these assumptions are valid, only twelve of the forty objects (mainly from this work) are potentially detectable. This is assuming an excitation temperature of 10 K at z=0, with the number decreasing with increasing temperatures (to zero detectable at 100 K).

We present the results of a mini-survey for associated HI 21-cm absorption at z < 0.42 with the Giant Metrewave Radio Telescope. Our targets are radio galaxies, selected on the basis that the 1216 Angstrom luminosities are below 10^23 W/Hz, above which there has never been a detection of 21-cm absorption. Of the three sources for which we obtained good data, two are unclassified active galactic nuclei (AGN) and one is type-2. Being a non-detection, the type-2 object is consistent with our previous result that 21-cm absorption in radio sources is not dictated by unified schemes of AGN. In the case of the detection, the absorption only occurs towards one of the two resolved radio lobes in PKS 1649-062. If the absorption is due to an another intervening galaxy, or cool HI gas in the intergalactic medium, covering only the south-west lobe, then, being at the same redshift, this is likely to be gravitationally bound to the optical object identified as PKS 1649-062. If the absorption is due to an inclined disk centred between the lobes, intervening the SW lobe while being located behind the NE lobe, by assuming that it covers the emission peak at 150 kpc from the nucleus, we estimate a dynamical mass of ~3 x 10^12 solar masses for the disk.

We report the detection of a third HI 21-cm absorber in the sight-line towards the z=2.64 quasar MG J0414+0534 (4C +05.19). In addition to the absorption at the host redshift and in the z=0.96 gravitational lens, we find, through a decimetre-wave spectral scan towards this source, strong absorption at z=0.38. We believe this may be associated with "Object X", an additional feature apparent in the field of the lensing galaxy and lensed images, on the basis of its close proximity to the quasar images and the possible detection of the [OIII] doublet in a published optical spectrum. If real, the strength of the [OIII] emission would suggest the presence of an active galactic nucleus, or a gas-rich galaxy undergoing rapid star formation, either of which is consistent with the strong outflows apparent in the 21-cm spectrum. Although this is the strongest intervening 21-cm absorber yet found, simultaneous observations failed to detect any of the 18-cm OH lines at the 21-cm redshift. This suggests that, as for the lensing galaxy, this is not the primary location of the intervening material responsible for the very red colour of MG J0414+0534.

From a survey for redshifted HI 21-cm and OH 18-cm absorption in the hosts of a sample of radio galaxies and quasars, we detect HI in three of the ten and OH in none of the fourteen sources for which useful data were obtained. As expected from our recent result, all of the 21-cm detections occur in sources with ultra-violet continuum luminosities of L < 10^23 W/Hz. At these "moderate" luminosities, we also obtain four non-detections, although, as confirmed by the equipartition of detections between the type-1 and type-2 objects, this near-50% detection rate cannot be attributed to unified schemes of active galactic nuclei (AGN). All of our detections are at redshifts of z < 0.67, which, in conjunction with our faint source selection, biases against UV luminous objects. The importance of ultra-violet luminosity (over AGN type) in the detection of 21-cm is further supported by the non-detections in the two high redshift z ~ 3.6 - 3.8 radio galaxies, which are both type-2 objects, while having L > 10^23 W/Hz. Our 21-cm detections in combination with those previously published, give a total of eight (associated and intervening) HI absorbing sources searched and undetected in OH. Using the detected 21-cm line strengths to normalise the limits, we find that only two of these eight may have been searched sufficiently deeply in OH, although even these are marginal.

Theories unifying gravity with other interactions suggest spatial and temporal variation of fundamental "constants" in the Universe. A change in the fine structure constant, alpha, could be detected via shifts in the frequencies of atomic transitions in quasar absorption systems. Recent studies using 140 absorption systems from the Keck telescope and 153 from the Very Large Telescope, suggest that alpha varies spatially. That is, in one direction on the sky alpha seems to have been smaller at the time of absorption, while in the opposite direction it seems to have been larger. To continue this study we need accurate laboratory measurements of atomic transition frequencies. The aim of this paper is to provide a compilation of transitions of importance to the search for alpha variation. They are E1 transitions to the ground state in several different atoms and ions, with wavelengths ranging from around 900 - 6000 A, and require an accuracy of better than 10^{-4} A. We discuss isotope shift measurements that are needed in order to resolve systematic effects in the study. The coefficients of sensitivity to alpha-variation (q) are also presented.

The detection of a spatial variation of the fine-structure constant, alpha, based on study of quasar absorption systems has recently been reported. The physics that causes this alpha-variation should have other observable manifestations, and this motivates us to look for complementary astrophysical effects. In this paper we propose a method to test whether spatial variation of fundamental constants existed during the epoch of big bang nucleosynthesis. Using existing measurements of primordial deuterium abundance we find very weak indications that such a signature might exist, but the paucity of measurements precludes any firm conclusion. We also examine existing quasar absorption spectra data that are sensitive to variation of the electron-to-proton mass ratio, mu, and x = (alpha^2 mu g_p) for spatial variation.

The Keck telescope's HIRES spectrograph has previously provided evidence for a smaller fine-structure constant, alpha, compared to the current laboratory value, in a sample of 143 quasar absorption systems: da/a=(-0.57+/-0.11)x10^{-5}. This was based on a variety of metal-ion transitions which, if alpha varies, experience different relative velocity shifts. This result is yet to be robustly contradicted, or confirmed, by measurements on other telescopes and spectrographs; it remains crucial to do so. It is also important to consider new possible instrumental systematic effects which may explain the Keck/HIRES results. Griest et al. (2009, arXiv:0904.4725v1) recently identified distortions in the echelle order wavelength scales of HIRES with typical amplitudes +/-250m/s. Here we investigate the effect such distortions may have had on the Keck/HIRES varying alpha results. We demonstrate that they cause a random effect on da/a from absorber to absorber because the systems are at different redshifts, placing the relevant absorption lines at different positions in different echelle orders. The typical magnitude of the effect on da/a is ~0.4x10^{-5} per absorber which, compared to the median error on da/a in the sample, ~1.9x10^{-5}, is relatively small. Consequently, the weighted mean value changes by less than 0.05x10^{-5} if the corrections we calculate are applied. Nevertheless, we urge caution, particularly for analyses aiming to achieve high precision da/a measurements on individual systems or small samples, that a much more detailed understanding of such intra-order distortions and their dependence on observational parameters is important if they are to be avoided or modelled reliably. [Abridged]

Spectroscopy of redshifted radio absorption of atomic and molecular species provide excellent probes of the cold component of the gas in the early Universe which can be used to address many important issues, such as measuring baryonic content, probing large-scale structure and galaxy evolution, as well as obtaining independent measurements of various combinations of fundamental constants at large look-back times. However, such systems are currently very rare with only 80 detected in HI 21-cm and five in OH and millimetre-band species. Here we summarise the main selection criteria responsible for this and how the next generation of radio telescopes are expected to circumvent these through their wide instantaneous bandwidths and fields-of-view.

Recent attempts to constrain cosmological variation in the fine structure constant, alpha, using quasar absorption lines have yielded two statistical samples which initially appear to be inconsistent. One of these samples was subsequently demonstrated to not pass consistency tests; it appears that the optimisation algorithm used to fit the model to the spectra failed. Nevertheless, the results of the other hinge on the robustness of the spectral fitting program VPFIT, which has been tested through simulation but not through direct exploration of the likelihood function. We present the application of Markov Chain Monte Carlo (MCMC) methods to this problem, and demonstrate that VPFIT produces similar values and uncertainties for (Delta alpha)/(alpha), the fractional change in the fine structure constant, as our MCMC algorithm, and thus that VPFIT is reliable.

We present a strong constraint on variation of the proton-to-electron mass ratio, mu, over cosmological time scales using molecular hydrogen transitions in optical quasar spectra. Using high quality spectra of quasars Q0405-443, Q0347-383 and Q0528-250, variation in mu relative to the present day value is limited to (Delta mu)/mu=(2.6 +/- 3.0) x 10^(-6). We reduce systematic errors compared to previous works by substantially improving the spectral wavelength calibration method and by fitting absorption profiles to the forest of hydrogen Lyman alpha transitions surrounding each H_2 transition. Our results are consistent with no variation, and inconsistent with a previous approx 4 sigma detection of mu variation involving Q0405-443 and Q0347-383. If the results of this work and those suggesting that alpha may be varying are both correct, then this would tend to disfavour certain grand unification models.

Previous surveys for HI 21-cm absorption in z > 0.1 radio galaxies and quasars yield a 40% detection rate, which is attributed to unified schemes of active galactic nuclei (AGN). In this paradigm absorption is only witnessed in (close to) type-2 objects, where the central obscuration is viewed (nearly) edge-on and thus absorbs the rest frame 1420 MHz emission along our sight-line. However, we find this mix of detections and non-detections to only apply at low redshift (z < 1): From a sensitive survey of eight z > 3 radio sources we find no 21-cm absorption, indicating a low abundance of cold neutral gas in (the sight-lines searched in) these objects. Analysing the spectral energy distributions of these sources, we find that our high redshift selection introduces a bias where our sample consists exclusively of quasars with ultra-violet luminosities in excess of 10e23 W/Hz. This may suggest that we have selected a class of particularly UV bright type-1 objects. Whatever the cause, it must also be invoked to explain the non-detections in an equal number of z < 0.7 sources, where we find, for the first time, the same exclusive non-detections at > 10e23 W/Hz. These objects also turn out to be quasars and, from these exclusive high UV luminosity--21-cm non-detections, it is apparent that orientation effects alone cannot account for the mix of 21-cm detections and non-detections at any redshift.

We present a new catalogue of variable stars compiled from data taken for the University of New South Wales Extrasolar Planet Search. From 2004 October to 2007 May, 25 target fields were each observed for 1-4 months, resulting in ~87000 high precision light curves with 1600-4400 data points. We have extracted a total of 850 variable light curves, 659 of which do not have a counterpart in either the General Catalog of Variable Stars, the New Suspected Variables catalogue or the All Sky Automated Survey southern variable star catalogue. The catalogue is detailed here, and includes 142 Algol-type eclipsing binaries, 23 beta Lyrae-type eclipsing binaries, 218 contact eclipsing binaries, 53 RR Lyrae stars, 26 Cepheid stars, 13 rotationally variable active stars, 153 uncategorised pulsating stars with periods <10 d, including delta Scuti stars, and 222 long period variableswith variability on timescales of >10 d. As a general application of variable stars discovered by extrasolar planet transit search projects, we discuss several astrophysical problems which could benefit from carefully selected samples of bright variables. These include: (i) the quest for contact binaries with the smallest mass ratio, which could be used to test theories of binary mergers; (ii) detached eclipsing binaries with pre-main-sequence components, which are important test objects for calibrating stellar evolutionary models; and (iii) RR Lyrae-type pulsating stars exhibiting the Blazhko-effect, which is one of the last great mysteries of pulsating star research.

In their Letter [Phys. Rev. Lett. 92, 121302 (2004)] (also [Astron. Astrophys. 417, 853 (2004)]), Srianand et al. analysed optical spectra of heavy-element species in 23 absorption systems along background quasar sight-lines, reporting limits on relative variations in the fine-structure constant: da/a=(-0.06+/-0.06) x 10^{-5}. Here we demonstrate basic flaws in their analysis, using the same data and absorption profile fits, which led to spurious values of da/a and significantly underestimated uncertainties. We conclude that these data and fits offer no stringent test of previous evidence for a varying alpha. In their Reply (arXiv:0711.1742) to this Comment, Srianand et al. state or argue several points regarding their original analysis and our new analysis. We discuss these points here, dismissing all of them because they are demonstrably incorrect or because they rely on a flawed application of simple statistical arguments.

We critically review the current null results on a varying fine-structure constant, alpha, derived from VLT/UVES quasar absorption spectra, focusing primarily on the many-multiplet analysis of 23 absorbers from which Chand et al. (2004) reported a weighted mean relative variation of da/a=(-0.06+/-0.06) x 10^{-5}. Our analysis of the same reduced data, using the same fits to the absorption profiles, yields very different individual da/a values with uncertainties typically larger by a factor of ~3. We attribute the discrepancies to flawed parameter estimation techniques in the original analysis and demonstrate that the original da/a values were strongly biased towards zero. Were those flaws not present, the input data and spectra should have given a weighted mean of da/a=(-0.44+/-0.16) x 10^{-5}. Although this new value does reflect the input spectra and fits (unchanged from the original work -- only our analysis is different), we do not claim that it supports previous Keck/HIRES evidence for a varying alpha: there remains significant scatter in the individual da/a values which may stem from the overly simplistic profile fits in the original work. Allowing for such additional, unknown random errors by increasing the uncertainties on da/a to match the scatter provides a more conservative weighted mean, da/a=(-0.64+/-0.36) x 10^{-5}. We highlight similar problems in other current UVES constraints on varying alpha and argue that comparison with previous Keck/HIRES results is premature.

We present the results of a survey for 21-cm absorption in four never previously searched damped Lyman-alpha absorption systems (DLAs) with the Westerbork Synthesis Radio Telescope. We add our results to other recent studies in order to address the important issues regarding the detection of cold gas, through 21-cm absorption, in DLAs: Although, due to the DLAs identified with spiral galaxies, there is a mix of spin temperature/covering factor ratios at low redshift, two recent high redshift end points confirm that this ratio does not generally rise over the whole redshift range searched (up to z = 3.39). That is, if the covering factors of many of these galaxies were a factor of >2 smaller than for the spirals, then no significant difference in the spin temperatures between these two classes would be required. Furthermore, although it is difficult to separate the relative contributions of the spin temperature and covering factor, the new results confirm that 21-cm detections tend to occur at low angular diameter distances, where the coverage of a given absorption cross section is maximised. Finally, we also find an apparent 21-cm line strength--MgII equivalent width correlation, which appears to be due to a coupling of the velocity structure between the components that each species traces. That is, the gas seen in 21-cm absorption could be the same as that seen in optical absorption. Combined with the known equivalent width--metallicity relation, this may be manifest as a spin temperature--metallicity anti-correlation, which is non-evolutionary in origin.

We report the detection of \HI 21-cm absorption in the $z=0.96$ early-type lensing galaxy towards MG J0414+0534 with the Green Bank Telescope. The absorption, with total $N_{\rm HI}=1.6 \times 10^{18} (T_{\rm s}/f) {\rm cm}^{-2}$, is resolved into two strong components, probably due to the two strongest lens components, which are separated by 0.4\arcsec. Unlike the other three lenses which have been detected in \HI, J0414+0534 does not exhibit strong OH absorption, giving a OH/\HI column density ratio of $N_{\rm OH}/N_{\rm HI}\lapp10^{-6}$ (for $T_{\rm s}=100$ K, $T_{\rm x}=10$ K and $f_{\rm HI}=f_{\rm OH}=1$). This underabundance of molecular gas may indicate that the extreme optical--near-IR colour ($V-K=10.26$) along the line-of-sight is not due to the lens. We therefore suggest that despite the strong upper limits on molecular absorption at the quasar redshift, as traced by millimetre lines, the extinction occurs primarily in the quasar host galaxy.

We report the discovery of the first high-amplitude delta Scuti star in an eclipsing binary, which we have designated UNSW-V-500. The system is an Algol-type semi-detached eclipsing binary of maximum brightness V = 12.52 mag. A best-fitting solution to the binary light curve and two radial velocity curves is derived using the Wilson-Devinney code. We identify a late A spectral type primary component of mass 1.49+/-0.02 M_sun and a late K spectral type secondary of mass 0.33+/-0.02 M_sun, with an inclination of 86.5+/-1.0 degrees, and a period of 5.3504751+/-0.0000006 d. A Fourier analysis of the residuals from this solution is performed using PERIOD04 to investigate the delta Scuti pulsations. We detect a single pulsation frequency of f_1 = 13.621+/-0.015 c/d, and it appears this is the first overtone radial mode frequency. This system provides the first opportunity to measure the dynamical mass for a star of this variable type; previously, masses have been derived from stellar evolution and pulsation models.

We report the detection of 21-cm absorption at z = 0.656 towards 1622+238 (3C336). The line is very broad with a Full-Width Half Maximum (FWHM) of 235 km/s, giving a velocity integrated optical depth of 2.2 km/s. The centroid of the line is offset from that of the known damped Lyman-alpha absorption (DLA) system by 50 km/s, and if the Lyman-alpha and 21-cm absorption are due to the same gas, we derive a spin temperature of < 60 K, which would be the lowest yet in a DLA. The wide profile, which is over four times wider than that of any other DLA, supports the hypothesis that the hydrogen absorption is occurring either in the disk of a large underluminous spiral or a group of dim unidentified galaxies, associated with the single object which has been optically identified at this redshift.

Echelle spectrographs currently provide some of the most precise and detailed spectra in astronomy, the interpretation of which sometimes depends on the wavelength calibration accuracy. In some applications, such as constraining cosmological variations in the fundamental constants from quasar absorption lines, the wavelength calibration is crucial. Here we detail an algorithm for selecting thorium-argon (ThAr) emission lines for wavelength calibration which incorporates the properties of both a new laboratory wavelength list and the spectrograph of interest. We apply the algorithm to the Very Large Telescope Ultraviolet and Visual Echelle Spectrograph (UVES) and demonstrate a factor of >3 improvement in the wavelength calibration residuals (i.e. random errors) alone. It is also found that UVES spectra calibrated using a previous, widely distributed line-list contain systematic +/-30-75 m/s distortions of the wavelength scale over both short and long wavelength ranges. These distortions have important implications for current UVES constraints on cosmological variations in the fine-structure constant. The induced systematic errors are most severe for Mg/FeII quasar absorbers in the redshift range 1.2 < z < 2.3, with individual absorbers studied by recent authors containing systematic errors up to 4 times larger than quoted statistical errors.

We present an analysis of a new, detached, double-lined eclipsing binary system with K7 Ve components, discovered as part of the University of New South Wales Extrasolar Planet Search. The object is significant in that only 6 other binary systems are known with comparable or lower mass. Such systems offer important tests of mass-radius theoretical models. Follow-up photometry and spectroscopy were obtained with the 40-inch and 2.3m telescopes at SSO respectively. An estimate of the radial velocity amplitude from spectral absorption features, combined with the orbital inclination (83.5 deg) estimated from lightcurve fitting, yielded a total mass of M=(1.041 +/- 0.06)M_sun and component masses of M_A=(0.529 +/- 0.035)M_sun and M_B=(0.512 +/- 0.035)M_sun. The radial velocity amplitude estimated from absorption features (167 +/- 3)kmps was found to be less than the estimate from the H_alpha emission lines (175 +/- 1.5)kmps. The lightcurve fit produced radii of R_A=(0.641 +/- 0.05)R_sun and R_B=(0.608 +/- 0.06)R_sun, and a temperature ratio of T_B/T_A=0.980 +/- 0.015. The apparent magnitude of the binary was estimated to be V=13.9 +/- 0.2. Combined with the spectral type, this gave the distance to the binary as 169 +/- 14 pc. The timing of the secondary eclipse gave a lower limit on the eccentricity of the binary system of 0.0025 +/- 0.0005. This is the most statistically significant non-zero eccentricity found for such a system, possibly suggesting the presence of a third companion.

A correlation between the metallicity, [M/H], and rest-frame MgII equivalent width, EW, is found from 49 DLAs and strong sub-DLAs drawn from the literature over the redshift range 0.2<z_abs<2.6. The correlation is significant at 4.2 sigma and improves to 4.7 sigma when the mild evolution of [M/H] with redshift is taken into account. Even when including only the 26 DLAs (i.e. excluding sub-DLAs) which have Zn metallicities and EW>0.7A, the correlation remains at >3 sigma significance. Since the MgII2796 transition is predominantly saturated in DLAs (which always have EW greater than 0.3A), EW is far more sensitive to the kinematic spread of the metal velocity components across the absorption profile than it is to [M/H]. Thus, the observed [M/H]--EW correlation points to a strong link between the absorber metallicity and the mechanism for producing and dispersing the velocity components. We also note that approximately half of the 13 known molecular hydrogen absorbers have very high EW and very broad velocity structures which show characteristics usually associated with outflows. Follow-up ultraviolet- and blue-sensitive high-resolution spectra of high-EW systems, initially identified in low-resolution spectra, may therefore yield a large number of new H_2 discoveries.

Current analyses of VLT/UVES quasar spectra disagree with the Keck/HIRES evidence for a varying fine-structure constant, alpha. To investigate this we introduce a simple method for calculating the minimum possible uncertainty on Delta(alpha)/alpha for a given quasar absorber. For many absorbers in Chand et al. (2004) and for the single-absorber constraint of Levshakov et al. (2006) the quoted uncertainties are smaller than the minimum allowed by the UVES data. Failure of this basic consistency test prevents reliable comparison of the UVES and HIRES results.

Nine quasar absorption spectra at 21-cm and UV rest-wavelengths are used to estimate possible variations in x=alpha^2 g_p mu, (alpha is the fine structure constant, g_p the proton g-factor and mu=me/mp the electron-to-proton mass ratio). We find <Delta x/x>^weighted_total(=Dxxwt)=(0.63+-0.99) 10^-5 over 0.23~<z_abs~<2.35 (2.7 to 10.5 Gyr, look-back time, t_lb). A linear best fit against t_lb, tied to Delta x/x=0 at z=0, gives (dot x)/x=(-0.6+-1.2) 10^-15 /yr. Our large sample demonstrates that intrinsic line-of-sight velocity differences between the 21-cm and UV absorption redshifts, (on average Delta_vlos~6km/s), with random sign and magnitude in each absorption system, limit our precision. Combining our Delta x/x measurement with absorption-line constraints on alpha-variation yields strong limits on the variation of mu. Our most conservative estimate, obtained by assuming no variations in alpha or g_p is Delta mu/mu(=Dmm)=Dxxwt. If we use only the four high-redshift absorbers in our sample, we obtain Dmm=(0.58+-1.95) 10^-5, which agrees (2sigma) with recent, more direct estimates from two absorption systems containing molecular hydrogen, also at high redshift, and which have hinted at a possible mu-variation, Dmm=(-2.0+-0.6) 10^-5. Our method of constraining Dmm is completely independent from the molecular hydrogen observations. If we include the low-redshift systems, our Dmm result differs significantly from the high-redshift molecular hydrogen results. We detect a dipole variation in mu across the sky, but this model is required by the data at only the 88 per cent confidence level. Clearly, much larger samples of 21-cm and molecular hydrogen absorbers are required to adequately resolve the issue of the variation of mu and x.(Abridged)

We are currently undertaking a large survey for redshifted atomic and molecular absorption ... only one clear and one tentative detection were obtained: HI absorption at z = 0.097 in PKS 1555-140 and OH absorption at z =0.126 in PKS 2300-189, respectively... In order to determine why no clear molecular absorption was detected in any of the 13 sources searched, we investigate the properties of the five redshifted systems currently known to exhibit OH absorption. In four of these, molecules were first detected via millimetre-wave transitions and the flat radio spectra indicate compact background continuum sources, which may suggest a high degree of coverage of the background source by the molecular clouds in the absorber. Furthermore, for these systems we find a relationship between the molecular line strength and red optical--near infrared (V-K) colours, thus supporting the notion that the reddening of these sources is due to dust, which provides an environment conducive to the formation of molecules. Upon comparison with the V-K colours of our sample, this relationship suggests that, presuming the reddening occurs at the host galaxy redshift at least in some of the targets, many of our observations still fall short of the sensitivityrequired to detect OH absorption, although a confirmation of the ``detection'' of OH in 2300-189 could contravene this.

In this paper we investigate the possible reasons why HI 21-cm absorption in damped Lyman-alpha systems (DLAs) has only been detected at low redshift... We suggest that the lack of 21-cm absorption detections at high redshift arises from the fact that these DLAs are at similar angular diameter distances to the background quasars (i.e. the distance ratios are always close to unity): Above z~1.6 the covering factor becomes largely independent of the DLA--QSO distance, making the high redshift absorbers much less effective at covering the background continuum emission. At low redshift, small distance ratios are strongly favoured by the 21-cm detections, whereas large ratios are favoured by the non-detections. This mix of distance ratios gives the observed mix of detections and non-detections at z<1.6.In addition to the predominance of large distance ratios and non-detections at high redshift, this strongly suggests that the observed distribution of 21-cm absorption in DLAs is dominated by geometric effects.

We report on the current status of the University of New South Wales Extrasolar Planet Search project, giving details of the methods we use to obtain millimagnitude precision photometry using the 0.5m Automated Patrol Telescope. We use a novel observing technique to optimally broaden the PSF and thus largely eliminate photometric noise due to intra-pixel sensitivity variations on the CCD. We have observed 8 crowded Galactic fields using this technique during 2003 and 2004. Our analysis of the first of these fields (centred on the open cluster NGC 6633) has yielded 49 variable stars and 4 shallow transit candidates. Follow-up observations of these candidates have identified them as eclipsing binary systems. We use a detailed simulation of our observations to estimate our sensitivity to short-period planets, and to select a new observing strategy to maximise the number of planets detected.

Quasar absorption spectra at 21-cm and UV rest-wavelengths are used to estimate the time variation of x=alpha^2 g_p mu, where alpha is the fine structure constant, g_p the proton g factor, and m_e/m_p=mu the electron/proton mass ratio. Over a redshift range 0.24 ~< z_abs ~< 2.04, <Delta x/x>_total_weighted=(1.17+-1.01) x 10^-5. A linear fit gives d/dt{x}/x=(-1.43+-1.27) x 10^-15 yr^-1. Two previous results on varying alpha yield the strong limits Delta mu/mu=(2.31+-1.03) x 10^-5 and Delta mu/mu=(1.29+-1.01) x 10^-5. Our sample, 8x larger than any previous, provides the first direct estimate of the intrinsic 21-cm and UV velocity differences ~6 km s^-1.

We investigate the practice of assigning high spin temperatures to damped Lyman-alpha absorption systems (DLAs) not detected in HI 21-cm absorption. In particular, Kanekar & Chengaulr (2003) have attributed the mix of 21-cm detections and non-detections in low redshift (z<2.04) DLAs to a mix of spin temperatures, while the non-detections at high redshift were attributed to high spin temperatures. Below z=0.9, where some of the DLA host galaxy morphologies are known, we find that 21-cm absorption is normally detected towards large radio sources when the absorber is known to be associated with a large intermediate (spiral) galaxy. Furthermore, at these redshifts, only one of the six 21-cm non-detections has an optical identification and these DLAs tend to lie along the sight-lines to the largest background radio continuum sources. For these and many of the high redshift DLAs occulting large radio continua, we therefore expect covering factors of less than the assumed/estimated value of unity. This would have the effect of introducing a range of spin temperatures considerably narrower than the current range of >9000 K, while still supporting the hypothesis that the high redshift DLA sample comprises a larger proportion of compact galaxies than the low redshift sample.

We have identified a new Lyman limit absorption system towards PKS1937-1009, with log N(HI)=18.25 +/- 0.02 at z=3.256. It is suitable for measuring D/H, and we find a 68.3% confidence range for D/H of 1.6^{+0.25}_{-0.30} \times 10^{-5}, and a 95.4% range of 1.6^{+0.5}_{-0.4} \times 10^{-5}. The metallicity of the cloud where D/H was measured is low, [Si/H] = -2.0 +/- 0.5. At these metallicities we expect that D/H will be close to the primordial value. Our D/H disagrees at a level of 99.4% with the predicted D/H using the baryon density calculated from the cosmic background radiation measured by WMAP, 2.60^{+0.19}_{-0.17} \times 10^{-5}. Our result also exacerbates the scatter in D/H values around the mean primordial D/H.

Strong constraints on possible variations in fundamental constants can be derived from HI 21-cm and molecular rotational absorption lines observed towards quasars. With the aim of forming a statistical sample of constraints we have begun a program of systematic searches for such absorption systems. Here we describe molecular rotational searches in 25 damped Lyman-alpha systems where, in many cases, we set optical depth limits an order of magnitude better than that required to detect the 4 known redshifted millimeter-wave absorbers. We also discuss the contributory factors in the detectability of HI 21-cm absorption, focusing on possible biases (e.g.low covering factors) in the currently known sample of absorbers and non-detections.

Spectroscopic observations of distant quasars have resulted in the detection of molecular hydrogen in intervening damped Lyman-alpha absorption clouds (DLAs). We use observations compiled from different experimental groups to show that the molecular hydrogen abundance exhibits a dramatic increase over a cosmological time period corresponding to 13% to 24% of the age of the universe. We also tentatively show that the heavy element abundances in the same gas clouds exhibit a faster and more well-defined cosmological evolution compared to the general DLA population over the same time baseline. We argue that this latter point is unsurprising, because the general DLA population arises in a wide variety of galaxy types and environments, and thus a spans broad range of ISM gas-phases and abundances at the same cosmic time. DLAs exhibiting H2 absorption may therefore circumvent this problem, efficiently identifying a narrower class of objects, and provide a more sensitive probe of cosmological chemical evolution.

We have used the Green Bank Telescope (GBT) and Berkeley-Illinois-Maryland Association (BIMA) array to search for redshifted millimetre absorption in a sample of damped Lyman-alpha absorption systems (DLAs). This brings the number of published systems searched from 18 to 30. In 17 cases we reach $3\sigma$ limits of $\tau\leq0.1$, which is a significant improvement over the previous searches and more than sufficient to detect the 4 known redshifted millimetre absorbers ($\tau\gapp1$). While the CO rotational (millimetre) column density limits obtained are weaker than the electronic (optical) limits, they may provide useful limits below the atmospheric cut-off for the Lyman and Werner \MOLH-bands in the UV ($z_{\rm abs}\lapp1.8$). Using a model for the DLA metallicity evolution combined with assumed HCO$^+$/\MOLH ~and CO/\MOLH \~conversion ratios, we use the molecular column density limits to calculate plausible \MOLH ~molecular fraction limits. Finally, we use these results to discuss the feasibility of detecting rotational CO transitions in DLAs with the next generation of large radio telescopes.

The chemical abundances in damped Lyman-alpha systems (DLAs) show more than 2 orders of magnitude variation at a given epoch, possibly because DLAs arise in a wide variety of host galaxies. This could significantly bias estimates of chemical evolution. We explore the possibility that DLAs in which H_2 absorption is detected may trace cosmological chemical evolution more reliably since they may comprise a narrower set of physical conditions. The 9 known H_2 absorption systems support this hypothesis: metallicity exhibits a faster, more well-defined evolution with redshift than in the general DLA population. The dust-depletion factor and, particularly, H_2 molecular fraction also show rapid increases with decreasing redshift. We comment on possible observational selection effects which may bias this evolution. Larger samples of H_2-bearing DLAs are clearly required and may constrain evolution of the UV background and DLA galaxy host type with redshift.

We present evidence for variations in the fine-structure constant from Keck/HIRES spectra of 143 quasar absorption systems over the redshift range 0.2 < z_abs < 4.2. This includes 15 new systems, mostly at high-z (z_abs > 1.8). Our most robust estimate is a weighted mean da/a=(-0.57 +/- 0.11)x10^{-5}. We respond to recent criticisms of the many-multiplet method used to extract these constraints. The most important potential systematic error at low-z is the possibility of very different Mg heavy isotope abundances in the absorption clouds and laboratory: {\it higher} abundances of {25,26}Mg in the absorbers may explain the low-z results. Approximately equal mixes of {24}Mg and {25,26}Mg are required. Observations of Galactic stars generally show {\it lower} {25,26}Mg isotope fractions at the low metallicities typifying the absorbers. Higher values can be achieved with an enhanced population of intermediate mass stars at high redshift, a possibility at odds with observed absorption system element abundances. At present, all observational evidence is consistent with the varying-alpha results. Another promising method to search for variation of fundamental constants involves comparing different atomic clocks. Here we calculate the dependence of nuclear magnetic moments on quark masses and obtain limits on the variation of alpha and m_q/Lambda_QCD from recent atomic clock experiments with hyperfine transitions in H, Rb, Cs, Hg+ and an optical transition in Hg+.

We use Big Bang Nucleosynthesis calculations and light element abundance data to constrain the relative variation of the deuteron binding energy since the universe was a few minutes old, $\delta Q = Q(BBN)-Q(present)$. Two approaches are used, first treating the baryon to photon ratio, $\eta$, as a free parameter, but with the additional freedom of varying $\delta Q$, and second using the WMAP value of $\eta$ and solving only for $\delta Q$. Including varying $Q$ yields a better fit to the observational data than imposing the present day value, rectifying the discrepancy between the $^4He$ abundance and the deuterium and $^7Li$ abundances, {\it and} yields good agreement with the independently determined $\eta_{WMAP}$. The minimal deviation consistent with the data is significant at about the 4-$\sigma$ level; $\delta Q/Q= -0.019 \pm 0.005$. If the primordial $^4$He abundance lies towards the low end of values in the literature, this deviation is even larger and more statistically significant. Taking the light element abundance data at face-value, our result may be interpreted as variation of the dimensionless ratio $X=m_s/\Lambda_{QCD}$ of the strange quark mass and strong scale: $\delta X/X=(1.1 \pm 0.3) \times 10^{-3}$. These results provide a strong motivation for a more thorough exploration of the potential systematic errors in the light element abundance data.

Millimetre-band scans of the frequency space towards optically dim quasars is potentially a highly efficient method for detecting new high redshift molecular absorption systems. Here we describe scans towards 7 quasars over wide bandwidths (up to 23 GHz) with sensitivity limits sufficient to detect the 4 redshifted absorbers already known. With wider frequency bands, highly efficient searches of large numbers of possibly obscured objects will yield many new molecular absorbers.

We report a further analysis of the ratio of deuterium to hydrogen (D/H) using HST spectra of the z=0.701 Lyman limit system towards the QSO PG1718+481. Initial analyses of this absorber found it gave a high D/H value, 1.8 - 3.1 \times 10^{-4} (Webb et al. 1998), inconsistent with several higher redshift measurements. It is thus important to critically examine this measurement. By analysing the velocity widths of the DI, HI and metal lines present in this system, Kirkman et al. (2001) report that the additional absorption in the blue wing of the lya line can not be DI, with a confidence level of 98%. Here we present a more detailed analysis, taking into account possible wavelength shifts between the three sets of HST spectra used in the analysis. We find that the constraints on this system are not as strong as those claimed by Kirkman et al. The discrepancy between the parameters of the blue wing absorption and the parameters expected for DI is marginally worse than 1 sigma. Tytler et al.(1999) commented on the first analysis of Webb et al.(1997,1998), reporting the presence of a contaminating lower redshift Lyman limit system, with log[N(HI)] = 16.7 at z=0.602, which biases the N(HI) estimate for the main system. Here we show that this absorber actually has log[N(HI)] < 14.6 and does not impact on the estimate of N(HI) in the system of interest at z = 0.701. The purpose of the present paper is to highlight important aspects of the analysis which were not explored in previous studies, and hence help refine the methods used in future analyses of D/H in quasar spectra.

[Abridged] We previously presented evidence for a varying fine-structure constant, alpha, in two independent samples of Keck/HIRES QSO spectra. Here we present a detailed many-multiplet analysis of a third Keck/HIRES sample containing 78 absorption systems. We also re-analyse the previous samples, providing a total of 128 absorption systems over the redshift range 0.2<z_abs<3.7. All three samples separately yield consistent, significant values of da/a. The analyses of low- and high-z systems rely on different ions/transitions with very different dependencies on alpha, yet they also give consistent results. We identify additional random errors in 22 high-z systems characterized by transitions with a large dynamic range in apparent optical depth. Increasing the statistical errors on da/a for these systems gives our fiducial result, a weighted mean da/a=(-0.543+/-0.116)x10^-5, representing 4.7-sigma evidence for a smaller weighted mean alpha in the absorption clouds. Assuming that da/a=0 at z_abs=0, the data marginally prefer a linear increase in alpha with time: dota/a=(6.40+/-1.35)x10^-16 yr^-1. The two-point correlation function for alpha is consistent with zero over 0.2-13 Gpc comoving scales and the angular distribution of da/a shows no significant dipolar anisotropy. We therefore have no evidence for spatial variations in da/a. We extend our previous searches for possible systematic errors, identifying atmospheric dispersion and isotopic structure effects as potentially the most significant. However, overall, known systematic errors do not explain the results. Future many-multiplet analyses of QSO spectra from different telescopes and spectrographs will provide a now crucial check on our Keck/HIRES results.

We use the dynamics of a galaxy, set up initially at a constant proper distance from an observer, to derive and illustrate two counter-intuitive general relativistic results. Although the galaxy does gradually join the expansion of the universe (Hubble flow), it does not necessarily recede from us. In particular, in the currently favored cosmological model, which includes a cosmological constant, the galaxy recedes from the observer as it joins the Hubble flow, but in the previously favored cold dark matter model, the galaxy approaches, passes through the observer, and joins the Hubble flow on the opposite side of the sky. We show that this behavior is consistent with the general relativistic idea that space is expanding and is determined by the acceleration of the expansion of the universe -- not a force or drag associated with the expansion itself. We also show that objects at a constant proper distance will have a nonzero redshift; receding galaxies can be blueshifted and approaching galaxies can be redshifted.

We describe initial results of a search for redshifted molecular absorption toward four millimetre-loud, optically faint quasars. A wide frequency bandwidth of up to 23 GHz per quasar was scanned using the Swedish-ESO Sub-millimetre Telescope at La Silla. Using a search list of commonly detected molecules, we obtained nearly complete redshift coverage up to z_abs=5. The sensitivity of our data is adequate to have revealed absorption systems with characteristics similar to those seen in the four known redshifted millimetre-band absorption systems, but none were found. Our frequency-scan technique nevertheless demonstrates the value of wide-band correlator instruments for searches such as these. We suggest that a somewhat larger sample of similar observations should lead to the discovery of new millimetre-band absorption systems.

We report a lambda > 1cm search for rotational molecular absorption towards quasars, now possible with the upgraded Australia Telescope Compact Array (ATCA). The targets were PKS 0201+113, PKS 0336-017 and Q 0537-286, where known damped Lyman-alpha absorption systems (DLAs) could cause redshifted molecular absorption in the 12mm band of the telescope. We place 3-sigma upper limits on any HCO+ 0-1 absorption features of < 30mJy per 3.4km/s channel. The non-detections could be attributed to the inherent low metallicities in DLAs leading to generally low H_2, and thus HCO+, column densities. In general, the detection of molecular rotational transitions in DLAs could be further hindered by a lower than expected CO-to-H_2 conversion ratio, whether due either to photoionization of carbon or its relative under-abundance at high redshift.

We have previously presented 5.7-sigma evidence for a smaller alpha at redshifts 0.2 < z_abs < 3.7 from a sample of 128 Keck/HIRES quasar absorption systems: da/a = (-0.57 +/- 0.10) x 10^{-5}. A non-zero da/a manifests itself as a distinct pattern of shifts in the measured absorption line wavelengths. The statistical error is now small: we do detect small line shifts in the HIRES data. Whether these shifts are due to systematic errors or due to real variation in alpha is now the central question. Here we summarize the two potentially most important systematic effects: atmospheric dispersion and isotopic abundance evolution. Previously, these have been difficult to quantify/model but here we find that neither of them can explain our results. Furthermore, the HIRES spectra themselves contain no evidence for these effects. Independent measurements of da/a with a different telescope and spectrograph are now crucial if we are to rule out or confirm the present evidence for a variable alpha.

We report preliminary results from a third sample of quasar absorption line spectra from the Keck telescope which has been studied to search for any possible variation of the fine structure constant, alpha. This third sample, which is larger than the sum of the two previously published samples, shows the same effect, and also gives, as do the previous two samples, a significant result. The combined sample yields a highly significant effect, da/a = (alpha_z - alpha_0)/alpha_0 = -0.57 +/- 0.10 x 10^{-5}, averaged over the redshift range 0.2 < z < 3.7. We include a brief discussion of small-scale kinematic structure in quasar absorbing clouds. However, kinematics are unlikely to impact significantly on the averaged non-zero da/a above, and we have so far been unable to identify any systematic effect which can explain it. New measurements of quasar spectra obtained using independent instrumentation and telescopes are required to properly check the Keck results.

We present a short review of the current quasar (QSO) absorption line constraints on possible variation of the fine structure constant, alpha = e^2/(hbar*c). Particular attention is paid to recent optical Keck/HIRES spectra of 49 absorption systems which indicate a smaller alpha in the past (Murphy et al. 2001, Webb et al. 2001). Here we present new preliminary results from 128 absorption systems: da/a = (-0.57 +/- 0.10) x 10^{-5} over the redshift range 0.2 < z < 3.7, in agreement with the previous results. Known potential systematic errors cannot explain these results. We compare them with strong `local' constraints and discuss other (radio and millimetre-wave) QSO absorption line constraints on variations in alpha^2 * g_p and alpha^2 * g_p * m_e/m_p (g_p is the proton g-factor and m_e/m_p is the electron/proton mass ratio). Finally, we discuss future efforts to rule out or confirm the current 5.7 sigma optical detection.

We are setting up a new search for transiting extrasolar planets using the 0.5m Automated Patrol Telescope at Siding Spring Observatory, Australia. We will begin regular observations in September 2002. We expect to find ~7 new planets per year.

We have used the SEST 15-metre and Onsala 20-metre telescopes to perform deep (r.m.s. >~ 30 mJy) integrations of various molecular rotational transitions towards damped Lyman-alpha absorption systems (DLAs) known to occult millimetre-loud quasars. We have observed 6 new systems and improved the existing limits for 11 transitions. These limits may be approaching the sensitivities required to detect new systems and we present a small number of candidate systems which we believe warrant further observation.

We present a catalogue of the 322 damped Lyman alpha absorbers taken from the literature. All damped Lyman alpha absorbers are included, with no selection on redshift or quasar magnitude. Of these, 123 are candidates and await confirmation using high resolution spectroscopy. For all 322 objects we catalogue the radio properties of the background quasars, where known. Around 60 quasars have radio flux densities above 0.1 Jy and approximately half of these have optical magnitudes brighter than V = 18. This compilation should prove useful in several areas of extragalactic/cosmological research.

Although the extra-solar planets discovered so far are of the giant, gaseous, type, the increased sensitivity of future surveys will result in the discovery of lower mass planets. The detection of O2 in the atmosphere of a rocky extra-solar planet would be a potential indicator of a life. In this paper we address the specific issue of whether we would be able to detect the O2 A-band absorption feature in the atmosphere of a planet similar to the Earth, if it were in orbit around a nearby star. Our method is empirical, in that we use observations of the Earth's O2 A-band, with a simple geometric modification for a transiting extra-solar planet, allowing for limb-darkening of the host star. We simulate the spectrum of the host star with the superposed O2 A-band absorption of the transiting planet, assuming a spectral resolution of 7 km/s (typical of current echelle spectrographs), for a range of spectral signal-to-noise ratios. The main result is that we could reliably detect the O2 A-band of the transiting planet for host stars with radii 0.3 solar or less. However, using existing instrumentation and 8m telescopes, this requires target M stars with m(V) of approximately 10 or brighter for integration times of about 10 hours or less. The number of such stars over the sky is small. Larger aperture telescopes and/or improved instrumentation efficiency would enable surveys of M stars down to m(V) = 13 and greatly improve the chances of discovering life elsewhere.

We describe the results of a search for time variability of the fine structure constant, alpha, using absorption systems in the spectra of distant quasars. Three large optical datasets and two 21cm/mm absorption systems provide four independent samples, spanning 23% to 87% of the age of the universe. Each sample yields a smaller alpha in the past and the optical sample shows a 4-sigma deviation: da/a = -0.72 +/- 0.18 x 10^{-5} over the redshift range 0.5 < z < 3.5. We find no systematic effects which can explain our results. The only potentially significant systematic effects push da/a towards positive values, i.e. our results would become more significant were we to correct for them.

An experimental search for variation in the fundamental coupling constants is strongly motivated by modern high-energy physics theories. Comparison of quasar absorption line spectra with laboratory spectra provides a sensitive probe for variability of the fine structure constant, alpha, over cosmological time-scales. We have previously developed and applied a new method providing an order of magnitude gain in precision over previous optical astrophysical constraints. Here we extend that work by including new quasar spectra of damped Lyman-alpha absorption systems. We also re-analyse our previous lower redshift data and confirm our initial results. The constraints on alpha come from simultaneous fitting of absorption lines of subsets of the following species: Mg I, Mg II, Al II, Al III, Si II, Cr II, Fe II, Ni II and Zn II. We present a detailed description of our methods and results based on an analysis of 49 quasar absorption systems (towards 28 QSOs) covering the redshift range 0.5 < z < 3.5. There is statistical evidence for a smaller alpha at earlier epochs: da/a = (-0.72 +/- 0.18) * 10^{-5}. The new and original samples are independent but separately yield consistent and significant non-zero values of da/a. We summarise the results of a thorough investigation of systematic effects published in a companion paper. The value we quote above is the raw value, not corrected for any of these systematic effects. The only significant systematic effects so far identified, if removed from our data, would lead to a more significant deviation of da/a from zero.

Comparison of quasar absorption spectra with laboratory spectra allow us to probe possible variations in the fundamental constants over cosmological time-scales. In a companion paper we present an analysis of Keck/HIRES spectra and report possible evidence suggesting that the fine structure constant, alpha, may have been smaller in the past: da/a = (-0.72 +/- 0.18) * 10^{-5} over the redshift range 0.5 < z < 3.5. In this paper we describe a comprehensive investigation into possible systematic effects. Most of these do not significantly influence our results. When we correct for those which do produce a significant systematic effect in the data, the deviation of da/a from zero becomes more significant. We are lead increasingly to the interpretation that alpha was slightly smaller in the past.

Comparison of quasar absorption line spectra with laboratory spectra provides a precise probe for variability of the fine structure constant, alpha, over cosmological time-scales. We constrain variation in alpha in 21 Keck/HIRES Si IV absorption systems using the alkali doublet (AD) method in which changes in alpha are related to changes in the doublet spacing. The precision obtained with the AD method has been increased by a factor of 3: da/a = (-0.5 +/- 1.3) * 10^{-5}. We also analyse potential systematic errors in this result. Finally, we compare the AD method with the many-multiplet method which has achieved an order of magnitude greater precision and we discuss the future of the AD method.

QSO absorption spectra provide an extremely useful probe of possible cosmological variation in various physical constants. Comparison of H I 21cm absorption with corresponding molecular (rotational) absorption spectra allows us to constrain variation in y=alpha^2*g_p where alpha is the fine structure constant and g_p is the proton g-factor. We analyse spectra of two QSOs, PKS 1413+135 and TXS 0218+357, and derive values of dy/y at absorption redshifts of z=0.2467 and 0.6847 by simultaneous fitting of the H I 21cm and molecular lines. We find dy/y=(-0.20 +/- 0.44)*10^{-5} and dy/y=(-0.16 +/- 0.54)*10^{-5} respectively, indicating an insignificantly smaller y in the past. We compare our results with other recent constraints from the same two QSOs (Drinkwater et al. 1998; Carilli et al. 2000) and with our recent optical constraints which indicated a smaller alpha at higher redshifts.

We present an analysis of galaxy and QSO absorption line pairs toward 24 QSOs at redshifts between z~0.2 and 1 in an effort to establish the relationship between galaxies and absorption lines in physical proximity to QSOs. We demonstrate the existence of a galaxy proximity effect, in that galaxies in the vicinities of QSOs do not show the same incidence and extent of gaseous envelopes as galaxies far from QSOs. We show that the galaxy proximity effect exists to galaxy-QSO velocity separations of ~ 3000 km/s, much larger than the size of a typical cluster (~ 1000 km/s), i.e. it is more comparable to the scale of the sphere of influence of QSO ionizing radiation rather than the scale of galaxy-QSO clustering. This indicates that the QSO ionizing radiation rather than some dynamical effect from the cluster environment is responsible for the galaxy proximity effect. We combine previous findings that (1) many or most Lya absorption lines arise in extended galaxy envelopes, and (2) galaxies cluster around QSOs to show that the magnitude of the Lya forest proximity effect is underestimated. Consequently, determinations of the UV ionizing background intensity using the proximity effect are likely overestimated. We use the galaxy-QSO cross-correlation function measured from our data to estimate the magnitude of this overestimate and find that it could be as high as a factor of 20 at z<1. This can have strong implications for models of the origin and evolution of the ionizing background, and may indicate that QSOs produce sufficient ionizing flux at all redshifts to account for the entire background radiation field.

(Abridged) We present results of a program to obtain and analyze HST WFPC2 images and ground-based images of galaxies identified in an imaging and spectroscopic survey of faint galaxies in fields of HST spectroscopic target QSOs. We confirm and improve the results presented by Lanzetta et al. (1995) and Chen et al. (1998) that (1) extended gaseous envelopes are a common and generic feature of galaxies of a wide range of luminosity and morphological type, (2) the extent of tenuous gas (N(H I) > 10^{14} cm^{-2} around galaxies scales with galaxy B-band luminosity as r\propto L_B^{0.39+/-0.09}, and (3) galaxy interactions do not play an important role in distributing tenuous gas around galaxies in most cases. We further demonstrate that (4) the gaseous extent of galaxies scales with galaxy K-band luminosity as r\propto L_K^{0.28+/-0.08}, and (5) tenuous gas around typical L_* galaxies is likely to be distributed in spherical halos of radius ~ 180 h^{-1} kpc of covering factor of nearly unity. Because galaxies of all morphological types possess extended gaseous halos and because the extent of tenuous gas around galaxies scales with galaxy K-band luminosity, we argue that galaxy mass---rather than recent star-formation activity---is likely to be the dominant factor that determines the extent of tenuous gas around galaxies. Nevertheless, applying the scaling relationship between the extent of Lya absorbing gas around galaxies and galaxy B-band luminosity, the results of our analysis also suggest that the number density evolution of Lya absorption systems may serve to constrain the evolution of the comoving galaxy B-band luminosity density (at least for the redshift interval between z~0 and z~1 that has been studied in our survey).

(Abridged) We report the discovery of extended CIV gaseous envelopes around galaxies of a wide range of luminosity and morphological type. First, we show that CIV absorption systems are strongly clustered around galaxies on velocity scales of v < 250 km/s and impact parameter scales of rho < 100 h^{-1} kpc but not on larger velocity or impact parameter scales. Next, adopting measurements of galaxy properties presented in previous papers, we examine how properties of the CIV absorption systems depend on properties of the galaxies. On the basis of 14 galaxy and absorber pairs and 36 galaxies that do not produce corresponding CIV absorption lines to within sensitive upper limits, we find that: (1) Galaxies of a range of morphological type and luminosity appear to possess extended CIV gaseous envelopes of radius R ~ 100 h^{-1} kpc, with abrupt boundaries between the CIV absorbing and non-absorbing regions. (2) The extent of CIV-absorbing gas around galaxies scales with galaxy B-band luminosity as R \propto L_B^{0.5 +/- 0.1} but does not depend strongly on galaxy surface brightness, redshift, or morphological type. And (3) the covering factor of CIV clouds within ~ 100 h^{-1} kpc of galaxies is nearly unity, but there is a large scatter in the mean number of clouds encountered along the line of sight. The most significant implication of the study is that galaxies of a wide range of luminosity and morphological type are surrounded by chemically enriched gas that extends for at least ~ 100 h^{-1} kpc. We consider various scenarios that may have produced metals at large galactic distance and conclude that accreting satellites are most likely to be responsible for chemically enriched gas at large galactic distances to regular looking galaxies.

Theories unifying gravity and other interactions suggest the possibility of spatial and temporal variation of physical ``constants'' in the Universe. Detection of high redshift absorption systems intersecting the sight lines towards distant quasars provide a powerful tool for measuring these variations. In the present paper we demonstrate that high sensitivity to variation of the fine structure constant alpha can be obtained by comparing cosmic and laboratory spectra of the Ni II ion. Relativistic effects in Ni II reveal many interesting features. The Ni II spectrum exhibits avoided level crossing phenomenon under variation of alpha and the intervals between the levels have strong nonlinear dependencies on relativistic corrections. The values of the transition frequency shifts, due to the change of alpha, vary significantly from state to state including change of the sign. This enhances the sensitivity to the variation of alpha and reduces possible systematic errors. The calculations of alpha-dependence of the nickel ion spectral lines that are detectable in quasar absorption spectra have been performed using a relativistic configuration interaction method.

The mechanism for gamma ray bursters and the detection of gravitational waves (GWs) are two outstanding problems facing modern physics. Many models of gamma ray bursters predict copious GW emission, so the assumption of an association between GWs and GRBs may be testable with existing bar GW detector data. We consider Weber bar data streams in the vicinity of known GRB times and present calculations of the expected signal after co-addition of 1000 GW/GRBs that have been shifted to a common zero time. Our calculations are based on assumptions concerning the GW spectrum and the redshift distribution of GW/GRB sources which are consistent with current GW/GRB models. We discuss further possibilities of GW detection associated with GRBs in light of future bar detector improvements and suggest that co-addition of data from several improved bar detectors may result in detection of GWs (if the GW/GRB assumption is correct) on a time scale comparable with the LIGO projects.

The spatial distribution of Ly-alpha forest absorption systems towards ten QSOs has been analysed to search for large-scale structure over the redshift range 2.2 < z < 3.4. The QSOs form a closely spaced group on the sky and are concentrated within a 1 deg^2 field. We have employed a technique based on the first and second moments of the transmission probability density function which is capable of identifying and assessing the significance of regions of over- or underdense Ly-alpha absorption. We find evidence for large-scale structure in the distribution of Ly-alpha forest absorption at the > 99 per cent confidence level. In individual spectra we find overdense Ly-alpha absorption on scales of up to 1200 km s^-1. There is also strong evidence for correlated absorption across line of sight pairs separated by < 3 h^-1 proper Mpc (q_0 = 0.5). For larger separations the cross-correlation signal becomes progressively less significant.

We present new GHRS observations of Lyman alpha absorption lines associated with a group of galaxies towards the QSO 1545+2101. We have identified eight distinct Lyman alpha absorption features in the spectrum of QSO 1545+2101 at a mean redshift of z=0.2648 with a velocity dispersion of 163 km/s. A group of galaxies is detected in the vicinity of this QSO at a mean redshift of z=0.2645 and velocity dispersion 239 km/s. The identification of discrete absorption systems indicates that they arise in clouds of neutral hydrogen rather than in a diffuse intragroup medium. Our analysis suggests that the Lyman alpha absorption lines are associated with individual galaxies in the group, although a one-to-one relationship between absorbers and galaxies is difficult to establish in such a dense environment.

A method offering an order of magnitude sensitivity gain is described for using quasar spectra to investigate possible time or space variation in the fine structure constant, alpha. Applying the technique to a sample of 30 absorption systems, spanning redshifts 0.5 < z< 1.6, obtained with the Keck I telescope, we derive limits on variations in alpha over a wide range of epochs. For the whole sample Delta(alpha)/alpha = -1.1 +/- 0.4 x 10^{-5}. This deviation is dominated by measurements at z > 1, where Delta(alpha)/alpha = -1.9 +/- 0.5 x 10^{-5}. For z < 1, Delta(alpha)/alpha = -0.2 +/- 0.4 x 10^{-5}, consistent with other known constraints. Whilst these results are consistent with a time-varying alpha, further work is required to explore possible systematic errors in the data, although careful searches have so far not revealed any.

Here we describe results of our attempt to determine what types of galaxies are responsible for the Lyman-alpha forest absorption systems, based on our ongoing imaging and spectroscopic survey of faint galaxies in fields of HST spectroscopic target QSOs. Our primary conclusions are that the bulk of the Lyman-alpha forest arises in more or less normal galaxies (that span the normal range of luminosity and morphology) and that any ``unseen'' low surface brightness galaxies are unlikely to contribute significantly to the luminosity density of the universe.

We present results of an ongoing program to study the gaseous extent of galaxies and the origin of QSO absorption line systems. For \lya absorption systems, we find that absorption equivalent width depends strongly on galaxy impact parameter and galaxy $B$-band luminosity, and that the gaseous extent of individual galaxies scales with galaxy $B$-band luminosity as $r\propto L_B^{0.40\pm0.09}$. Applying the results to galaxies in the Hubble Deep Field to calculate the predicted number density of \lya absorption lines as a function of redshift and comparing it with observations, we find that at least 50% and perhaps as much as 100% of \lya absorption systems with $W\apg0.32$ \AA can be explained by the extended gaseous envelops of normal galaxies. The anti-correlation analysis has also been performed on \civ absorption line systems, and the results show that the ionized gas cross section scales with galaxy $B$-band luminosity as $r\propto L_B^{0.76\pm0.26}$.

We present new statistical study on the origin of Lyman-alpha absorption systems that are detected in QSO fields. The primary objective is to investigate environmental dependence of the relationship between galaxies and Lyman-alpha absorption systems. We find that galaxies exhibit no environmental dependence in giving rise to Lyman-alpha absorbers in their halos. Conversely, Lyman-alpha absorbers do not exhibit preference as to the environment where they reside.

We present a new technique for detecting structure on Mpc scales in the Lyman-alpha forest. The technique is easy to apply in practice since it does not involve absorption line fitting but is rather based on the statistics of the transmitted flux. It identifies and assesses the statistical significance of regions of over- or underdense Lyman-alpha absorption and is fairly insensitive to the quality of the data. Using extensive simulations we demonstrate that the new method is significantly more sensitive to the detection of large-scale structure in the Lyman-alpha forest than a traditional two-point correlation function analysis of fitted absorption lines.

We present high spectral resolution HST-GHRS observations of the white dwarf G191-B2B, and derive the interstellar D/H ratio on the line of sight. We have observed and analysed simultaneously the interstellar lines of HI, DI, NI, OI, SiII, and SiIII. We detect three absorbing clouds, and we measure an average (D/H) ratio: N(DI)_Total / N(HI)_Total = 1.12 +/- 0.08 x 10**-5, in disagreement with the value of the D/H ratio as measured toward Capella by Linsky et al. (1993, 1995). We re-analyze the GHRS data of the Capella line of sight, and confirm their estimate, as we find: (D/H)_Capella = 1.56 +/- 0.1 x 10**-5, in the Local Interstellar Cloud in which the Sun is embedded. This shows that the D/H ratio varies by at least 30% within the local interstellar medium (not to be confused with the Local Interstellar Cloud). Furthermore, the Local Interstellar Cloud is also detected toward our target G191-B2B, and we show that the D/H ratio in this component, toward G191-B2B, can be made compatible with that derived toward Capella. However, this comes at the expense of a much smaller value for the D/H ratio as averaged over the other two components, of order 0.9 x 10**-5, and in such a way that the D/H ratio, as averaged over all three components, remains at the above value, i.e. 1.12 x 10**-5. We thus conclude that, either the D/H ratio varies from cloud to cloud, and/or the D/H ratio varies within the Local Interstellar Cloud in which the Sun is embedded, although we note that our observations neither prove nor disprove this latter possibility.

The Hubble Deep Field-South (HDF-S) was chosen to have a QSO (RA 22:33:37.6 Dec -60:33:29 J2000 and B=17.5) in the field to allow for studies of absorption systems intersecting the sight line to the QSO. To assist in the planning of HDF-S observations we present here a ground-based spectrum of the QSO. We measure a redshift of z=2.24 for the QSO and find associated absorption in the spectrum at z=2.204 as well as additional absorption features.

The Hubble Deep Field South (HDF-S) HST observations are expected to begin in October 1998. We present a composite spectrum of the QSO in the HDF-S field covering UV/optical/near IR wavelengths, obtained by combining data from the ANU 2.3m Telescope with STIS on the HST. This intermediate resolution spectrum covers the range 1600-10000A and allows us to derive some basic information on the intervening absorption systems which will be important in planning future higher resolution studies of this QSO.

Recent detections of high-redshift absorption by both atomic hydrogen and molecular gas in the radio spectra of quasars have provided a powerful tool to measure possible temporal and spatial variations of physical `constants' in the Universe. We have compared the frequency of high-redshift Hydrogen 21 cm absorption with that of associated molecular absorption in two quasars to place new (1 sigma) upper limits on any variation in y = g_p.alpha^2 (where alpha is the fine structure constant and g_p is the proton g-factor) of |Delta(y)/y| < 5.10^{-6} at redshifts z=0.25 and z=0.68. These quasars are separated by a comoving distance of 3000 Mpc (H_0=75 km/s/Mpc and q_0=0). We also derive limits on the time rates of change of |(d.g_p/d.t)/g_p| < 10^{-15}/y and |(d.alpha/d.t)/\alpha| < 5.10^{-16}/y between the present epoch and z=0.68. These limits are more than an order of magnitude smaller than previous results derived from high-redshift measurements.

We present initial results of a program to obtain and analyze HST WFPC2 images of galaxies identified in an imaging and spectroscopic survey of faint galaxies in fields of HST spectroscopic target QSOs. We measure properties of 87 galaxies, of which 33 are associated with corresponding \lya absorption systems and 24 do not produce corresponding \lya absorption lines to within sensitive upper limits. Considering only galaxy and absorber pairs that are likely to be physically associated and excluding galaxy and absorber pairs within 3000 \kms of the background QSOs leaves 26 galaxy and absorber pairs and seven galaxies that do not produce corresponding \lya absorption lines to within sensitive upper limits. Redshifts of the galaxy and absorber pairs range from 0.0750 to 0.8912 with a median of 0.3718, and impact parameter separations of the galaxy and absorber pairs range from 12.4 to $157.4 h^{-1}$ kpc with a median of $62.4 h^{-1}$ kpc. The primary result of the analysis is that the amount of gas encountered along the line of sight depends on the galaxy impact parameter and B-band luminosity but does not depend strongly on the galaxy average surface brightness, disk-to-bulge ratio, or redshift. This result confirms and improves upon the anti-correlation between \lya absorption equivalent width and galaxy impact parameter found previously by Lanzetta et al. (1995). There is no evidence that galaxy interactions play an important role in distributing tenuous gas around galaxies in most cases. Galaxies might account for all \lya absorption systems with $W > 0.3$ \AA, but this depends on the unknown luminosity function and gaseous cross sections of low-luminosity galaxies as well as on the uncertainties of the observed number density of \lya absorption systems.

Recent HST/GHRS observations of the z = 0.7010 absorber toward the QSO 1718+4807 (selected because of apparently ideal characteristics for measuring D/H) yield a D/H significantly higher D/H than some recent high-redshift measurements. Our analysis indicates D/H}= 1.8 - 2.5 x 10^{-4}. This may indicate a cosmological inhomogeneity in the deuterium abundance of at least a factor of ten.

We have compared the frequency of high-redshift Hydrogen HI 21 cm absorption with that of associated molecular absorption in two quasars to place new (1 sigma) upper limits on any variation in y = g_p.alpha^2 (where alpha is the fine structure constant and g_p is the proton g-factor) of |Delta(y)/y| < 5.10^{-6} at redshifts z=0.25 and z=0.68.

We describe an analysis of the galaxy-absorber cross-correlation function as a function of velocity separation, impact parameter separation, and absorber Lyman alpha equivalent width on impact parameter scales that extend up to 1 Mpc.

We present initial results of a program to obtain and analyze HST WFPC2 images of galaxies in fields of HST spectroscopic target QSOs. The goal of the program is to investigate how the properties of \lya absorption systems observed in the spectra of background QSOs vary with the properties of intervening galaxies. We found that \lya absorption equivalent width depends strongly on galaxy impact parameter and galaxy B-band luminosity, and that the gaseous extent of individual galaxies scales with galaxy B-band luminosity by $r\propto L_B^{0.35\pm0.10}$.

We present here new GHRS observations of Lyman-alpha absorption lines associated with groups or clusters toward the QSOs 1545+2101 and 0850+4400. In the first case we have identified at least eight distinct Lyman-alpha absorption features, with a mean redshift of <z> = 0.2648 \pm 0.0002 and a velocity dispersion of 163 \pm 57 km/s. We have also identified a group or cluster of galaxies in the vicinity of this QSO with a mean redshift of <z> = 0.2643 \pm 0.0004 and velocity dispersion 223 \pm 91 km/s . The spectrum of QSO 0850+4400, of poorer quality, reveals two Lyman-alpha absorption features at z=0.0909500 \pm 0.0000070 (which is just resolved) and z=0.0948215 \pm 0.0000090, separated by $\sim 1060 km/s. A group or cluster of galaxies is also present in the vicinity of the QSO line-of-sight with a mean redshift <z>=0.0901 \pm 0.0007 and velocity dispersion of 530 \pm 200 km/s. The results of this work establish that Lyman-alpha absorption can occur in denser than average galaxy environments, and that it arises in discrete objects spanning a velocity range similar to that of the cluster galaxies. Although a one-to-one relationship between absorbers and galaxies is difficult to establish in such a dense environment, the results obtained here are indeed consistent with the Lyman-alpha absorption lines being associated with individual galaxies also in groups and clusters. Moreover, the data shows clearly that line clustering takes place in the Lyman-alpha forest.

The spatial distribution of Ly-alpha forest absorption systems toward a group of 8, closely spaced QSOs has been analysed and evidence for large scale structure has been found at <z> = 2.8. Our technique is based on the first and second moments of the transmission probability density function which is capable of identifying and assessing the significance of regions of over- or underdense Ly-alpha absorption. The data has revealed at least two interesting features. 1. An overdense structure at z = 2.27 which extends at least over ~8 h^{-1} comoving Mpc (q_0 = 0.5) in the plane of the sky. Metal absorption lines have been found at the same redshift and thus a cluster or proto-cluster of galaxies seems to have been discovered. 2. A void at z = 2.97, extending over ~20 h^{-1} comoving Mpc in the plane of the sky, possibly caused by a locally increased UV ionising flux due to a foreground QSO.

Of the light elements, the primordial abundance of deuterium, (D/H)_p, provides the most sensitive diagnostic for the cosmological mass density parameter Omega_B. Recent high redshift (D/H) measurements are highly discrepant, although this may reflect observational uncertainties. The larger (D/H) values, which imply a low Omega_B and require the Universe to be dominated by non-baryonic matter (dynamical studies indicate a higher total density parameter), cause problems for galactic chemical evolution models since they have difficulty in reproducing the large decline down to the lower present-day (D/H). Conversely, low (D/H) values imply an Omega_B greater than derived from ^7Li and ^4He abundance measurements, and may require a deuterium abundance evolution that is too low to easily explain. Here we report the first measurement at intermediate redshift, where the observational difficulties are smaller, of a gas cloud with ideal characteristics for this experiment. Our analysis of the z = 0.7010 absorber toward 1718+4807 indicates (D/H) = 2.0 +/- 0.5 x 10^{-4} which is in the high range. This and other independent observations suggests there may be a cosmological inhomogeneity in (D/H)_p of at least a factor of ten.

We report new HST and ground-based observations of a damped Lyman-alpha absorption system toward the QSO 0850+4400. The redshift of the absorption system is z = 0.163770 and the neutral hydrogen column density of the absorption system is log N = 19.81 cm**-2. The absorption system is by far the lowest redshift confirmed damped Lyman-alpha absorption system yet identified, which provides an unprecedented opportunity to examine the nature, impact geometry, and kinematics of the absorbing galaxy in great detail. The observations indicate that the absorption system is remarkable in three respects: First, the absorption system is characterized by weak metal absorption lines and a low metal abundance, possibly less than 4% of the solar metal abundance. This cannot be explained as a consequence of dust, because the neutral hydrogen column density of the absorption system is far too low for obscuration by dust to introduce any significant selection effects. Second, the absorption system is associated with a moderate-luminosity early-type S0 galaxy, although the absorption may actually arise in one of several very faint galaxies detected very close to the QSO line of sight. Third, the absorbing material moves counter to the rotating galaxy disk, which rules out the possibility that the absorption arises in a thin or thick co-rotating gaseous disk. These results run contrary to the expectation that low-redshift damped Lyman-alpha absorption systems generally arise in the gas- and metal-rich inner parts of late-type spiral galaxies. We suggest instead that mounting evidence indicates that low-redshift galaxies of a variety of morphological types may contain significant quantities of low metal abundance gas at large galactocentric distances.

We use new observations of very weak CIV absorption lines associated with high-redshift Lyman-alpha absorption systems to measure the high-redshift Lyman-alpha line two-point correlation function (TPCF). These very weak CIV absorption lines trace small-scale velocity structure that cannot be resolved by Lyman-alpha absorption lines. We find that (1) high-redshift Lyman-alpha absorption systems with N(HI) > 3.10^14 cm^{-2} are strongly clustered in redshift, (2) previous measurements of the Lyman-alpha line TPCF underestimated the actual clustering of the absorbers due to unresolved blending of overlapping velocity components, (3) the present observations are consistent with the hypothesis that clustering of Lyman-alpha absorption systems extends to lower column densities, but maybe with smaller amplitude in the correlation function, and (4) the observed clustering is broadly compatible with that expected for galaxies at z \sim 2-3. We interpret these results as suggesting that many or most Lyman-alpha absorbers may arise in galaxies even at high redshifts, and, therefore, that the Lyman-alpha forest probes processes of galaxy formation and evolution for redshifts \lesssim 5.

The influence of a foreground QSO on the lyman a forest of another QSO with higher redshift has been investigated by analyzing the spectra of three such objects at redshifts $z=2 - 2.7$. This influence is not contaminated by any projection effects, as opposed to the inverse effect along the line of sight, where incomplete coverage of the QSO continuum emitting region by the lyman a clouds could contribute to the relative lack of lines. Our results are consistent with the existence of a proximity effect due to the foreground QSO, but due to its weakness we can only reject the absence of such effect at 1 sigma level. By modelling this proximity effect assuming that lyman a clouds are low-density highly ionized objects we find that the best value for the UV photoinizing intensity at those redshifts is 10^{-20.5}(...) at the Lyman limit which is consistent with previous estimates of the background based on the inverse effect. We also find an absolute lower limit (at 95% confidence) to the UV intensity at the level of 10^{-21.8}(...) which means the rejection of a number of models for the UV background.

We present the first moderate resolution (approx. 40--120 km/s) spectroscopic observations of the bright (V<16.7) high-redshift QSOs Q1107+487 (z_{em} = 2.965) and Q1442+295 (z_{em} = 2.669) (Sanduleak & Pesch 1989). The relatively high signal to noise reached in the spectra along with an extensive wavelength coverage of the Lyman alpha and the Lyman beta forest allowed us to obtain, through profile fitting, column densities and Doppler parameters of the Lyman clouds towards these QSOs. The spectral coverage of regions longward of the Lyman alpha emission line of the QSOs at the expected wavelengths of CIV 1548,1550 at the redshifts of the Lyman alpha forest allowed us to identify some heavy element absorption systems (hereafter HEASs) towards these QSOs. We have found no Lyman absorption system or HEAS towards these QSOs for which our data allow a deuterium measurement or to provide an interesting upper limit for the D/H ratio. The reason for this is that the Lyman lines with high column density detected towards these QSOs belong to absorption systems showing velocity structure.

We present the results of our study on three close pairs of QSOs. Our results are consistent with the existence of a proximity effect due to the foreground QSO, but due to its weakness we can only reject the absence of such effect at approx. 1 sigma level. By modelling this proximity effect in terms of a simple photoionisation model, we find the best value for the UV ionizing background to be approx. 10**(-20.5) erg cm**(-2) s**(-1) Hz**(-1) srad**(-1) at the Lyman limit, and an absolute lower limit (95\% confidence) of 10**(-21.8) in the same units. This lower limit rejects a number of models for the UV background where it is mostly contributed by QSOs and absorption by Lyman limit systems is taken into account.

We have re--examined the chemical composition of \lya clouds using the composite--cloud technique, in which each \lya line in a spectrum is shifted to its rest frame wavelength and all rest frame spectra are co--added to form an `averaged' \lya cloud spectrum. We illustrate how various estimates of the spectrum and redshift evolution of the background ionizing UV flux lead to very different predictions for the relative strengths of heavy element lines in the forest clouds. We also show how the potential abundance limits depend on the various observational quantities and how different procedures may be required for OVI and CIV. Preliminary results are presented from an analysis of two high redshift QSO echelle spectra.

We present initial results of an imaging and spectroscopic survey of faint galaxies in fields of Hubble Space Telescope spectroscopic target QSOs. The primary objectives of the survey are (a) to determine the incidence, extent, and covering factor of extended gaseous envelopes of luminous galaxies and (b) to determine the fraction of Lyman alpha absorption systems that arise in luminous galaxies. The goal of the survey is to identify in each field under consideration all objects with apparent r-band magnitudes satisfying r< 21.5 within angular distances to the QSOs satisfying theta < 1.3 arcmin. We find that at z \apl 1 most luminous galaxies are surrounded by extended gaseous envelopes of approx 160 h^{-1} kpc radius and roughly unit covering factor and that at z \apl 1 the fraction of Ly alpha absorption systems that arise in luminous galaxies is at least 0.35 \pm 0.10 and is likely 0.65 \pm 0.18.