\section{Introduction}

Galaxies show distinctive distributions on the plot of stellar mass ($M_{\rm star}$) and star-formation rate (SFR):
some galaxies are on the sequence showing a positive correlation between the two values, and some locate far below the correlation in the SFR direction.
The former are called the main sequence of star-forming (SF) galaxies \citep[hereafter SF main sequence, e.g.,][]{Noeske:2007lr,Elbaz:2007ub,Daddi:2007un} and the latter are called quenched or passive galaxies.
Investigating how the galaxy properties change along and across the SF main sequence leads to essential understanding of galaxy evolution \citep[e.g.,][]{Wuyts:2011uq}.
Provided that stars are formed in molecular clouds in the local Universe, it is vital to study the molecular gas properties along and across the SF main sequence.

It is reported that there are positive correlations between specific SFR (sSFR $=$ SFR/$M_{\rm star}$, or an offset from the SF main sequence, $\Delta$(MS)) with molecular gas mass ($M_{\rm mol}$) to $M_{\rm star}$ ratio ($\mu_{\rm mol}=M_{\rm mol}/M_{\rm star}$) and star-formation efficiency (SFE $=$ SFR/$M_{\rm mol}$), i.e., galaxies with higher sSFR tend to have more abundant molecular gas and form stars more efficiently from the molecular gas \citep[e.g.,][]{Saintonge:2012nj,Genzel:2015gn,Scoville:2017jw,Tacconi:2018pv}.
These correlations still hold for spatially resolved data (\citealt{Lin:2017ee,Ellison:2020dm,Ellison:2020kf}; Kajikawa in prep.).
In addition, it is known that these correlations are irrespective of the environment \citep{Koyama:2017lf} or mass concentration of stellar components of galaxies \citep{Koyama:2019vy}, indicating their essentiality for galaxy evolution.
Furthermore, recent studies claimed that the scatter of the SF main sequence is primarily due to the SFE variety with a secondary role of $\mu_{\rm mol}$ using spatially-resolved molecular-gas data and optical integral-field-unit data \citep{Ellison:2020kf}.

Molecular rotational transitions have been used to explore the physical properties of molecular gas in galaxies.
Especially, the lowest transition lines ($J=1-0$) of carbon monoxide, $^{12}$CO and its isotopes are ideal tracers of cold ($\sim10$~K) molecular gas which is raw material of star formation, since the corresponding energy gaps ($h \nu_{\rm rest}$ where $h$ is the Planck constant and $\nu_{\rm rest}$ is rest frequency of the line) divided by the Boltzmann constant are $\sim5.53$~K for $^{12}$CO($J$=1-0) ($\nu_{\rm rest}=115.271$~GHz, hereafter $^{12}$CO) and $\sim5.29$~K for $^{13}$CO($J$=1-0) ($\nu_{\rm rest}=110.201$~GHz, hereafter $^{13}$CO), respectively.
$^{12}$CO emission is generally optically thick but strong even in extragalactic objects, thus it is widely used to measure the molecular gas mass of galaxies by assuming CO-to-H$_2$ conversion factors \citep[e.g.,][and references therein]{Rickard:1975if,Solomon:1975jr,Young:1982sd,Bolatto:2013vn}.
On the other hand, $^{13}$CO is used as an optically thin tracer of molecular gas and known to be well correlated with dust extinction in molecular clouds in the Milky Way \citep[e.g.,][]{Dickman:1978nn,Frerking:1982zq,Lada:1994ya}.

Detailed studies of the line ratio of $^{12}$CO/$^{13}$CO (hereafter $R_{1213}$) have revealed a wide variety of $R_{1213}$ within the Milky Way.
In a giant molecular cloud (GMC) where stars are formed, $R_{1213}$ is reported to be $\sim3$ at the center and $\sim5$ for the whole cloud \citep{Gordon:1976ly,Solomon:1979zd,Polk:1988ly}, which are much smaller than the abundance ratio of [$^{12}$CO]/[$^{13}$CO] inferred from the $^{12}$C/$^{13}$C ratio which has a radial gradient, where $\sim20$ at the galaxy center and $\sim60$ in the Solar neighbourhood \citep{Milam:2005yb,Halfen:2017sd}.
This suggests that bulk of $^{12}$CO emission is optically thick in GMCs.
On the other hand, larger $R_{1213}$ values are reported in the peripheries of GMCs \citep[$R_{1213}\sim20$,][]{Sakamoto:1994di}, high-latitude molecular clouds \citep[$\sim10.5$,][]{Blitz:1984ky}, and small molecular clouds with a size of $<1.0$~pc \citep[$\sim40$,][]{Knapp:1988gf}, suggesting lower optical depth of $^{12}$CO line in those components than in GMCs.
\cite{Polk:1988ly} obtained a relatively high $R_{1213}$ value of $\sim7$ when they calculate an average $R_{1213}$ value of seven $0.5^\circ \times 0.5^\circ$ regions in the Galactic plane and suggested a non-negligible contribution of the lower optical depth (diffuse) molecular gas to large-scale $^{12}$CO emission of the Milky Way.

Since the early studies of CO isotopes of local galaxies in 1970s (\citealt{Encrenaz:1979dw}, upper limits by \citealt{Rickard:1977ox}), various and fairly larger $R_{1213}$ values than those of GMCs have been reported in nearby galaxies, for example, a typical $R_{1213}$ of $\sim10$ in spiral galaxies \citep[e.g.,][]{Young:1982sd,Young:1986hh,Weliachew:1988tx,Sandqvist:1988fh,Sage:1991rj,Wright:1993ic,Braine:1993iw,Aalto:1994yg,Xie:1994do,Matsushita:1998si,Paglione:2001if,Hirota:2010af,Watanabe:2011na,Davis:2014tq,Vila-Vilaro:2015ip,Morokuma-Matsui:2015cr,Cormier:2018lt,Lee:2018rt}, $\sim10$ in early-type galaxies \citep[ETGs, e.g.,][]{Eckart:1990rp,Sage:1990yz,Wild:1997dm,Crocker:2012vs,Davis:2014tq,Alatalo:2015ca}, $\sim15$ in irregular galaxies \citep[e.g.,][]{Becker:1991da}, $\sim30$ in merging galaxies \citep[e.g.,][]{Aalto:1991bn,Casoli:1992zy,Aalto:1997wk,Taniguchi:1998yi,Aalto:2010yx}, starburst galaxies \citep[e.g.,][]{Stark:1984av,Young:1984sl,Kikumoto:1998xw}, and infrared-bright galaxies \citep[e.g.,][]{Garay:1993wo,Aalto:1995ir,Cao:2017bf,Sliwa:2017aq,Herrero-Illana:2019lc}, $\sim10$ in Seyfert galaxies \citep[e.g.,][]{Papadopoulos:1998fh,Papadopoulos:1999wa}, and $\sim30$ in $z\sim3$ galaxies with $J>1$ transitions \citep[e.g.,][]{Henkel:2010av,Danielson:2013yl,Spilker:2014lc,Bethermin:2018jq}.
The large $R_{1213}$ values would reflect a presence of substantial amount of low-optical depth (diffuse) components in those galaxies \citep[e.g.,][]{Polk:1988ly,Aalto:1995ir,Garcia-Burillo:1992as,Sakamoto:1997gr,Hirota:2010af,Morokuma-Matsui:2015cr}.
These studies are suggestive of a dependence of $R_{1213}$ on SF activities in galaxies but the relationship between $R_{1213}$ and sSFR or $\Delta$(MS) has not been systematically explored so far.

