Research Article | Open Access

Shuhong Chen, Zhong Tan, "Solution Theory of Ginzburg-Landau Theory on BCS-BEC Crossover", *Advances in Mathematical Physics*, vol. 2014, Article ID 215672, 11 pages, 2014. https://doi.org/10.1155/2014/215672

# Solution Theory of Ginzburg-Landau Theory on BCS-BEC Crossover

**Academic Editor:**Salvador E. Venegas-Andraca

#### Abstract

We establish strong solution theory of time-dependent Ginzburg-Landau (TDGL) systems on BCS-BEC crossover. By the properties of Besov, Sobolev spaces, and Fourier functions and the method of bootstrapping argument, we deduce that the global existence of strong solutions to time-dependent Ginzburg-Landau systems on BCS-BEC crossover in various spatial dimensions.

#### 1. Introduction

In this paper, we consider strong solutions to time-dependent Ginzburg-Landau systems, where the definition domain of Laplacian is , and is a bounded domain in with Lipschitzian boundary. The parameters and are all coupling coefficients, , , and is generally complex.

The main purpose of this paper is to establish strong solutions theory of (1)â€“(4). The systems are the general form of Ginzburg-Landau theory for superfluid atomic Fermi gases describing the BCS-BEC crossover near the Feshbach resonance from the Fermion-Boson model (double-channel model) [1]. The BCS-BEC crossover phenomena have been found as early as in 1992 [2, 3]. Due to the strange feature taken on by the quantum phenomena, it attracts many scientistsâ€™ attention and interest [4â€“6]. Though so many results have been found, there are a few works got by mathematical framework.

Even though the Ginzburg-Landau theory can capture almost every unique feature that the superfluid exhibited macroscopically [7, 8], this leads to the fact that it played an important role in the history of superfluid atomic Fermi gases research. In fact, the Ginzburg-Landau equation (single-channel model) has proved fruitful for illustration of the connection between infinite dimensional dynamics and finite dimensional dynamical systems [7â€“11]. Thus, Machida-Koyama constructed a time-dependent Ginzburg-Landau theory for BCS-BEC crossover from the Fermion-Boson model [1]. Then, Chen and Guo found weak solutions [12] and classical solutions [13] theories of the time-dependent Ginzburg-Landau theory. In present paper, we would further show the global strong solutions of initial value problem (1)â€“(4) under some suitable conditions.

In the strong solution theory, the admissible parameter values of and the dimension are interrelated. This arbitrariness contrasts sharply with the theory of weak solutions [10]. Luckily, we can overcome this problem by Fourier transform and using the properties of Fourier functions. However, the higher the dimension is, the harder the problem is [14]. The difficulty lies in obtaining a priori estimate of higher derivatives of solutions. Furthermore, the approaches in both [15, 16] could not be applied to the current paper, since the crossover between the states BCS and BEC can bring some technical difficulties. In order to handle these difficulties, we have to use the properties of Besov and Sobolev spaces and the bootstrapping argument. But what make things worse is that a priori estimate even in the case of nonlinear wave [17â€“20] requires some restrictions on the nonlinear term, which cause the fact that one just can obtain the desired result in the case of . Thus, combining all techniques of Fourier transform, Besov and Sobolev spaces, and the energy method and bootstrapping argument, we overcome these difficulties and establish the strong solutions theory of (1)â€“(4) in various spatial dimensions.

Theorem 1. *Let be a bounded domain in , with Lipschitzian boundary. Assume that are positive coefficients and are real numbers. For , assume that
**
or let , provided
**
Then the problem (1)â€“(4) with initial conditions has a unique global strong solution.*

Theorem 2. *Under the coefficientsâ€™ conditions of Theorem 1, for and , the problem (1)â€“(4) has a unique global strong solution if the dispersive parameters lie in the region of the -plane, bounded by hyperbolae, which satisfies
*

*Remark 3. *The relations between physics and mathematics are closed and they can complement each other. The physical phenomenon can validate the mathematical results; on the other hand, the mathematical results would provide the theoretical support for physical phenomena. For example, Theoremâ€‰â€‰1 means that, under suitable conditions, the transformation process of the two kinds of particles in BCS-BEC crossover region exists a smooth state. Therefore, in practice, we can control whether the smooth state should take on by adjusting the related conditions. Theoremâ€‰â€‰2 shows that, under certain conditions, we can decide that which region should take on the smooth state by adjusting parameters.

