In this work we present new, high-to-intermediate spectral resolution evolutionary population synthesismodels, complementing and extending the widely used Maraston models. The new models are based on four popular libraries of empirical stellar spectra, which has necessitated some modifications to the original code, while keeping much of the original ingredients – such as stellar energetics, treatment of the thermally pulsating asymptotic giant branch, and mass loss recipe – intact. In addition, we have computed models at very high resolution (R = 20000) based on the theoretical MARCS library, that extends far into the infrared region of the electromagnetic spectrum. A library-dependent, but modelindependent comparison is made, where both photometric and spectroscopic similarities and discrepancies are highlighted. We find that stellar population models employing empirical stellar spectra exhibit considerably bluer (B-V) colours compared to models adopting theoretical spectra synthesised from the Kurucz model atmospheres (such as the BaSeL library), but that some differences arise between the empirical libraries due to, in particular, the adopted temperature scale. Furthermore, the results obtained with the theoretical MARCS library are fully consistent with the empirical libraries in this respect. The same effect can be found also in other EPS models that are based on empirical stars. We show that this discovery, whose origin can be traced mainly to cool stars (Teff < 5000 K), leads to improved photometric agreement with both galaxy and Milky Way globular cluster data. Spectral energy distributions of the latter are also used for testing the models concerning their ability to reproduce, through full SED-fitting, the cluster ages and metallicities as derived through independent fitting in colour-magnitude diagrams. In general, the agreement is very promising, although the higher resolution of the new models cannot alleviate the age-metallicity degeneracy in the optical in any significant way. A comparison with models of absorption line indices with variable abundance ratios is also made, both for the full SED-fitting procedure and when measuring indices directly on the SED. We obtain satisfying agreement in the first case, but in the second case only when a subsample of indices are used, the combination of which is little sensitive to abundance ratio effects. As a side track we exploit for the first time for population synthesis purposes the vast stellar database of the Sloan Digital Sky Survey, but find that the metallicity and age range of the stars is currently too narrow for computing sensible population models. Finally, we re-invigorate the Ca K line for the purpose of absorption feature diagnostics, presenting a new version of the index which will help in separating solutions for a key case of the age-metallicity degeneracy.

Cosmology aims to study the origin, composition and evolution of the entire Universe. The standard model for cosmology, called ΛCDM , represents a good fit to most of the observations we have, but it is a phenomenological model with no strong theoretical foundation, so one of the biggest challenges in cosmology (but important for the entire physics) will be to understand if this is the correct model (and so try to find a theoretical framework for it) or if a model with some sort of “new” physics will take place as the standard one. From the theoretical point of view there are several attempts to solve open problems in cosmology, such as the origin of the Universe and the nature of dark energy; their solution could shed some light on profound and interesting questions potentially revolutionising our understanding of nature. Important data revealing the nature of dark energy will be provided by forthcoming and planned galaxy surveys, that will reach a high precision in their measurements. Data available in the next years will allow us to constrain much better the cosmic expansion history, the geometry of the Universe and the growth of structures within it. For this reason, in this thesis we focused on observational tests of one of the key aspects of a cosmological model, the growth of structures; this allowed us to perform tests of cosmological models and General Relativity. We performed studies of the evolution of growth and clustering of cosmological structures and the evolution of the gravitational potential, comparing effects that depend on them against observations coming from various datasets. In particular, in Chapter 2 we test the growth of structures and their clustering using Redshift-Space Distortions (RSD), developing a new methodology to carefully analyse large scale spectroscopic galaxy surveys; we implement and test a practical application of the wide-angle formalism and then we investigate the significance of different systematics that affect measurements of large scale RSD. In Chapter 3 we use the Integrated Sachs-Wolfe (ISW) effect to test cosmological models to search for possible deviations from the ΛCDM model and then to test a model for the evolution of low frequency radio sources. In Chapter 4 we forecast cosmological measurements it will be possible to obtain using forthcoming radio surveys, using different probes such as the auto-correlation of radio sources, the ISW effect, the Cosmic Magnification and a joint analysis, in order to show how they can be used to test deviations from the standard cosmological constant and General Relativity models.