In this study, we compare the $R_{1213}$ and sSFR of galaxies targeted in the CO Multi-line Imaging of Nearby Galaxies project \citep[COMING,][]{Sorai:2019hs}, a legacy project making use of the 45-m radio telescope at the Nobeyama Radio Observatory\footnote{The Nobeyama 45-m radio telescope is operated by Nobeyama Radio Observatory, a branch of National Astronomical Observatory of Japan.} in order to investigate the molecular gas properties across the SF main sequence of galaxies.
The COMING project {\it simultaneously} obtained $^{12}$CO, $^{13}$CO, and C$^{18}$O($J=1-0$) mapping data of 147~galaxies with various sSFR located at the distances of $\lesssim40$~Mpc, which makes it one of the largest mapping  surveys of galaxies in multiple CO isotopes.
Provided that the $^{12}$CO and $^{13}$CO data in the COMING project are obtained in the same observing conditions including the accuracy of telescope pointings, $R_{1213}$ measurements do not suffer from the relevant uncertainties.

The structure of this paper is as follows:
we first briefly describe the data we used and data analysis, specifically stacking analysis, in section~\ref{sec:data_analysis},
we present the basic properties and the relationship between $R_{1213}$ and sSFR of our sample galaxies in section~\ref{sec:results},
based on the results, we discuss the origin of the $R_{1213}$ variation and its implication for the consequent variety in CO-to-H$_2$ conversion factor ($\alpha_{\rm CO}$) in local galaxies as well as high-redshift galaxies in section~\ref{sec:discussions},
and then we summarize this study in section~\ref{sec:summary}.
Throughout this study, we assume the Kroupa initial mass function (IMF) \citep{Kroupa:2001xf,Kroupa:2003do}.
For a correlation check, we refer to three values, Pearson's correlation coefficient $r$, Spearman's rank correlation coefficient $\rho$, and Kendall's rank correlation coefficient $\tau$.
The first one is a parametric measure and the latter two are non-parametric measures.
Although the obtained values are different from method to method, overall trends are almost consistent in the three methods.
We basically adopt the widely-used Spearman's $\rho$ to discuss the correlation of two galaxy parameters in this paper.

\section{Data \& Analysis}
\label{sec:data_analysis}

\begin{figure*}[h]
\begin{center}
\includegraphics[width=\textwidth, bb=0 0 1431 502]{basic_properties_20200520.pdf}
\end{center}
\caption{
Basic properties of the COMING galaxies:
(a) $M_{\rm star}$-SFR relation,
(b) $M_{\rm mol}$-SFR relation (``Kennicutt-Schmidt relation''), and
(c) sSFR-$\mu_{\rm mol}$ relation.
COMING galaxies are indicated as circle symbols, and among them, those with $R_{1213}$ measurements are indicated as filled black circles.
In the panel (a), the grey-filled contour indicates galaxies observed in the xCOLD GASS project and the thick contours with different shades of blue indicates xCOLD GASS galaxies with CO detection \citep{Saintonge:2011hl,Saintonge:2017ve}.
The SF main sequence at $z=0$ is indicated as cyan solid line and the 0.2-dex scatter range is indicated as cyan shade \citep{Speagle:2014by}.
In the panels (b) and (c), the blue-filled contour indicates xCOLD GASS galaxies with CO detection, which are same galaxies indicated as thick contours in the panel (a).
}
\label{fig:basic}
\end{figure*}

In this study, we used physical parameters of the COMING galaxies such as $M_{\rm star}$, SFR, sSFR, $M_{\rm mol}$, SFE, $\mu_{\rm mol}$, $R_{1213}$, and velocity dispersion of stacked {\bf $^{12}$CO} spectra.
Here, we briefly mention the derivation of these parameters from the observed data in this section.

\subsection{Molecular lines: COMING $^{12}$CO and $^{13}$CO data}
\label{sec:mmol}

We used $^{12}$CO and $^{13}$CO data obtained in the COMING project \citep{Sorai:2019hs} to calculate the integrated intensity ratio of:
\begin{equation}
R_{1213} = I_{\rm ^{12}CO}/I_{\rm ^{13}CO},
\end{equation}
where $I_{\rm ^{12}CO}$ and $I_{\rm ^{13}CO}$ are integrated intensities in units of K~km~s$^{-1}$ of $^{12}$CO and $^{13}$CO emission lines, respectively.
The angular and velocity resolutions of the COMING data are $\sim17$ arcsec ($\sim1.6$~kpc at 20~Mpc) and 10~km~s$^{-1}$, respectively.
We stacked all the $^{12}$CO and $^{13}$CO spectra of pixels where $^{12}$CO emission is detected with a signal-to-noise ratio of larger than three for each galaxy after aligning the line center expected from the first-moment map of $^{12}$CO line \citep[hereafter, ``VA'' stacking,][]{Schruba:2011zr,Caldu-Primo:2013lb,Morokuma-Matsui:2015cr}.
The number of stacked pixels differs from galaxy to galaxy and ranges from 20 to 2431 (the median, first and third quartiles are 274.5, 154.5, and 447.5, respectively)\footnote{Since the number of galaxies with $^{13}$CO detection is even (80), the median, first and third quartiles are calculated as means of the central two values of the whole sample, the first quarter sample, and the last quarter sample, respectively when the galaxies are sorted according to the number of the stacked pixels.}.
With the VA stacking analysis, $^{13}$CO emission line is detected in 80 galaxies with a signal-to-noise ratio of $>3$.
It should be noted that the number of galaxies with a ``secure'' detection of $^{13}$CO emission in the original data (more than three pixels with $\geq4$ sigma detection in the $^{13}$CO integrated intensity and $\geq3$ sigma detection in the $^{12}$CO integrated intensity) is only 52 \citep{Sorai:2019hs}.

Velocity dispersions of stacked $^{12}$CO spectra, $\sigma_{{\rm ^{12}CO}}$~are derived by the single Gaussian fitting.
Since this value encompasses both random motion among clouds and/or volume-filling gas and velocity gradient due to galactic rotation within the $\sim1$~kpc beam, we correct the effect of inclination ($i$) by dividing the value by sin($i$) as $\sigma_{{\rm ^{12}CO},i}$=$\sigma_{{\rm ^{12}CO}}$/sin($i$).
$^{12}$CO data is also used to calculate molecular gas mass of galaxies.
We adopt $M_{\rm mol}$ from \cite{Sorai:2019hs} which assumed the standard CO-to-H$_2$ conversion factor of the Milky Way of $2.0\times10^{20}$~cm$^{-2}$~(K km s$^{-1}$)$^{-1}$ \citep{Bolatto:2013vn}.

\subsection{Stellar mass and SFR: WISE and GALEX}
\label{sec:mstarsfr}

