\section{Introduction}

The first decades of the XXth century witnessed the ascension of the general relativity (GR) as the revolutionary paradigm for the gravitational interaction. Among all possible adjectives that can be assigned to GR, it is clearly a simple theory. Simple, in the sense that its field equations are obtained from the standard variational principle where the gravitational Lagrangian equals the Ricci scalar $\cal{L}\propto R$ which is the simplest Lagrangian in four dimensions (up to a cosmological constant term) leading to second order differential equations.

Although GR remains being considered the standard gravitational theory the observations of galaxy rotation curves and the inference of the accelerated expansion of the universe lead to the concept of dark matter and dark energy phenomena, respectively. The conclusion from these facts is that either the energy content of the universe is composed by strange forms of particles/fields or GR fails in describing the dynamics of galaxies and other cosmological observables. The latter assumption has led to the construction of several alternatives to GR. Most of them based on the fact that the gravitational Lagrangian has a non trivial dependence on geometrical quantities or even matter fields. There are also theories which relax some of the fundamental pillars over which GR has been built.

In this work we focus on the cosmological aspects of a recent gravitational theory proposed in \cite{lazo} in which dissipative processes are incorporated in the gravity by means of a generalization of the least action principle. This new theory of gravity is discussed in the next section. In section III we explore the cosmological scenario emerging from the modified field equations and derive the dynamical equations at background and perturbative levels. The results are compared either with the usual viscous fluid formulation and with the $\Lambda$CDM model.
In the section IV we bring a discussion on the perturbative aspect of such cosmology using the synchronous gauge. Lastly, the section V is dedicated to our concluding remarks.

\section{General equations}

A geometrical viscous gravity model can be obtained from first principles by an action-dependent Lagrangian formulation \cite{lazo}. A generalization of the action Principle for Action-dependent Lagrangians was introduced for the first time in the 30's by Herglotz in order to give a variational principle to dissipative phenomena \cite{her1,her2}. More precisely, the original Herglotz variational problem consists in the problem of determining the path $x(t)$ that extremize an Action of the form
\begin{eqnarray}
\label{H}
S = \int {\cal L}(x, \dot x, S) d t.
\end{eqnarray}
Herglotz proved \cite{her1,her2} that a necessary condition for a path $x(t)$ to be an extremizer of the variational problem (\ref{H}) is given by the generalized Euler-Lagrange equation
\begin{eqnarray}
\label{HEL}
\frac{d}{dt}\frac{\partial {\cal L}}{\partial \dot x} - \frac{\partial {\cal L}}{\partial x}  - \frac{\partial{\cal L}}{\partial S}\frac{\partial {\cal L}}{\partial \dot x} = 0.
\end{eqnarray}
The application of Herglotz problem to non-conservative systems is evident even in the simplest case where the dependence of the Lagrangian function on the Action is linear. For example, the Lagrangian ${\cal L}=\frac{m\dot{x}^2}{2}-U(x)-\frac{\gamma}{m} S$ describes a particle under viscous forces and, from (\ref{HEL}), the resulting equation of motion includes the well known dissipative force $\gamma \dot{x}$. In this context, the term linear on $S$ in the Lagrangian can be interpreted as a potential function for the non-conservative force.

However, despite the Herglotz problem was introduced in 1930, a covariant generalization of (\ref{H}) for several variables was obtained only recently \cite{lazo}. The Lagrangian including the gravity sector considered in \cite{lazo} is given by
\begin{eqnarray}
{\cal L} = \sqrt{-g}(R - \lambda_{\mu}s^\mu) + {\cal L}_m,
\end{eqnarray}
where $s^\mu$ is an action-density field and $\lambda_\mu$ is coupling term which may depend on the space-time coordinate. The development of this proposal is given in details in Ref. \cite{lazo}. The inclusion of term $\lambda_{\mu}s^{\mu}$ can be seen as a covariant implementation of the linear dependence on the action approached by Herglotz in the classical mechanics context. So that, it is expected that such a term carries the dissipative nature of the theory. In this sense, the introduction of the two four-vectors $\lambda^{\mu}$ and $s^{\mu}$ revealed the most simple way to perform such generalization, and even a natural choice as a careful look at \cite{lazo} can show us. Notice that in that reference the proposal of the authors consisted in providing a covariant version for the set of classical equations (1). However, it implies in generalizing the time derivative of the action appearing in the first of the equations (1). A simple way to implement it could be by extending this time derivative of $S$ to a divergence of a certain auxiliary four-vector $s^{\mu}$, which was introduced by hand. For $s^{\mu}$ the classical action would be associated somehow with its component $s^{0}$, as the authors briefly discuss in the paragraph just below the equations (8). On the other hand $\lambda^{\mu}$ is just a backgroung four-vector, playing the role of a coupling parameter associated with the dependence of the gravitational lagrangian upon the action. As a first attempt, it is assumed constant, although it could be coordinate-dependent in a more general scenario. One expects that the symmetry of the problem which is being studied helps us to define a specific form for this four-vetor. Besides, different from $s^{\mu}$, which disappears during the variation of the action, $\lambda^{\mu}$ survives and appears in the field equations (see below) as the actual parameter of this model, encoding the dissipative properties that this gravitational theory will manifest. Therefore, any observable quantity of this theory shall be given in terms of $\lambda^{\mu}$.

The resulting field equations are given by \cite{lazo},
\begin{eqnarray}
\label{ei}
R_{\mu\nu} + K_{\mu\nu} - \frac{1}{2}g_{\mu\nu}(R + K) = 8\pi G T_{\mu\nu},
\end{eqnarray}
where $R_{\mu\nu}$ and $T_{\mu\nu}$ are the Ricci and the Hilbert stress-energy tensors, respectively, and
\begin{eqnarray}
K_{\mu\nu} = \lambda_\rho \Gamma^\rho_{\mu\nu} - \frac{1}{2}\biggr(\lambda_\mu \Gamma^\rho_{\nu\rho} + \lambda_\nu \Gamma^\rho_{\mu\rho}\biggl)
\end{eqnarray}
is a tensor related to viscous geometric dissipations.

Since the Bianchi identities are still valid, when the divergence of (\ref{ei}) is taken, some relations involving the tensor $K_{\mu\nu}$, its trace and matter sector are obtained. One possibility is to suppose that the gravitational coupling is not constant \cite{lazo}. However, this implies to consider another field responsible for the evolution of $G$, as in the Brans-Dicke theory \cite{BD}. Another possibility is to consider that the divergence of the energy-momentum tensor is not zero, somehow as in the Rastall theory \cite{Rastall}, or as in the Brans-Dicke theory reformulated in the Einstein frame through a conformal transformation.
In what follows we will consider that $G$ is constant. In this case, the usual conservation equations are replaced by,
\begin{eqnarray}
{K^{\mu\nu}}_{;\mu} - \frac{K^{;\nu}}{2} = 8\pi G {T^{\mu\nu}}_{;\mu}.
\end{eqnarray}
This relation complement the field equations (\ref{ei}). In contrast to the Rastall theory, however, now the non-conservation of the energy-momentum tensor has a geometrical origin, with a basis on a variational principle.

Finally, it is important to note that the introduction of a preferential direction given by the coupling vector $\lambda_\mu$ breaks the space-time symmetry resulting in the existence of preferred reference frames, while it still remains compatible with the main principles of the general theory of relativity. Moreover, this feature should be expected in a dissipative theory since the non-conservation of energy-momentum tensor is directly related to the space-time symmetry broken. The relation between energy-momentum non-conservation and space-time symmetry broken by the introduction of a preferred referential frame is also found in the gravitational aether scenario introduced in the context of the cosmological constant problem \cite{aether,Wondrak}. However, while in the gravitational aether theory the space-time symmetry is broken by the introduction of a vector field, the "aether", on the matter side of Einstein field equations, in the present work the symmetry is broken by the geometric coupling $\lambda_\mu$.

\section{Background equations}

Choose the flat Friedmann-Lema\^{\i}tre-Robertson-Walker metric:
\begin{eqnarray}
ds^2 = dt^2 - a(t)^2[dx^2 + dy^2 + dz^2].
\end{eqnarray}
Moreover, let us choose an {\it Ansatz} for $\lambda_\mu$:
\begin{eqnarray}
\lambda_0 &=& const. \neq 0,\\
\lambda_i &=& 0.
\end{eqnarray}

