We present numerical simulations of the evolution of low-mass, isothermal, molecular cores which are subjected to an increase in external pressure $P\xt$. If $P\xt$ increases very slowly, the core approaches instability quite quasistatically. However, for larger (but still quite modest) $dP\xt/dt$ a compression wave is driven into the core, thereby triggering collapse from the outside in. If collapse of a core is induced by increasing $P\xt$, this has a number of interesting consequences. (i) The density profile is approximately flat in the centre during the prestellar phase (i.e. before the compression wave converges on the centre creating the central protostar). (ii) During the prestellar phase there are (subsonic) inward velocities in the outer layers of the core, whilst the inner parts are still approximately at rest. (iii) There is an initial short phase of rapid accretion (notionally the Class 0 phase), followed by a longer phase of slower accretion (the Class I phase). All these features accord well with observation, but are at variance with the predictions of the standard theory of star formation based on the inside-out collapse of a singular isothermal sphere. We note that the setting up of a coherent inward velocity field appears to be a generic feature of compression waves; and we speculate that interactions and interference between such velocity fields may play a crucial r\^ole in initiating the fragmentation of cores and the genesis of multiple star systems.

We analyse and compare the clustering of young stars in Chamaeleon I and Taurus. We compute the mean surface-density of companion stars \bar{N} as a function of angular displacement \theta from each star. We then fit \bar{N}(\theta) with two simultaneous power laws. For Chamaeleon I, the exponents of the power laws are 1.97 and 0.28, with the elbow at ~0.011 degrees. For Taurus, we obtain 2.02 and 0.87, with the elbow at ~0.013 deg. For both star clusters the observational data make quite large systematic excursions from the best fitting curve in the binary regime. These excursions may be attributable to evolutionary effects of the types discussed recently by Nakajima et al. and Bate et al. In the clustering regime the data conform to the best fitting curve very well. We also calculate the box-dimensions for the two star clusters. However, the limited dynamic range makes these estimates simply descriptors of the large-scale clustering, and not admissible evidence for fractality. We also propose two algorithms for objectively generating maps of constant stellar surface-density in young star clusters. Such maps are useful for comparison with molecular-line and dust-continuum maps of star-forming clouds, and with the results of numerical simulations. These maps retain information which is suppressed in the evaluation of \bar{N}(\theta). Analyzing such maps can be used to discriminate between fractal structure and single-level clustering, and to determine the degree of central condensation in small-N clusters. Algorithm I uses a universal smoothing length, and therefore has a restricted dynamic range, but it is implicitly normalized. Algorithm II uses a local smoothing length, which gives it much greater dynamic range, but it has to be normalized explicitly.