We adopt $M_{\rm star}$ values derived in \cite{Sorai:2019hs}, which used $3.4$~$\mu$m data obtained with the Wide-field Infrared Survey Explorer \citep[WISE,][]{Wright:2010oi} and an empirical relation between $3.4$~$\mu$m luminosity and $M_{\rm star}$ presented in \cite{Wen:2013wt} assuming the Kroupa IMF.
\cite{Sorai:2019hs} reported that their $M_{\rm star}$ values are $\sim20$~\% lower than those calculated in the Spitzer Survey of Stellar Structure in Galaxies project \citep[S$^4$G][]{Sheth:2010mz} due to their conservative subtractions of foreground stars.
SFR map of the COMING galaxies is generated using WISE $22$~$\mu$m and far-ultraviolet (FUV) imaging data obtained with the Galaxy Evolution Explorer \citep[GALEX,][]{Martin:2005wd} as
\begin{equation}
\Sigma_{\rm SFR} = 3.2\times10^{-3} \times I_{22\mu m} + 8.1 \times 10^{-2} \times I_{\rm FUV},
\end{equation}
where $\Sigma_{\rm SFR}$ is SFR surface densities in units of $M_\odot$~yr$^{-1}$~kpc$^{-2}$ and $ I_{22\mu m}$ and $I_{\rm FUV}$ are WISE/$22\mu$m and GALEX/FUV intensities, respectively, in units of MJy~sr$^{-1}$ \citep{Casasola:2017zv}.
The SFR calibration method adopted in \cite{Casasola:2017zv} assumes the default IMF in {\tt STATBURST99} \citep{Leitherer:1999ye}.
To convert from the SFR with the {\tt STATBURST99} IMF to SFR with the Salpeter IMF \citep{Salpeter:1955bz}, one should multiply the value by a factor of 1.59.
Furthermore, to convert from the SFR with Salpeter IMF to SFR with Kroupa IMF, one should multiply the value by a factor of 0.67.
Details in estimating SFR and star-formation properties will be presented in the forthcoming paper (Takeuchi et al. in prep.).
In this study, we limit our sample to the COMING galaxies with both the WISE $22$~$\mu$m and the GALEX FUV detections, resulting in 107 galaxies.
The number of galaxies with $M_{\rm star}$ and SFR measurements and $^{13}$CO detections reduces to $63$.
The $63$ galaxies consist of spiral galaxies \citep[S, SA, SAB, or SB classifications in the Third Reference Catalogue of Bright Galaxies, RC3,][]{de-Vaucouleurs:1991pr} and no interacting galaxies, irregular galaxies, nor ETGs are included.

\subsection{Dust emission: IRAS $60~\mu$m and $100~\mu$m data}
\label{sec:iras}
We crossmatch the 147 COMING galaxies with the IRAS Revised Bright Galaxy Sample \citep[RBGS,][]{Sanders:2003ci}, which lists total fluxes covering extended emission of galaxies in IRAS bands.
The IRAS RBGS consists of 629 galaxies with $60~\mu$m flux of $>5.24$~Jy.
The matching results in 111 galaxies, and 62 galaxies have all the measurements of $M_{\rm star}$, SFR, stacked $^{12}$CO and $^{13}$CO emissions, and IRAS fluxes at $60~\mu$m and $100~\mu$m.

\section{Results}
\label{sec:results}

\begin{figure*}[h]
\begin{center}
\includegraphics[width=\textwidth, bb=0 0 1397 500]{r1213_ssfr_20200627.pdf}
\end{center}
\caption{
The obtained relations between (a) $R_{1213}$ and sSFR and (b) $R_{1213}$ and $\Delta$(MS) of the 63 COMING galaxies.
Panel (c) shows the histogram of the $R_{1213}$ values for the 80 galaxies with $R_{1213}$ measurements and the 63 galaxies with both $R_{1213}$ and sSFR measurements.
The fitting functions with ODR are presented at the bottom of the panels (a) and (b).
Grey-dashed line and grey-shaded regions in panel (c) indicates an error-weighted mean and a standard deviation of the 63 COMING galaxies.
}
\label{fig:r1213_ssfr}
\end{figure*}

In figure~\ref{fig:basic}, we present basic properties of the COMING galaxies and compare them with galaxies in the extended CO Legacy Database for GASS project \citep[xCOLD GASS,][]{Saintonge:2011hl,Saintonge:2017ve}, which is the largest $^{12}$CO survey targeting 532 galaxies at $0.01<z<0.05$.
The SF main sequence at the local Universe is also indicated as a solid cyan line as a reference in the panel (a) \citep{Speagle:2014by}.
We can see that the most COMING sample is above or on the SF main sequence and some are located on the so-called ``green valley'' regime.
It should be noted that our sSFR estimations of the COMING galaxies are expected to be higher than those in literatures since the adopted $M_{\rm star}$ values are known to be $\sim20$~\% lower than the previous studies as mentioned in section~\ref{sec:mstarsfr}.
Overall, the COMING and the CO-detected xCOLD-GASS galaxies share similar star-formation properties such as sSFR, SFE and $\mu_{\rm mol}$.
In panel (c) of the same figure, it is also seen that galaxies without $^{13}$CO detection indicated as open black circles tend to have low $\mu_{\rm mol}$ and/or high sSFR values.

An error-weighted mean and a standard deviation of $R_{1213}$ for the 80 COMING galaxies with $^{12}$CO and $^{13}$CO measurements are $10.9$ and $7.0$, respectively.
In case we limit sample to the 63 galaxies that are plotted on the figure~\ref{fig:r1213_ssfr}a, these value becomes $11.4$ (mean) and $6.7$ (standard deviation).
This is consistent with previous studies on $R_{1213}$ of local star-forming galaxies \citep[$\sim10-15$, e.g.,][]{Young:1986hh,Paglione:2001if,Vila-Vilaro:2015ip}.
Using $R_{1213}$ and sSFR of the 63 galaxies, we find a moderate correlation between these two values with a Spearman rank correlation coefficient $\rho$ of $0.46$ (see figure~\ref{fig:r1213_ssfr}a).
The possibility of a fake correlation for the $R_{1213}$-sSFR relation is discussed in the Appendix section.
The fitting function by the orthogonal distance regression (ODR) is shown in this figure.