#### 2. Local Existence of Strong Solutions

In this section, we would establish local strong solutions to (1)â€“(4). At first, the initial boundary value problem (1)-(2) may be rewritten as follows:

For with , (1)â€“(4) become where Obviously, is a nonnegative defined matrix.

Consider the evolution of in a Banach space to be governed by the abstract initial-value problem (10) with and the perturbation is often a nonlinear and noncontinuous map over . Then, via the corresponding integral equation assume that is a locally Lipschitz continuous map from into itself; that is, (i) for every ,(ii) for every ,

where and are nondecreasing functions. Then, employing the contraction mapping theorem, one can get the following basic result.

Lemma 4 (local existence theorem). *For every there exists a time such that, for every initial data , with , there exists a unique satisfying the mild formulation (12). Moreover, is a locally Lipschitz continuous function of .*

*Definition 5. *A function satisfying the mild formulation (12) is called a mild solution for the initial value problem (10).

*Definition 6. *A function satisfying (10) with in space and in time is called a strong solution for the initial value problem (10).

Lemma 7 (-local strong solutions). *For every there exists a time such that for every initial data with there exists a unique
**
satisfying the initial value problem (10). Moreover, for every initial data , one has .*

*Proof. *By Fourier transform and the properties of Fourier functions, the associated operator of (12) acting on can be written as a convolution, , with its Green function given by

The integral equation (12) recast in terms of this Green function takes the form

We first appeal to . The perturbation as a map from into itself is clearly locally Lipschitz continuous. Lemma 4 yields a unique mild solution of (12) over a time interval that depends only on the norm of . This solution is the limit of a sequence of successive iterates of (15), say, that defined by
that converge in for some chosen sufficiently small such that the sequence contracts.

In order to elevate these mild solutions to strong solutions, we use standard bootstrapping argument. First, evoking the regularity of , the gradient of (15) yields
Then, use the -estimate of with
where is a constant depending only on the dimension , applied to the successive iterates (16) to show that each lies in and that the sequence converges. It follows that with

Then, if it was the case that , then the singularity in this estimate at disappears and one sees that and that is a solution of
where represents the derivative of with respect to , where in this case is given explicitly by

A repetition of the above regularity argument starting from (20) shows that is in . Moreover, because (8) and (9) relate the first time derivative to the second space derivative, the solution must also be in and is therefore a strong solution of systems (1)â€“(4) so long as it is a mild solution.

Lemma 8 (local solutions). *Let for some positive integer . Then for every there exists a time such that for every initial data with there exists a unique
**
satisfying (1)â€“(4). Moreover, for every initial data , one has
*

*Proof. *From (21), we observe that the lowest degree of homogeneity for the factors and appearing in a term of the derivative of will be , and this can be controlled whenever . In that case, the bootstrapping argument will gain additional spatial derivatives, showing that the solution is in . This is Lemma 8.

Refinements of the basic existence argument given above greatly enlarge the class of initial data that evolve into strong solutions for finite times.

Lemma 9 (-local mild solutions). *If satisfies
**
then for every there exists a time such that for every initial data with and , with , , and , there exists a unique
**
satisfying the mild formulation (15) of the initial value problem (1)â€“(4).*

*Proof. *Noting that function (14) satisfies the -estimate

Controlling the norm of follows from the bound
for some positive absolute constants and . By an interpolation argument, the norm of is estimated by

Thus, one sees that

Then, a direct estimate of the mild formulation (15) for every satisfying
yields
where and satisfy
The idea is to recast (32) as an inequality for the ratio
whence obtaining

It is easily seen from (30) that whenever
the condition
holds. The iterates (16) contract in the space that is the completion of for some sufficiently small, in the norm

At last, we show that Lemma 9 holds just as the condition (25) is satisfied. By the second equality of (33) and then the first one,
thus, we have , which means that the second inequality of (36) is seen to be equivalent to the condition
Again, using (33), the first inequality of (36) becomes
which implies that
By a direct calculation and using the condition that , we deduce
then, one sees that (42) is satisfied by choosing ; thus, condition (37) is met and the contraction mapping argument yields a unique solution in for .

