
In this paper, we study the phase field models with fractionalorder in time.
The phase field models have been widely used to study coarsening dynamics of
material systems with microstructures. It is known that phase field models are
usually derived from energy variation so that they obey some energy dissipation
laws intrinsically. Recently, many works have been published on investigating
fractionalorder phase field models, but little is known of the corresponding
energy dissipation laws. We focus on the timefractional phase field models and
report that the effective free energy and roughness obey a universal powerlaw
scaling dynamics during coarsening. Mainly, the effective free energy and
roughness in the timefractional phase field models scale by following a
similar power law as the integer phase field models, where the power is
linearly proportional to the fractional order. This universal scaling law is
verified numerically against several phase field models, including the
CahnHilliard equations with different variable mobilities and molecular beam
epitaxy models. This new finding sheds light on potential applications of time
fractional phase field models in studying coarsening dynamics and crystal
growths.

In this paper, we propose a timefractional molecular beam epitaxy (MBE)
model with slope selection and its efficient, accurate, full discrete, linear
numerical approximation. The numerical scheme utilizes the fast algorithm for
the Caputo fractional derivative operator in time discretization and Fourier
spectral method in spatial discretization. Refinement tests are conducted to
verify the $2\alpha$ order of time convergence, with $\alpha \in (0, 1]$ the
fractional order of derivative. Several numerical simulations are presented to
demonstrate the accuracy and efficiency of our newly proposed scheme. By
exploring the fast algorithm calculating the Caputo fractional derivative, our
numerical scheme makes it practice for long time simulation of MBE coarsening,
which is essential for MBE model in practice. With the proposed fractional MBE
model, we observe that the scaling law for the energy decays as $
O(t^{\frac{\alpha}{3}})$ and the roughness increases as
$O(t^{\frac{\alpha}{3}})$, during the coarsening dynamics with random initial
condition. That is to say, the coarsening rate of MBE model could be
manipulated by the fractional order $\alpha$, and it is linearly proportional
to $\alpha$. This is the first time in literature to report/discover such
scaling correlation. It provides a potential application field for fractional
differential equations. Besides, the numerical approximation strategy proposed
in this paper can be readily applied to study many classes of timefractional
and high dimensional phase field models.

In this paper, we consider the numerical approximations for the fourth order
CahnHilliard equation with concentration dependent mobility, and the
logarithmic FloryHuggins potential. One challenge in solving such a diffusive
system numerically is how to develop proper temporal discretization for
nonlinear terms in order to preserve the energy stability at the timediscrete
level. We resolve this issue by developing a set of the first and second order
time marching schemes based on a novel, called "Invariant Energy
Quadratization" approach. Its novelty is that the proposed scheme is linear and
symmetric positive definite because all nonlinear terms are treated
semiexplicitly. We further prove all proposed schemes are unconditionally
energy stable rigorously. Various 2D and 3D numerical simulations are presented
to demonstrate the stability, accuracy and efficiency of the proposed schemes
thereafter.

How to develop efficient numerical schemes while preserving the energy
stability at the discrete level is a challenging issue for the three component
CahnHilliard phasefield model. In this paper, we develop first and second
order temporal approximation schemes based on the "Invariant Energy
Quadratization" approach, where all nonlinear terms are treated
semiexplicitly. Consequently, the resulting numerical schemes lead to a
wellposed linear system with the symmetric positive definite operator to be
solved at each time step. We rigorously prove that the proposed schemes are
unconditionally energy stable. Various 2D and 3D numerical simulations are
presented to demonstrate the stability and the accuracy of the schemes.

In this paper, we consider the numerical approximations for the fourth order
viscous CahnHilliard equation with the hyperbolic relaxation. The main
challenge in solving such a diffusive system numerically is how to develop high
order temporal discretization for the hyperbolic and nonlinear terms that
allows large time step while preserving the unconditional energy stability,
i.e., the energy dissipative structure at the timediscrete level. We resolve
this issue by developing two second order time marching schemes using the
recently developed "Invariant Energy Quadratization" approach where all
nonlinear terms are treated semiexplicitly. In each time step, one only needs
to solve a linear system that is symmetric positive definite. We rigorously
prove all proposed schemes are unconditionally energy stable, which is further
verified by time step refinement test numerically. Moreover, a number of 2D and
3D numerical simulations are presented to demonstrate the stability, accuracy
and efficiency of the proposed schemes.

Let $G=(V,E)$ and $H$ be two graphs. Packing problem is to find in $G$ the
largest number of independent subgraphs each of which is isomorphic to $H$. Let
$U\subset{V}$. If the graph $GU$ has no subgraph isomorphic to $H$, $U$ is a
cover of $G$. Covering problem is to find the smallest set $U$. The
vertexdisjoint tree packing was not sufficiently discussed in literature but
has its applications in data encryption and in communication networks such as
multicast routing protocol design. In this paper, we give the kind of
$(k+1)$connected graph $G'$ into which we can pack independently the subgraphs
that are each isomorphic to the $(2^{k+1}1)$order perfect binary tree $T_k$.
We prove that in $G'$ the largest number of vertexdisjoint subgraphs
isomorphic to $T_k$ is equal to the smallest number of vertices that cover all
subgraphs isomorphic to $T_k$. Then, we propose that $T_k$ does not have the
\emph{Erd\H{o}sP\'{o}sa} property. We also prove that the $T_k$ packing
problem in an arbitrary graph is NPhard, and propose the distributed
approximation algorithms.

We say that a cyclotomic polynomial \Phi_{n}(x) has order three if n is the
product of three distinct primes, p<q<r. Let A(n) be the largest absolute value
of a coefficient of \Phi_{n}(x) and M(p) be the maximum of A(pqr). In 1968,
Sister Marion Beiter conjectured that A(pqr)<=(p+1)/2. In 2008, Yves Gallot and
Pieter Moree showed that the conjecture is false for every p>=11, and they
proposed the Corrected Beiter conjecture: A(pqr)<=2p/3. Here we will give a
proof of this conjecture.

We say that a cyclotomic polynomial Phi_{n}(x) has order three if n is the
product of three distinct primes, p<q<r. Let A(n) be the largest absolute value
of a coefficient of Phi_{n}(x) and M(p) be the maximum of A(pqr). In 1968,
Sister Marion Beiter conjectured that A(pqr)<=(p+1)/2. In 2008, Yves Gallot and
Pieter Moree showed that the conjecture is false for every p>=11, and they
proposed the Corrected Beiter conjecture: M(p)<=2p/3. Here we will give a
sufficient condition for the Corrected Beiter conjecture and prove it when p=7.