We also examine the relationship between $R_{1213}$ and the SFR offset from the SF main sequence, $\Delta$(MS) that is calculated as
\begin{equation}
\Delta({\rm MS}) = \log({\rm SFR}) - \log({\rm SFR}_{\rm MS}(M_{\rm star}, t)),
\end{equation}
where SFR$_{\rm MS}$($M_{\rm star}$, $t$) is SFR of the SF main-sequence galaxies with mass of $M_{\rm star}$ [$M_\odot$] at the age of the Universe of $t$ [Gyr], which is calculated as
\begin{equation}
\log({\rm SFR}_{\rm MS}(M_{\rm star}, t)=(0.84-0.026) \log(M_{\rm star}) - (6.51-0.11t)
\end{equation}
\citep{Speagle:2014by}\footnote{\cite{Speagle:2014by} compiled 64 observations of main-sequence galaxies at $0<z<6$ from 25 studies and derived the SF main sequence as a function of redshift after converting them to the same calibrations.}.
Here we adopt the same cosmology parameters as \cite{Speagle:2014by} of $(h, \Omega_{M}, \Omega_{\Lambda})=(0.7, 0.3, 0.7)$ to calculate the age of the Universe at $z=0$.
$R_{1213}$ seems to similarly but rather weakly correlate with $\Delta$(MS) compared to sSFR.
It should be noted that this tendency does not change even if another widely-used definition of the SF main sequence \citep{Whitaker:2012ok} is adopted, where the obtained Pearson's $r$, Spearman's $\rho$, and Kendall's $\tau$ are 0.41, 0.39, and 0.27, respectively.
Since $\Delta$(MS) depends on the definition of the main sequence of SF galaxies, we prefer to use sSFR, a more assumption free value, rather than $\Delta$(MS) to see possible variations in physical properties of molecular gas {\it across} the main sequence of SF galaxies in the following sections.

\section{Discussions}
\label{sec:discussions}

\begin{figure*}[h]
\begin{center}
\includegraphics[width=\textwidth, bb=0 0 2807 2821]{comparison_r1213_20200520.pdf}
\end{center}
\caption{
Correlation plots among six key parameters (from top to bottom and right to left: $R_{1213}$, $f_{60}/f_{100}$, $\sigma_{{\rm ^{12}CO},i}$, sSFR, SFR, and $M_{\rm star}$).
All the galaxies with measurements of both abscissa and ordinate values are indicated as open circle, and filled circle indicates galaxies with $R_{1213}$ and sSFR measurements among them.
The colors of the symbols indicate the $R_{1213}$ value.
The grey-shaded region indicates $\sigma_{{\rm ^{12}CO},i}$~$<10$~km~s$^{-1}$ that is the velocity resolution of the COMING data.
The number of galaxies with measurements of both abscissa and ordinate values is also indicated in each panel.
}
\label{fig:comp}
\end{figure*}

We found that $R_{1213}$ somehow relates to sSFR (and $\Delta$(MS)) suggesting that physical properties of molecular gas traced by $^{12}$CO can be different across the main-sequence of SF galaxies.
Here, we investigate relationship of $R_{1213}$ with the other physical parameters of galaxies in section~\ref{sec:corrplots} and discuss causes of the $R_{1213}$ variety in sections~\ref{sec:causes} and \ref{sec:morph}, and its implications for variety in CO-to-H$_2$ conversion factor ($\alpha_{\rm CO}$) in section~\ref{sec:implication}.
We also discuss possible effects of the $\alpha_{\rm CO}$ variation across the SF main sequence of galaxies on the known relations of sSFR with SFE and $\mu_{\rm mol}$ in section~\ref{sec:galevo}.

\subsection{Correlation plots for other parameters}
\label{sec:corrplots}

\begin{figure*}[h]
\begin{center}
\includegraphics[width=\textwidth, bb=0 0 1685 519]{corrcoeffs_r1213_20200520.pdf}
\end{center}
\caption{
Correlation coefficients (upper row) and $p$-values (lower row) of the plots in each panel in figure~\ref{fig:comp} for galaxies with the $R_{1213}$ and sSFR measurements.
The correlation coefficients are measured in three difference methods:
(a) Pearson's correlation coefficient, $r$,
(b) Spearman's rank correlation coefficient, $\rho$,
(c) Kendall's $\tau$ coefficient.
The pixel color indicates the correlation coefficient value.
Note that (a) is parametric and (b) and (c) are non-parametric methods.
(c) is normally more robust than (b) and generally (c)$<$(b) (see Appendix).
}
\label{fig:corr}
\end{figure*}

In order to investigate the $R_{1213}$-sSFR relation, we further compare $R_{1213}$ and sSFR with other four parameters including SFR, $M_{\rm star}$, $\sigma_{{\rm ^{12}CO},i}$, and $f_{60}/f_{100}$ in figure~\ref{fig:comp}.
We focus on these values because SFR and $M_{\rm star}$ are respectively a numerator and a denominator of sSFR, and $\sigma_{{\rm ^{12}CO},i}$ and $f_{60}/f_{100}$ are measures of physical properties of molecular gas such as turbulence and temperature.
Figure~\ref{fig:corr} shows correlation coefficients of each relation shown in figure~\ref{fig:comp}, which are calculated for galaxies with $R_{1213}$ and sSFR measurements on each panel, i.e., the 63 galaxies indicated as filled circles (or the 62 galaxies for the panels that are related to $f_{60}/f_{100}$).

We find moderate correlations (correlation coefficients of $0.3-0.5$) between $R_{1213}$ and two parameters: $\sigma_{{\rm ^{12}CO},i}$, and $f_{60}/f_{100}$ (figures~\ref{fig:comp} and \ref{fig:corr}).
Thus the $R_{1213}$ variation across the main sequence of SF galaxies is likely to be related to the varieties in physical properties of molecular gas.
Note that the correlation coefficients do not significantly change even if we simply use $\sigma_{{\rm ^{12}CO}}$~instead of $\sigma_{{\rm ^{12}CO},i}$.
On the other hand, no correlation is found for the $R_{1213}$-$M_{\rm star}$ and the $R_{1213}$-SFR relations.
This suggests the correlation between $R_{1213}$ and sSFR is not a consequence of neither a positive correlation between $R_{1213}$ and SFR nor a negative correlation between $R_{1213}$ and $M_{\rm star}$.

The strongest correlation is found in the $f_{60}/f_{100}$-sSFR relation.
Considering that $f_{60}/f_{100}$ is a measure of dust temperature, this is consistent with previous studies showing higher dust temperature for galaxies with higher sSFR \citep[e.g.,][]{Magnelli:2014zr,Matsuki:2017nl}.
The stronger correlation of the $f_{60}/f_{100}$-sSFR relation than the one of the $R_{1213}$-sSFR relation suggests that the former is more fundamental than the latter.
Thus, we consider that the $R_{1213}$ correlates with sSFR primarily through the correlations of the sSFR-$f_{60}/f_{100}$ and $R_{1213}$-$f_{60}/f_{100}$ relations.

$\sigma_{{\rm ^{12}CO},i}$~may also play a role in generating the $R_{1213}$-sSFR relation.
As aforementioned, $R_{1213}$ positively correlates with $\sigma_{{\rm ^{12}CO},i}$.
In panels (o) or (v) of figure~\ref{fig:comp} (a comparison of $\sigma_{{\rm ^{12}CO},i}$~and sSFR), we can see that the bulk of galaxies including those with $R_{1213}$ measurements are located on the upper left or lower right regions showing a moderate positive correlation with a Spearman's $\rho$ of $0.46$.
There also exist some outliers with $\sigma_{{\rm ^{12}CO},i}$~$\lesssim10$~km~s$^{-1}$ showing a wide range of sSFR.
Consequently, the overall correlation becomes quite weak when including those outliers.
Note that these trends are still seen even when we used the stacked spectra of only disk region of galaxies.
It should be also noted that the $\sigma_{{\rm ^{12}CO},i}$~values of the galaxies with $\sigma_{{\rm ^{12}CO},i}$~$\lesssim10$~km~s$^{-1}$ are upper limits, considering the velocity resolution of $10$~km~s$^{-1}$ of the COMING data.

In summary, figures~\ref{fig:comp} and \ref{fig:corr} show that
(1) $R_{1213}$ correlates with $\sigma_{{\rm ^{12}CO},i}$~and $f_{60}/f_{100}$,
(2) the correlation coefficients for the $\sigma_{{\rm ^{12}CO},i}$-sSFR (0.46) and $f_{60}/f_{100}$-sSFR (0.69) relations are comparable to or higher than the one for the $R_{1213}$-sSFR relation (0.46), and
(3) the correlation coefficient for the $f_{60}/f_{100}$-sSFR relation is higher than the one for the $\sigma_{{\rm ^{12}CO},i}$-sSFR relation.
Therefore, the observed trend of $R_{1213}$ along sSFR is considered to be mainly caused by $f_{60}/f_{100}$ with a secondary role of $\sigma_{{\rm ^{12}CO},i}$, while accurate measurements of $\sigma_{{\rm ^{12}CO},i}$ and $R_{1213}$ of those with $\sigma_{{\rm ^{12}CO},i}$~$\lesssim10$~km~s$^{-1}$ are required for a definite conclusion.

\subsection{Causes of the $R_{1213}$ variety}
\label{sec:causes}

As a cause of the $R_{1213}$ variation, four candidates have been mainly considered since the early studies of $R_{1213}$ in nearby galaxies: changes in
(1) physical conditions of molecular gas (density, temperature, and opacity),
(2) [$^{12}$CO]/[$^{13}$CO] abundance due to selective photodissociation,
(3) [$^{12}$CO]/[$^{13}$CO] abundance due to chemical fractionation, and
(4) $^{12}$C/$^{13}$C isotope abundance ratio due to stellar nucleosynthesis.
We consider that the scenario (1) is plausible for the observed $R_{1213}$ variation at least for our sample galaxies and the scenario (4) may play a secondary role.
In this section, we explain each scenario and reasons why the scenarios (2)-(4) are not expected to play a major role.

$R_{1213}$ can be approximated as a function of opacity of $^{12}$CO ($\tau_{^{12}{\rm CO}}$) and the abundance ratio of [$^{12}$CO]/[$^{13}$CO] under the local thermodynamic equilibrium (LTE) assumption as\footnote{
Some studies approximates as $R_{1213}\propto \frac{1-\exp{(-\tau_{^{12}{\rm CO}})}}{1-\exp{(-\tau_{^{13}{\rm CO}})}} \sim \frac{1}{\tau_{^{13}{\rm CO}}}$ by assuming that $\tau_{\rm ^{12}CO}$ is large enough \citep{Paglione:2001if,Hirota:2010af,Cao:2017bf}.}:
\begin{equation}
R_{1213}\propto
\frac{1-\exp{(-\tau_{^{12}{\rm CO}})}}{1-\exp{(-\tau_{^{13}{\rm CO}})}}
\sim \frac{1-\exp{(-\tau_{^{12}{\rm CO}})}}{1-\exp{\left (-\frac{\tau_{^{12}{\rm CO}}}{\rm [^{12}CO]/[^{13}CO]} \right)}},
\end{equation}
where $\tau_{^{13}{\rm CO}}$ is the optical depth of $^{13}$CO.
If $h\nu \ll k_{\rm B} T_{\rm k}$ (where $k_{\rm B}$, $\nu$, and $T_{\rm k}$ are the Boltzman constant, line frequency and kinetic temperature), $\tau_{^{12}{\rm CO}}$ scales as $\propto$~$\frac{N_{\rm ^{12}CO}}{\Delta v T^2_{\rm k}}$, where $N_{\rm ^{12}CO}$ and $\Delta v$ are $^{12}$CO column density and velocity dispersion, respectively.
Note that non-LTE analysis also predicts that $R_{1213}$ is a strong function of $\tau_{^{12}{\rm CO}}$ \citep{Cormier:2018lt}.
Thus, the $R_{1213}$ variation can be considered to reflect the $\tau_{^{12}{\rm CO}}$ variation as a first approximation, and $\tau_{^{12}{\rm CO}}$ is affected by temperature and turbulence of molecular gas.

Early studies have already noticed a correlation between $R_{1213}$ and $f_{60}/f_{100}$ (e.g., \citealt{Young:1986hh,Aalto:1995ir,Crocker:2012vs,Herrero-Illana:2019lc}, but see also \citealt{Vila-Vilaro:2015ip,Cao:2017bf}).
They discussed that high gas temperatures and turbulence (as traced by dust temperature) work against large optical depths in the $^{12}$CO.
However, no clear observational evidence has been found for a positive correlation between $R_{1213}$ with the velocity dispersion of $^{12}$CO line.
Although there are relevant studies to the $R_{1213}$-$\sigma_{{\rm ^{12}CO},i}$~relation, such as the relationships of $R_{1213}$ with velocity dispersion of $^{13}$CO line \citep{Meier:2004cy} and with inclinations of galaxies \citep{Sage:1991rj}, both studies reported null results.
In this study, we re-examined and confirmed the correlation between $R_{1213}$ and $f_{60}/f_{100}$ with the COMING data and additionally found the correlation between $R_{1213}$ and $\sigma_{{\rm ^{12}CO},i}$ when we limit the sample that are used to investigate the $R_{1213}$-sSFR correlation.
These correlations also support the scenario (1), when we assume $\Delta v \sim $$\sigma_{{\rm ^{12}CO},i}$, .

The abundance ratio of [$^{12}$CO]/[$^{13}$CO] can also affect $R_{1213}$.
The scenarios (2)-(4) are processes changing [$^{12}$CO]/[$^{13}$CO] ratio.
The selective photodissociation is expected to be dominant basically in low density ($n <10^2$~cm$^{-3}$) domains or in dense regions with very strong radiation fields \citep{Rollig:2013xw}.
UV radiation dissociates CO molecules and $^{13}$CO is less effective at self-shielding against UV photons than $^{12}$CO due to its lower abundance.
The chemical fractionation is expected to become important at deeper region of molecular clouds where the dust-shielding is effective for both $^{12}$CO and $^{13}$CO.
In this region with lower temperature of $<35$~K, $^{13}$CO is preferentially formed via the isotopic charge exchange reaction of:
\begin{equation}
^{12}{\rm CO} + ^{13}{\rm C}^+ \rightarrow ^{12}{\rm C}^+ + ^{13}{\rm CO} + \Delta {\rm E}.
\end{equation}

Assuming [$^{12}$CO]/[$^{13}$CO]$\sim^{12}$C/$^{13}$C, the $^{12}$C/$^{13}$C ratio is also potentially an important factor controlling $R_{1213}$.
The $^{12}$C/$^{13}$C ratio is a measure of the ratio of primary to secondary nucleosynthetic components,
thus the ratio is expected to decrease as stellar populations evolve, and to depend on the star-formation history \citep[e.g.,][]{Audouze:1975jo}, the slope of IMF \citep[e.g.,][]{Romano:2017gk}, and the treatments of stellar rotation which affects mixing in stars \citep[e.g.,][]{Limongi:2018fr,Romano:2019ke}.
It is expected that, in the early evolutionary phase of stellar population, $^{12}$C/$^{13}$C is high due to efficient supply of $^{12}$C from core-collapse supernovae (SNe) of massive stars, and decrease due to $^{13}$C production from intermediate-mass asymptotic giant branch (AGB) stars, then increase due to $^{12}$C production from low-mass AGB stars, and decrease again due to type-Ia SNe \cite[e.g.,][]{Kobayashi:2011gg}.

There are theoretical and observational studies of $R_{1213}$ in normal galaxies, which contradict with scenarios (2)-(4) as a dominant factor changing $R_{1213}$.
\cite{Szucs:2014qe} self-consistently calculated $R_{1213}$ of isolated molecular clouds with hydrodynamical simulations taking into account chemical network, various interstellar radiation fields, and cooling and heating processes.
They concluded that the selective photodissociation has minimal effect on $R_{1213}$ and the chemical fractionation can cause a factor of $2-3$ decreases.
Other theoretical studies also predicted that the photodissociation of CO molecule is not selective since self-shielding is found to be not important compared to dust- and H$_2$-shielding \citep{Krumholz:2014oe,Safranek-Shrader:2017hm}.

The observed anti-correlations between $R_{1213}$ and $^{13}$CO/C$^{18}$O reported in various galaxies also disfavor the selective photodissociation scenario \citep[e.g.,][]{Paglione:2001if,Tan:2011dy}, because the rarer C$^{18}$O is expected to be more reduced than $^{13}$CO if the selective photodissociation is working.
It is reported that the radial gradient of $R_{1213}$ in the Milky Way is in the opposite sense of the $^{12}$C/$^{13}$C abundance gradient \citep{Paglione:2001if}.
\cite{Crocker:2012vs} compared several line ratios using six molecular lines and various galaxy properties of $^{12}$CO-rich ETGs, and found that correlations between $R_{1213}$ and the galaxy properties were similarly seen even when they used $^{12}$CO/HCN ratio (no difference in carbon isotope) instead.
This suggests that the chemical fractionation and $^{12}$C/$^{13}$C variation due to stellar nucleosynthesis are not the dominant processes generating the $R_{1213}$ variation in these ETGs.

However, the scenario (4) is likely to become important for extreme dusty SF galaxies such as local ULIRGs, starburst galaxies \citep{Sliwa:2017aq,Brown:2019ge}, and high-$z$ sub-milimetre galaxies (SMGs) \citep{Danielson:2013yl,Zhang:2018mc}.
In these galaxies, high $R_{1213}$ and extremely low $^{13}$CO/C$^{18}$O ($\sim1$) values have been found.
Given that $^{18}$O is produced in massive stars ($\gtrsim8$~M$_\odot$) and assuming $^{13}$CO/C$^{18}$O~$\approx$~$^{13}$C/$^{18}$O, the $^{13}$CO/C$^{18}$O is expected to be low for young systems or for top-heavy IMF.
Indeed, chemical evolution models of galaxies show that such extreme line ratio can be explained by top-heavy IMF \citep{Romano:2017gk,Zhang:2018mc,Romano:2019ke}.
The observed correlation between $R_{1213}$ and sSFR in this study may suggest that the stellar nucleosynthesis is also important since sSFR can be a measure of ratio of massive/young and less massive/old stellar populations.
The dominant process determining $R_{1213}$ may be different across the SF main sequence.
In order to investigate the effect of $^{12}$C/$^{13}$C on the $R_{1213}$ trend across the main sequence, further deep observations in optically-thin lines are required.

\subsection{Dependence of the six parameters on galaxy morphologies}
\label{sec:morph}

\begin{table*}[t]
    \caption{Median and 1st and 3rd quartiles of six key parameters in figure~\ref{fig:morph} of galaxy subsamples with different morphology\label{tab:morph}}
        \begin{center}
            \begin{tabular}{ccccccc}\hline
            RC3 morphology\footnotemark[$*$] & SA & SAB+SB & SAB & SB & w/ pec & w/o pec\\
            \# of galaxies & 27 & 33 & 20 & 13 & 10 & 53\\
            \hline
            $M_{\rm star}$ [$10^{10}$~M$_\odot$] & 3.8\footnotemark[$**$] & 3.8 & 3.9 & 3.1 & 2.8 & 4.1\\
             &  (1.8 -- 7.1)\footnotemark[$***$] &  (2.4 -- 7.3) &  (2.3 -- 6.3) &  (2.6 -- 11) &  (1.7 -- 6.9) &  (2.4 -- 7.3)\\
            SFR [M$_\odot$ yr$^{-1}$] & 1.9 & 3.1 & 2.7 & 4.7 & 2.0 & 2.4\\
             &  (1.2 -- 2.7) &  (1.8 -- 5.2) &  (1.8 -- 3.8) &  (2.0 -- 8.3) &  (1.3 -- 3.4) &  (1.4 -- 4.9)\\
            sSFR [$10^{-11}$~yr$^{-1}$] & 6.2 & 7.6 & 7.3 & 8.0 & 6.8 & 7.3\\
             &  (3.8 -- 8.3) &  (6.2 -- 9.9) &  (6.2 -- 8.9) &  (6.7 -- 16) &  (4.3 -- 12) &  (5.2 -- 8.8)\\
            $\sigma_{{\rm ^{12}CO},i}$~[km/s] & 31 & 40 & 37 & 44 & 52 & 38\\
             &  (18 -- 47) &  (27 -- 53) &  (26 -- 45) &  (37 -- 58) &  (13 -- 59) &  (22 -- 46)\\
            $f_{60}/f_{100}$ & 0.41 & 0.47 & 0.43 & 0.48 & 0.49 & 0.43\\
             &  (0.37 -- 1.00) &  (0.40 -- 0.56) &  (0.39 -- 0.49) &  (0.46 -- 0.57) &  (0.41 -- 0.59) &  (0.39 -- 0.56)\\
            $R_{1213}$ & 12.2 & 15.3 & 16.9 & 14.1 & 17.1 & 14.1\\
             &  (10.5 -- 14.9) &  (13.2 -- 18.6) &  (14.2 -- 21.4) &  (12.9 -- 16.6) &  (13.2 -- 22.8) &  (10.8 -- 17.2)\\
            \hline
        \end{tabular}
    \end{center}
    \begin{tabnote}
    \footnotemark[$*$] There are three ``S'' galaxies among the 63 galaxies, i.e., no classification for bar structures.;
    \footnotemark[$**$] Median;
    \footnotemark[$***$] (1st quartile -- 3rd quartile)
    \end{tabnote}
\end{table*}

The turbulent and warm properties of interstellar medium (ISM) are a natural consequence for galaxies with high sSFR if star-formation feedbacks plays a major role in determining the ISM properties.
Theoretical studies showed that the molecular clouds are expanded due to feedback by \mbox{H\,{\sc ii}}~regions and are destroyed or displaced by SNe explosion \citep[e.g.,][]{Baba:2017fk,Grisdale:2018qm}.
It is also shown that the observed scaling relations of molecular clouds are reproduced when these feedbacks are taken into account \citep[e.g.,][]{Grisdale:2018qm,Fujimoto:2019dt}.
However, both the spatial and velocity resolutions of our data are not high enough to resolve a single molecular cloud and to discriminate the effects of star-formation feedback from the other larger-scale processes.

The shorter dynamical (or free-fall) timescale than the star-formation timescale \citep[$1-3$~\%,][]{Kennicutt:1998fk,Leroy:2008fb} suggests that secular processes are important for star formation in galaxies.
It is widely known that the bar-structures are expected to drive gas inflow by removing angular momentum of the gas and consequently induce active star formation at the central regions of galaxies \citep[e.g.,][]{Wada:1995ly,Sakamoto:1999rt,Kuno:2007uq}.
Galaxy mergers are also a key process inducing gas condensation and triggering star formation in galaxies \citep[e.g.,][]{Sanders:1996ae,Hopkins:2005tv}.
While enhancement of star formation due to mergers or stellar bars may not be so significant \citep[$10-15$~\%,][]{Rodighiero:2011rv,Saintonge:2012nj}, both processes help to enhance turbulence and broaden the linewidth of CO emission in galaxies \citep{Sorai:2012ix,Maeda:2018gw,Sun:2018wz,Yajima:2019xz}.

Figure~\ref{fig:morph} shows histograms of the six key values that are compared in figure~\ref{fig:comp} for subsamples with different types of bar structures \citep[SA/SAB/SB,][]{de-Vaucouleurs:1991pr} and table~\ref{tab:morph} summarizes medians and quartiles of these values for each subsample.
The galaxies categorized as ``pec'' in \cite{de-Vaucouleurs:1991pr} are also separately presented in the same figure and table.
According to the Kolmogorov-Smirnov (KS) test, there are no significant differences in $M_{\rm star}$ (p-value of $0.9$), sSFR ($0.07$), $\sigma_{{\rm ^{12}CO},i}$~($0.3$), and $f_{60}/f_{100}$ ($0.3$) between barred (SAB+SB) and non-barred (SA) galaxies, whereas SFR and $R_{1213}$ tend to be higher for barred galaxies than non-barred galaxies, whose $p$-values are $4\times10^{-3}$ and $3\times10^{-2}$, respectively.
Although the $p$-values of the KS test are not small enough, the median values of sSFR, $\sigma_{{\rm ^{12}CO},i}$, and $f_{60}/f_{100}$ are higher for barred galaxies than non-barred galaxies.
Therefore, the obtained tendency where galaxies with higher sSFR have broader $\sigma_{{\rm ^{12}CO},i}$~in this study might be {\it partially} attributed to those barred galaxies.

In terms of the morphological peculiarities indicated with the suffix of ``pec'', on the other hand, we do not find a clear differences in any of the six values (figure~\ref{fig:morph} and table~\ref{tab:morph}).
Although no closely-interacting galaxies is included in our sample, the morphological peculiarities are generally indicative of galaxy interaction.
This suggests that the galaxy interaction is unlikely to be a major cause of the wider $\sigma_{{\rm ^{12}CO},i}$~for galaxies with higher sSFR at least in our sample.
Therefore, high turbulence and temperature are likely to be primarily due to active star formation, and stellar bar might partially play some roles in driving the high turbulence of the molecular gas and the active star formation via efficient gas inflow.

\begin{figure*}[t]
\begin{center}
\includegraphics[width=150mm, bb=0 0 1395 723]{hist_r1213_morph_20200520.pdf}
\end{center}
\caption{
Histograms of the six key parameters for SA (white), SAB (light grey), SB (dark grey) galaxies and those with the ``pec'' suffix (black).
}
\label{fig:morph}
\end{figure*}

\subsection{Implications of the variety of CO-to-H$_2$ conversion factor across the main sequence}
\label{sec:implication}

We adopted a constant $\alpha_{\rm CO}$ to estimate molecular gas mass of the COMING galaxies in this study, while the $R_{1213}$ variety suggests different $\alpha_{\rm CO}$ among the COMING galaxies.
\cite{Cormier:2018lt} showed that the relationship between $R_{1213}$ and $\alpha_{\rm CO}$ is different from galaxies to galaxies, while their non-LTE radiative transfer calculations with {\tt RADEX} \citep{van-der-Tak:2007qb} predicts a clear anti-correlation between these two values, which is ascribed to $\tau_{\rm ^{12}CO}$ variation \citep{Sandstrom:2013qp}.
Indeed, \cite{Accurso:2017il} found that $\alpha_{\rm CO}$ depends not only on metallicity but also on $\Delta$(MS).
In this subsection, we discuss a possible variety in $\alpha_{\rm CO}$ and its impact on our results.

Our data show that galaxies with high sSFR tend to have high $R_{1213}$, suggesting low $\tau_{\rm ^{12}CO}$ owing to high gas temperature and large linewidths.
According to magnetohydrodynamic simulations of molecular clouds, it is claimed that $\alpha_{\rm CO}$ scales with temperature and $^{12}$CO linewidth as $\propto T^{-1/2}$ and $\propto \Delta v^{-1/2}$, respectively \citep{Shetty:2011ai}.
The dust temperatures of the COMING galaxies inferred from $f_{60}/f_{100}$ and a conventional modified black-body model (opacity of $\kappa_0=0.009$~cm$^{2}$~g$^{-1}$ at $\nu_0$=200~GHz, and emissivity index of $\beta=1.8$) ranges $25-40$~K.
The $\sigma_{{\rm ^{12}CO},i}$~ranges $\sim10-100$~km s$^{-1}$, although there also exists low-$\sigma_{{\rm ^{12}CO},i}$~($<10$~km s$^{-1}$) galaxies in our sample.
Assuming gas temperature equals to dust temperature and $\Delta v \sim$~$\sigma_{{\rm ^{12}CO},i}$, these varieties could impose an $\alpha_{\rm CO}$ variety by a factor of $\sim4$ {\it at most} for our sample.
Thus, the correlations between $R_{1213}$ and $\mu_{\rm mol}$ in figure~\ref{fig:basic}c would be actually weaker than the one obtained when using a constant $\alpha_{\rm CO}$.
With a sophisticated technique to derive $\alpha_{\rm CO}$ using CO, \mbox{H\,{\sc i}}~and dust maps of galaxies \citep{Leroy:2011rt,Sandstrom:2013qp}, $\alpha_{\rm CO}$ variations among the COMING galaxies are investigated in a forthcoming paper (Yasuda et al.~in prep.).

There are three caveats here:
(1) the resolutions of our data are not spatially or spectroscopically high enough to resolve a single molecular cloud,
(2) the ``luminosity-weighted'' dust temperature does not necessarily reflect the temperature of the bulk of molecular gas, and
(3) contribution from diffuse molecular gas.
The kpc-scale beam encompasses multiple molecular clouds and the observed velocity dispersion is likely to reflect inter-cloud velocity dispersion rather than that of a single cloud.
Additionally, the dust temperature derived by SED fitting is ``luminosity-weighted'' temperature rather than mass-weighted \citep{Scoville:2014rc,Scoville:2016de}.
But, the ranges of velocity dispersion and temperature we assumed here are not significantly different from the observed values for molecular clouds.
It is known that the velocity dispersions of molecular clouds including those at the central region of the Milky Way range $1-100$~km~s$^{-1}$ \citep[e.g.,][]{Oka:2001yv}.
The multi-transitions CO observations revealed that temperature of molecular clouds is reported to range from $10-20$~K for quiescent clouds to $15-100$~K for clouds with \mbox{H\,{\sc ii}}~regions in M~33 \citep{Wilson:1997mb}.
The assumed varieties of $\Delta v$ and $T$ here are not so inadequate to consider the possible maximum variety in the $\alpha_{\rm CO}$ value.

For the third point, several studies of both the Milky-Way clouds and extragalactic objects suggest that non-negligible fraction of $^{12}$CO comes from diffuse molecular component \citep[$\sim30$~\%; e.g.,][]{Wilson:1994tn,Rosolowsky:2007is,Sawada:2012fv,Pety:2013rr,Morokuma-Matsui:2015cr,Maeda:2020gd,Liszt:2020jj}.
This diffuse molecular component is considered to consist of volume-filling diffuse molecular gas and/or envelope of molecular clouds \citep{Wilson:1994tn,Meier:2000wc,Pety:2013rr}.
Considering that the diffuse component tends to have large velocity dispersion \citep{Garcia-Burillo:1992as, Pety:2013rr, Caldu-Primo:2013lb}, the positive correlations of the sSFR-$\sigma_{{\rm ^{12}CO},i}$ and $R_{1213}$-$\sigma_{{\rm ^{12}CO},i}$ relations indicates a larger fraction of the diffuse component for galaxies with high sSFR and $R_{1213}$.
The CO-to-H$_2$ conversion factor is reported to be similar in diffuse and dense molecular gas {\it on average} with a large scatter in the Milky Way because the small abundance ratio of [CO]/[H$_2$] in diffuse gas is compensated by a much higher brightness per CO molecule \citep{Liszt:2010te,Liszt:2012iz}.
However, it is important to investigate a fraction of the diffuse molecular-gas component of galaxies as a function of sSFR to understand the $\alpha_{\rm CO}$ variation across the main-sequence of SF galaxies.

\subsection{Implications for galaxy evolution on the SFR-Mass parameters space and high redshift galaxies}
\label{sec:galevo}

Our results suggest that there might be an $\alpha_{\rm CO}$ variety across the SF main sequence as a consequence of the differences in temperature and degree of turbulence of molecular gas in galaxies at least for the local Universe.
This further suggests that one may need to pay special attention when discussing either of which $\mu_{\rm mol}$ or SFE is the dominant factor controlling sSFR (or $\Delta$(MS)).
Previous studies found that $\Delta$(MS) and sSFR are correlated with both $\mu_{\rm mol}$ and SFE, i.e., galaxies with higher SFR at fixed $M_{\rm star}$ have more abundant molecular gas and are forming stars more efficiently from the molecular gas \citep{Saintonge:2012nj,Genzel:2015gn,Scoville:2017jw,Tacconi:2018pv}.
The possible $\alpha_{\rm CO}$ variation by a factor of $\sim4$ is unlikely to completely extinguish the reported correlations between sSFR and $\mu_{\rm mol}$ or SFE since the correlations are found over more than 2 orders of magnitude of the both values.
However, a magnitude relationship of the strength of the known positive correlations between sSFR and $\mu_{\rm mol}$ or SFE can be affected, since $M_{\rm mol}$ is numerator and denominator of $\mu_{\rm mol}$ and SFE, respectively.
If galaxies with higher sSFR tend to have lower $\alpha_{\rm CO}$, SFE would become relatively more important for the sSFR (or $\Delta$(MS)) variation than $\mu_{\rm mol}$, as claimed in \cite{Ellison:2020kf}.

The redshift evolution of $\alpha_{\rm CO}$ can be discussed in the same context combined with redshift evolution of metallicity.
sSFR of the main-sequence galaxies at fixed $M_{\rm star}$ is known to increase as redshift increases \citep[e.g.,][]{Whitaker:2012ok,Speagle:2014by}, making one expect higher turbulence and temperature of ISM in the higher-redshift galaxies from this study.
Indeed, high-redshift galaxies tend to have more highly turbulent molecular gas \citep{Genzel:2013pa,Ubler:2018ye,Girard:2019gz}, higher dust temperature \citep{Magdis:2012fh,Magnelli:2014zr,Schreiber:2018fe} and lower metallicity \citep[e.g.,][]{Savaglio:2005kq,Shapley:2005jk,Erb:2006jf,Liu:2008ry,Maiolino:2008ln,Mannucci:2009lp,Yuan:2013vg,Steidel:2014rw} than local galaxies.
With the Milky-Way type $\alpha_{\rm CO}$, a much larger $\mu_{\rm mol}$ value is reported in galaxies at $z\sim1-2$ \cite[$\sim50$~\%, e.g.,][and references therein]{Daddi:2010qq,Tacconi:2013qd,Genzel:2013pa,Riechers:2020nn,Tacconi:2020rs} than local galaxies \cite[$\lesssim5-10$~\%, e.g., ][]{Saintonge:2011hl,Saintonge:2017ve,Bolatto:2017jn,Sorai:2019hs}.
Since $\mu_{\rm mol}$ value is one of the key parameters to constraint on both star-formation and feedback models in cosmological simulations and semi-analytic models of galaxy formation \citep[e.g.,][]{Popping:2012gd,Popping:2014zm,Popping:2015fu,Morokuma-Matsui:2015vn,Saintonge:2017ve}, accurate $\alpha_{\rm CO}$ estimation at higher redshifts is essential.

Taking the three parameters (metallicity, velocity dispersion, and temperature) into consideration, we roughly estimate $\alpha_{\rm CO}$ of galaxies at high redshifts.
For galaxies with $M_{\rm star}\sim10^{10.5}$~M$_\odot$ at $z\sim1.5$, the metallicity, velocity dispersion of molecular gas, and dust temperature can be respectively $\sim0.2$-dex lower \citep{Yabe:2014dh,Zahid:2014om,Kashino:2017ar}\footnote{At the redshift, massive systems with $M_{\rm star}\sim10^{11}$~M$_\odot$ are known to be already enriched to the level observed in local counterparts \citep[][]{Zahid:2014om,Kashino:2017ar}.}, $\sim5$ times higher \citep{Genzel:2013pa}, and only $\sim1.2$~times higher than local counterparts.
In this case, the effects of metallicity and temperature are negligible \citep{Bolatto:2013vn} and $\alpha_{\rm CO}$ could become lower than the local value by a factor of 2.
Indeed, \cite{Carleton:2017jy} obtained comparable $\alpha_{\rm CO}$ for galaxies at $1<z<1.5$ to the local galaxies when they assume several models for redshift evolution of SFE.
Future high spatial- and velocity-resolution observations in $^{12}$CO of high redshift galaxies would shed light on the $\alpha_{\rm CO}$ variety across the SF main sequence of galaxies as a function of redshift.

\section{Summary}
\label{sec:summary}
Using the $^{12}$CO and $^{13}$CO mapping data obtained in the COMING project, we investigate the relationship between $R_{1213}$ and sSFR combined with other four parameters ($M_{\rm star}$, SFR, $f_{60}/f_{100}$, and $\sigma_{{\rm ^{12}CO},i}$) of local star-forming galaxies to understand the molecular gas properties traced with $^{12}$CO line on the SFR-Mass relation of galaxies.
Owing to the stacking analysis, $^{13}$CO emission is successfully detected in 80 galaxies out of the 147 COMING galaxies.
The obtained results and implications are as follows:
\begin{itemize}
\item The error-weighted mean of $R_{1213}$ for the 80 COMING galaxies with measurements of $^{12}$CO and $^{13}$CO is $10.9$ with a standard deviation of $7.0$, which is consistent with previous studies on $R_{1213}$ of local SF galaxies (section~\ref{sec:results}, figure~\ref{fig:r1213_ssfr}).

\item $R_{1213}$ moderately correlates with sSFR with a Spearman's rank correlation coefficient of $0.46$ (section~\ref{sec:results}, figure~\ref{fig:r1213_ssfr}).

\item $R_{1213}$ variation is likely to be attributed to $\tau_{\rm ^{12}CO}$ variation, which is supported by the positive correlations of the $R_{1213}$-$f_{60}/f_{100}$ and $R_{1213}$-$\sigma_{{\rm ^{12}CO},i}$~relations for galaxies with $R_{1213}$ measurements (section~\ref{sec:corrplots}, figures~\ref{fig:comp} and \ref{fig:corr}).

\item The correlations of the sSFR-$f_{60}/f_{100}$ and sSFR-$\sigma_{{\rm ^{12}CO},i}$~relations for galaxies with $R_{1213}$ measurements suggest that high turbulence and temperature of ISM are related to active star formation (section~\ref{sec:corrplots}, figures~\ref{fig:comp} and \ref{fig:corr}).
Stellar bars might play some roles in enhancing the turbulence and in inducing the active star formation at the central regions of galaxies via gas inflow, which is supported by slightly larger median values of $\sigma_{{\rm ^{12}CO},i}$~and sSFR for barred galaxies than non-barred galaxies in our sample although there is a large overlap between them (section~\ref{sec:morph}, figure~\ref{fig:morph} and table~\ref{tab:morph}).

\item The varieties in the dust temperature inferred from $f_{60}/f_{100}$ and the observed velocity dispersion of molecular gas in the COMING galaxies can impose $\alpha_{\rm CO}$ variation by a factor of $\sim4$ {\it at most} according to the relation of $\alpha_{\rm CO} \propto T^{-1/2} \Delta v^{-1/2}$ predicted in the numerical simulations of molecular clouds (sections~\ref{sec:causes} and \ref{sec:implication}).

\item The inferred $\alpha_{\rm CO}$ variation is not likely to completely extinguish the correlations between sSFR and $\mu_{\rm mol}$ or SFE reported in literatures.
However, a magnitude relationship of the strength of the positive correlations between sSFR and $\mu_{\rm mol}$ or SFE can be affected, since $M_{\rm mol}$ is numerator and denominator for $\mu_{\rm mol}$ and SFE, respectively (section~\ref{sec:galevo}).

\end{itemize}

\begin{ack}
We thank the anonymous referee for his/her comments which improve our paper.
This work was supported by JSPS KAKENHI Grant Numbers of 19J40004, 19H01931, 17H01110 and 19H05076.
This work has also been supported in part by the Sumitomo Foundation Fiscal 2018 Grant for Basic Science Research Projects (180923), and the Collaboration Funding of the Institute of Statistical Mathematics ``New Development of the Studies on Galaxy Evolution with a Method of Data Science''.
This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration.
\end{ack}

\bibliographystyle{apj}