This thesis is composed of two parts; the first part deals with observations of the molecular gas towards an unlensed, obscured quasar AMS12, and the second part investigates radio undetected, optically selected quasi-stellar objects (QSOs) to determine the nature of the radio flux density distributions of these objects. AMS12 is an unlensed, obscured, z = 2.767 quasar which we observed with the Plateau de Bure Interferometer to detect carbon monoxide rotational transitions and atomic carbon fine structure lines in the molecular gas. We present new detections of the CO(5-4), CO(7-6), [CI]( ³P₁- ³P₀) and [CI](³P₂− ³P₁) molecular and atomic line transitions in this thesis. AMS12 is the first unlensed, high redshift source to have both atomic carbon ([CI]) transitions detected. The highly excited molecular gas probed by CO(3-2), (5-4) and (7-6), is modelled with large velocity gradient models. The gas kinetic temperature TG, density n(H₂), and the characteristic size r₀, are determined using the dust temperature from the far-infrared spectral energy distribution which had the following best-fitting parameters log₁₀[LFIR/L☉] = 13.5, dust temperature TD = 88 K and emissivity index β=0.6, as a prior for the gas temperature. The best fitting parameters are TG = 89.6 K, n(H₂) = 10 3.9 cm⁻³ and r₀ = 0.8 kpc. The ratio of the [CI] lines gives a [CI] excitation temperature of 43 ± 10 K, indicating the [CI] and the high-excitation CO are not in thermal equilibrium. The [CI] excitation temperature is below that of the dust temperature and the gas kinetic temperature of the high excitation CO, perhaps because [CI] lies at a larger radius where there may also be a large reservoir of CO at a cooler temperature, which may be detectable through the CO(1-0). Using the [CI]( ³P₁− ³P₀) line we can estimate the strength of the CO(1-0) line and hence the gas mass. This suggests that a significant fraction (~30%) of the molecular gas is missed from the high-excitation line analysis, giving a gas mass higher than that inferred from the assumption that the high-excitation gas is a good tracer of the low-excitation gas. The stellar mass was estimated from the mid-/near-infrared spectral energy distribution to be M* ~ 3 × 10¹¹M☉. The Eddington limited black hole mass is found from the bolometric luminosity to be M• ≳ 1.5×10⁹M☉. These give a black hole - bulge mass ratio of M•/M* ≳ 0.005. This is in agreement with studies on the evolution of the M•/M* relationship at high redshifts, which find a departure from the local value ~ 0.002. In the second half of the thesis we investigate the possible existence of a lower envelope in the radio luminosity versus optical luminosity plane. We select a population of QSOs from the Sloan Digital Sky Survey photometric quasar catalogue from Richards et al. The QSOs are within a narrow redshift band 0.3 < zphot < 0.5 and cross-matched with the 1.4 GHz National Radio Astronomy Observatory Very Large Array Sky Survey. The radio images extracted from the positions of the optical QSOs are retained if the flux integrated over the beam size of the radio survey is less than 3σIrms ≤ 1.35 mJy. The radio-undectected QSO population is split into eight samples depending on their optical magnitudes and stacked to determine the mean flux in each sample. The stacked mean flux is detected in all but the faintest optical magnitude sample. The radio versus optical luminosity relation from the stacked samples hint at a lower envelope in the radio luminosity which may be interpreted as there being a minimum radio jet power for a given accretion rate. Stacking assumes the underlying distribution of the property being measured is fairly represented by the stacked result. We investigate the underlying distribution of the radio flux density from the QSOs taking the noise of the sample into account. We find the distribution of the QSO flux density is modelled by a power-law with a negative index in all eight optical magnitude samples. This implies the mean stacked result is not a good representation of the distribution of the flux density of the QSOs and that there is no lower envelope. This highlights the danger of interpreting results from stacking without verifying the distribution is characterised by the mean stacked value. We appear to recover the quasar optical luminosity function when we model the distribution of radio loudness parameters suggesting that, since we are essentially holding the radio flux density fixed, the radio loudness is a function of the optical luminosity. This suggest that the radio loudness is not a fundamental property of the QSO but rather the ratio of two independent properties, the radio and optical luminosities. We convert the radio loudness parameter to jet efficiencies and find a minimum jet efficiency of ηmin = 4 × 10⁻⁴. We find there is no sign of a minimum jet efficiency as far as our data’s sensitivity limit allows, so we expect η<ηmin. Hence we provide an observational constraint for theoretical models of jet production in the minimum jet efficiency.