We find:
\begin{eqnarray}
R_{00} &=& - 3\frac{\ddot a}{a} = -3(\dot H + H^2), \\
R_{ij} &=& (a\ddot a + \dot a^2)\delta_{ij} = (\dot H + 3H^2)a^2 \delta_{ij},\\
R &=& - 6\biggr(\frac{\ddot a}{a} + \frac{\dot a^2}{a^2}\biggl) = - 6(\dot H + 2H^2),\\
K_{00} &=& - 3\lambda_0 \frac{\dot a}{a} = - 3\lambda_0H,\\
K_{ij} &=& \lambda_0a\dot a\delta_{ij} = \lambda_0 Ha^2\delta_{ij},\\
K &=& - 6\lambda_0\frac{\dot a}{a} = - 6 \lambda_0 H.
\end{eqnarray}

The energy-momentum tensor $T^{\mu\nu}$ are that of a perfect fluid,
\begin{eqnarray}
T^{\mu\nu} = (\rho + p)u^\mu u^\nu - pg^{\mu\nu},
\end{eqnarray}
with the non null components,
\begin{eqnarray}
T^{00} &=& \rho,\\
T^{ij} &=& p a^{-2}\delta^{ij}.
\end{eqnarray}

The equations of motion are now,
\begin{eqnarray}
3H^2 &=& 8\pi G\rho,\\
2\dot H + 3H^2 + 2\lambda_0H &=& - 8\pi Gp,
\end{eqnarray}
where we have defined,
\begin{eqnarray}
H = \frac{\dot a}{a}.
\end{eqnarray}

On the other hand, we will consider that the Bianchi identities imply,
\begin{eqnarray}
K^{\mu\nu}_{;\mu} - \frac{K^{;\nu}}{2} = 8\pi GT^{\mu\nu}_{;\mu} &=& 0.
\end{eqnarray}
This implies, using the previous component, the equation,
\begin{eqnarray}
8\pi G\{\dot\rho + 3H(\rho + p)\} = - 6\lambda_0H^2.
\end{eqnarray}

Defining,
\begin{eqnarray}
\frac{2}{3}\lambda_0 = - 8\pi G\xi_0,
\end{eqnarray}
the ensemble of equations take the following form:
\begin{eqnarray}
3H^2 &=& 8\pi G\rho,\\
2\dot H + 3H^2 &=& - 8\pi G(p - 3\xi_0H),\\
\dot\rho + 3H(\rho + p - 3\xi_0H) &=& 0.
\end{eqnarray}

