Using N-body+hydro simulations we study relations between the local environments of galaxies on 0.5 Mpc scale and properties of the luminous components of galaxies. Our numerical simulations include effects of star formation and supernova feedback in different cosmological scenarios: the standard Cold Dark Matter model, the Broken Scale Invariance model (BSI), and a model with cosmological constant (LCDM). In this paper, we concentrate on the effects of environment on colors and morphologies of galaxies, on the star formation rate and on the relation between the total luminosity of a galaxy and its circular velocity. We demonstrate a statistically significant theoretical relationship between morphology and environment. In particular, there is a strong tendency for high-mass galaxies and for elliptical galaxies to form in denser environments, in agreement with observations. We find that in models with denser environments (CDM scenario) ~ 13 % of the galactic halos can be identified as field ellipticals, according to their colors. In simulations with less clustering (BSI and LCDM), the fraction of ellipticals is considerably lower (~ 2-3 %). The strong sensitivity of morphological type to environment is rather remarkable because our results are applicable to ``field'' galaxies and small groups. If all galaxies in our simulations are included, we find a statistically significant dependence of the galaxy luminosity - circular velocity relation on dark matter overdensity within spheres of radius 0.5 Mpc, for the CDM simulations. But if we remove ``elliptical'' galaxies from our analysis to mimic the Tully-Fisher relation for spirals, then no dependence is found in any model.

Using 3D hydrodynamical simulations of galaxy formation with supernova feedback and a multiphase medium, we derive theoretical relations analogous to the observed Tully-Fisher (TF) relations in various photometric bands. This paper examines the influence of self-regulation mechanisms including supernova feedback on galaxy luminosities and the TF relation in three cosmological scenarios (CDM, \Lambda CDM and BSI (broken scale invariance)). The galaxy catalogs derived from our hydrodynamical simulations lead to an acceptably small scatter in the theoretical TF relation amounting to \Delta M =0.2-0.4 in the I band, and increasing by 0.1 magnitude from the I-band to the B-band. Our results give strong evidence that the tightness of the TF relation cannot be attributed to supernova feedback alone. However, although eliminating supernova feedback hardly affects the scatter, it does influence the slope of the TF relation quite sensitively. With supernova feedback, L \propto V_c^{3-3.5} (depending on the strength of supernova feedback). Without it, L \propto V_c^{2} as predicted by the virial theorem with constant M/L. The luminosity functions in the B and K bands are quite sensitive to supernova feedback at the faint end studied here. We find that the faint end of the B-band luminosity function (-18 \leq M_B \leq -15) has a slope that is steeper than the Stromlo-APM estimate, but in rough agreement with the recent ESO Slice Project estimates.

The time evolution of thermal and thermo-reactive instabilities of primordial gas in the presence of ionizing UV radiation is studied. We obtain conditions (depending on density, temperature, and UV radiation intensity) favorable for the formation of a multi-component medium. Nonlinear effects, especially those attributable to opacity, can play an important role. In certain parameter regimes, dramatic, rapid evolution of ionization states away from ionization equilibrium may occur even if system control parameters evolve slowly and the system starts in or near ionization equilibrium. Mathematically, these rapid changes may best be understood as manifestations of bifurcations in the solution surface corresponding to ionization equilibrium. The astrophysical significance of these instabilities is that they constitute a nongravitational mechanism (thermo-reactive instability) for decoupling isolated spatial regions with moderate initial fluctuations from their surroundings, even under the influence of heating processes due to ionizing UV radiation. In particular, isolated, decoupled regions with relatively high neutral hydrogen density may result, and these could be associated with Lyman alpha clouds of high column density or with Lyman limit systems.

