Observations of the molecular gas in galaxies are vital to understanding the evolution and star-forming histories of galaxies. However, galaxies with molecular gas maps of their whole discs having sufficient resolution to distinguish galactic structures are severely lacking. Millimeter wavelength studies at a high angular resolution across multiple lines and transitions are particularly needed, severely limiting our ability to infer the universal properties of molecular gas in galaxies. Hence, we conducted a legacy project with the 45 m telescope of the Nobeyama Radio Observatory, called the CO Multi-line Imaging of Nearby Galaxies (COMING), which simultaneously observed 147 galaxies with high far-infrared flux in $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=1-0$ lines. The total molecular gas mass was derived using the standard CO-to-H$_2$ conversion factor and found to be positively correlated with the total stellar mass derived from the WISE $3.4 \mu$m band data. The fraction of the total molecular gas mass to the total stellar mass in galaxies does not depend on their Hubble types nor the existence of a galactic bar, although when galaxies in individual morphological types are investigated separately, the fraction seems to decrease with the total stellar mass in early-type galaxies and vice versa in late-type galaxies. No differences in the distribution of the total molecular gas mass, stellar mass, and the total molecular gas to stellar mass ratio was observed between barred and non-barred galaxies, which is likely the result of our sample selection criteria, in that we prioritized observing FIR bright (and thus molecular gas-rich) galaxies.

We investigate the quasi-equilibrium state of one-dimensional self-gravitating systems. If the null virial condition is satisfied at initial time, it is found that the number density around the center of the system at the quasi-equilibrium state has the universality similar to two- and three-dimensional self-gravitating systems reported in \cite{Tashiro16,Tashiro10}. The reason why the null virial condition is sufficient for the universality is unveiled by the envelope equation. We present a phenomenological model to describe the universal structure by using a special Langevin equation with a distinctive random noise to self-gravitating systems. Additionally, we unveil a mechanism which decides the radius of the system.

We study a universal structure of two and three dimensional self-gravitating systems in the quasi-equilibrium state. It is shown numerically that the two dimensional self-gravitating system in the quasi-equilibrium state has the same kind of density profile as the three dimensional one. We develop a phenomenological model to describe this universal structure by using a special Langevin equation with a distinctive random noise to self-gravitating systems. We find that the density profile derived theoretically is consistent well with results of observations and simulations.

We clarify how magnetic reconnection can be derived from magnetohydrodynamics (MHD) equations in a way that is easily understandable to university students. The essential mechanism governing the time evolution of the magnetic field is diffusion dynamics. The magnetic field is represented by two components. It is clarified that the diffusion of a component causes a generation of another component that is initially zero and, accordingly, that the magnetic force lines are reconnected. For this reconnection to occur correctly, the initial magnetic field must be directed oppositely in the two regions, e.g., $y>0$ and $y<0$; must be concave (convex) for $y>0$ ($y<0$); and must be saturated for $y$ far from the x axis, which would indicate the existence of the current sheet. It will be clear that our comprehension based on diffusion runs parallel to the common qualitative explanation about the magnetic reconnection.

We will construct a theory which can explain the dynamics toward the steady state self-gravitating systems (SGSs) where many particles interact via the gravitational force. Real examples of SGS in the universe are globular clusters and galaxies. The idea is to represent an interaction by which a particle of the system is affected from the others by a special random force. That is, we will use a special Langevin equation, just as the normal Langevin equation can unveil the dynamics toward the steady state described by the Maxwell-Boltzmann distribution. However, we cannot introduce the randomness into the system without any evidence. Then, we must confirm that each orbit is random indeed. Of course, it is impossible to understand orbits of stars in globular clusters from observations. Thus we use numerical simulations. From the numerical simulations of SGS, grounds that we use the random noise become clear. The special Langevin equation includes the additive and the multiplicative noise. By using the random process, we derive the non-Maxwellian distribution of SGS especially around the core. The number density can be obtained through the steady state solution of the Fokker-Planck equation corresponding to the random process. We exhibit that the number density becomes equal to the density profiles around the core by adjusting the friction coefficient and the intensity of the multiplicative noise. Moreover, we also show that our model can be applied in the system which has a heavier particle, corresponding to the black hole in a globular cluster.

We derive the non-Maxwellian distribution of self-gravitating $N$-body systems around the core by a model based on the random process with the additive and the multiplicative noise. The number density can be obtained through the steady state solution of the Fokker-Planck equation corresponding to the random process. We exhibit that the number density becomes equal to that of the King model around the core by adjusting the friction coefficient and the intensity of the multiplicative noise. We also show that our model can be applied in the system which has a heavier particle. Moreover, we confirm the validity of our model by comparing with our numerical simulation.

We propose the self-organized relaxation process which drives a collisionless self-gravitating system (SGS) to the equilibrium state satisfying local virial (LV) relation. During the violent relaxation process, particles can move widely within the time interval as short as a few free fall times, because of the effective potential oscillations. Since such particle movement causes further potential oscillations, it is expected that the system approaches the critical state where such particle activities, which we call gravitational fugacity, is independent of the local position as much as possible. Here we demonstrate that gravitational fugacity can be described as the functional of the LV ratio, which means that the LV ratio is a key ingredient estimating the particle activities against gravitational potential. We also demonstrate that LV relation is attained if the LV ratio exceeds the crticaial value $b=1$ everywhere in the bound region during the violent relaxation process. The local region which does not meet this criterion can be trapped into the pre-saturated state. However, small phase-space perturbation can bring the inactive part into the LV critical state.