The recent observational data in cosmology seem to indicate that the universe is currently expanding in an accelerated way. This unexpected conclusion can be explained assuming the presence of a non-vanishing yet extremely fine tuned cosmological constant, or invoking the existence of an exotic source of energy, dark energy, which is not observed in laboratory experiments yet seems to dominate the energy budget of the Universe. On the other hand, it may be that these observations are just signalling the fact that Einstein's General Relativity is not the correct description of gravity when we consider distances of the order of the present horizon of the universe. In order to study if the latter explanation is correct, we have to formulate new theories of the gravitational interaction, and see if they admit cosmological solutions which fit the observational data in a satisfactory way. Quite generally, modifying General Relativity introduces new degrees of freedom, which are responsible for the different large distance behaviour. On one hand, often these new degrees of freedom have negative kinetic energy, which implies that the theory is plagued by ghost instabilities. On the other hand, for a modified gravity theory to be phenomenologically viable it is necessary that the extra degrees of freedom are efficiently screened on terrestrial and astrophysical scales. One of the known mechanisms which can screen the extra degrees of freedom is the Vainshtein mechanism, which involves derivative self-interaction terms for these degrees of freedom. In this thesis, we consider two different models, the Cascading DGP and the dRGT massive gravity, which are candidates for viable models to modify gravity at very large distances. Regarding the Cascading DGP model, we consider the minimal (6D) set-up and we perform a perturbative analysis at first order of the behaviour of the gravitational field and of the branes position around background solutions where pure tension is localized on the 4D brane. We consider a specific realization of this set-up where the 5D brane can be considered thin with respect to the 4D one. We show that the thin limit of the 4D brane inside the (already thin) 5D brane is well defined, at least for the configurations that we consider, and confirm that the gravitational field on the 4D brane is finite for a general choice of the energymomentum tensor. We also confirm that there exists a critical tension which separates background configurations which possess a ghost among the perturbation modes, and background configurations which are ghost-free. We find a value for the critical tension which is different from the value which has been obtained in the literature; we comment on the difference between these two results, and perform a numeric calculation in a particular case where the exact solution is known to support the validity of our analysis. Regarding the dRGT massive gravity, we consider the static and spherically symmetric solutions of these theories, and we investigate the effectiveness of the Vainshtein screening mechanism. We focus on the branch of solutions in which the Vainshtein mechanism can occur, and we truncate the analysis to scales below the gravitational Compton wavelength, and consider the weak field limit for the gravitational potentials, while keeping all non-linearities of the mode which is involved in the screening. We determine analytically the number and properties of local solutions which exist asymptotically on large scales, and of local (inner) solutions which exist on small scales. Moreover, we analyze in detail in which cases the solutions match in an intermediate region. We show that asymptotically flat solutions connect only to inner configurations displaying the Vainshtein mechanism, while non asymptotically flat solutions can connect both with inner solutions which display the Vainshtein mechanism, or with solutions which display a self-shielding behaviour of the gravitational field. We show furthermore that there are some regions in the parameter space of the theory where global solutions do not exist, and characterize precisely in which regions the Vainshtein mechanism takes place.

