Marcinkiewicz-type multiplier theorem for q-variation (q > 1)

Not long ago we discussed one of the main direct applications of the Littlewood-Paley theory, namely the Marcinkiewicz multiplier theorem. Recall that the single-variable version of this theorem can be formulated as follows:

Theorem 1 [Marcinkiewicz multiplier theorem]: Let {m} be a function on \mathbb{R} such that

  1. m \in L^\infty
  2. for every Littlewood-Paley dyadic interval L := [2^k, 2^{k+1}] \cup [-2^{k+1},-2^k] with k \in \mathbb{Z}

    \displaystyle \|m\|_{V(L)} \leq C,

    where \|m\|_{V(L)} denotes the total variation of {m} over the interval L .

Then for any {1 < p < \infty} the multiplier {T_m} defined by \widehat{T_m f} = m \widehat{f} for functions f \in L^2(\mathbb{R}) extends to an L^p \to L^p bounded operator,

\displaystyle \|Tf\|_{L^p} \lesssim_p (\|m\|_{L^\infty} + C) \|f\|_{L^p}.

You should also recall that the total variation V(I) above is defined as

\displaystyle \sup_{N}\sup_{\substack{t_0, \ldots, t_N \in I : \\ t_0 < \ldots < t_N}} \sum_{j=1}^{N} |m(t_j) - m(t_{j-1})|,

and if {m} is absolutely continuous then {m'} exists as a measurable function and the total variation over interval I is given equivalently by \int_{I} |m'(\xi)|d\xi . We have seen that the “dyadic total variation condition” 2.) above is to be seen as a generalisation of the pointwise condition |m'(\xi)|\lesssim |\xi|^{-1} , which in dimension 1 happens to coincide with the classical differential Hörmander condition (in higher dimensions the pointwise Marcinkiewicz conditions are of product type, while the pointwise Hörmander(-Mihklin) conditions are of radial type; see the relevant post). Thus the Marcinkiewicz multiplier theorem in dimension 1 can deal with multipliers whose symbol is somewhat rougher than being differentiable. It is an interesting question to wonder how much rougher the symbols can get while still preserving their L^p mapping properties (or maybe giving up some range – recall though that the range of boundedness for multipliers must be symmetric around 2 because multipliers are self-adjoint).

Coifman, Rubio de Francia and Semmes came up with an answer to this question that is very interesting. They generalise the Marcinkiewicz multiplier theorem (in dimension 1) to multipliers that have bounded {q} -variation with {q} > 1. Let us define this quantity rigorously.

Definition: Let q \geq 1 and let I be an interval. Given a function f : \mathbb{R} \to \mathbb{R}, its {q} -variation over the interval {I} is

\displaystyle \|f\|_{V_q(I)} := \sup_{N} \sup_{\substack{t_0, \ldots t_N \in I : \\ t_0 < \ldots < t_N}} \Big(\sum_{j=1}^{N} |f(t_j) - f(t_{j-1})|^q\Big)^{1/q}

Notice that, with respect to the notation above, we have \|m\|_{V(I)} = \|m\|_{V_1(I)} . From the fact that \|\cdot\|_{\ell^q} \leq \|\cdot \|_{\ell^p} when p \leq q we see that we have always \|f\|_{V_q (I)} \leq \|f\|_{V_p(I)} , and therefore the higher the {q} the less stringent the condition of having bounded {q} -variation becomes (this is linked to the Hölder regularity of the function getting worse). In particular, if we wanted to weaken hypothesis 2.) in the Marcinkiewicz multiplier theorem above, we could simply replace it with the condition that for any Littlewood-Paley dyadic interval L we have instead \|m\|_{V_q(L)} \leq C . This is indeed what Coifman, Rubio de Francia and Semmes do, and they were able to show the following:

Theorem 2 [Coifman-Rubio de Francia-Semmes, ’88]: Let q\geq 1 and let {m} be a function on \mathbb{R} such that

  1. m \in L^\infty
  2. for every Littlewood-Paley dyadic interval L := [2^k, 2^{k+1}] \cup [-2^{k+1},-2^k] with k \in \mathbb{Z}

    \displaystyle \|m\|_{V_q(L)} \leq C.

