The results of a near-infrared (J H K LP) imaging linear polarimetry survey of 20 young stellar objects (YSOs) in rho Ophiuchi are presented. The majority of the sources are unresolved, with K-band polarizations, P_K < 6 per cent. Several objects are associated with extended reflection nebulae. These objects have centrosymmetric vector patterns with polarization discs over their cores; maximum polarizations of P_K > 20 per cent are seen over their envelopes. Correlations are observed between the degree of core polarization and the evolutionary status inferred from the spectral energy distribution. K-band core polarizations >6 per cent are only observed in Class I YSOs. A 3D Monte Carlo model with oblate grains aligned with a magnetic field is used to investigate the flux distributions and polarization structures of three of the rho Oph YSOs with extended nebulae. A rho proportional to r^(-1.5) power law for the density is applied throughout the envelopes. The large-scale centrosymmetric polarization structures are due to scattering. However, the polarization structure in the bright core of the nebula appears to require dichroic extinction by aligned non-spherical dust grains. The position angle indicates a toroidal magnetic field in the inner part of the envelope. Since the measured polarizations attributed to dichroic extinction are usually <10 per cent, the grains must either be nearly spherical or very weakly aligned. The higher polarizations observed in the outer parts of the reflection nebulae require that the dust grains responsible for scattering have maximum grain sizes <=1.05 microns.

We present high quality near infrared imaging polarimetry of HL Tau at 0.4 to 0.6 arcsec resolution, obtained with Subaru/CIAO and UKIRT/IRCAM. 3-D Monte Carlo modelling with aligned oblate grains is used to probe the structure of the circumstellar envelope and the magnetic field, as well as the dust properties. At J band the source shows a centrosymmetric pattern dominated by scattered light. In the H and K bands the central source becomes visible and its polarisation appears to be dominated by dichroic extinction, with a position angle inclined by ~40 degrees to the disc axis. The polarisation pattern of the environs on scales up to 200 AU is consistent with the same dichroic extinction signature superimposed on the centrosymmetric scattering pattern. These data can be modelled with a magnetic field which is twisted on scales from tens to hundreds of AU, or alternatively by a field which is globally misaligned with the disc axis. A unique solution to the field structure will require spatially resolved circular polarisation data. The best fit Monte Carlo model indicates a shallow near infrared extinction law. When combined with the observed high polarisation and non-negligible albedo these constraints can be fitted with a grain model involving dirty water ice mantles in which the largest particles have radii slightly in excess of 1 um. The best fit model has an envelope structure which is slightly flattened on scales up to several hundred AU. Both lobes of the bipolar outflow cavity contain a substantial optical depth of dust (not just within the cavity walls). Curved, approximately parabolic, cavity walls fit the data better than a conical cavity. The small inner accretion disc observed at millimetre wavelengths is not seen at this spatial resolution.