\section{\label{sec:level1}First-level heading}
\noindent\textit{\textbf{Introduction}.}
Axions are ultra-light pseudo-scalar particles that are generically predicted in the Peccei-Quinn (PQ) mechanism \cite{Peccei:1977hh,Peccei:1977ur,Wilczek:1977pj}, a compelling solution for the Strong CP problem in the Standard Model of particle physics. In recent years, axions have also attracted substantial interest as a leading dark matter (DM) candidate alternative to WIMPs, and the scope of related studies has extended to axion-like particles (ALPs) beyond the originally proposed QCD axions \cite{Weinberg:1977ma,Abbott:1982af,Dine:1982ah,Preskill:1982cy}.\\
 \indent Understanding the production mechanism and the relic abundance prediction of axions is critical for determining their potential as a viable DM candidate and related phenomenology \cite{Irastorza:2018dyq,Marsh:2015xka}. Despite extensive literature on this subject, our understanding are not yet complete. For instance, for post-inflationary PQ symmetry breaking, axion topological defects (cosmic strings, domain walls) necessarily form and through their subsequent decays may contribute to axion relic abundance ($\Omega_a$) in significant ways \cite{Sikivie:1982qv,Vilenkin:1984ib,Davis:1986xc,Vincent:1996rb,Kawasaki:2014sqa,Vilenkin:1986ku}. Due to the challenges in simulating and analyzing the dynamics of axion topological defects, the prediction for such a potentially important contribution still bears large uncertainties while a growing effort has recently been made \cite{Klaer:2017qhr,Gorghetto:2018myk,Kawasaki:2018bzv,Martins:2018dqg,Buschmann:2019icd,Hindmarsh:2019csc,Martins:2018,Hook:2018dlk}. Meanwhile, our understanding of possible outcomes of the misalignment production may not be complete either. According to the conventional misalignment mechanism, axion field starts with an initial value, $\theta_i$ ($\theta\equiv a/Nf_a$), away from the true vacuum, then begins to oscillate around the minimum when its mass $m_a\sim3H$ ($H$: Hubble expansion rate), and behaves like cold dark matter since then. In order to solve the equation of motion for axion evolution, the initial velocity $\dot{\theta}_i$ also needs to be specified, which is implicitly assumed to be zero in the conventional misalignment, and directly affects the $\Omega_a$ prediction. However, this simple assumption of $\dot{\theta}_i=0$ is only one of the possible initial conditions for $\dot{\theta}_i$. A nonzero $\dot{\theta}_i$ can arise from an explicit breaking of the PQ symmetry effective in the early Universe which is generally expected for an approximate symmetry. \textit{Then how would a nonzero initial velocity of the axion field influence the axion relic abundance and phenomenology?}  
 
In this \textit{Letter} we propose and systematically investigate an alternative misalignment mechanism with $\dot{\theta}_i\neq0$ initial condition. With benchmark examples we demonstrate that in some cases axion relic density prediction can be significantly different from the conventional, with potentially dramatic enhancement or suppression depending on specifics with $\dot{\theta}_i$ and whether the PQ symmetry breaks before/during or after inflation. \\