This solution is in . Indeed, it is clear from the definition of the -norm (38) that . To see that is also in , one only needs to check the continuity of at . First, subtracting from both sides of the mild formulation (15), a direct estimate gives
where and are defined in (33). The strong continuity of the linear semigroup implies that the first term on the right side vanishes as tends to zero. The boundedness of in and the bound (27) of imply that the second term on the right side will also vanish as tends to zero whenever the singularity of at is integrable. By (33), this will be the case if and only if
which was already controlled in (36) by the assumption that . Therefore, is continuous at and hence is in . This is Lemma 9.

Lemma 10 (-local strong solutions). *If satisfies
**
then for every there exists a time such that for every initial data with there exists a unique
**
satisfying the initial value problem (1)â€“(4).*

*Proof. *In order to verify that a mild solution in is a strong solution in , we follow from Lemma 4 upon showing it is in . Thus, if then every satisfies these criteria, and in particular . If, on the other hand, then the solution is at least in for . Then, when , a second application of the above argument then shows that . More generally, a simple argument shows that whenever
after applying the above argument times, it follows that , and one more application then shows that .

Combining the above argument with Lemmas 4 and 8 yields the desired result.

#### 3. Global Existence of Strong Solutions

From Lemma 10 in the last section, one can find that it suffices to obtain global control of any norm where satisfies (46). Precisely, we have the following lemma.

Lemma 11. *For and assume that coefficients parameters satisfy the conditions in Theorem 1, there is
**
And let , provided
**
then
*

*Proof. *A direct calculation and an application of HÃ¶lderâ€™s inequality yield

For this differential inequality shows

By Gronwallâ€™s inequality, we have
which implies that the norm of is uniformly bounded in time.

For , multiplying (8) by and taking the real part,

Notice that
then,

Noting that , (55) becomes

Using Fourier transform, (9) can be reformed into
where .

Let and ; then, from the properties of Besov space [21, 22], for ,

The integral equation (59) gives [23]
where . Thus,

Noting that , then, by HÃ¶lderâ€™s and Youngâ€™s inequalities, we have

And (58) implies

The former two inequalities mean that

By Gronwallâ€™s inequality, we obtain

This completes the proof of Lemma 11.

*Proof of Theorem 1. *From Lemma 11, we can see that the norm of is uniformly bounded in time. Then whenever
Lemma 10 then implies that the problem (1)â€“(4) has global strong solutions.

When (67) is violated, one must control more than the norm. Thus when in addition to (50) one also has
Lemma 10 then implies that the problem (1)â€“(4) has global strong solutions. We can find a satisfying both (50) and (68) provided, which means that

Thus, by Lemma 10, we get the global strong solutions for (1)â€“(4) with the initial conditions. This completes the Proof of Theorem 1.

*Proof of Theorem 2. *An alternation to directly controlling the norm is to rather first directly control the norm and then control through Sobolev estimate. As the norm is already controlled by (49), it suffices to control norm of . Assuming that and utilizing the decomposition identity
then a direct calculation shows that
Similarly, as the estimate of (57),
Provided that
the last term of (71) will be nonpositive. And then it becomes

Now, multiplying (9) with and taking the real part,

Combining (74) and (75) gives

Via Gronwallâ€™s inequality, we gain
This differential inequality shows that the norm of for is uniformly bounded in time. This then gives a priori control of the norm if condition (73) is satisfied.

Control of the norm yields, by a Sobolev inequality, the norm for every satisfying

When, in addition to (73) and (78), one also has

Lemma 10 implies that the problem (1)â€“(4) has global strong solutions. Given that , , and satisfy (73), we can find a satisfying both (78) and (79), provided that

This encompasses all subcritical and critical nonlinearities in every spatial dimension (), the supercritical cubic nonlinearity () in dimension , and all nonlinearities in dimensions .

Equations (65) and (77) may be combined to gain control on a set that is actually much larger than the union of the strips (50) and (73) about the axes of the -plane. This is reflected in the condition of (7) in Theorem 2.

We will produce a bounds on a functional of the form
where the parameter is to be chosen later.

From (70), we can obtain the identities

Taking a linear combination of (55), (62), (74), and (75), while using (82), and noting that give
where