In this thesis, we present a comparison of the evolution of the massive galaxies in the 7.8Gyr since redshift z=1 to the evolution predicted from galaxy formation models. Observing the most massive galaxies in the Universe at high redshift is challenging due to their red colours, owing to both their intrinsically red Spectral Energy Distributions (SEDs) and their redshift. In Chapter 1, We produce a method using catalogue-level data to produce matched aperture photometry for the SDSS and UKIDSS surveys in order to extend the wavelength coverage of a sample of galaxies in order to improve the precision with which models can be fitted to photometric data for these high redshift galaxies. Our matched photometry has consistent colours with those of the full processing of SDSS+UKIDSS images performed by the GAMA survey, and produces magnitudes within ∼0.1 magnitudes of the GAMA photometry for all galaxies. This is reduced to within 0.04 magnitudes when all blended sources are excluded. We compute stellar masses by fitting a Maraston et al. (2009) LRG model to both our derived photometry and that of the GAMA processing, and find that our photometry’s best fit stellar masses are within ∼0.2 dex of that which comes from the GAMA photometry, demonstrating that the method is consistent with that of a full processing, and that it is possible to quickly compute matched photometry for large area surveys of complimentary wavelength coverage. This is of vital importance for upcoming surveys e.g. DES, VISTA, EUCLID etc. Fitting Stellar Population Models to galaxy photometry is a widely used technique in order to convert from observables (colours, magnitudes) to physical properties (mass, absolute magnitude, age). In spite of their widespread use, the optical and Near Infrared (NIR) properties of stellar population models are still subject to debate. Two of the most commonly used models are those of (Maraston, 2005) (M05) and (Bruzual & Charlot, 2003) (BC03), which can differ greatly in the NIR due to the M05 models’ inclusion of the TP-AGB phase, which was neglected for BC03 models. We explore the ability of these models to reproduce measured optical+NIR properties of galaxies in Chapter 3. We produce matched optical+NIR photometry for the subsample of the galaxies surveyed by Zibetti et al. (2013) (Z13) which lie within the UKIDSS imaging area in an attempt to reproduce the findings of Z13, who conclude that their optical and NIR spectroscopy is better fit by models from Bruzual & Charlot (2003) than similar models from Maraston et al (2005). We compare the observed optical+NIR Spectral Energy Distributions (SEDs) to those of BC03 and M05 models, as well as the approximate Z13 NIR fluxes. Z13 found that M05 models fitted to the optical data and extrapolated into the NIR displayed excess flux in the NIR relative to the data, and BC03 models are better at reproducing the data. However, we show that our data is consistent with both sets of models, and on average brighter in the NIR than that of Z13. We also compare the strength of spectral features in the optical to rest frame optical and optical-NIR colours, and show that our set of Composite Stellar Population (CSP) models agree well with data, with a preference for the M05 models, showing the validity of using these models on massive galaxies. A measurement of the Stellar Mass Function (SMF) of galaxies is a powerful tool in detecting evolution of the galaxy population. With a statistically complete sample of a galaxy population down to a given stellar mass, it is possible to calculate a statistically complete SMF down to this mass. Comparison of the shape of this SMF to that of a similar sample over a different redshift interval allows the evolution of galaxies over this redshift interval to be calculated, in order to determine whether these galaxies are forming stars, merging or simply passively evolving. For this purpose, in 4 compute matched SDSS+UKIDSS photometry for the AA omega KIDSS SDSS (AUS) survey. This is a 145.416 deg² area survey of Luminous Red Galaxies (LRGs) from redshift z∼0.5 to z∼1 located within Stripe 82. We fit this photometry to a Maraston et al. (2009) Luminous Red Galaxy (LRG) template to give stellar masses, and scale masses according to the magnitude difference between the matched photometry and the SDSS model photometry in order to produce “total” stellar masses. We produce a volume-weighted SMF for the survey, and find that our SMF is consistent with the Maraston et al. (2013) SMF from the BOSS survey, meaning that the most massive galaxies in the universe are evolving passively from z=1 to the present day, which is a challenge to hierarchical models of galaxy formation. Comparison of observed SMFs to those produced by galaxy formation models is a method of testing the ability of the models to reproduce the evolution displayed by the real galaxy population. This is therefore a test of the physics included within the models, with the level of agreement between the simulation and the real galaxy SMF being indicative of whether the modelling has incorporated all the processes in action in the real universe. In order to test the ability of the state of the art semi analytical models of Henriques et al. (2013) (H13 hereafter), we compare SMFs of the simulated galaxies to those of the AUS and BOSS surveys in Chapter 5. The H13 galaxies were tailored via the application of both the AUS and BOSS colour and magnitude cuts, and SMFs calculated within lightcones of the same area as the surveys in order to compare equal volumes. Our findings extend the conclusions of Maraston et al. (2013), namely that the most massive galaxies in the simulations are not sufficiently massive to agree with the observed galaxy population at this redshift. By extending this analysis to redshift z∼1, we can confirm that the discrepancy is larger at higher redshift, with the difference between the most massive galaxies in the simulations and those observed being log(ΔM/M⊙) ≃0.2 at z≃0.6–0.7, whereas going beyond this to the range z≃0.7–1 the difference becomes log(ΔM/M⊙) ≃0.25, as can be seen in Figure 5.6, which demonstrates that the simulations are failing to either form, or assemble, the mass quickly enough to reproduce the observations. Instead, the simulations continue to assemble mass through to low redshift at a higher rate than is seen in the galaxy SMF. These discrepancies may indicate that the physics of the simulations is not fully accounting for the real processes in the Universe,and that we do not yet have a model capable of reproducing the galaxy population in the real universe. Clearly semi analytical galaxy simulations need to be modified in order to reproduce the observations, before being further challenged by upcoming spectroscopic surveys of galaxies at redshifts as high as z=2 eg. eBoss, DESI.

