We have performed optical imaging observations of the dark cloud L1251 at multiple wavelengths, B, V, R, and I, using the 105 cm Schmidt telescope at the Kiso Observatory, Japan. The cloud has a cometary shape with a dense "head" showing star formation activity and a relatively diffuse "tail" without any signs of star formation. We derived extinction maps of A_B and A_V with a star count method, and also revealed the color excess (E_{B-V}, E_{V-R}, and E_{V-I}) distributions. On the basis of the color excess measurements we derived the distribution of the ratio of total to selective extinction R_V over the cloud using an empirical relation between R_V and A_lambda/A_V reported by Cardelli et al. In the tail of the cloud, R_V has a uniform value of ~3.2, close to that often found in the diffuse interstellar medium (~3.1), while higher values of R_V=4-6 are found in the dense head. Since R_V is closely related to the size of dust grains, the high R_V-values are most likely to represent the growth of dust grains in the dense star-forming head of the cloud.

We investigate the thermal and dynamical evolution of primordial gas clouds in the universe after decoupling. Comparing the time-scale of dynamical evolution with that of fragmentation, we can estimate the typical fragmentation scale. We propose the following scenario of the formation process of first luminous objects consisting of large number stars. First, by pancake collapse of the overdensity regions in the expanding universe or collision between clouds in potential wells, quasi-plane shocks form. If the shock-heated temperature is higher than about $10^4$ K, the postshock gas cools down to several hundred K by \Hbun line cooling, and the shock-compressed layer fragments into filamentary clouds. The filamentary cloud collapses dynamically once more and fragments into cloud cores. Finally, a primordial star forms in a cloud core. We show that the minimum mass of the first star is essentially determined by the Chandrasekhar mass. Also, we investigate the dynamical collapse of cloud cores by numerical simulation and show that the evolution paths of the central regions of the cores depend only very weakly on the total core mass. After mass accretion, a massive star may be formed in a core, since the estimated mass accretion rate is very large. In such a case, it may be possible for many massive stars form almost simultaneously. Then the clouds can be luminous objects. On the other hand, if the shock-heated temperature is lower, effective star formation is delayed significantly.

The contraction processes of metal-free molecular clouds of starlike mass (or cloud cores) are investigated. We calculate radiative transfer of the H_2 lines and examine quasi-static contraction with radiative cooling. Comparing two time-scales, the free-fall time t_ff and the time-scale of quasi-static contraction t_qsc (nearly equal to t_cool, the cooling time) of these cores, we find that the ratio of the two time-scales t_ff/t_qsc, i.e., the efficiency of cooling, becomes larger with contraction even under the existence of cold and opaque envelopes. In particular, for fragments of primordial filamentary clouds, for which t_ff is nearly equal to t_qsc at the fragmentation epoch, they collapse dynamically in the free-fall time-scale. This efficiency of cooling is unique to line cooling.

We re-examine the thermal evolution of the postshock layer in primordial gas clouds. Comparing the time scales, we find that the evolutionary paths of postshock regions in primordial gas clouds can be basically understood in terms of the diagram drawn in the ionization degree vs temperature plane. The results obtained from the diagram are independent of the density in the case that we do not include photodissociation and photoionization. We also argue that the diagram is not only relevant to the case of the steady postshock flow, but also to the isochorically cooling gas.

The lower limit of the mass of the first stars suggested recently may imply the formation of massive stars of mass greater than 8 solar mass irrespective of the details of the initial mass function. The production of heavy metals from the first stars will ensure a requisite for the existence of life without the anthropic principle.

We discuss the fragmentation of primordial gas clouds in the universe after decoupling. Comparing the time scale of collapse with that of fragmentation, we obtain the typical mass of a fragment both numerically and analytically. It is shown that the estimated mass gives the minimum mass of a fragment which is formed from the primordial gas cloud and is essentially determined by the Chandrasekhar mass.

We investigated the thermal evolution of non-spherical primordial gas clouds of mass $ M=10^6M_\odot $. We studied two cases: 1) prolate and oblate clouds without angular momentum, 2) initially spherical, oblate and prolate clouds with angular momentum. In the spherical case, the magnitude of the angular momentum is the key quantity which determines the fragment mass. The fragment mass is found to be $M \gsim 100 M_\odot$ for $10^{-2} \lsim \lambda \lsim 1$, where $\lambda$ is the cosmological spin parameter. For the oblate shape of initial gas clouds, the angular momentum is almost always not important, In the case of prolate clouds with some angular momentum, collapse proceeds in the same way as spherical case, and the bounce occurs to form a disk due to the angular momentum. In any case, the estimated ``fragment mass'', must be above $100 M_\odot$ for the primordial gas clouds with typical angular momentum and typical oblateness or prolateness.