Then for any {1 < p < \infty} such that {\Big|\frac{1}{2} - \frac{1}{p}\Big| < \frac{1}{q} } the multiplier {T_m} defined by \widehat{T_m f} = m \widehat{f} extends to an L^p \to L^p bounded operator,

\displaystyle \|Tf\|_{L^p} \lesssim_p (\|m\|_{L^\infty} + C) \|f\|_{L^p}.

The statement is essentially the same as before, except that now we are imposing control of the {q} -variation instead and as a consequence we have the restriction that our Lebesgue exponent {p} satisfy {\Big|\frac{1}{2} - \frac{1}{p}\Big| < \frac{1}{q} }. Taking a closer look at this condition, we see that when the variation parameter is 1 \leq q \leq 2 the condition is empty, that is there is no restriction on the range of boundedness of T_m : it is still the full range {1} < {p} < \infty , and as {q} grows larger and larger the range of boundedness restricts itself to be smaller and smaller around the exponent p=2 (for which the multiplier is always necessarily bounded, by Plancherel). This is a very interesting behaviour, which points to the fact that there is a certain dichotomy between variation in the range below 2 and the range above 2, with 2 -variation being the critical case. This is not an isolated case: for example, the Variation Norm Carleson theorem is false for {q} -variation with {q \leq 2} ; similarly, the Lépingle inequality is false for 2-variation and below (and this is related to the properties of Brownian motion).

Today, as a natural continuation to the posts on Littlewood-Paley theory and its applications, I am going to present the interesting proof of this very nice theorem. I found it impossible to get a hold of the original paper without taking a trip to the library, so being lazy I am going to follow instead Lacey’s nice presentation (which I suppose is very close to the original one anyway). The reasons why I believe the proof to be interesting are several: for starters, the proof is surprisingly simple; moreover, it relies on a generalisation of the Littlewood-Paley theorem due to Rubio de Francia which is worth seeing once you have studied the basic theory of square functions; finally, one very important ingredient in the proof is a partition of \mathbb{R} dependent on the symbol {m} which is very much alike that beautiful idea which is used to prove the maximal Hausdorff-Young inequality of Christ-Kiselev.
To summarise the proof before we start: we will build sub-partitions of the partition of \mathbb{R} given by the Littlewood-Paley dyadic intervals, in such a manner that on each subpartition we have a uniform control of the L^\infty norm of (part of) the multiplier symbol {m} ; then we will use square functions adapted to these partitions to control the resulting multipliers, somewhat analogously to what is done in the Marcinkiewicz case (but conceptually even simpler). Finally we will combine everything together to conclude. In the next section we will present the necessary inequalities of Rubio de Francia (without proof) before presenting the proof.

1. Rubio de Francia square functions

Recall that the Littlewood-Paley theorem says the following: if \Delta_I  denotes the frequency projection given by \widehat{\Delta_I f} = \mathbf{1}_I \widehat{f} and L = [2^k, 2^{k+1}] \cup [-2^{k+1},-2^k] are the Littlewood-Paley dyadic intervals, the collection of which we denote by \mathbb{L} , then we have for any {1 < p < \infty}

\displaystyle \|f\|_{L^p (\mathbb{R})} \sim_p \Big\|\Big(\sum_{L \in \mathbb{L}} |\Delta_{L} f|^2\Big)^{1/2}\Big\|_{L^p (\mathbb{R})}.

The first person to go beyond this statement was Lennart Carleson, who investigated the square function

\displaystyle f \mapsto \Big(\sum_{n \in \mathbb{Z}} |\Delta_{[n,n+1]} f|^2\Big)^{1/2}.