The properties of galaxies, such as the galaxy red fraction and galaxy stellar mass function, have been shown to depend upon their environment in the local Universe. Large scale photometric surveys such as the DES and in the future Euclid, will be vital to gain insight into the evolution of galaxy properties and the role of environment through cosmic time. Large samples come at the cost of redshift precision and this affects the measurement of galaxy environment. In this thesis an analysis pipeline is constructed to derive galaxy parameters including absolute magnitudes, stellar masses and galaxy environments. The analysis pipeline consists of well established components, such as HYPERZ, that performs SED fitting and components that I have developed and tested, including codes to compute galaxy environment. Five methods to compute galaxy environment are implemented, including three fixed aperture methods, based on spheres, cylinders and cones, the Nth nearest neighbour method and the adaptive Gaussian method. The codes are optimized and parallelized and are executed on Portsmouth’s high performance computer cluster. The codes are thoroughly tested using mock data. Further testing is conducted employing GAMA data, with an external collaborator. The pipeline is applied to two datasets and the results lead to two scientific papers: Etherington & Thomas (2015) and Etherington et al. (in DES collaboration review). The first study is based on a low redshift sample drawn from the SDSS. Spectroscopic and photometric redshifts and also simulated photometric redshifts with a range of uncertainties are employed to study the impact of photometric redshift uncertainty on measures of environment as a function of the aperture parameters. The photometric environments are found to have a smaller dynamic range compared to the spectroscopic measurements because uncertain redshifts scatter galaxies from dense environments into less dense environments. With the optimal aperture parameter values, even for large redshift uncertainties, ∼ 0.1, there is a Spearman Rank Correlation Coefficient of ∼ 0.4 between the photometric measurements and the spectroscopic benchmark environments. This is sufficient to extract an environment signal from large scale photometric surveys. The second study in this thesis is based on the science verification data from the DES. This is the first set of observations from the survey. This study uses ∼3.2 million galaxies from the SPT-East (South Pole Telescope) field that covers approximately 100 sq. deg. of the sky. From the grizY photometry the analysis pipeline is used to derive galaxy stellar masses and absolute magnitudes. The errors on these properties are assessed using Monte-Carlo simulations sampled from the full photometric redshift probability distributions. Galaxy environments are computed using a fixed conical aperture method, for a range of scales. Galaxy environment probability distribution functions are constructed and the dependence of the environment errors on the aperture parameters is investigated. The environment components of the galaxy stellar mass function for the redshift range: 0.15 < z < 1.05 are calculated. For z < 0.75 it is found that the fraction of massive galaxies is larger in high density environment than low density environments. The low and high density components converge with increasing redshift to z ∼ 1.0 where the shapes of the mass function components are indistinguishable. This redshift is important because it marks the transition between an earlier epoch where the mass distribution of galaxies is independent of environment and a later epoch where the mass distribution does depend on galaxy environment. This study shows the build up of high density structures around massive galaxies, through cosmic time. The results in this thesis demonstrate that large scale photometric surveys can produce competitive galaxy evolution science, enabling further investigations of the role of galaxy environment. This is hugely encouraging for current and future experiments.

With the substantial improvement in instrumentation and our ability to now probe ever greater regions of space, the study of the Universe in its totality has moved towards the regime of a precision discipline. Several probes are now used by cosmologists to study the underlying cosmological model and understand its constituents. Modern cosmology has honed in on a concordance model that tells us that the Universe is predominantly composed of ‘dark’ components which still remain elusive to discovery. Weak lensing has emerged as one such powerful tool in probing the cosmological model. Its clean application of General Relativity, as well as its insensitivity in distinguishing between luminous and ‘dark matter' make it an attractive probe of the large scale structure in the Universe. To date, almost all weak lensing studies have been conducted using optical data. This is due primarily to the constraints required for a weak lensing study, i.e. high angular resolution and a high number density of distant sources, being most readily met at these wavelengths. The primary goal of this thesis is to address the feasibility of conducting such weak lensing experiments at the much longer, radio wavelengths. Many of the existing radio facilities are either undergoing (e.g. Extended Multi-Element Radio Linked Interferometer (eMERLIN1), Expanded Very Large Array (EVLA2)) or are scheduled to undergo major upgrades resulting in them being able to provide high resolution and high sensitivity data over a large field-of-view; aiding greatly in the detection of many more galaxies, a primary goal for weak lensing. Coupled with this, new, large interferometric arrays (e.g. Low Frequency Array (LOFAR3), and eventually the Square Kilometre Array (SKA4)) are in the process of being built and they too will provide the necessary quality of data for weak lensing experiments.