\noindent\textit{\textbf{The Origin of a Nonzero Initial Velocity}} 
Axion originates from the phase of a complex scalar $\Phi$ whose vacuum expectation value $f_a/\sqrt{2}$ leads to the spontaneous breaking of the $U(1)_{\hbox{\scriptsize{PQ}}}$ symmetry (or a generic global $U(1)$ for ALPs) \footnote{Without loss of generality we focus on the simplest scenario where domain wall number $N=1$.}:
\begin{align}
\label{Eq:PQscalar}
\Phi \equiv \frac{1}{\sqrt{2}} \left( f_a + \phi \right)e^{i a/f_a}
\end{align}
where $\phi$ and $a$ are the radial and angular (axion) modes, respectively. The conserved Noether charge associated with the PQ symmetry is $Q=R^3 f_a^2 \dot{\theta}$ where $R$ is the cosmic scale factor. A $\dot{\theta}_i\neq0$ condition thus corresponds to the rotation of the $\Phi$ field and an asymmetry of the global PQ charge. Such a charge asymmetry can result from higher dimensional operators that explicitly breaks $U(1)_{\hbox{\scriptsize{PQ}}}$ in the early Universe, in analogy to the Affleck-Dine (AD) mechanism for baryogenesis \cite{Affleck:1984fy,Dine:1995uk,Dine:1995kz}. This effect in fact can be generic since global symmetry violating higher dimensional operators is generally expected to be induced by Planck-scale physics \cite{Kamionkowski:1992mf,Dine:1995uk,Dine:1995kz}. Alternatively $\dot{\theta}_i\neq0$ may originate from axion models with a small dimensionful symmetry-breaking term which introduces a slope in axion potential \cite{Graham:2015cka}. Such PQ-breaking effects should be absent today in order not to undermine the solution to the Strong CP problem, which can be realized by tying its strength to the Hubble rate or a dynamical field that has a larger value in the early Universe. Although the specifics of $\dot{\theta}_i$ is model-dependent, important phenomenological insights can be obtained by studying the axion evolution with the benchmark examples of $\dot{\theta}_i$ that we will demonstrate in the paper.\\

\begin{figure}[h]
\includegraphics[width=0.48\textwidth]{Plot-Sin-v3.pdf} 
\caption{\label{F1} Cartoon illustration of the axion field evolution for the two representative possibilities with $\dot{\theta}_i\neq0$, as explained in the text. The green trajectory represents the sequence of motion. The arrows at $t_i$ indicate the direction of initial velocity, $t_o$ is the onset time of oscillation. In the conventional misalignment, the field starts with $\dot{\theta}_i=0$. }
\end{figure}

\noindent\textit{\textbf{Axion Misalignment Mechanism with a Nonzero Initial Velocity.}}
We first present Fig.~\ref{F1} as a cartoon illustration for two representative possibilities of $\dot{\theta}_i \neq0$ initial condition (IC), with related technical details elaborated later in the text. In conventional misalignment, axion field starts at rest with the rescaled field value $\theta_i$, then roll down the periodic potential well, and start oscillating around the true vacuum when $m_a\sim3H$ (at $t_o^{con}$) and behaves like matter. In both cases we show, the field starts with $\theta_i$ at time $t_i$. The lower panel demonstrates the possibility where the axion field has a negative moderate initial velocity to allow it to roll down further in the potential well so that the field value becomes smaller than $\theta_i$ when oscillation begins. This would lead to a suppression of the axion relic density. The upper panel shows the possibility with a high initial velocity, which may delay the onset of oscillation (at $t_o$) to be later than $t^{con}_o$, thus reduce the entropy $s(t_o)$ then and enhance the relic density of the axion ($\Omega_a\propto Y_a\equiv n_a/s$).