We numerically simulate some of the most critical physical processes in galaxy formation: The supernova feedback, in conjunction with gasdynamics and gravity, plays a crucial role in determining how galaxies arise within the context of a model for large-scale structure. Our treatment incorporates a multi-phase model of the interstellar medium and includes the effects of cooling, heating and metal enrichment by supernovae, and evaporation of cold clouds. The star formation happens inside the clouds of cold gas, which are produced via thermal instability. We simulate the galaxy formation in standard biased CDM model for a variety of parameters and for several resolutions in the range 2--20$h^{-1}$kpc. In our picture, supernova feedback regulates the evolution of the gas components and star formation. The efficiency of cloud evaporation by supernova strongly influences star formation rates. This feedback results in a steady rate of star formation in large galaxies (mass larger than $2-3x10^{11}\Msun$) at a level of $(1-10)\Msun\yr$ for $z<3$. Supernova feedback has an even stronger effect on the evolution of dwarf galaxies, most of which have a small fraction of stars and extremely low luminosities: $M_R>-15$. In the case of both large and small galaxies, the distribution of luminous matter (stars) is strongly BIASED with respect to the dark matter. We find an approximate biasing measure of the form $\rho_{lum}= (\rho_{dm}/133)^{1.7}$ for z=0 and overdensities exceeding 1000. Deviations from this relation (a factor 2-3) depend on the environment. For halo masses exceeding $2x10^{10}\Msun$, the dependence of the absolute magnitude on the total mass can be approximated as $M_V=-18.5-4\log(M_{tot}/10^{11}\Msun)$, with a scatter of less than 0.5mag.

We study the evolution with redshift of the number density of Ly$\alpha$ clouds in a CDM model using numerical simulations including photo-ionization and cooling of the baryonic component. The photo-ionizing flux is consistently taken as proportional to the rate at which material cools below 5000~K in the simulation. The number density of lines with $N$(HI)~$>$~10$^{14}$~cm$^{-2}$ given by the simulation can be approximated as d$n$/d$z$~$\sim$~16.8(1+$z$)$^{0.8}$+0.687(1+z)$^{2.5}$ which fits well the data at any redshift. The number of weaker lines is reproduced within a factor of two and does not decrease with time as fast as the number of strong lines. For $N$(HI)~$>$~10$^{13}$~cm$^{-2}$, we find d$n$/d$z$~$\propto$~(1+$z$)$^{0.7}$. \par At $z$~=~0, the two-point correlation function for Lyman$\alpha$ clouds with $N$(HI)~$>$~10$^{13}$ and $>$~10$^{14}$~cm$^{-2}$ has a slope $\sim$~-1.75 comparable to the slope of the two point correlation function for the galaxy distribution. The signal is higher for stronger lines for separations smaller than 1~Mpc. Part of the Lyman$\alpha$ forest is indeed associated with galaxies however the overall picture is that Ly$\alpha$ absorption lines originate in the warm gas that traces the potential wells of filamentary structures defined by the dark matter.

We have simulated the formation of large-scale structure arising from COBE-normalized spectra computed by convolving a primordial double-inflation perturbation spectrum with the CDM transfer function. Due to the broken scale invariance ('BSI') characterizing the primordial perturbation spectrum, this model has less small-scale power than the (COBE-normalized) standard CDM model. The particle-mesh code (with $512^3$ cells and $256^3$ particles) includes a model for thermodynamic evolution of baryons in addition to the usual gravitational dynamics of dark matter. It provides an estimate of the local gas temperature. In particular, our galaxy-finding procedure seeks peaks in the distribution of gas that has cooled. It exploits the fact that ``cold" particles trace visible matter better than average and thus provides a natural biasing mechanism. The basic picture of large-scale structure formation in the BSI model is the familiar hierarchical clustering scenario. We obtain particle in cell statistics, the galaxy correlation function, the cluster abundance and the cluster-cluster correlation function and statistics for large and small scale velocity fields. We also report here on a semi-quantitative study of the distribution of gas in different temperature ranges. Based on confrontation with observations and comparison with standard CDM, we conclude that the BSI scenario could represent a promising modification of the CDM picture capable of describing many details of large-scale structure formation.

We study the distribution of low-redshift Ly$\alpha$ clouds in a CDM model using numerical simulations including photoionization and cooling of the baryonic component. The ionizing background is found to be efficient enough to keep most of the gas warm ($T$~=~1--5~10$^4$~K) and to prevent most of it from collapsing. In this scenario, the numerous low redshift ($z$~$<$~0.5) neutral hydrogen Ly$\alpha$ lines recently revealed by HST observations of QSOs originate in the warm gas that traces the potential wells of filamentary structures defined by the dark matter. The predicted line number density is consistent with observations. About 25\% of the Ly$\alpha$ lines are found to arise in extended envelopes of luminous galaxies or groups of galaxies ($M$~$>$~3~10$^{11}$~M$_{\odot}$) of radius of the order of 1.5~Mpc. The exact number, however, is model dependent. The remainder arises in gas independent of any galaxy and at a distance of 2 to 7~Mpc from the first neighbour galaxy.