This dissertation examines two fundamentally different types of isolated systems in general relativity. In part 1, an exact solution of the Einstein-Maxwell equations representing the exterior field of two arbitrary charged essentially spherically symmetric (Reissner-Nordström) bodies in equilibrium is studied. Approximate solutions representing the gravitational- electrostatic balance of two arbitrary point sources in general relativity have led to contradictory arguments in the literature with respect to the condition of balance. Up to the present time, the only known exact solutions which can be interpreted as the nonlinear superposition of two Reissner-Nordström bodies without an intervening strut has been for critically charged masses, [special characters omitted]. In this dissertation . the invariant physical charge for each source is found by direct integration of Maxwell's equations. The physical mass for each source is invariantly defined in a manner similar to which the charge was found. It is shown that balance without tension or strut can occur for non-critically charged bodies. It is demonstrated that other authors have not identified the correct physical parameters for the masses and charges of the sources. Examination of the fundamental parameters of the space-time suggests a refinement of the nomenclature used to describe the physical properties is necessary. Such a refinement is introduced. Further properties of the solution, including the multipole structure and comparison with other parameterizations, are examined. Part 2 investigates the viability of constructing gravitational and electromagnetic geons: zero-rest-mass field concentrations, consisting of gravitational or electromagnetic waves, held together for long periods of time by their gravitational attraction. In contrast to an exact solution, the method studied involves solving the Einstein or Einstein-Maxwell equations for perturbations on a static background metric in a self-consistent manner. The Brill-Hartle gravitational geon construct as a spherical shell of small amplitude, high-frequency gravitational waves is reviewed and critically analyzed. The spherical shell in the proposed Brill-Hartle geon cannot be regarded as an adequate geon construct because it does not meet the regularity conditions required for a non-singular source. An attempt is made to build a non- singular solution to meet the requirements of a gravitational geon. Construction of a geon requires gravitational waves of high-frequency and the field equations are decomposed accordingly. A geon must also possess the property of quasi-stability on a time-scale longer than the period of the comprising waves. It is found that only unstable equilibrium solutions to the gravitational and electromagnetic geon problem exist. A perturbation analysis to test the requirement of quasi-stability resulted in a contradiction. Thus it could not be concluded that either electromagnetic or gravitational geons meet all the requirements for existence. The broader implications of the result are discussed with particular reference to the problem of with particular reference to the problem of gravitational energy. / Graduate

The work reported in this thesis is based on the vierbein field formulation of gravitational theory, used in conjunction with the method of the compensating field. It is shown that the most general linear equations of second order for a tensor field, which are invariant under orientations of the local inertial frame and under gauge transformations of the vierbein field components are identical with Einstein's field equation written down in the weak field approximation. An attempt is made to take into account any possibly existing gravitating effect of gravitation by applying the method of the compensating field to the weak field Lagrangian, resulting in a set of nonlinear field equations. The invariance properties of the modified field equations are examined, and some special solutions are exhibited. / Science, Faculty of / Physics and Astronomy, Department of / Graduate

Students' beliefs about free-fall motion were explored using structured interviews. The sample of 24 students was composed of 6 students (3 boys, 3 girls) selected from grades 6, 8, 10 and 11 respectively. Three sets of tasks involving one actual experiment and a number of simulated thought problems were used to investigate (a) the student's beliefs about the motion of a single object, (b) the relevance of the variables of height difference, initial velocity difference, frame of reference difference and weight difference for two objects and (c) the combined action of these factors when more than one was present. It was found that the interview methodology and tasks used were effective for collecting the data required in an exploratory study of this type. It was possible to categorize the mode of action of each variable in terms of: (a) not operating, (b) as operating in a short impulse only, (c) as operating but slowly dissipating, or (d) as operating with a continuous action. Examples of most response categories occurred at all levels but a number of possible developmental trends by grade were evident. Also there was a possible indication of the resistance of certain intuitive beliefs to standard kinematics instruction for the grade 11 physics group. The results of this study could be useful to the classroom teacher as well as to the designer of a science curriculum. It suggests that students are able to explore some problems of motion beginning at the grade 6 level; that they should be allowed to explore the relevance of related variables; and that they should be encouraged to express and explore their own beliefs which they bring to the classroom about motion. / Education, Faculty of / Curriculum and Pedagogy (EDCP), Department of / Graduate