As you can see, the intervals on which we are taking the frequency projections are no longer dyadic – rather, they all have the same (unit) length. If you recall, the heuristic motivation behind the Littlewood-Paley theory is that, since f = \sum_{L \in \mathbb{L}} \Delta_{L} f and the different frequency pieces are dyadically separated, we should expect a random cancellation between the different terms of the frequency decomposition – which means that most of the time the sum should have magnitude approximately \big(\sum_{L \in \mathbb{L}} |\Delta_{L} f|^2\big)^{1/2}. This heuristic no longer works in the case where the frequency intervals are the [n,n+1] ‘s, because there are now many many terms in the frequency decomposition that have comparable frequencies and will therefore be “aligned” for longer periods of time, or “correlated” if you prefer – there will not be much cancellation between them. This leads us to suspect that the analogue of the Littlewood-Paley theorem for Carleson’s square function should fail, at least partially. This is indeed the case. In fact, Carleson showed that the inequality \Big\|\Big(\sum_{n \in \mathbb{Z}} |\Delta_{[n,n+1]} f|^2\Big)^{1/2}\Big\|_{L^p (\mathbb{R})} \lesssim_p \|f\|_{L^p(\mathbb{R})} is FALSE when p < 2. This is actually very simple to check: just take {f} such that \widehat{f} = \mathbf{1}_{[0,N]} for {N} a large integer ({f} is a continuous version of the Dirichlet kernel). It is easy to estimate that \|f\|_{L^p} = \|\check{\mathbf{1}}_{[0,N]}\|_{L^p} \sim N^{1/{p'}} (the contribution from the peak around the origin dominates). As for the square function, we simply have

\displaystyle \Big(\sum_{n \in \mathbb{Z}} |\Delta_{[n,n+1]} \check{\mathbf{1}}_{[0,N]}|^2\Big)^{1/2} = |\check{\mathbf{1}}_{[0,1]}| N^{1/2},

and therefore the inequality can only be true if N^{1/2} \lesssim N^{1/{p'}} for any large {N} , which in turn is only possible if p \geq 2 . Carleson then went on to prove that in this range the inequality is actually true. His paper, titled “On the Littlewood-Paley theorem” and published in 1967 on the Report of the Mittag-Leffler Institute seems to have vanished – not even MathSciNet has it in their records. Indeed, I think the result went largely unnoticed, since more than a decade later (1981) Córdoba re-proved the theorem in his work on Bochner-Riesz multipliers.

Nevertheless, this work was later subsumed by work of Rubio de Francia, who generalised it to arbitrary collections of disjoint intervals. Indeed, he proved the following:

Theorem 3 [Rubio de Francia, ’85]: Let \mathcal{I} be a collection of disjoint intervals and let S_{\mathcal{I}} denote the associated square function

\displaystyle S_{\mathcal{I}}f := \Big(\sum_{I \in \mathcal{I}} |\Delta_I f|^2 \Big)^{1/2}.

For any {2 \leq p < \infty} we have

\displaystyle \|S_{\mathcal{I}}f \|_{L^p} \lesssim_p \|f\|_{L^p}. \ \ \ \ \ (1)

Importantly, the constant is independent of the collection \mathcal{I} .

Some remarks are in order:

  • The sharp reader will have already noticed that, since the intervals are disjoint, the p = 2 case of Theorem 3 is actually a trivial consequence of Plancherel.
  • As seen in the proof of the Littlewood-Paley theorem, if the collection \mathcal{I} partitions \mathbb{R} then we have by duality and Cauchy-Schwarz that inequality (1) for exponent {p \geq 2} has the consequence that

    \displaystyle \|f\|_{L^{p'}} \lesssim \|S_{\mathcal{I}}f \|_{L^{p'}},

    where now {1 < p' \leq 2} . Thus we can see Rubio de Francia's result as saying that the heuristic above still applies to arbitrary partitions of \mathbb{R} , but only in the range p \in (1,2] – the heuristic being that most of the time {f} \text{ ``}\lesssim  \text{'' }S_{\mathcal{I}}f .

  • Of course, depending on the collection of intervals, the condition p\geq 2 might not be sharp, as is the case when we take the Littlewood-Paley dyadic intervals and can thus go below 2. It is not currently known (despite the result being nearly 40 years old) when the condition is sharp, that is, there is no characterisation of the collections for which the theorem fails in the range p < 2. It is conjectured that the condition is sharp essentially for all collections that are not "lacunary" in some sense, though it is even hard to understand what the correct notion of lacunarity to state a conjecture should be.

There are a number of proofs of Rubio de Francia’s theorem, though we will not see a single one in here. Rubio de Francia’s original proof worked roughly as follows: first of all, using classical Littlewood-Paley theory and Whitney decompositions of the intervals, one reduces to the case where the intervals are all well-separated, meaning that if I,J \in \mathcal{I} are disjoint then 5I,5J are also disjoint; then, the theorem is reduced to proving the same statement for a square function G with the same frequency information but smoother frequency projections (much like in the proof of the Littlewood-Paley theorem); finally, the boundedness of G is deduced by interpolation between the trivial p=2 case and the endpoint inequality \|Gf\|_{\mathrm{BMO}}\lesssim \|f\|_{L^\infty} – which is not hard to prove as a consequence of some simple vector-valued kernel estimates.
Bourgain reproved the theorem by proving the endpoint of the dual inequality, thus extending the result somewhat; that is, he proved that when \mathcal{I} is a partition of the real line then

\displaystyle \|f\|_{H^1} \lesssim \|S_{\mathcal{I}}f\|_{L^1},

where \|\cdot\|_{H^1} denotes the quasi-norm of the (real) Hardy space H^1(\mathbb{R}) . Bourgain’s paper “On square functions on the trigonometric system” is very hard to get, having been published in a journal that does not exist in that form anymore. The proof he gives is beautiful (not a surprise) and I will probably talk about it in the future.
Another vividly distinct proof of Theorem 3 is given by Lacey in the aforementioned notes, in which the theorem is reformulated in time-frequency language (using wavepackets) and a time-frequency proof is given, with plenty of details. These notes include a discussion of many issues related to the Rubio de Francia inequalities and are a must-read for whoever is interested in such inequalities.
Finally, there is also a not-yet-published paper of Benea and Muscalu in which they re-prove Rubio de Francia’s theorem in yet another time-frequency way (distinct from Lacey’s).
Other inputs, variations and extensions to higher dimensions have been given by Soria, Sjölin, Journé, Sato, Zhu and maybe others I am forgetting at the moment; I will not discuss these here.

2. Proof of Theorem 2
With the Rubio de Francia inequalities at hand, we are ready to prove Theorem 2.

Proof: The idea is to reduce to the simpler case of multipliers which are simple functions on each Littlewood-Paley dyadic interval. In particular, assume that {m} is of the form \sum_{L \in \mathbb{L}} \sum_{I \in \mathcal{I}_L} m_I \mathbf{1}_I, where each collection \mathcal{I}_L is a partition of the Littlewood-Paley interval {L} and m_I is a complex coefficient. Assume furthermore that these subpartitions of Littlewood-Paley intervals are bounded in cardinality, that is assume that for every L \in \mathbb{L} we have \# \mathcal{I}_L \leq N . Then we can argue as follows, following the footsteps of the proof of the Marcinkiewicz multiplier theorem (Theorem 1 above): if we let {S} denote the Littlewood-Paley square function, we have by Littlewood-Paley theorem that

\displaystyle \|T_m f\|_{L^p} \sim_p \|ST_m f\|_{L^p},


\displaystyle ST_m f = \Big( \sum_{L \in \mathbb{L}} | \Delta_L T_m f|^2 \Big)^{1/2}.

Due to the simple form the symbol {m} has, we see that

\displaystyle |\Delta_L T_m f| = \Big|\sum_{I \in \mathcal{I}_L} m_I \Delta_I f \Big| \leq \|m\|_{L^\infty} \sum_{I \in \mathcal{I}_L} |\Delta_I f|;

having an \ell^1 sum is inconvenient in this context, and therefore we apply Cauchy-Schwarz to the latter to get a nice square function instead,

\displaystyle \begin{aligned} |\Delta_L T_m f| \leq & \|m\|_{L^\infty} (\# \mathcal{I}_L)^{1/2} \Big(\sum_{I \in \mathcal{I}_L} |\Delta_I f|^2 \Big)^{1/2} \\ \leq & \|m\|_{L^\infty} N^{1/2} S_{\mathcal{I}_L}f. \end{aligned}

Performing the \ell^2 -summation in L we see that we have shown

\displaystyle |S T_m f| \leq \|m\|_{L^\infty} N^{1/2} S_{\mathcal{I}}f,

where \mathcal{I} is the collection of all intervals, that is \mathcal{I} = \bigcup_{L} \mathcal{I}_L . Now the object on the RHS is a Rubio de Francia square function, which we know is bounded at least when p \geq 2 . Assume therefore that this is the case (which is not a limitation, because multipliers have range of boundedness symmetric about exponent 2), and as a consequence of (1) we have therefore that

\displaystyle \|T_m f\|_{L^p} \lesssim \|m\|_{L^\infty} N^{1/2} \|f\|_{L^p}

for all p\geq 2 . This bound is not so great because there is a large loss in the parameter {N} , but at this point we should observe something: when p=2 this factor is not there! Indeed, in that case we have simply \|T_m f\|_{L^2} \leq \|m\|_{L^\infty} \|f\|_{L^2} by Plancherel; but this means that for any p\geq 2 we can (complex) interpolate between all these estimates and lower the exponent 1/2 somewhat. Indeed, for a fixed exponent p\geq 2, we can write {p} as an interpolation exponent between 2 and any extremely large exponent {r} ; in practice, the result will be the same as if we had assumed r=\infty (once we take a limit), although obviously we are not allowed to use precisely this exponent. The result is the following: if \theta \in (0,1) is such that

\displaystyle \frac{1}{p} = \frac{1 - \theta}{2} + \frac{\theta}{\infty} = \frac{1 - \theta}{2}

then we have by interpolation that \|T_m f\|_{L^p} \lesssim (N^{1/2})^{\theta} \|f\|_{L^p} , and a simple computation shows that (N^{1/2})^{\theta} = N^{1/2 - 1/p} (recall that p\geq 2 in this argument), that is we have shown

\displaystyle \|T_m f\|_{L^p} \lesssim \|m\|_{L^\infty} N^{|1/2 - 1/p|} \|f\|_{L^p} \ \ \ \ \ (2)

for all p \in (1,\infty) . We have improved the constant a little! This small improvement will go a long way though.

Now that we have some partial result, how can we exploit it? Can we reduce the multiplier symbol in Theorem 2 to a symbol of the type just considered above? It turns out that it is not at all hard to reduce the symbol {m} to a sum of symbols of the above form, in such a way that (2) will give a summable contribution. The argument is an ingenious decomposition of {m} in martingale differences where the martingale is dictated from {m} itself (by its {q} -variation, precisely). Let us see how.
We assume for simplicity that the constant C is 1, that is we normalise the symbol so that for any Littlewood-Paley interval {L} we have \|m\|_{V_q(L)} \leq 1 . Fix such an L \in \mathbb{L} and let j \in \mathbb{N} ; we want to partition {L} into intervals that carry uniform “{q}-variation-mass”, so to speak. This is easy to achieve in the following way: let \mu_L \, : \, L \to [0,1] denote the function

\displaystyle \mu_L(\xi) := (\| m \|_{V_q (L \cap (-\infty,\xi])})^q;

that is, \mu_L(\xi) is the {q} -variation of {m} in the interval from the left endpoint of {L} to the point \xi \in L , raised to the power {q} (so that we have additivity). Function \mu_L is clearly a monotone increasing function, and therefore has a well-defined inverse function. Split therefore the interval [0,1] into 2^j equal disjoint intervals, that is,

\displaystyle [0,1] = J_1 \sqcup \ldots \sqcup J_{2^j}

with J_k := [2^{-j}(k-1), 2^{-j}k) ; we define then for any k = 1, \ldots, 2^j

\displaystyle I_{k, L} := \mu_L^{-1}(J_k),

which is an interval, by the monotonicity of \mu_L . With this definition, we have \|m\|_{V_q(I_{k,L})} \leq 2^{-j/q} by construction.
Define the collection \mathcal{J}_j to be the collection of all the intervals I_{k,L} resulting from this procedure:

\displaystyle \mathcal{J}_j := \{ I_{k,L} \; : \; L \in \mathbb{L}, k = 1,\ldots, 2^j\};

we have that:

  1. each \mathcal{J}_j collection is a partition of \mathbb{R} ;
  2. each \mathcal{J}_j collection is a refinement1 of the Littlewood-Paley partition and each L \in \mathbb{L} is partitioned by \mathcal{J}_j into at most 2^j intervals;
  3. the collections \mathcal{J}_j are refinements of each other: if j' > j then \mathcal{J}_{j'} is a refinement of \mathcal{J}_j . Notice that each interval of \mathcal{J}_{j} is split into (at most) 2 subintervals by \mathcal{J}_{j+1} .

Observe that property 2.) above points in the direction of (2), in that we have uniform control on the cardinality of the (sub-)partitions, but we still have to decompose the symbol. Using properties 1.) and 3.), we do a martingale decomposition of {m} adapted to these collections: let \mathcal{F}_j = \sigma(\mathcal{J}_j) (the sigma-algebra generated by the collection \mathcal{J}_j ) and notice that (\mathcal{F}_j)_{j\in\mathbb{Z}} is an increasing sequence of sigma-algebras; therefore, if we introduce the martingale differences

\displaystyle \mathbf{D}_j f := \mathbf{E}[f | \mathcal{F}_{j+1}] - \mathbf{E}[f | \mathcal{F}_{j}],

we can decompose every function {f} into

\displaystyle f = \mathbf{E}[f | \mathbb{L}] + \sum_{j \in \mathbb{N}} \mathbf{D}_j f.

Now it is worth observing what happens when we apply the martingale decomposition to {m} itself. Consider a fixed j and a fixed \xi_0 and observe that

\displaystyle \mathbf{D}_j m(\xi_0) = \frac{1}{|I|} \int_{I} m \,d\xi - \frac{1}{|\widehat{I}|} \int_{\widehat{I}} m \,d\xi,

where I is the unique interval in \mathcal{J}_{j+1} such that \xi_0 \in I and \widehat{I} is the unique interval in \mathcal{J}_{j} that contains I . Observe that there is a unique affine map \varphi that maps I to \widehat{I} while preserving the ordering; using this map as a change of variables we can thus write

\displaystyle \mathbf{D}_j m(\xi_0) = \frac{1}{|I|} \int_{I} m(\xi) - m(\varphi(\xi)) \,d\xi.

Notice though that for any \xi \in I we have that trivially
|m(\xi) - m(\varphi(\xi))| \leq \|m\|_{V_q(\widehat{I})} ; but by construction this quantity is controlled by 2^{-j/q} ! Therefore, since the above expression is an average, we have the pointwise bound

\displaystyle \| \mathbf{D}_j m\|_{L^\infty} \leq 2^{-j/q}.

Combining this bound with property 2.) of the collections \mathcal{J}_j we see that for the symbols m_j := \mathbf{D}_j m we have by inequality (2) that

\displaystyle \| T_{m_j} f\|_{L^p} \lesssim \|m_j\|_{L^\infty} (\sup_{L \in \mathbb{L}}\#\{I \subset L \, : \, I \in \mathcal{J}_j\})^{|1/2 - 1/p|} \|f\|_{L^p} \leq 2^{-j/q} 2^{j|1/2 - 1/p|} \|f\|_{L^p}.

If |1/2 - 1/p| < 1/q the overall exponent of 2^j above is negative and therefore the quantity is summable; by triangle inequality we have therefore that

\displaystyle \| T_{m} f\|_{L^p} \lesssim (\|m\|_{L^\infty} + C)\|f\|_{L^p}

provided the condition |1/2 - 1/p| < 1/q is satisfied (the \|m\|_{L^\infty} term comes from bounding \mathbf{E}[f | \mathbb{L}]  using the standard Littlewood-Paley square function). This concludes the nice proof of Theorem 2. \Box

We close this post with a remark and a question:

Remark: If we had not improved the exponent of N in inequality (2), we would have only concluded that \| T_{m_j} f\|_{L^p} \lesssim 2^{-j/q} 2^{j/2} \|f\|_{L^p} . This is still summable in {j} when q \in [1,2) , so we would still have concluded something, but we would have missed the q\geq 2 part of Theorem 2.

A question: Is the range \Big| \frac{1}{2} - \frac{1}{p}\Big| 2, can we find a symbol {m} with {\sup_{L \in \mathbb{L}} \|m\|_{V_q(L)} } finite but such that the multiplier T_m is unbounded for {p} such that { \Big| \frac{1}{2} - \frac{1}{p}\Big| \geq \frac{1}{q} } ?

Presumably the answer is yes, but I have not thought about the answer, honestly.

1: In the interest of clarity, by refinement we mean the following: for every I \in \mathcal{J}_j there exists an L \in \mathbb{L} such that I \subset L, and for each L \in \mathbb{L} there exist I_1, \ldots, I_n \in \mathcal{J}_j such that L = \bigcup_{\ell} I_\ell. [go back]

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google photo

You are commenting using your Google account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s