We now study the dynamics of axion evolution in details. The equation of motion (e.o.m) of axion field with rescaled $\theta(t) \equiv a(t)/f_a$ (mod $2\pi$) in FRW cosmology is \footnote{We assume that by $t_i$ of our consideration the radial mode has settled into its true vacuum, therefore does not influence axion evolution through coupled e.o.m's. }
\begin{align}
\label{Eq:EoM}
\ddot{\theta} + 3 H \dot{\theta} + m_a^2(T) \theta=0. 
\end{align}
In the conventional misalignment mechanism, initial condition is set as $\dot{\theta}_i = 0$, and the axion field freezes at a random initial field value $\theta_i$ from PQ breaking to QCD phase transition. For post-inflation PQ breaking, axion relic abundance $\Omega_a$ is obtained by averaging over the randomly distributed $\theta_i$ over all causal patches. For QCD axion, we will assume $m_a(T)\propto T^{-4}$ as found by instanton calculation \cite{diCortona:2015ldu, Bonati:2018blm, Petreczky:2016vrs, Burger:2018fvb, Gorghetto:2018ocs}, while a constant $m_a$ may be realistic for general ALPs. For each case in our studies with $\dot{\theta} (t_i)\neq 0$ we will first analyze the simple scenario with constant $m_a$, then consider the QCD axion case.

\begin{figure}[h]
\includegraphics[width=0.48\textwidth]{Axion_Misalignment_v9.png}
\caption{\label{postInf_evol} The time evolution of the axion field with Type-I IC in post-inflation PQ-breaking scenario (for an individual $\theta_i$, before averaging). Green/orange: case (i), blue: case (ii), black: case (iii), red: case (iv). Details given in the text.}
\end{figure}

We start by investigating the axion evolution at early times well before oscillation starts. The starting time $t_i$ is generally assumed to be at PQ breaking scale, but can be later times when the axion picks up a nonzero $\dot{\theta}_i$. In this regime the potential term in Eq.~\ref{Eq:EoM} is negligible. Allowing a nonzero initial velocity $\dot{\theta} (t_i) \neq0$, and dropping the potential term  in Eq.~\ref{Eq:EoM} which is negligible for this early regime, we find the following solution in general cosmology:
\begin{align} 
\label{eq: theta_sol}
\theta(t) =&\,\left\{           
\begin{aligned}
 &\,\theta_i + \frac{\dot{\theta}_i}{H_i} \left( \frac{2}{6-n}\right) \left[1-\left( \frac{R(t)}{R(t_i)} \right)^{n/2-3} \right], \; (n\neq 6)\\
&\, \theta_i + \frac{\dot{\theta}_i}{3H_i} \hbox{ln}\left[ \frac{t}{t_i} \right] , \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\, (n= 6)
\end{aligned}
\right.\\
\dot{\theta}(t) =&\, \dot{\theta_i} \left( \frac{R(t_i)}{R(t)} \right)^3,
\end{align}
where the background energy density $\rho\propto R^{-n}$. We will assume standard cosmology (i.e. radiation domination, $n=4$, thus $\theta(t)\propto1/\sqrt{t}$) except when considering inflationary effects. The energy density of the axion evolves as $\rho_a(t)=\frac{1}{2}\dot{\theta}^2(t)f_a^2+\frac{1}{2}m_a(T)^2\theta(t)^2f_a^2$, and relic density can be estimated as 
\begin{equation}
\Omega_a= m_a(T_o)m_a(T=0)\theta_o^2f_a^2\frac{s_0}{{s(t_o)}\rho_c},
\end{equation}
where $\theta_o\equiv\theta(t_o)$, $s_0$ and $\rho_c$ are the current-day entropy and critical density, respectively.

We now specify two benchmark types of initial condition of $\dot{\theta}_i$ to find concrete form of solutions given by Eqs.~\ref{eq: theta_sol}.\\
\underline{\textbf{Type-I IC: $\dot{\theta}_i=-\delta H_i$,}} where $\delta$ is a constant parameter independent of ${\theta}_i$. To simplify discussions we choose the convention of $\delta\geq0$ without loss of generality \footnote{The late time evolution of the field with $\delta<0$ is nearly the same as their $\delta>0$ counterpart.}. While the detailed UV physics leading to such initial conditions is not our focus here, they may arise from an AD-like scenario \cite{Dine:2003ax}, which readily gives $|\dot{\theta}(t_i)| \sim H(t_i)$ upon PQ-breaking at $t_i \sim m_{\hbox{\scriptsize{pl}}}/f_a^2$ \footnote{The specifics of initial velocity depends on the radial mode in the UV model. The AD mechanism can generate a initial velocity with a varying ratio to the Hubble rate, see e.g. \cite{Dine:2003ax, Akita:2017ecc,Akita:2018zma}.}. Applying this IC to the Eq.~\ref{eq: theta_sol} in RD we find
\begin{align}\label{eq: sol_IC1}
\theta(t)= \theta_i-\delta +\delta \left( \frac{R(t_i)}{R(t)}\right),~~{\rm Mod}[2\pi]
\end{align}
We define the field value at oscillation time as $\theta_o \equiv \theta(t_o)$, where $t_o$ implicitly depends on $\dot{\theta}_i$. In the conventional misalignment, oscillation starts when $m_a(T)\sim 3H$, i.e. $t_o^{con}\sim 1/m_a$. Provided a small/moderate $\dot{\theta}(t_i)$, the oscillation onset $t_o$ in the new scenario is also close to $t_o^{con}$. However, with sufficiently large $\dot{\theta}_i$ the kinetic energy (KE) could be larger than the potential energy $V$ at $t_o^{con}$, and oscillation can only start later when KE becomes comparable to $V$. Such a delayed $t_o$ may enhance axion relic abundance since  $\Omega_a\sim \rho_a(t_o)/s(t_o)$. For high $\dot{\theta_i}$ we can estimate $t_o$ as the earliest time when KE and V become equal ($ \theta_o\approx2\pi$): 
\begin{equation}\label{eq: tosc}
\dot{\theta}_i\left( t_i/t_o \right)^{3/2} \approx \, 2 \pi \, m_a (t_o).
\end{equation}
Combining Eq.~\ref{eq: tosc} and $m_a(t_o)\approx3H(t_o)$, we then can find the critical $\delta_c$: for $\delta>\delta_c$ a notable delay of $t_o$ relative to $t_o^{con}$ would occur: 
\begin{align}
\label{Eq: delta_c}
\delta_c \approx 6\pi  \left({t_o^{con}}/{t_i} \right)^{1/2} \simeq 2 \times 10^{11} \left( \frac{f_a}{10^{11}\,\hbox{GeV}} \right)^{7/6},
\end{align}
where we assumed $t_i$ around the PQ-breaking time $\sim M_{P}/f_a^2$. It is worth noting that $\delta_c$ can be much smaller if $t_i$ (when axion picks up a nonzero velocity) is much later than PQ-breaking time. We are now ready to discuss the evolution patterns in different $\delta$ parameter regions based on Eqs.~\ref{eq: theta_sol}, \ref{eq: sol_IC1}. The two distinct scenarios of PQ breaking before and during/after inflation will be considered in order.\\
\twocolumngrid
\begin{figure}[h] 
\includegraphics[width=0.45\textwidth]{Omega_vs_delta_QCDaxion_v20191124.pdf}
\caption{\label{Fig:OmegaDelta_IC1_post} The dependence of axion relic abundance on initial velocity ($\delta$) for Type-I IC, post-inflation PQ breaking ($f_a=10^{11}\,$GeV). The kink around $\delta\sim10^{20}$ is due to the change in the number of relativistic degrees of freedom, $g_*$, which is accounted for in our numerical calculation.} 
\end{figure}

\noindent\textbf{\textit{ IC-I: Post-inflationary $U(1)_{PQ}$ breaking.}}\\
In each of the four following cases we will first consider the field evolution with a random $\theta_i$, then discuss the impact on $\Omega_a$ after averaging over causal patches after inflation, with the averaged field value at oscillation onset $\langle \theta_o^2\rangle=\langle \theta_i^2 \rangle \sim 2.67\pi^2/3$ \cite{Visinelli:2009zm,Diez-Tejedor:2017ivd}.\\ 
{\underline{(i) $0<\delta<\theta_i$}}, i.e., with generally small initial KE. For generic $\theta_i\sim 1$, although the $\frac{R(t)}{R(t_i)}$ term in Eq.~\ref{eq: theta_sol} slightly reduces $\theta(t)$ over time, the effect is negligible relative to the constant term and $\theta_o\sim\theta_i$. Therefore the axion evolution and the prediction $\Omega_a\propto\theta_o^2$ is very similar to the conventional case, and remains so after averaging over causal patches. This is illustrated with green dashed lines ($\delta \sim 0$, overlapping with conventional case) in Fig.~\ref{postInf_evol} for constant $m_a$, while for QCD axion the evolution before $t_{o}^{con}$ follows the orange line (overlapping with green for $w(t)$ and $\theta(t)$).\\

\noindent{\underline{(ii) $\theta_i<\delta<\delta_c$}}, i.e., with moderate initial KE. In this case KE becomes important at early stages of the field evolution, but is not yet sufficient to cause a notable delay of the $t_o$ relative to $t_o^{con}$. Notice that in this regime as $\delta$ increases $\theta_o$ and thus $\Omega_a$ become more sensitive to $\delta$ (oscillatory dependence) due to the periodic nature of the field potential. In particular, $\Omega_a$ can be much suppressed for particular $\delta$ values that causes cancelation among the $\theta_i$ and $\delta$-dependent terms in Eqs.~\ref{eq: sol_IC1}, i.e., when 
\begin{align}
\label{Eq: thetai-delta}
\theta_o = 2\pi k, \;\;\;\;\; k \in \mathds{Z}.
\end{align}
 However, with a constant $\delta$ such an accidental cancelation only occurs for certain $\theta_i$'s, and its effect is washed out after averaging. Consequently, the $\Omega_a$ prediction is comparable to $\Omega_a^{con}$. The evolution for this case is illustrated with blue lines  in Fig.~\ref{postInf_evol} (for $\theta_i$ away from the cancelation region Eq.~\ref{Eq: thetai-delta}): due to the dominance of KE over V in early times, both the constant $m_a$ and the QCD cases essentially follow same evolution (this also applies to the case (iii) and (iv) below). We can see that due to the large $\theta_i$ the field's equation of state in early stage is kination-like. \\
\begin{figure}[h]
\includegraphics[width=0.48\textwidth]{Axion_Misalignment_inflation_v9.png}
\caption{\label{preInf_evol} The time evolution of the axion field with with Type-I IC, in the pre-inflationary PQ-breaking scenario. Color codes are the same as in Fig.~\ref{postInf_evol}, details given in the text.} 
\end{figure}

\noindent{\underline{(iii) $\delta>\delta_c$}}, i.e., with high initial KE. This case is similar to the above (ii), yet the difference is that KE is reshifted to $\sim V$ after $t_o^{con}$, so the oscillation is delayed. This case is illustrated with black lines in Fig.~\ref{postInf_evol}. Although $\rho(t_o)$ is the same as in case (ii) at the onset of axion oscillation, the entropy then is diluted as $s(t_o)\sim s(t^{con}_o)\left(t_o^{con}/t_o\right)^{3/2}\propto \delta^{-1}$ (using Eqs.~\ref{eq: sol_IC1}, \ref{eq: tosc}). Like in (ii), accidental cancelation Eq.~\ref{Eq: thetai-delta} can happen but becomes irrelevant after averaging. Therefore the relic abundance is enhanced relative to $\Omega_a^{con}$ by a factor of $\sim\delta/\delta_c$.\\

\noindent{\underline{(iv) $\delta=\theta_i$ or $\delta\approx\theta_i$,}} the special regime where $\delta$ is equal to or in the close vicinity of $\theta_i$. This case is qualitatively different from the previous ones. Due to the (almost) cancelation between $\theta_i$ and $\delta$, $\delta \left( \frac{R(t_i)}{R(t)}\right)$ term in Eq.~\ref{eq: sol_IC1} dominates the evolution which causes a potentially dramatic decrease in $\theta(t)$ until oscillation starts around $t_o\sim t_o^{con}$, with $\theta(t_o)=\theta_o\ll\theta_i$. This case is demonstrated in red in Fig.~\ref{postInf_evol} for $\delta=\theta_i$, in yellow for $\delta\approx\theta_i$. $\rho_a(t_o)$ would be expected to be suppressed by $\sim\theta_o^2/\theta_i^2$ relative to the conventional. For post-inflationary PQ breaking as considered here though, with a constant $\delta$ as assumed in IC-II, the cancelation between $\delta$ and $\theta_i$ only occurs for patches with peculiar $\theta_i$, and the effect disappears after averaging. 

\begin{figure}[h] 
\includegraphics[width=0.45\textwidth]{Omega_vs_delta_QCDaxion_v20191126-Pre-Inf_v2.pdf}
\caption{ \label{Fig:OmegaDelta_IC1_pre}The dependence of axion relic abundance on initial velocity ($\delta$) for Type-I IC, pre/during-inflation PQ breaking ($f_a=10^{11}\,$GeV, $t_I = 10^3 t_i$). The oscillatory dependence in large $\delta$ region is shown, and the lowest value of $\Omega_a$ can be zero for $\theta(t_I) = 0$.}
\end{figure}


We show the numerical results for $\Omega_a$-$\delta$ dependence in Fig.~\ref{Fig:OmegaDelta_IC1_post} (with $m(T)\propto T^{-4}$, and the example of $f_a=10^{11}$ GeV). The analytical estimate of $\Omega_a$ is summarized later in Eq.~\ref{Eq: Omega_a}.
\textit{The upshot for the post-inflation scenario is: $\Omega_a$ can be similar to or enhanced relative to the conventional misalignment due to $\dot{\theta}(t_i)\neq0$, while the potential suppression effect in specific patches would be resolved after averaging over post-inflationary causal patches.}

\noindent\textbf{\textit{{IC-I: Before or during inflation $U(1)_{PQ}$ breaking.}}} During inflation the KE in the axion field is rapidly diluted, while $\theta_i$ or the potential energy freezes in as its value at the onset of inflation (or at PQ breaking time if PQ breaks during inflation). Many of the discussions for the post-inflationary scenario apply here, but there are key differences due to inflationary effects. We briefly summarize the results for the same four cases as follows:\\
{\underline{(i) $0<\delta<\theta_i$:}} similar to the conventional pre-inflationary case, for the same reasons as in the post-inflationary scenario.\\
{\underline{(ii) $\theta_i<\delta<\delta_c$:}} in general similar to the conventional case, but unlike the post-inflationary scenario, the accidental cancelation/suppression on $\Omega_a$ ($\theta_o\approx 0$) persists without the averaging. \\
{\underline{(iii) $\delta>\delta_c$:}} the situation with $\rho_a$ evolution is similar to the above (iii), but the enhancement due to the diluted $s(t_o)$ is absent due to the intervention of inflation which cuts short the KE domination time (unless inflation happens after QCD phase transition). Therefore the result is in general similar to the conventional, but accidental suppression is possible for certain $\delta$.\\
{\underline{(iv) $\delta=\theta_i$ or $\delta\approx\theta_i$.}} For PQ breaking during inflation, the situation is similar to the conventional, since the initial KE is quickly depleted by inflation before it can drive down $\theta_i$ value. However, if the PQ breaks before inflation with a moderate/large separation in their scales, the $\dot{\theta}(t_i)\neq0$ initial condition can leave a trace despite inflation: the field value is already driven down to $\theta_I \equiv\theta(t_I)\approx\frac{R(t_i)}{R(t_I)}\theta_i$ for $\theta_i-\delta\rightarrow0$, where $t_I$ is when inflation starts.   $\theta(t_I)$ then freezes in as the effective new initial condition for axion misalignment after inflation. \\
\indent We illustrate the representative field evolutions in this pre/during inflation scenario in Fig.~\ref{preInf_evol}, with the same color codes as in their post-inflationary counterparts. \textit{The upshot for this scenario is: $\Omega_a$ is similar to or suppressed relative to the conventional due to $\dot{\theta}(t_i)\neq0$ , while the potential enhancement effective in post-inflationary case is generally wiped out by inflation.} 
The $\Omega_a$-$\delta$ relation for this scenario is illustrated in Fig.~\ref{Fig:OmegaDelta_IC1_pre} based on numerical solution with $m(T)\propto T^{-4}$. 
\begin{figure}[h] 
\includegraphics[width=0.46\textwidth]{Omega_vs_delta_QCDaxion_v20191127-Fine-Tuned_v5.pdf}
\caption{\label{Fig: Omega_IC2_postI}. The dependence of axion relic abundance on initial velocity ($\gamma$) for Type-II IC ($f_a=10^{11}$ GeV), post-inflationary PQ-breaking. }
\end{figure}

We now summarize the prediction for $\Omega_a$ with IC-I as follows for post- and pre-inflation cases in order:
\begin{align}
\label{Eq: Omega_a}
\Omega_a^{\rm post-I} \simeq \left\{            
\begin{aligned}
&\,\Omega_a^{con} \simeq 0.02\, \langle\theta_o^2\rangle\left( \frac{f_a}{10^{11}\,\hbox{GeV}} \right)^{7/6}, \hspace{0.1cm} \delta < \delta_c\\*[-0.1\normalbaselineskip]
&\, \Omega_a^{con} \frac{\delta}{\delta_c} \simeq 0.01\,  \langle\theta_o^2\rangle \left( \frac{\delta}{10^{11}} \right),\hspace{0.9cm}  \delta \geq \delta_c,
\end{aligned}
\right.
\end{align}
and
\begin{align}
\Omega_a^{\rm pre-I} = \,\Omega_a^{con} \frac{\theta_I^2}{\langle\theta_o^2\rangle} \simeq \,\Omega_a^{con} \frac{(\theta_i-\delta)^2}{\langle\theta_o^2\rangle},
\end{align}
 Note that for the last equality in the 2nd line, we used Eq.~\ref{Eq: delta_c} ($\delta_c\propto f_a^{7/6}$) which assumes $t_i$ at PQ breaking scale. Other $t_i$ choices would change the numerics. 
\\
Notice that \textit{case (iv) of this scenario provides a new possibility for $\Omega_a$ to account for $\Omega_{DM}$ with $f_a\gtrsim10^{11}$ GeV} due to specific relations between $\theta_i$ and $\dot{\theta}_i$ even with a natural $\theta_i\sim1$. On the other hand, one could argue that this is trading one type of fine-tuning to another type. Although challenging it is curious to see if it is possible to distinguish the two types of fine-tuning combining various avenues of observational data. \\

\textbf{\underline{Type-II IC: $\dot{\theta}_i=-(1-\gamma)\theta_iH_i$}}, where $\gamma \geq 0$. This IC is inspired by case (iv) with Type-I IC. Most results in IC-I apply, with the correspondence of $\delta\rightarrow (1-\gamma)\theta_i$. The key difference though is that here the relation $\dot{\theta}_i\propto\theta_i$ is assumed to be valid for any $\theta_i$, \textit{therefore the aforementioned suppression effect in case (iv) (here $\gamma\rightarrow 0$) is robust and survives even after averaging patches for post-inflationary PQ breaking.} The realization of such an IC generally requires an explicit PQ breaking term in the axion potential effective in the early Universe, which is beyond the scope of this work yet is interesting to explore. \textit{We highlight this possibility since it provides a novel dynamic way to accommodate large $f_a\gtrsim10^{11}$ GeV QCD axion as DM candidate for post-inflation PQ breaking.} This solution is intriguing despite requiring a special relation between $\theta_i$ and $\dot{\theta}_i$, especially considering that for post-inflationary PQ breaking in the conventional misalignment there is no way to even fine-tune to accommodate $f_a\gtrsim10^{11}$ GeV, since $\theta_i$ is averaged to $O(1)$!  We show the $\Omega_a$-$\gamma$ relation for IC-II in Figs.~\ref{Fig: Omega_IC2_postI} and \ref{Fig: Omega_IC2_preI}, for post- and pre-inflation PQ-breaking, respectively. \\
\begin{figure}[h] 
\includegraphics[width=0.46\textwidth]{Omega_vs_delta_QCDaxion_v20191126-Fine-Tuned-Pre-Inflation.pdf}
\caption{\label{Fig: Omega_IC2_preI}. The dependence of axion relic abundance on initial velocity ($\gamma$) for Type-II IC ($f_a=10^{11}$ GeV and $t_I = 10^3 t_i$), pre/during-inflationary PQ-breaking. }
\end{figure}
The difference from IC-I in $\Omega_a$ prediction for case (iv) (small $\gamma$) in the post-inflation scenario is demonstrated with the following formula:
\begin{align}
\Omega_a^{\rm post-I} \simeq \left\{            
\begin{aligned}& \,\Omega^{con}_a \gamma^2, \hspace{1.3cm} \frac{R(t_i)}{R(t_{o})}  \leq \gamma \leq 1\\
&\, \Omega^{con}_a  \left(\frac{R(t_i)}{R(t_{o})} \right)^{2}  , \hspace{0.45cm} 0 \leq \gamma \leq \frac{R(t_i)}{R(t_{o})},  
\end{aligned}
\right.
\end{align}
where we can clearly see the suppression relative to the conventional by a factor of $\gamma^2$ or $\left( \frac{R(t_i)}{R(t_{o})}\right)^2$.\\

\noindent\textit{\textbf{Conclusion.}} 
In this \textit{Letter}, we propose a new misalignment mechanism where the axion initial velocity is non-zero and demonstrate its impact on axion relic abundance for a variety of benchmark cases. While in many cases $\Omega_a$ remains similar to the conventional prediction, it may be significantly enhanced with a large initial velocity $\dot{\theta}_i$ or suppressed when $\dot{\theta}_i$ and $\theta_i$ satisfy peculiar relations. Such $\dot{\theta}_i\neq0$ initial conditions may arise from explicitly PQ-breaking operators active in the early Universe which are well-motivated from the understanding that the PQ symmetry is an accidental symmetry. We will investigate the detailed realization of these initial conditions in future work. As an outcome of this new scenario, new viable parameter space for the QCD axion DM opens up, allowing $f_a$ both much above and much below the conventional $\sim10^{11}$ GeV scale. Further phenomenological consequences for both QCD axion and general ALPs are worth further exploration. Meanwhile a caveat for the scenario of post-inflationary PQ-breaking is that, as in conventional misalignment, topological defects' contribution to the axion relic density may dominate over the misalignment contribution, and is worth revisiting to advance and complete our understanding of the axion production mechanism \cite{Kawasaki:2014sqa, Chang:2019mza, furture, Gorghetto:2018myk, Buschmann:2019icd}. \\

\noindent\textit{{Side Notes: During the final stage of completing this manuscript, \cite{Co:2019jts} came out which has overlap with our results. The two papers are complementary to each other with different emphases/approaches, of note we demonstrate more possible ways that axion relic density may be affected by the new misalignment mechanism.}}\\
\begin{acknowledgments}
\noindent \textit{{\textbf{Acknowledgments.}}}
We thank Anson Hook for discussion. The authors are supported in part by the US Department of Energy grant DE-SC0008541. CC thanks Yin Chin Foundation of U.S.A. for its support. YC thanks the Kavli Institute for Theoretical Physics (supported by the National Science Foundation under Grant No. NSF PHY-1748958), Erwin Schr\"odinger International Institute, and Galileo Galilei Institute for hospitality while the work was being completed.
\end{acknowledgments}



\bibliographystyle{apsrev4-1}
