Advances in Mathematical Physics
Volume 2014 (2014), Article ID 215672, 11 pages
http://dx.doi.org/10.1155/2014/215672

Research Article

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

Shuhong Chen¹ and Zhong Tan²

¹School of Mathematics and Statistics, Minnan Normal University, Zhangzhou, Fujian 363000, China
²School of Mathematical Science, Xiamen University, Xiamen, Fujian 361005, China

Received 28 June 2014; Accepted 22 September 2014; Published 2 November 2014

Academic Editor: Salvador E. Venegas-Andraca

Copyright © 2014 Shuhong Chen and Zhong Tan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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

Then, from (84), we have

for some to be chosen. The last term above will be nonpositive provided the matrix

is nonnegative defined, that is, whenever

In that case, neglecting the last three terms on the right of (85) gives

This differential inequality shows that is uniformly bounded in time when and is exponentially bounded in time when . This then gives a priori control if , and are such that (88) is satisfied.
Control of the obviously gives control of the norm of . Assume that (80) and (88) are satisfied. The following choices for and will maximize the upper bounds (88) on . If , then for any and we may choose a large enough to satisfy (88) for any fixed value of , say, , and there is no restriction, where we have used (69). If , the result is obvious. If and , then choose and there is again no restriction on . If and and have the opposite signs, then choose for any value of . The maximum of the upper bounds (88) on is then obtained by setting and yields (69).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work is supported by National Natural Science Foundation of China (Grants nos. 11201415, 11271305); Program for New Century Excellent Talents in Fujian Province University (Grant no. JA14191).

References

M. Machida and T. Koyama, “Time-dependent Ginzburg-Landau theory for atomic Fermi gases near the BCS-BEC crossover,” Physical Review A, vol. 74, no. 3, Article ID 033603, 2006.

M. Drechsler and W. Zwerger, “Crossover from BCS-superconductivity to Bose-condensation,” Annals of Physics, vol. 504, no. 1, pp. 15–23, 1992.

C. A. R. Sá de Melo, M. Randeria, and J. R. Engelbrecht, “Crossover from BCS to Bose superconductivity: transition temperature and time-dependent Ginzburg-Landau theory,” Physical Review Letters, vol. 71, no. 19, pp. 3202–3205, 1993.

Y. He, Q. Chen, C. C. Chien, and K. Levin, “First- and second-sound-like modes at finite temperature in trapped Fermi gases from BCS to BEC,” Physical Review A, vol. 76, no. 5, Article ID 051602, 2007.

T. N. de Silva, “Breathing mode frequencies of a rotating Fermi gas in the BCS-BEC crossover region,” Physical Review A, vol. 78, no. 2, Article ID 023623, 2008.

Z. Q. Yu, K. Huang, and L. Yin, “Induced interaction in a Fermi gas with a BEC-BCS crossover,” Physical Review A, vol. 79, no. 5, Article ID 053636, 2009.

P. G. de Gennes, Superconductivity of Metals and Alloys, Westview Press, 1999.

A. A. Abrikosov, Fundamentals of the Theory of Metals, North-Holland, 1988.

C. R. Doering, J. D. Gibbon, D. D. Holm, and B. Nicolaenko, “Low-dimensional behaviour in the complex Ginzburg-Landau equation,” Nonlinearity, vol. 1, no. 2, pp. 279–309, 1988.

C. R. Doering, J. D. Gibbon, and C. D. Levermore, “Weak and strong solutions of the complex Ginzburg-Landau equation,” Physica D: Nonlinear Phenomena, vol. 71, no. 3, pp. 285–318, 1994.

L. Sirovich and J. D. Rodriguez, “Coherent structures and chaos: a model problem,” Physics Letters. A, vol. 120, no. 5, pp. 211–214, 1987.

S. Chen and B. Guo, “Existence of the weak solution of coupled time-dependent Ginzburg-Landau equations,” Journal of Mathematical Physics, vol. 51, no. 3, Article ID 033507, 2010.

S. Chen and B. Guo, “Classical solutions of time-dependent Ginzburg-Landau theory for atomic Fermi gases near the BCS-BEC crossover,” Journal of Differential Equations, vol. 251, no. 6, pp. 1415–1427, 2011.

P. Brenner and W. von Wahl, “Global classical solutions of nonlinear wave equations,” Mathematische Zeitschrift, vol. 176, no. 1, pp. 87–121, 1981.

M. Reed, Abstract Non-Linear Wave Equations, vol. 507 of Lecture Notes in Mathematics, Springer, Berlin, Germany, 1976.

T. Peter and S. Huang, “Global smooth solutions of the complex Ginzburg-Landau equation and their dynamical properties,” Discrete and Continuous Dynamical Systems, vol. 5, no. 4, pp. 825–848, 1999.

N. Hayashi, “Classical solutions of nonlinear Schrödinger equations,” Manuscripta Mathematica, vol. 55, no. 2, pp. 171–190, 1986.

H. Pecher, “-Abschätzungen und klassische Lösungen für nichtlineare Wellengleichungen I,” Mathematische Zeitschrift, vol. 150, pp. 159–183, 1976.

M. Tsutsumi and N. Hayashi, “Classical solutions of nonlinear Schrödinger equations in higher dimensions,” Mathematische Zeitschrift, vol. 177, no. 2, pp. 217–234, 1981.

H. Pecher, “-Abschätzungen und klassische Lösungen für nichtlineare Wellengleichungen II,” Manuscripta Mathematica, vol. 20, pp. 227–244, 1977.

J. Bergh and J. Löfström, Interpolation Spaces, Springer, New York, NY, USA, 1976.

H. Triebel, Interpolation Theory, Function Spaces, Differential Operators, North Holland, Amsterdam, The Netherlands, 1978.

H. Pecher and W. von Wahl, “Time dependent nonlinear Schrödinger equations,” Manuscripta Mathematica, vol. 27, no. 2, pp. 125–157, 1979.