Remember that the bulk viscosity (Eckart's theory) leads to a pressure,
\begin{eqnarray}
p^* &=& p - \xi(\rho)u^\mu_{;\mu} \nonumber\\
&=& p - 3\xi(\rho)H.
\end{eqnarray}
Hence, the previous construction corresponds to a constant bulk viscosity coefficient:
\begin{eqnarray}
\xi(\rho) = \xi_0.
\end{eqnarray}
The bulk viscosity coefficient $\xi_0$ must be positive, implying that $\lambda_0 < 0$ in order to retain the analogy. In fact there is a more fundamental reason to impose a negative sign for $\lambda_0$: the entropy production observed in this scenario, an extra aspect which reinforce the resemblance to the non-causal viscous model.
Following the procedure properly detailed in \cite{winfried,roy} one finds the time evolution of the specific entropy
\begin{equation}
\label{entropy}
\dot{s}=-\frac{2\lambda_{0}\rho}{nT},
\end{equation}
where $T$ and $n$ are the temperarature and the particle number density, respectively. So, $\lambda_0$ is obliged to have necessarily negative sign, in order to ensure a non-negative entropy rate production predicted by the sencond law of thermodynamics.

\section{Perturbed equations}

Now we shall perturb the model described in the previous section, by introducing small fluctuations around the background metric:
\begin{eqnarray}
\tilde g_{\mu\nu} = g_{\mu\nu} + \delta g_{\mu\nu},
\end{eqnarray}
where $\tilde g_{\mu\nu}$ is the inhomogeneous, metric $g_{\mu\nu}$ is the background metric and $\delta g_{\mu\nu}$ represents the fluctuation around it.
From now on, we will note,
\begin{eqnarray}
h_{\mu\nu} \equiv \delta g_{\mu\nu}.
\end{eqnarray}
Due to the inverse metric relation,
\begin{eqnarray}
g^{\mu\rho}g_{\rho\nu} = \delta^\mu_\nu,
\end{eqnarray}
we have,
\begin{eqnarray}
\delta g^{\mu\nu} = - h^{\mu\nu}, \quad h^{\mu\nu} = g^{\mu\rho}g^{\mu\sigma}h_{\rho\sigma}.
\end{eqnarray}

We will work in the synchronous coordinate condition:
\begin{eqnarray}
h_{\mu0} = 0.
\end{eqnarray}
It comes out more convenient, using the synchronous coordinate condition, to rewrite the field equations as,
\begin{eqnarray}
R_{\mu\nu} + K_{\mu\nu} &=& 8\pi G\biggr[T_{\mu\nu} - \frac{1}{2}g_{\mu\nu}T\biggl],\\
K^{\mu\nu}_{;\mu} - \frac{K^{;\nu}}{2} &=& 8\pi G T^{\mu\nu}_{;\mu}.
\end{eqnarray}

The perturbed field equations are:
\begin{eqnarray}
\delta R_{\mu\nu} + \delta K_{\mu\nu} &=& 8\pi G\biggr[\delta T_{\mu\nu} - \frac{1}{2}h_{\mu\nu}T - \frac{1}{2}g_{\mu\nu}\delta T\biggl],\\
\delta (K^{\mu\nu}_{;\mu}) - \delta\biggr(\frac{K^{;\nu}}{2}\biggr) &=& 8\pi G \delta (T^{\mu\nu}_{;\mu}).
\end{eqnarray}

The pertubation of the Ricci and $K_{\mu\nu}$ tensors read,
\begin{eqnarray}
\delta R_{\mu\nu} &=& \partial \delta\Gamma^\rho_{\mu\nu} -\partial_\nu\delta\Gamma^\rho_{\mu\rho} + \Gamma^\rho_{\rho\sigma}\delta\Gamma^\sigma_{\mu\nu}\nonumber\\
&-& \Gamma^\sigma_{\rho\mu}\delta\Gamma^\rho_{\sigma\nu} - \Gamma^\sigma_{\nu\rho}\delta\Gamma^\rho_{\mu\sigma} + \Gamma^\sigma_{\mu\nu}\delta\Gamma^\rho_{\sigma\rho},\\
\delta K_{\mu\nu} &=& \lambda_\alpha \delta\Gamma^\alpha_{\mu\nu} - \frac{1}{2}\biggr(\lambda_\mu\delta\Gamma^\alpha_{\nu\alpha} + \lambda \delta\Gamma^\alpha_{\mu\alpha}\biggl),
\end{eqnarray}
where we have supposed that $\lambda_\alpha$ is constant. In the above expressions, the perturbation of the Christoffel symbol reads,
\begin{eqnarray}
\delta\Gamma^\rho_{\mu\nu} = \frac{1}{2}g^{\rho\sigma}\biggr(\partial_\mu h_{\nu\sigma} + \partial_\nu h_{\mu\sigma} - \partial_\sigma h_{\mu\nu} - 2\Gamma^\lambda_{\mu\nu}h_{\sigma\lambda}\biggl).
\end{eqnarray}

The relevant components of the perturbed Ricci tensor for the study of the scalar modes are:
\begin{eqnarray}
\delta R_{00} &=& \frac{\ddot h}{2} + H\dot h, \\
\delta R_{0i} &=& \frac{1}{2}\biggr(\partial_i \dot h - \partial_k\bar h_{ki}\biggl),
\end{eqnarray}
where,
\begin{eqnarray}
h \equiv \frac{h_{kk}}{a^2}, \quad \bar h_{ij} \equiv \frac{h_{ij}}{a^2}.
\end{eqnarray}

The non-null components of the perturbed $K_{\mu\nu}$ tensor are:
\begin{eqnarray}
\delta K_{00} &=& \lambda_0\frac{\dot h}{2}, \\
\delta K_{0i} &=& \frac{\lambda_0}{4}\partial_i h,\\
\delta K_{ij} &=& - \lambda_0 \frac{\dot h_{ij}}{2},
\end{eqnarray}
implying,
\begin{eqnarray}
\delta K &=& \lambda_0\dot h.
\end{eqnarray}
It is useful also to write the perturbations of $K^{\mu\nu}$ (the contravariant form):
\begin{eqnarray}
\delta K^{00} &=& \lambda_0\frac{\dot h}{2}, \\
\delta K^{0i} &=& - \frac{\lambda_0}{4a^2}\partial_i h,\\
\delta K^{ij} &=& 2\frac{\lambda_0 H}{a^4}h_{ij}- \frac{\lambda_0}{a^4} \frac{\dot h_{ij}}{2}.
\end{eqnarray}

The non-null components of the energy-momentum tensor are:
\begin{eqnarray}
\delta T^{00} &=& \delta\rho, \\
\delta T^{i0} &=& (\rho + p)\delta u^i, \\
\delta T^{ij} &=& \frac{1}{a^4}\biggr(ph_{ij} + \delta p\,a^2\delta_{ij}\biggl),
\end{eqnarray}
implying,
\begin{eqnarray}
\delta T = \delta\rho - 3\delta p.
\end{eqnarray}

In computing all these expressions we have used, of course, the synchronous coordinate condition and the fact that $\lambda_0$ is constant.

The final set of perturbed equations are:
\begin{eqnarray}
&&\ddot h + (2H + \lambda_0)\dot h = 3H^2(1 + 3v_s^2)\delta, \\
&&\partial_i\dot h + \frac{\lambda_0}{2}\partial_i h - \partial_k\dot{\bar h}_{ki} = - 6H^2(1 + \omega)\delta u^i a^2,\\
&&2H\lambda_0\dot h - \frac{\lambda_0}{4a^2}\nabla^2h = 3H^2\biggr\{\dot\delta + [3H(v_s^2 - \omega) - 2\lambda_0]\delta
+ (1 + \omega)\biggr(\theta - \frac{\dot h}{2}\biggl)\biggl\},\\
\frac{\lambda_0}{4a^2}\partial_i\dot h &-& \frac{3}{4} H\frac{\lambda_0}{a^2}\partial_i h - \frac{\lambda_0}{2a^2}\partial_k\dot{\bar h}_{ki}=
3H^2\biggr\{(1 + \omega)\delta\dot u^i + [(1 + \omega)\frac{\dot\rho}{\rho}+5(1 + \omega)H]\delta u^i + \frac{v_s^2}{a^2}\partial_i\delta\biggl\}.
\end{eqnarray}

In these expressions, we have defined,
\begin{eqnarray}
\delta &=& \frac{\delta\rho}{\rho},\\
\theta &=& \partial_k\delta u^k,\\
\omega &=& \frac{p}{\rho}, \\
v_s^2 &=& \frac{\delta p}{\delta\rho}.
\end{eqnarray}

Using the background relations,
\begin{eqnarray}
&&3H^2 = 8\pi G\rho,\\
\dot\rho &=& - 2\lambda_0 - 3H(1 + \omega),
\end{eqnarray}
defining
\begin{eqnarray}
\biggr(\frac{\partial_k\partial_l h_{kl}}{a^2}\biggl)^. = g,
\end{eqnarray}
and performing a Fourier mode decomposition, we have the following set of equations:
\begin{eqnarray}
&&\ddot h + (2H + \lambda_0)\dot h = 3H^2(1 + 3v_s^2)\delta, \\
&&k^2(\dot h + \frac{\lambda_0}{2}h) + g = 6H^2(1 + \omega)\theta a^2,\\
&&2H\lambda_0\dot h + \frac{\lambda_0}{4a^2}k^2h = 3H^2\biggr\{\dot\delta + [3H(v_s^2 - \omega) - 2\lambda_0]\delta + (1 + \omega)\biggr(\theta - \frac{\dot h}{2}\biggl)\biggl\},\\
&-& \frac{\lambda_0}{4a^2}\biggr\{k^2[\dot h - 3 Hh] + 2g\biggl\} =
3H^2\biggr\{(1 + \omega)\dot\theta+(1 + \omega)[- 2\lambda_0 + (2 - 3\omega)H]\theta - k^2\frac{v_s^2}{a^2}\delta\biggl\}.
\end{eqnarray}

On the other hand the background equations admit (for the one fluid case) the analytical solution,
\begin{eqnarray}
a = \biggr\{ - c\biggr[e^{- \lambda_0 (t - t_0)} - 1\biggl] + 1\biggl\}^\frac{2}{3(1 + \omega)},
\end{eqnarray}
where $c$ is a constant and $t_0$ is the present time.

For the zero pressure case ($\omega = v_s^2 = 0$) the system of perturbed equations reduces to,

\begin{eqnarray}
&&\ddot h + (2H + \lambda_0)\dot h = 3H^2\delta, \\
&&\dot\delta - 2\lambda_0\delta + \biggr(\theta - \frac{\dot h}{2}\biggl)=\frac{2\lambda_0}{3H}\dot h + \frac{\lambda_0}{12a^2}\frac{k^2}{H^2}h,\\
&&\dot\theta + [ 2H- \lambda_0]\theta =\frac{\lambda_0 k^2}{12 H^2 a^2}[\dot h + (3H + \lambda_0)h].
\end{eqnarray}

We can promote a direct comparison between the above set of equations with the case of a single viscous fluid

\begin{eqnarray}
&&\ddot h + \biggr(2H - \frac{H_0\bar\xi_0}{2}\biggl)\dot h = 3H^2\delta - H_0\bar\xi_0\theta,\\
&&\dot\delta + H_0 \bar\xi_0\delta + \biggr(1-\frac{2 H_0 \bar\xi_0}{3 H}\biggl)\biggr(\theta - \frac{\dot h}{2}\biggl) = 0, \\
&&\dot\theta + \biggr(2 H + \frac{H_0\bar\xi_0}{2}\biggl)\theta = \frac{H_0\bar\xi_0 k^2}{6\dot H a^2}\biggr(\theta - \frac{\dot h}{2}\biggl).
\end{eqnarray}

In the above equations $(76)-(78)$ the bulk viscous parameter has been redefined as the dimensionless parameter $\bar\xi_0 = 24 \pi G \xi_0 / H_0$. It is worth noting that making $\lambda_0 =0$ the set $(73)-(75)$ coincides with the equations $(76)-(78)$ for $\xi_0=0$. This corresponds to the pressureless Cold Dark Matter (CDM) case.

We solve numerically the above systems of equations in order to obtain the evolution of the linear density contrast $\delta$. In Fig.1 we show the behavior for the CDM model ($\delta\sim a$) in the red line. For the geometrical model we fix the parameter $\lambda_0 =-0.001$ (in $H_0$ units) in the black dashed line of the left panel. In the right panel we adopt $\lambda_0=-0.01$. The equivalent (at background level) viscous models have the bulk viscous parameter $\bar\xi_0=+0.002$ (left panel) and $\bar\xi_0=+0.02$ (right panel). The growth behavior for such bulk viscous matter is seen in the blue line in both panels. Both curves have the same initial conditions which is equivalent to a $k=0.2 h Mpc^{-1}$ mode deep in the matter dominated epoch. As expected such scale just recently entered the nonlinear regime $\delta\sim 1$.

The curves for the viscous model show the expected pathological behavior as already shown in Refs. \cite{Velten:2011bg} (see also \cite{Li:2009mf}
). Structure growth is highly suppressed in pure viscous cosmologies. On the other hand, the geometrical model follows the CDM behavior. Indeed by increasing the magnitude of $\lambda_0$ a suppression is expected. Let us recall that a slight suppression of the matter clustering on small scales can be important to alleviate the tensions of the $\Lambda$CDM model due to an excess of power verified in numerical simulations \cite{sats}.

In the same spirit of several proposals presented in the last decades, this model is an attempt to abandon cosmological constant, by replacing it by an another gravitational theory. In the present case the cosmic acceleration process would be in charge of the dissipative effects naturally emerging in this new gravity. Dissipative processes are commonly found in the nature, and there is no reason to discard them at all as a possible ingredient in the dynamics of the universe. However, it is reasonable to include them in the cosmology through a viable cosmological model, consistent with the observational data. The triumph of $\Lambda$CDM in this aspect place it as a fiducial model to be compared with. This duty forces us to look at some of the main existing datasets, both at background and at perturbative level, in order to pursue a suitable comparison with the standard scenario and also check the consistency of our model. Such data analysis is going to be performed in a future investigation.

\begin{figure}[t]
\includegraphics[width=0.4\textwidth]{deltas.pdf}
\includegraphics[width=0.4\textwidth]{deltas2.pdf}
\caption{Evolution of the matter density contrast as a function of the scale factor $a$. Red line represents the CDM model. For the non-conservative model black-dashed lines are used with values $\lambda_0=-0.001$ (left panel) and $\lambda_0=-0.01$ (right panel). The blue line represents the case in which dark matter possesses bulk viscosity.}
\end{figure}

\section{Final Remarks}

Ref. \cite{lazo} developed the gravitational field equations for a class of theories based on action-dependent Lagrangians. We have studied in this work the flat FLRW cosmology of the resulting theory. Interestingly, we found a deep connection between such formalism and the unified bulk viscous cosmologies (in the Eckart formalism) which have been widely studied in the literature \cite{Fabris:2005ts, Colistete:2007xi,
	HipolitoRicaldi:2009je, HipolitoRicaldi:2010mf, Fabris:2011wk, Velten:2011bg}. This should be seem as the main result of this contribution.

We have also presented the perturbative dynamics of this model focusing on the scalar matter density fluctuations in the synchronous gauge. Our preliminary analysis shows that the geometrical viscous model does not present the same pathological behavior as the fluid viscous one concerning the growth of cosmic structures. This allows us to promote a proper comparison with matter clustering data in a future work. Such analysis will set the viability of cosmological scenarios based on the new class of action-dependent gravitational theories.

{\bf Acknowledgements:} The authors are grateful to CNPq (Brazil) and FAPES (Brazil) for financial support. We also thank the anonymous referees for their helpful comments and suggestions.


\section{introduction}

A classical Theorem of Dirichlet~\cite[Theorem 185]{HW} states that if $\alpha$ is irrational then there exist infinitely many pairs of integers $(m,n)$ satisfying the inequality
\begin{equation} \label{dis1}
\left | \alpha - \frac{m}{n} \right | < \frac{1}{n^{2}}.
\end{equation}
Such problems fall in the area of Diophantine approximation.
Equivalently by defining
 $\lVert x \rVert = \min_{v \in \mathbb{Z}} |x-v|$ to be the distance from $x$ to the nearest integer, (\ref{dis1}) implies
$$
\lVert \alpha n \rVert < n^{-1}
$$
for infinitely many positive integers $n$.
A natural extension to this problem is to consider if there exist infinitely many  solutions to
\begin{equation} \label{dio-ineq}
\lVert \alpha n + \beta \rVert < n^{-\kappa}.
\end{equation}
Here $\beta \in \mathbb{R}$, $\kappa >0$, $n \in {\mathcal N} \subseteq \mathbb{N}$ and ${\mathcal N}$ is some set of arithmetic interest.

Various results have been obtained when ${\mathcal N}$ is the set of prime numbers. These include the results of Harman~\cite{H1,H2}, Heath-Brown \& Jia~\cite{Heath} and Jia~\cite{Jia}, which involve sieve methods and bounds of exponential sums, while Vaughan~\cite{V} obtained his result by applying what is now known as the {\it Vaughan Identity}~\cite{V2}, together with bounds for exponential sums.

We remark that Harman's sieve method first appeared when  Harman studied this exact problem~\cite{H1} to get the exponent $\kappa = 3/10$. For more details on the Harman sieve, see the monograph~\cite{H}. The best result to date is by Matom\"aki~\cite{M} with $\kappa < 1/3$ under the condition $\beta =0$ and employs Harman's sieve method where the arithmetic information comes from bounds for averages of Kloosterman sums.

Let $k>1$ be a fixed positive integer and ${\mathcal N}$ the set of $k$th powers of primes. Baker \& Harman~\cite{BH} showed that we can take $\kappa < 3/20$ if $k=2$ and $\kappa < (3 \cdot 2^{k-1})^{-1}$ if $k \ge 3$. In particular, when $k=2$ this improves a result of Ghosh~\cite{Chosh}. Later, Wong~\cite{W} provided an improvement of Baker \& Harman~\cite{BH} in the range $3 \le k \le 12$.

When ${\mathcal N}$ is the set of square-free numbers, the best result is due to Heath-Brown~\cite{Heath1} with $\kappa <2/3$ using an essentially elementary method. The result improves the previous work by Harman~\cite{H3} and Balog \& Perelli~\cite{BP}, who independently showed we can essentially take the exponent $\kappa < 1/2$.

In this paper, we consider the problem of establishing an asymptotic formula for the number of $[y,z]$-smooth $n \le x$ solutions  to~(\ref{dio-ineq}) (numbers with prime factors in the interval $[y,z]$) with $\kappa =1/4 -\varepsilon$, where $\varepsilon >0$. We also consider a hybrid problem that interpolates between square-free and $[y,z]$-smooth integers.

We note that in the special case of smooth numbers we are able to obtain a non-trivial lower bound immediately. Indeed, for any fixed $\varepsilon >0$, consider the set
$$
\{n=ab \le x : p|ab \implies p < x^{\varepsilon}, x^{1/3} \le a < 2 x^{1/3} \}.
$$
By applying Lemma~\ref{expansion} and our Type II information (Lemma~\ref{TypeI,II-bound} with~(\ref{le 1/2})), we can show immediately that for infinitely many $x$ chosen correctly there are at least $x^{2/3+\varepsilon-o(1)}$ integers up to $x$ which are $x^{\varepsilon}$-smooth  and satisfies
$$
\lVert \alpha n + \beta \rVert < x^{-1/3 +\varepsilon}.
$$
The statement remains valid if we additionally require $n$ to be square-free.

\section{Notation}
For complex valued functions $f$ and $g$, we use the notation $f =O(g)$ and $f \ll g$ to mean there exists an absolute constant $C > 0$ such that $ |f(x)|\le C |g(x)|$ for all sufficiently large $x$.

Write $m \sim M$ when $ M \le m < 2M$, $e(x)= \exp(2 \pi i x)$ and let $\mu$ be the M\"obius function.
Given a positive integer $n \ge 2$, we denote by $P^-(n)$ and $P^+(n)$ the smallest and largest prime factor of $n$ respectively. We always assume $p$ with or without sub/superscript are primes.

\section{Main Results}
For any real numbers $0<y \le z \le x$, we denote by
$$
S(x;y,z)= \{a \in [1,x] \cap \mathbb{N}: p|a \implies p \in [y,z] \}
$$
 the set of all integers in $[1,x]$ all of whose prime factors are in $[y,z]$. Moreover, we denote the cardinality of $S(x;y,z)$ by $\Psi(x;y,z)$.

For any positive integer $A > 0$, we write
\begin{align*}
{\mathcal I}(A) = \{ \alpha \in \mathbb{R} \backslash \mathbb{Q} : \alpha =[a_0; a_1 , \ldots ], |a_j| \le A, j \ge 0 \}.
\end{align*}
We note that $\cup_{A \in \mathbb{N}} {\mathcal I}(A)$ contains the set of all quadratic irrationals. We can now state our first result.

\begin{theorem} \label{smooth}
There exists an increasing sequence  $( x_{k} )_{k \in \mathbb{N}}$ of positive integers such that if $2 \le y < x_k^{1/2 }$, $y < z \le x_k$, $\varepsilon > 0$ and $\delta = x_k^{-1/4 +\varepsilon}$, then
\begin{equation} \label{y-z smooth eqn}
\sum_{\substack{n \in S(x_k;y,z) \\ \lVert \alpha n + \beta \rVert < \delta \\ }} 1   =  2\delta \Psi(x_k;y,z)  + O(x_k^{ 3/4+ \varepsilon /2 +o(1) }).
\end{equation}
Moreover,~$\normalfont{(\ref{y-z smooth eqn})}$ \textit{holds for the sequence $( x_{k})_{k \in \mathbb{N}} = (k)_{k \in \mathbb{N}}$  uniformly for all $\alpha \in {\mathcal I}(A)$ and any fixed positive integer $A$}.
\end{theorem}

For any real numbers $0 <y \le z \le x$, we denote by
$$
S^*(x;y,z) =\{ a \in [1,x] \cap \mathbb{N} : \mu^2(a)=1, p|a \implies p \in [y,z] \}
$$
the set of all square-free integers in $[1,x]$ all of whose prime factors are in $[y,z]$. We also denote the cardinality of $ S^*(x;y,z)$ by $\Psi^*(x;y,z)$. The next theorem is essentially Theorem~\ref{smooth}, where we also require the integers we are counting to be square-free.

\begin{theorem} \label{y-z smooth}
There exists an increasing sequence $( x_{k} )_{k \in \mathbb{N}}$  of positive integers such that if $2 \le y < x_k^{1/2 }$, $y < z \le x_k$, $\varepsilon > 0$ and $\delta = x_k^{-1/4 +\varepsilon}$, then we have
\begin{equation} \label{s.f y-eqn}
\sum_{\substack{n \in S^*(x_k;y,z) \\ \lVert \alpha n + \beta \rVert < \delta \\ }} 1   =  2\delta \Psi^*(x_k;y,z)  + O(x_k^{ 3/4+ \varepsilon/2 +o(1) }).
\end{equation}
Moreover,~$\normalfont{(\ref{s.f y-eqn})}$ \textit{holds for the sequence $( x_{k})_{k \in \mathbb{N}} = (k)_{k \in \mathbb{N}}$  uniformly for all $\alpha \in {\mathcal I}(A)$ and any fixed positive integer $A$}.
\end{theorem}

If we fix $0 < u_2 < u_1$ with $2 < u_1 $, $u_2 < \lfloor u_1 \rfloor$ and set $y = x_k^{1/u_1 }$ and $z=x_k^{1/u_2 }$ then both  $\Psi(x_k;y,z)$, $\Psi^*(x_k;y,z)$ in the main term of~(\ref{y-z smooth eqn}) and~(\ref{s.f y-eqn}) respectively are bounded below by the number of integers $n=p_1 \ldots p_j$ that are products of $j = \lfloor u_1 \rfloor$ distinct primes  with $p_i \in [y,z]$. This gives the lower bound
$$
x_k^{1-o(1)} < \Psi(x_k;y,z),\Psi^*(x_k;y,z).
$$
It follows that our Theorem~\ref{smooth} and~\ref{y-z smooth} are non-trivial in this region.

We note that by~\cite[Theorem 1]{F} of Friedlander, we can obtain an asymptotic formula for $\Psi(x;y,z)$ in certain regions. We also mention that Saias~\cite{S1,S2,S3} has  extensively studied this quantity. In particular, our  $\Psi(x_k;y,z)$ is $\theta(x_k,z,y)$ or $\Theta(x_k,z,y)$ in the notation of Friedlander and Saias respectively.

We remark that we assume $y < x_k^{1/2}$ in Theorem~\ref{smooth} and~\ref{y-z smooth} since if $y \ge x_k^{1/2}$ and $z \ge x_k$ then both $S(x_k;y,z)$ and $S^*(x_k;y,z)$ essentially count primes in the interval $[y,x_k]$ and the result of Harman~\cite[Theorem 3.2]{H} covers this case with $\delta = x_k^{-1/4+\varepsilon}$.

It is easy to see that the proof of Theorem~\ref{smooth} can be adjusted to prove Theorem~\ref{y-z smooth}, so we will only give full details for the proof of Theorem~\ref{smooth}.

\section{Preparations}

For any $\delta >0$, we define
$$
\chi(r) =
\begin{cases}
1 &\mbox{ if } \lVert r \rVert < \delta, \\
0 &\mbox{ otherwise.}
\end{cases}
$$
We recall a Lemma from~\cite[Chapter 2]{B}, which provides a finite Fourier approximation to $\chi$. This converts the problem of detecting solutions to~(\ref{dio-ineq}) to a problem about estimates for exponential sums.
\begin{lemma} \label{expansion}
For any positive integer $L$, there exist complex sequences $( c_{\ell}^{-} )_{|\ell |\le L}$ and $ ( c_{\ell}^{+} )_{|\ell |\le L} $ such that
$$
2 \delta - \frac{1}{L+1} + \sum_{0 < |\ell | \le L} c_{\ell}^{-} e(\ell \theta) \le \chi ( \theta) \le 2 \delta + \frac{1}{L+1} + \sum_{0 < | \ell| \le L} c_{\ell}^{+} e(\ell \theta)
$$
where
\begin{equation*} \label{fourier-expansion-coeff}
| c_{\ell}^{\pm} |  \le \min \left \{ 2 \delta + \frac{1}{L+1}, \frac{3}{2 \ell} \right \}.
\end{equation*}
\end{lemma}

Let $\alpha ,x \in \mathbb{R}$ be such that $|\alpha - a/q| < 1/q^2$ where $1 \le a \le q \ll x^{O(1)}$ and $(a,q)=1$. Let  $(a_{m})_{m \in \mathbb{N}}$, $(b_{n})_{n \in \mathbb{N}}$, $(c_{\ell})_{\ell \in \mathbb{N} }$ be complex sequences satisfying
$$
|a_{m}| \le \tau(m), \quad |b_{n}| \le \tau(n), \quad
|c_{\ell}| \le \min \Big \{ 2 \delta + \frac{1}{\lfloor x \rfloor +1}, \frac{3}{2\ell} \Big \}
$$
whenever $m,n, \ell \in \mathbb{N}$. We will also use the classical divisor bound $\tau(r) \ll x^{o(1)}$ for $r \le x$. Lastly for $\varepsilon >0$, we denote $\delta = x^{-\frac{1}{4} + \varepsilon}$.

The next two Lemmas provide our Type I and II information respectively. They can be obtained by following the method of~\cite[Section 2.3]{H}.

\begin{lemma} \label{typeI exp-bound}
We have
\begin{align*}
 \sum_{\ell \le x}  c_{\ell}  & \sum_{\substack{mn \le x \\ m \sim  M }}  a_{m}  e(\alpha \ell mn) \ll  (M +  xq^{-1} + \delta q)x^{o(1)}  .
\end{align*}
\end{lemma}

\begin{lemma} \label{TypeI,II-bound}
We have
\begin{equation} \label{le 1/2}
\sum_{\ell \le x}  c_{\ell}  \sum_{\substack{mn \le x \\ m \sim M }} a_{m} b_{n} e(\alpha \ell mn) \ll
x^{1+o(1)}\left( \frac{\delta}{M}  + \frac{M}{x}  + \frac{1}{q} + \frac{q \delta }{x} \right)^{ \frac{1}{2} }
\end{equation}
for $M \ll x^{1/2 }$ and
\begin{equation} \label{ge 1/2}
\sum_{\ell \le x}  c_{\ell}  \sum_{\substack{mn \le x \\ m \sim M }} a_{m} b_{n} e(\alpha \ell mn) \ll x^{1+o(1)} \left(\frac{M \delta}{x} + \frac{1}{M} + \frac{1}{q} + \frac{q\delta }{x} \right)^{ \frac{1}{2} }
\end{equation}
for $M \gg x^{1/2 }$.
\end{lemma}

\subsection{Estimates for Type I \& II sum} \label{sec:type I, II estimate}

Denote
$$
{\mathcal B} = \{n \in \mathbb{N} : 2 \le n \le x \} \qquad \mbox{and} \qquad
  {\mathcal A} = \{ n \in {\mathcal B}: \lVert \alpha n + \beta \rVert < \delta \}.
$$

We will first state our Type I estimate.

\begin{lemma}[Type I estimate] \label{type one estimate}
For $q = x^{2/3}$ we have
$$
\sum_{\substack{mn \in {\mathcal A} \\ m \sim M}} a_m = 2 \delta \sum_{\substack{mn \in {\mathcal B} \\ m \sim M}} a_m + O(x^{3/4 +\varepsilon/2 })
$$
whenever $M \ll x^{ 3/4 }$.
\end{lemma}

\begin{proof}
By Lemma~\ref{expansion} with $L = \lfloor  x \rfloor$, we get
$$
\sum_{\substack{mn \in {\mathcal A} \\ m \sim M}} a_m = \sum_{\substack{mn \in {\mathcal B} \\ m \sim M}} a_m \chi(\alpha mn + \beta ) = 2\delta \sum_{\substack{mn \in {\mathcal B} \\ m \sim M}} a_m +O(E_1 + E_2)
$$
where
\begin{align*}
E_1 & = \frac{1}{L+1} \sum_{\substack{mn \le x \\ m \sim M}} a_m  \ll x^{o(1)}, \\
E_2 & = \sum_{0 < |\ell| \le L} |c_{\ell}^{\pm}| \cdot \Bigg |\sum_{\substack{mn \le x \\ m \sim m}} a_m e(\ell(\alpha mn + \beta)) \Bigg |.
\end{align*}
Clearly there exists $\xi_{\ell} \in \mathbb{C}$ with $|\xi_{\ell}|=1$ such that
\begin{align*}
E_{2}
& = \sum_{0 < |\ell| \le L} \xi_{\ell}c_{\ell}^{\pm} \sum_{\substack{mn \le x \\ m \sim M  }} a_{m} e( \alpha \ell  mn).
\end{align*}
 Applying Lemma~\ref{typeI exp-bound} immediately gives the bound
$$
E_2 \ll (M + xq^{-1} + \delta q)x^{o(1)}  \ll x^{ \frac{3}{4} + \frac{\varepsilon}{2} }
$$
whenever $M \ll x^{3/4 }$.
\end{proof}

Next, we state our Type II estimate.

\begin{lemma}[Type II estimate] \label{type two estimate}
For $q = x^{ 2/3 }$ we have
$$
\sum_{\substack{mn \in {\mathcal A} \\ x^{\gamma} \le m \le x^{\gamma+\tau} }} a_m b_n = 2\delta \sum_{\substack{mn \in {\mathcal B} \\ x^{\gamma} \le m \le x^{\gamma+\tau}}} a_m b_n + O(x^{ 3/4 + \varepsilon/2 + o(1) })
$$
uniformly for $1/4 \le \gamma$, $\gamma +\tau \le 3/4$.
\end{lemma}

\begin{proof}
We follow the method of Lemma~\ref{type one estimate}. Partition the summation over $m$ into dyadic intervals and apply Lemma~\ref{TypeI,II-bound}.
\end{proof}

\subsection{Sieve estimates}

For positive real numbers $w, u, v$, we denote
$$
P(w) = \prod_{p < w}p \qquad \mbox{and} \qquad P(u,v] = \prod_{u < p\le v} p.
$$
We also set
$$
Y=x^{ 3/4+ \varepsilon/2 + o(1)}.
$$

For any set ${\mathcal A} \subseteq [2,x]$ of integers and any positive integer $s$, we denote
$$
{\mathcal A}_s = \{ n : ns \in {\mathcal A}  \}.
$$
For real numbers $0<y\le z \le x$, we denote
$$
{\mathcal S}({\mathcal A};y,z) = \# \{  a \in {\mathcal A} : p|a \implies p \in [y,z]     \}
$$
and in the special case $z=x$, we denote
$$
{\mathcal S}({\mathcal A};y,z) = {\mathcal S}( {\mathcal A}; y).
$$
We state a variant of the Buchstab identity for $[y,z]$-smooth numbers which is based on taking out the largest prime factor of the integers we are counting.

\begin{lemma} \label{lem:Buchstab}
For any $2 \le  y \le z$ and any set ${\mathcal A} \subseteq [2,x] $ of integers, we have
$$
{\mathcal S}({\mathcal A};y,z) = \sum_{y \le p \le z} {\mathcal S}({\mathcal A}_p;y,p) + {\mathcal S}({\mathcal C}; x^{1/2}) +O(x^{ 1/2 })
$$
where ${\mathcal C} =\{a \in {\mathcal A}: y \le a \le z \}$.
\end{lemma}

\begin{proof}
Take $a \in \{ a \in {\mathcal A} : p|a \implies p \in [y,z] \}$. If $a$ is a prime then $a$ is counted in $\{a \in {\mathcal A}: y \le a \le z , a \mbox{ is prime} \}$. Otherwise $a$ has at least two prime factors and we can write $a = P^+(a)n$ where $n>1$ and the prime factors of $n$ lie in the interval $[y,P^+(a)]$; the result follows immediately.
\end{proof}

Next we state three lemmas which give sieve estimates for different regions. In particular, the proofs will rely on ingredients coming from the Harman sieve~\cite[Lemma 2]{H2}.

Our first sieve estimate is essentially based on an application of our Type II estimate (Lemma~\ref{type two estimate}).

\begin{lemma} \label{1/4 3/4}
For $x^{ 1/4 } \le y < z \le x^{ 3/4 }$, we have
\begin{equation*}
\sum_{y  \le p \le z} {\mathcal S} ({\mathcal A}_{p}; y,p)  = 2\delta \sum_{y  \le p \le z} {\mathcal S} ({\mathcal B}_{p}; y,p) +  O(Y).
\end{equation*}
\end{lemma}

\begin{proof}
We have
\begin{align*}
 \sum_{y  \le p \le z} {\mathcal S} ({\mathcal A}_{p}; y,p)
& = \sum_{y \le  p \le z } \sum_{\substack{np \in {\mathcal A} \\ (n,P(y)P(p,x])=1}} 1  = \sum_{ \substack{mn \in {\mathcal A} \\ y \le m \le z }}     a(m,n).
\end{align*}
Here
\begin{align*}
a(m,n) & =   \mathbbm{1}_{\mathbb{P}}(m) \sum_{(n,P(y)P(m,x])=1} 1  =    \mathbbm{1}_{\mathbb{P}}(m) \sum_{\substack{d|n \\ d|P(y)P(m,x]}} \mu(d)
\end{align*}
where $\mathbbm{1}_{\mathbb{P}}(\cdot)$ is the characteristic function for the primes. We may assume that $d$ is square-free, as otherwise it does not contribute to the sum $a(m,n)$.

Write
$$
a(m,n) =  \mathbbm{1}_{\mathbb{P}}(m) \sum_{\substack{d_1 d_2|n \\ d_2|P(y) \\ P^-(d_1) > m}} \mu(d_1 d_2).
$$
We may assume $d_1 >1$; the case $d_1=1$ follows immediately.
We now get rid of the condition $P^-(d_1) > m$ by appealing to the truncated Perron formula
$$
\frac{1}{\pi} \int_{-T}^{T} e^{i \gamma t} \frac{\sin \rho t}{t} dt =
\begin{cases}
1 + O(T^{-1}(\rho - |\gamma|)^{-1} ) &  \mbox{ if $|\gamma | \le \rho$,}\\
O(T^{-1}(|\gamma| - \rho)^{-1} ) & \mbox{ if $|\gamma| > \rho$}.
\end{cases}
$$
Applying the above formula with $\rho = \log (P^-(d_1) -\frac{1}{2})$, $\gamma = \log m$ and $T = x^2 \delta^{-1}$, we have
\begin{align*}
 \sum_{ \substack{mn \in {\mathcal A} \\ y \le m \le z }}    a(m,n) = M({\mathcal A}) + O(R).
\end{align*}
Here
\begin{align*}
M ({\mathcal A}) & =  \frac{1}{\pi} \int_{-T}^{T} \sum_{\substack{mn \in {\mathcal A} \\ y \le m \le z}} \mathbbm{1}_{\mathbb{P}}(m) e^{i\gamma t} \Bigg (\sum_{\substack{d_1 d_2|n \\ d_2|P(y) }} \mu(d_1 d_2) \sin  ( \rho t ) \Bigg ) \frac{dt}{t}\\
R & =T^{-1} \sum_{\substack{mn \in {\mathcal A} \\ y \le m \le z}} \mathbbm{1}_{\mathbb{P}}(m) \sum_{\substack{d_1 d_2|n \\ d_2|P(y) }}  \frac{\mu(d_1d_2)}{|\log (P^-(d)- \frac{1}{2}) - \log m|}.
\end{align*}
We will consider the remainder term first. By the mean value theorem,
$$
\frac{1}{\log(P^-(d)-\frac{1}{2}) - \log m} = \frac{\eta}{P^-(d)-\frac{1}{2}-m}
$$
where $\eta \in [m,P^-(d)-\frac{1}{2}]$ or $[P^-(d)-\frac{1}{2},m]$. In any case, we have $\eta \le \max \{m, n- \frac{1}{2} \}$, so we bound
$$
\frac{\eta}{P^-(d) -\frac{1}{2}-m} \ll \max \{ m,n  \} \le x.
$$ Therefore $R \ll T^{-1}x^{2+o(1)}\ll \delta x^{o(1)}$.

It remains to estimate the main term. Note that the integral in the main term between $-1/T$ and $1/T$ can be trivially bounded by $\ll T^{-1} x^{1+o(1)} \ll \delta x^{-1}$. Applying our Type II estimate (Lemma~\ref{type two estimate}) in the integral over $\mathfrak{R}(T) = (-T,-1/T) \cup (1/T, T)$, we get
\begin{align*}
& M({\mathcal A}) \\
 & =  \frac{2\delta}{\pi}  \int_{\mathfrak{R}(T)}  \sum_{\substack{mn \in {\mathcal B} \\ y \le m \le z}} \mathbbm{1}_{\mathbb{P}}(m) m^{it} \sum_{\substack{d_1 d_2|n \\ d_2|P(y) }} \mu(d_1 d_2) \sin \Big (t \log \Big (P^-(d_1)-\frac{1}{2} \Big ) \Big )  \frac{dt}{t}\\
&  \quad + O \Big (Y \int_{\mathfrak{R}(T)} \frac{dt}{t}  \Big )   + O(\delta x^{-1})\\
& = \frac{2\delta}{\pi} \int_{-T}^{T}  \sum_{\substack{mn \in {\mathcal B} \\ y \le m \le z}} \mathbbm{1}_{\mathbb{P}}(m) m^{it} \sum_{\substack{d_1 d_2|n \\ d_2|P(y)}} \mu(d_1 d_2) \sin \Big (t \log \Big (P^-(d_1)-\frac{1}{2} \Big ) \Big )  \frac{dt}{t}\\
&  \quad + O(Y) \\
&  = 2\delta M({\mathcal B})+ O(Y)
\end{align*}
where the last line follows from the truncated Perron formula once again. The result follows immediately.
\end{proof}

The next sieve estimate is based on an idea which dates back to Vinogradov who applied it to estimate sum over primes. The idea is to systematically take out the largest prime factor until we get sums which we can estimate by our Type II estimate (Lemma~\ref{type two estimate}).

\begin{lemma}   \label{1/4}
For $2 \le y <z < x^{1/4 }$, we have
\begin{equation*}
\sum_{ y \le p \le z} {\mathcal S} ({\mathcal A}_{p}; y,p) =2\delta \sum_{ y \le p \le z} {\mathcal S} ({\mathcal B}_{p}; y,p) + O(Y).
\end{equation*}
\end{lemma}

\begin{proof}
By the method of Lemma~\ref{lem:Buchstab} we have
\begin{align*}
\sum_{y \le p \le z} {\mathcal S}({\mathcal A}_p; y,p) & = \sum_{y \le p_1 \le p \le z} {\mathcal S}({\mathcal A}_{pq_1}; y, p_1) +O( x^{1/2} ).
\end{align*}
We will first consider the sum on the right with $p p_1 >x^{1/4}$. Clearly we have
$$
\sum_{\substack{y \le p_1 \le p \le z \\ pp_1 >x^{1/4} }} {\mathcal S}({\mathcal A}_{pp_1}; y, p_1) = \sum_{\substack{y \le p_1 \le p \le z \\ pp_1 >x^{1/4 }}} \sum_{\substack{npp_1 \in {\mathcal A} \\ (n, P(y)P(p_1,x] )=1 }} 1.
$$
Set $m = pp_1$ and note that $m < x^{1/2}$. It follows that
$$
\sum_{\substack{y \le p_1 \le p \le z \\ pp_1 >x^{1/4} }} {\mathcal S}({\mathcal A}_{pp_1}; y, p_1) = \sum_{\substack{mn \in {\mathcal A} \\ x^{1/4} < m < x^{1/2} }} a_m \sum_{\substack{d_1d_2|n \\ d_1|P(y) \\ P^-(d_2) > P^-(m)}} \mu(d_1 d_2)
$$
where $a_m$ is 1 if $m=pp_1$ with $p,p_1 \in [y,z]$ and 0 otherwise. Applying the method of Lemma~\ref{1/4 3/4}, we get
$$
\sum_{\substack{y \le p_1 \le p \le z \\ pp_1 >x^{1/4} }} {\mathcal S}({\mathcal A}_{pp_1}; y, p_1) = 2\delta \sum_{\substack{y \le p_1 \le p \le z \\ pp_1 >x^{1/4} }} {\mathcal S}({\mathcal B}_{pp_1}; y, p_1) +O( Y ).
$$
The only part left to consider is the sum
$$
\sum_{ \substack{ y \le p_1 \le p \le z \\ pp_1 \le x^{1/4} }} {\mathcal S}({\mathcal A}_{pp_1}; y, p_1).
$$
If it is zero then we are done, otherwise we take out the next largest prime factor to get
$$
\sum_{ \substack{ y \le p_1 \le p \le z \\ pp_1 \le x^{1/4} }} {\mathcal S}({\mathcal A}_{pp_1}; y, p_1) = \sum_{\substack{ y \le p_2 \le p_1 \le p \le z \\ pp_1 \le x^{1/4} }}  {\mathcal S}({\mathcal A}_{pp_1 p_2}; y , p_2) +O(x^{1/2 }).
$$
The sum on the right with $pp_1 p_2 >x^{1/4 }$ can be dealt with again by the method of Lemma~\ref{1/4 3/4}. By induction this can go on for at most $O(\log x)$ steps. Since we have an asymptotic formula for every sum, the result follows.
\end{proof}

For our next sieve estimate, we note that our Type II estimates are not sufficient in this region. We bypass this complication by a role reversal that minimises the length of summation in exchange for sifting primes.

\begin{lemma}   \label{3/4}
For $x^{{3/4}} \le z \le x$, we have
\begin{equation*}
\sum_{ x^{3/4} < p \le z} {\mathcal S} ({\mathcal A}_{p}; y,p) = 2\delta \sum_{ x^{3/4} < p \le z} {\mathcal S} ({\mathcal B}_{p}; y,p)   +O(Y).
\end{equation*}
\end{lemma}

\begin{proof}
Take $x^{3/4} < p \le z$ and an element $np \in {\mathcal A}$ such that if $q|n$ where $q$ is prime then $q \in [y,p]$. This gives
\begin{equation} \label{switching roles}
\sum_{ x^{3/4} < p \le z} {\mathcal S} ({\mathcal A}_{p}; y,p) = \sum_{n < x^{1/4} } c_{n} {\mathcal S}({\mathcal A}_n' ; (x/n)^{1/2 }) + O(x^{o(1)}).
\end{equation}
Here
\begin{align*}
{\mathcal A}'_n & = \{ \max \{ y+\Delta(y),x^{3/4} \} < m \le z : mn \in {\mathcal A} \} \\
& = \{ \max \{ y+\Delta(y),x^{3/4} \} < m \le \min \{ z, x/n\} : \lVert \alpha mn + \beta \rVert < \delta \}
\end{align*}
where $c_n$ is 1 if all prime factors of $n$ are at least $y$ and 0 otherwise and $\Delta(y)$ is -1/2 if $y$ is an integer and zero otherwise. We also recall
$$
{\mathcal S}({\mathcal W}; k)= \# \{ w \in {\mathcal W} : (w,P(k))=1 \}.
$$

We note that the sum in~(\ref{switching roles}) may be empty depending on the choice of $y$. To show~(\ref{switching roles}), take $m \in S({\mathcal A}_n' ; (x/n)^{\frac{1}{2}})$ and suppose that $m$ has two prime factors say $q_1, q_2$ then we must have
$$
m \ge q_1 q_2 \ge (x/n)^{1/2}(x/n)^{1/2} = x/n.
$$
This is a contradiction unless $mn =x$ and this occurs  at most $O(x^{o(1)})$ times. Moreover, it is evident that
\begin{align*}
\sum_{n < x^{1/4} } c_n {\mathcal S}({\mathcal A}_n' ; (x/n)^{1/2 }) &= \sum_{n < x^{1/4}} c_n ({\mathcal S}({\mathcal A}'_n ; x^{1/2 }) +O(x^{1/2 })) \\
& = \sum_{n < x^{1/4}} c_n {\mathcal S}({\mathcal A}'_n ; x^{1/2}) +O(x^{3/4 }).
\end{align*}

Let $n \sim N$ where $N \ll x^{1/4}$ and $\max \{ y + \Delta(y),x^{3/4}\} n \ll M(n)=M \ll \min \{zn,x\}$ and consider the set
$$
{\mathcal A}^{(M)}= \{ m \sim M : \lVert \alpha m +\beta \rVert < \delta  \}.
$$
Then we have
$$
{\mathcal A}_n^{(M)} = \{ m \sim M/n : \lVert \alpha mn + \beta \rVert < \delta \}
$$
and by the Harman sieve~\cite[Lemma 2]{H2} using the  Type I and II estimate (Lemma~\ref{type one estimate} and~\ref{type two estimate}), we get
$$
\sum_{n \sim N} c_n  {\mathcal S}({\mathcal A}_n^{(M)}; x^{1/2}) = 2 \delta  \sum_{n \sim N} c_n {\mathcal S}({\mathcal B}_n^{(M)}; x^{1/2 }) + O(Y).
$$
Therefore by summing over $N,M$ we obtain
$$
\sum_{n < x^{1/4}} c_n {\mathcal S}({\mathcal A}'_n; x^{1/2}) = 2\delta  \sum_{n < x^{1/4}} c_n  {\mathcal S}({\mathcal B}'_n; x^{1/2}) + O(Y),
$$
which is what we needed to show.
\end{proof}

For a set ${\mathcal C} \subseteq [2,x]$ of integers and for integers $2 \le y < z \le x$, we denote
$$
{\mathcal T}({\mathcal C};y,z) = \# \{ c \in {\mathcal C} : \mu^2(c)=1, p|c \implies p\in [y,z] \}.
$$

The next three results are square-free analogue of Lemma~\ref{1/4 3/4}-\ref{3/4}.

\begin{lemma} \label{1/4 3/4 sf}
For $x^{ 1/4 } \le y < z \le x^{3/4 }$, we have
\begin{equation*}
\sum_{y  \le p \le z} {\mathcal T} ({\mathcal A}_{p}; y,p-1)  = 2\delta \sum_{y  < p \le z} {\mathcal T} ({\mathcal B}_{p}; y,p-1) +  O(Y).
\end{equation*}
\end{lemma}

\begin{proof}
By writing
\begin{align*}
\sum_{y \le p \le z} {\mathcal T}({\mathcal A}_p;y,p-1)   &=\sum_{y \le p \le z} \sum_{ \substack{np \in {\mathcal A} \\ (n,P(y),P(p-1,x])=1} } \mu^2(np) \\
& = \sum_{\substack{mn \in {\mathcal A} \\ y \le m \le z}} \mathbbm{1}_{\mathbb{P}}(m) \mu^2(n) \sum_{\substack{d_1d_2|n \\d_1|P(y) \\ P^-(d_2) \ge m}} \mu(d_1 d_2).
\end{align*}
 We see that the method of Lemma~\ref{1/4 3/4} gives the result.
\end{proof}

\begin{lemma} \label{1/4 sf}
For $2 \le y <z < x^{1/4 }$, we have
$$
\sum_{y \le p \le z} {\mathcal T}({\mathcal A}_p;y,p-1) = 2\delta \sum_{y \le p \le z} {\mathcal T}({\mathcal B}_{p};y,p-1) +O(Y).
$$
\end{lemma}

\begin{proof}
We follow as in Lemma~\ref{1/4} but we successively take out the largest distinct prime factors.
\end{proof}

\begin{lemma}   \label{3/4 sf}
For $x^{3/4} \le z \le x$, we have
\begin{equation*}
\sum_{ x^{3/4} < p \le z} {\mathcal T} ({\mathcal A}_{p}; y,p-1) = 2\delta \sum_{ x^{3/4} < p \le z} {\mathcal T} ({\mathcal B}_{p}; y,p-1)   +O(Y).
\end{equation*}
\end{lemma}

\begin{proof}
This follows immediately from the proof of Lemma~\ref{3/4} by taking $c_n$ there to be 1 if $n$ is square-free and all prime factors of $n$ are at least $y$ and 0 otherwise.
\end{proof}

\section{Proof of Theorem~\ref{smooth}}

By Dirichlet's theorem, we have $|\alpha - a/q|<1/q^2$ for infinitely many $q$, let $(q_k)_{k \in \mathbb{N}}$ be an increasing sequence of such denominators. Then we can choose an increasing sequence $(x_k)_{k \in \mathbb{N}}$ so that $q_k=x_k^{2/3 }$. Now we take $x_k$ to be a sufficiently large element in the sequence $(x_k)_{k \in \mathbb{N}}$.

By Lemma~\ref{lem:Buchstab}, we assert
$$
{\mathcal S}({\mathcal A};y,z) = \sum_{y \le p \le z} {\mathcal S}({\mathcal A}_p;y,p) + {\mathcal S}({\mathcal C}; x_k^{ 1/2 }) +O(x_k^{ 1/2 }).
$$
By~\cite[Theorem 2]{V} and partial summation, we have
$$
{\mathcal S}({\mathcal C}; x_k^{1/2}) = 2\delta {\mathcal S}({\mathcal D};  x_k^{ 1/2 }) +O(Y)
$$
where ${\mathcal D}= \{ n \in \mathbb{N} : y\le n \le z\}$ and $Y =x_k^{3/4 + \varepsilon/2 + o(1) }$. Therefore by substituting this into the above, we obtain
$$
{\mathcal S}({\mathcal A};y,z) = \sum_{y \le p \le z} {\mathcal S}({\mathcal A}_p;y,p) + 2\delta {\mathcal S}({\mathcal D}; x_k^{ 1/2 }) +O(Y).
$$
Suppose $y < x_k^{1/4 }$. Then we write
\begin{align*}
{\mathcal S}({\mathcal A};y,z) &= \Bigg (\sum_{y \le p \le z <x_k^{1/4}} + \sum_{y \le p < x_k^{1/4}} + \sum_{x_k^{1/4}\le p \le z \le x_k^{3/4}} +  \sum_{x_k^{1/4}\le p \le x_k^{3/4} \le z } \\
&\quad + \sum_{x_k^{3/4} < p \le z} \Bigg ) {\mathcal S}({\mathcal A}_p;y,p)  + 2\delta {\mathcal S}({\mathcal D}; x_k^{\frac{1}{2} }) +O(Y).
\end{align*}
We note that some of the sums above may be empty depending on the choice of $z$. Applying Lemma~\ref{1/4 3/4}-\ref{3/4} then Lemma~\ref{lem:Buchstab} to the above we get
$$
{\mathcal S}({\mathcal A};y,z) = 2\delta {\mathcal S}({\mathcal B} ;y,z) +O(Y).
$$
The case $x_k^{1/4 } \le y < x_k^{1/2 }$ is similar to the case above and the first result follows.

Now suppose $\alpha \in {\mathcal I}(A)$ then we can assume $\alpha \in (0,1) \backslash \mathbb{Q}$. Indeed, we can write $\alpha = b + r$ where $b \in \mathbb{Z}$ and $r \in (0,1) \backslash \mathbb{Q}$ so that $\lVert \alpha n + \beta \rVert = \lVert r n + \beta \rVert$. Therefore we have the continued fraction expansion
$$
\alpha = [0; a_{1}, a_2, \ldots]
$$
with $0 < a_{i} \le A$ for some absolute constant $A$. The convergents $p_{s}/q_{s}$ to $\alpha$ satisfy the inequality
$|\alpha -p_{s}/q_{s}| < 1/q_s^2$.
Clearly $( q_{s} )_{s \ge 1}$ is a strictly increasing sequence. Therefore, for any sufficiently large $x$ there exist $s \in \mathbb{N}$ such that $q_{s-1} \le x^{2/3} \le q_{s}$. Since $q_{s}=a_{s}q_{s-1}+q_{s-2}$ for $s \ge 2$, we get $q_{s} \le  (A+1)q_{s-1} \le  2A x^{2/3} \ll x^{2/3}$, because $A$ is fixed. Hence  $x^{2/3} \ll q_{s} \ll x^{2/3}$ and the result follows by the argument above.

\section{Proof of Theorem~\ref{y-z smooth}}

We proceed as in the proof of Theorem~\ref{smooth} but instead we apply Lemma~\ref{1/4 3/4 sf}-\ref{3/4 sf}.

\section*{Acknowledgement}
The author is grateful to the anonymous referee for their excellent comments which resulted in a complete rewrite of this paper, and to B. Kerr and I. E. Shparlinski for helpful comments on a previous draft. Lastly the author thanks T. Evans and N. Poznanovi\'c for stimulating conversations on a related topic at the beach.

 The author was supported by an Australian Government Research Training Program (RTP) Scholarship.

\normalsize