Review Special Issues

Recent progress on mathematical analysis and numerical simulations for Maxwell's equations in perfectly matched layers and complex media: a review

  • Received: 18 December 2023 Revised: 02 February 2024 Accepted: 18 February 2024 Published: 04 March 2024
  • In this paper, we presented a review on some recent progress achieved for simulating Maxwell's equations in perfectly matched layers and complex media such as metamaterials and graphene. We mainly focused on the stability analysis of the modeling equations and development and analysis of the numerical schemes. Some open issues were pointed out, too.

    Citation: Jichun Li. Recent progress on mathematical analysis and numerical simulations for Maxwell's equations in perfectly matched layers and complex media: a review[J]. Electronic Research Archive, 2024, 32(3): 1901-1922. doi: 10.3934/era.2024087

    Related Papers:

    [1] Yuri Prokhorov . Conic bundle structures on Q-Fano threefolds. Electronic Research Archive, 2022, 30(5): 1881-1897. doi: 10.3934/era.2022095
    [2] János Kollár . Relative mmp without Q-factoriality. Electronic Research Archive, 2021, 29(5): 3193-3203. doi: 10.3934/era.2021033
    [3] Meng Chen, Yong Hu, Matteo Penegini . On projective threefolds of general type with small positive geometric genus. Electronic Research Archive, 2021, 29(3): 2293-2323. doi: 10.3934/era.2020117
    [4] Jiange Liu, Yu Chen, Xin Dai, Li Cao, Qingwu Li . MFCEN: A lightweight multi-scale feature cooperative enhancement network for single-image super-resolution. Electronic Research Archive, 2024, 32(10): 5783-5803. doi: 10.3934/era.2024267
    [5] Haijun Wang, Wenli Zheng, Yaowei Wang, Tengfei Yang, Kaibing Zhang, Youlin Shang . Single hyperspectral image super-resolution using a progressive upsampling deep prior network. Electronic Research Archive, 2024, 32(7): 4517-4542. doi: 10.3934/era.2024205
    [6] Zhengyu Liu, Yufei Bao, Changhai Wang, Xiaoxiao Chen, Qing Liu . A fast matrix completion method based on truncated L2,1 norm minimization. Electronic Research Archive, 2024, 32(3): 2099-2119. doi: 10.3934/era.2024095
    [7] Zhiming Cai, Liping Zhuang, Jin Chen, Jinhua Jiang . Lightweight high-performance pose recognition network: HR-LiteNet. Electronic Research Archive, 2024, 32(2): 1145-1159. doi: 10.3934/era.2024055
    [8] Jongho Kim, Woosuk Kim, Eunjeong Ko, Yong-Shin Kang, Hyungjoo Kim . Estimation of spatiotemporal travel speed based on probe vehicles in mixed traffic flow. Electronic Research Archive, 2024, 32(1): 317-331. doi: 10.3934/era.2024015
    [9] Limin Guo, Weihua Wang, Cheng Li, Jingbo Zhao, Dandan Min . Existence results for a class of nonlinear singular p-Laplacian Hadamard fractional differential equations. Electronic Research Archive, 2024, 32(2): 928-944. doi: 10.3934/era.2024045
    [10] Li Du, Xiaoqin Yuan . The minimality of biharmonic hypersurfaces in pseudo-Euclidean spaces. Electronic Research Archive, 2023, 31(3): 1587-1595. doi: 10.3934/era.2023081
  • In this paper, we presented a review on some recent progress achieved for simulating Maxwell's equations in perfectly matched layers and complex media such as metamaterials and graphene. We mainly focused on the stability analysis of the modeling equations and development and analysis of the numerical schemes. Some open issues were pointed out, too.



    Smooth varieties have certain nice properties, and both algebraic and analytic methods can be applied to them. However, when studying problems in birational geometry, particularly those related to the minimal model program, it becomes necessary to investigate varieties with singularities. Fortunately, as Hironaka's famous theorem states, every variety in characteristic zero has a resolution of singularity. The existence of resolutions of singularities provides a way to study singular varieties. For instance, by comparing a variety to its resolution of singularities, one can measure the complexity of a singularity. This is a fundamental technique in higher-dimensional birational geometry.

    Since we want to understand a singularity through its resolution, it is natural to inquire about the difference between two distinct resolutions of singularities. For an algebraic surface S, there exists a smooth surface known as the minimal resolution of S. This is a resolution of singularities ˉSS such that ρ(ˉS/S) is minimal. The minimal resolution ˉS is unique, and any birational morphism SS from a smooth surface S to S factors through ˉSS.

    In this paper, we want to find a higher-dimensional analog of the minimal resolution for surfaces. It is not reasonable to assume the existence of a unique minimal resolution for higher-dimensional singularities. For instance, if XX is a smooth flop over W, then X and X are two different resolutions of singularities for W. Since flops are symmetric (at least in dimension three), it appears that X and X are both minimal. Thus, we need to consider the following issues:

    (1) To define "minimal resolutions", which ideally should be resolutions of singularities with some minimal geometric invariants.

    (2) To compare two different minimal resolutions. We need some symmetry between them, so that even if minimal resolutions are not unique, it is not necessary to distinguish them.

    (3) To compare a minimal resolution with an arbitrary resolution of singularities.

    Inspired by the two-dimensional case, it is natural to consider resolutions of singularities with minimal Picard number (we will call such resolutions P-minimal resolutions, see Section 7 for a more precise definition). We know that a fixed singularity may have more than one P-minimal resolution, and two different P-minimal resolutions can differ by a smooth flop. It is also possible that two different P-minimal resolutions "differ by a singular flop": consider XX as a possibly singular flop over W. Let ˜XX and ~XX be P-minimal resolutions of X and X, respectively. Then, because of the symmetry between flops, one may expect that ˜X and ~X are two different P-minimal resolutions of W. We call the birational map ˜X~X a P-desingularization of the flop XX (a precise definition can be found in Section 7). If we consider P-desingularizations of flops as elementary birational maps, then in dimension three, P-minimal resolutions have nice properties.

    Theorem 1.1. Assume that X is a projective threefold over the complex numbers and ˜X1, ˜X2 are two different P-minimal resolutions of X. Then ˜X1 and ˜X2 are connected by P-desingularizations of terminal and Q-factorial flops.

    Moreover, if X has terminal and Q-factorial singularities, then the birational map ˜X1˜X2 has an Ω-type factorization.

    Please see Section 6 for the definition of Ω-type factorizations.

    Theorem 1.2. Assume that X is a projective threefold over the complex numbers and WX is a birational morphism from a smooth threefold W to X. Then, for any P-minimal resolution ˜X of X, one has a factorization

    W=˜Xk...˜X1˜X0=˜X

    such that ˜Xi+1˜Xi is either a smooth blow-down or a P-desingularization of a terminal Q-factorial flop.

    Since three-dimensional terminal flops are topologically symmetric, some topological invariants like Betti numbers will not change after P-desingularizations of terminal flops. Hence, it is easy to see that P-minimal resolutions are the resolution of singularities with minimal Betti numbers.

    Corollary 1.3. Assume that X is a projective threefold over the complex numbers and WX is a birational morphism from a smooth threefold W to X. Then, for any P-minimal resolution ˜X of X, one has that bi(˜X)bi(W) for all i=0, ..., 6.

    Although in dimension three P-minimal resolutions behave well, for singularities of dimension greater than three, P-minimal resolutions may not be truly "minimal". A simple example is a smooth flip. If XX is a smooth flip over W, then both X and X are P-minimal resolutions of W, but X is better than X. Notice that the only known smooth flips are standard flips [1, Section 11.3], and if XX is a standard flip, then it is easy to see that bi(X)bi(X) for all i and the inequality is strict for some i. Thus, the resolution of singularities with minimal Betti numbers may be the right minimal resolution for higher-dimensional singularities. Because of Corollary 1.3, in dimension three P-minimal resolutions are exactly those smooth resolutions which have minimal Betti numbers. Therefore, this new definition of minimal resolutions is compatible with our three-dimensional theorems.

    We now return to the proof of our main theorems. Let X be a threefold and WX be a resolution of singularities. One can run KW-MMP over X as

    W=X0X1...Xk=X.

    Let ˜Xi be a P-minimal resolution of Xi. Then ˜X0=W and it is easy to see that ˜Xk is also a P-minimal resolution of X. Thus, our main theorems can be easily proved if we know the relation between ˜Xi and ˜Xi+1. Since Xi has only terminal and Q-factorial singularities, studying P-minimal resolutions of Xi becomes simpler.

    In [2], Chen introduced feasible resolutions for terminal threefolds, which is a resolution of singularities consisting of a sequence of divisorial contractions to points with minimal discrepancies (see Section 2.3.3 for more detail). Given a terminal threefold X and a feasible resolution ˉX of X, one can define the generalized depth of X to be the integer ρ(ˉX/X). The generalized depth is a very useful geometric invariant of a terminal threefold. In our application, the crucial factor is that one can test whether a resolution of singularities WX is a feasible resolution or not by comparing ρ(W/X) and the generalized depth of X. We need to understand how generalized depths change after steps of the minimal model program. After that, we can prove that for terminal and Q-factorial threefolds, P-minimal resolutions and feasible resolutions coincide.

    Now we only need to figure out the following two things: how generalized depths change after a step of minimal model program (MMP), and how P-minimal resolutions change after a step of MMP. To answer those questions, we have to factorize a step of MMP into more simpler birational maps. In [3], Chen and Hacon proved that three-dimensional terminal flips and divisorial contractions to curves can be factorized into a composition of (inverses of) divisorial contractions and flops. In this paper, we construct a similar factorization for divisorial contractions to points. After knowing the factorization, we are able to answer the two questions above and prove our main theorems.

    In addition to the above, we introduce the notion of Gorenstein depth for terminal threefolds. The basic idea is as follows: given a sequence of steps of MMP of terminal threefolds

    X0X1...Xk,

    one can show that the generalized depth of Xk is bounded above by the integer k and the generalized depth of X0. That is to say, the number of steps of MMP bounds the singularities on the minimal model. One may ask whether there is an opposite bound. Specifically, if we know the singularities of the minimal model Xk, can we bound singularities of X0? In this paper, we define the Gorenstein depth of terminal threefolds, which roughly speaking measures only Gorenstein singularities. One can show that the Gorenstein depth is always non-decreasing when running three-dimensional terminal MMP. Our result on Gorenstein depth will have important applications in [4].

    This paper is structured as follows: Section 2 is a preliminary section. In Section 3, we develop some useful tools to construct relations between divisorial contractions to points. Those tools, as well as the explicit classification of divisorial contractions, will be used in Section 4 to construct links of different divisorial contractions to points. In Section 5, we prove the property of the generalized depth. The construction of diagrams in Theorem 1.1 will be given in Section 6. All our main theorems will be proved in Section 7. In the last section, we discuss possible higher-dimensional generalizations of the notion of minimal resolutions, and possible applications of our main theorems.

    In this paper we only consider varieties over complex numbers.

    Let X and Y be two algebraic varieties. We say that X and Y are birational if there exists Zariski open sets UX and VY such that U and V are isomorphic. If X and Y are birational, we say that ϕ:XY is a birational map, and we will denote ϕ|U to be the isomorphism UV. If ϕ:XY is a morphism between X and Y and there exists a Zariski open set UX such that ϕ|U is an isomorphism, then we say that ϕ is a birational morphism.

    For a divisorial contraction, we mean a birational morphism YX which contracts an irreducible divisor E to a locus of codimension at least two, such that KY is Q-Cartier and is anti-ample over X. We will denote by vE the valuation that corresponds to E.

    Let G be a cyclic group of order r generated by τ. For any Z-valued n-tuple (a1,...,an), one can define a G-action on An(x1,...,xn) by τ(xi)=ξaixi, where ξ=e2πir. We will denote the quotient space An/G by An(x1,...,xn)/1r(a1,...,an).

    We say that w is a weight on W/G=An(x1,...,xn)/G defined by w(x1,...,xn)=1r(b1,...,bn) if w is a map OW1rZ0 such that

    w((i1,...,in)Zn0c(i1,...,in)xi11...xinn)=min{1r(b1i1+...+bnin) |c(i1,...,in)0}.

    Assume that ϕ:XY is a birational map. Let UX be the largest open set such that ϕ|U is an isomorphism and ZX be an irreducible subset such that Z intersects U non-trivially. We will denote by ZY the closure of ϕ|U(Z|U).

    Let W=An and G be a finite cyclic group, such that ˉW=W/GAn(x1,...,xn)/1r(a1,...,an). There is an elementary way to construct a birational morphism WˉW, so called the weighted blow-up, defined as follows.

    We write everything in the language of toric varieties. Let N be the lattice e1,...,en,vZ, where e1, ..., en is the standard basis of Rn and v=1r(a1,...,an). Let σ=e1,...,enR0. We have ˉWSpecC[Nσ].

    Let w=1r(b1,...,bn) be a vector such that bi=λai+kir for λN and kiZ with bi0. We define a weighted blow-up of ˉW with weight w to be the toric variety defined by the fan consisting of the cones

    σi=e1,...,ei1,w,ei+1,...,en.

    Let Ui be the toric variety defined by the cone σi and lattice N, namely

    Ui=SpecC[Nσi].

    Lemma 2.1. One has that

    UiAn/τ,τ

    where τ is the action given by

    xiξrbixi,xjξbjbixj,ji

    and τ is the action given by

    xiξaibixi,xjξajbiaibjrbixj,ji.

    Here, ξk denotes a k-th roots of unity for any positive integer k.

    In particular, the exceptional divisor of WˉW is P(b1,...,bn)/G where G is a cyclic group of order m where m is an integer that divides λ.

    Proof. Let Ti be a linear transformation such that Tiej=ej if ji and Tiw=ei. One can see that

    Tiei=rbi(eijibjrej)

    and

    Tiv=jiajrej+airrbi(eijibjrej)=aibiei+jiajbiaibjrbjej.

    Under this linear transformation, σi becomes the standard cone e1,...,enR0. Note that

    kiTiei+λTiv=kir+λaibiei+jiλ(ajbiaibj)kibjrrbiej=ei+jiλajbibibjrbiej=eijikjej.

    Hence, eiTiN and TiN=e1,...,en,Tiei,TivZ. Now Tiei corresponds to the action τ and Tiv corresponds to the action τ. This means that UiAn/τ,τ.

    The computation above shows that τki=τλ. If we glue (xi=0)An/τ together, then we get a weight projective space P(b1,...,bn). The relation τki=τλ implies that (xi=0)Ui can be viewed as P(b1,...,bn)/G where G is a cyclic group of order m for some factor m of λ.

    Corollary 2.2. Let x1, ..., xn be the local coordinates of W and let y1, ..., yn be the local coordinates of Ui. The change of coordinates of the morphism UiˉW are given by xj=yjybjri and xi=ybiri.

    Proof. The change of coordinates is defined by Tti, where Ti is defined as in Lemma 2.1.

    Corollary 2.3. Assume that

    S=(f1(x1,...,xn)=...=fk(x1,...,xn)=0)ˉW

    is a complete intersection and S is the proper transform of S on W. Assume that the exceptional locus E of SS is irreducible and reduced. Then

    a(E,S)=b1+...+bnrki=1w(fk)1.

    Proof. Assume first that k=0. Denote ϕ:WˉW. Then, on Ui, we have

    ϕdx1...dxn=birybir1i(jiybjri)dy1...dyn,

    hence KW=ϕKˉW+(b1+...+bnr1)F where F=exc(WˉW).

    Now the statement follows from the adjunction formula.

    Corollary 2.4. Let F=exc(WˉW). Then

    Fn=(1)n1rn1b1...bnm.

    Here, m is the integer in Lemma 2.1.

    Proof. From the change of coordinate formula in Corollary 2.2, one can see that F|F=OP(b1,...,bn)/G(r). It follows that

    Fn=(F|F)n1=(1)n1rn1b1...bnm.

    Definition 2.5. Let ϕi:UiˉW be the morphism in Corollary 2.2. For any G-semi-invariant function uOW, we can define the strict transform of u on Ui by (ϕ1i)(u)=ϕ(u)/yw(u)i.

    In this paper, we will consider terminal threefolds which are embedded into a cyclic quotient of A4 or A5

    XA4(x,y,z,u)/1r(a,b,c,d)orXA5(x,y,z,u,t)/1r(a,b,c,d,e).

    We say that YX is a weighted blow-up with weight w if Y is the proper transform of X inside the weighted blow-up of A4(x,y,z,u)/1r(a,b,c,d) or A5(x,y,z,u,t)/1r(a,b,c,d,e) with weight w.

    Notation 2.6. Assume that X is of the above form and let YX be a weighted blow-up. The notation Ux, Uy, Uz, Uu and Ut will stand for U1, ..., U5 in Lemma 2.1.

    Notation 2.7. Assume that w is a weight on An(x1,...,xn) determined by w(x1,...,xn)=(a1,...,an) and

    f(x1,...,xn)=(i1,...,in)Zn0λi1,...,inxi11...xinn

    is a regular function on An. We denote

    fw=a1i1+...+anin=w(f)λi1,...,inxi11...xinn.

    The local classification of terminal threefolds was done by Reid [5] for Gorenstein cases and Mori [6] for non-Gorenstein cases.

    Definition 2.8. A compound Du Val point PX is a hypersurface singularity which is defined by f(x,y,z)+tg(x,y,z,t)=0, where f(x,y,z) is an analytic function which defines a Du Val singularity.

    Theorem 2.9 ([5,Theorem 1.1]). Let PX be a point of threefold. Then PX is an isolated compound Du Val point if and only if PX is terminal and KX is Cartier near P.

    Theorem 2.10 ([6, cf. [7,Theorem 6.1]). Let PX be a germ of three-dimensional terminal singularity such that KX has Cartier index r>1. Then

    X(f(x,y,z,u)=0)A4(x,y,z,u)/1r(a1,...,a4)

    such that f, r and ai are given by Table 1.

    Table 1.  Classification of terminal threefolds.
    Type f(x,y,z,u) r ai condition
    cA/r xy+g(zr,u) any (α,α,1,r) gm2P α and r are coprime
    cAx/4 xy+z2+g(u)x2+z2+g(y,u) 4 (1,1,3,2) gm3P
    cAx/2 xy+g(z,u) 2 (0,1,1,1) gm4P
    cD/3 x2+y3+z3+u3x2+y3+z2u+yg(z,u)+h(z,u)x2+y3+z3+yg(z,u)+h(z,u) 3 (0,2,1,1) gm4Phm6P
    cD/2 x2+y3+yzu+g(z,u)x2+yzu+yn+g(z,u)x2+yz2+yn+g(z,u) 2 (1,0,1,1) gm4P,n4n3
    cE/2 x2+y3+yg(z,u)+h(z,u) 2 (1,0,1,1) g,hm4Ph40

     | Show Table
    DownLoad: CSV

    Assume that PX is a three-dimensional terminal singularity. Then there exists a section H|KX| which has Du Val singularities (referred to as a general elephant). Please see [7,(6.4)] for details.

    Divisorial contractions to points between terminal threefolds are well-classified by Kawamata [8], Hayakawa [9,10,11], Kawakita [12,13,14] and Yamamoto [15].

    Theorem 2.11. Assume that YX is a divisorial contraction to a point between terminal threefolds. Then there exists an embedding XW with W=A4(x,y,z,u) or A5(x,y,z,u,t) and a weight w(x,y,z,u)=1r(a1,...,a4) or w(x,y,z,u,t)=1r(a1,...,a5), respectively, such that YX is a weighted blow-up with respect to w.

    The defining equation of XW and the weight are given in Table 2 to Table 11.

    Table 2.  Divisorial contractions to cAx/r points.
    No. defining equations (r;ai)weight typea(X,E) condition
    A1 xy+zrk+gka(z,u) (r;β,β,1,r)1r(b,c,a,r) cA/ra/r baβ(modr),b+c=rka
    A2 x2y2+z3+xu2+g6(x,y,z,u) (1;)(4,3,2,1) cA23 xzg(x,y,z,u)

     | Show Table
    DownLoad: CSV
    Table 3.  Divisorial contractions to cAx/r points.
    No. defining equations (r;ai)weight typea(X,E) condition
    Ax1 x2+y2+g2k+12(z,u) (4;1,3,1,2)14(b,c,1,2) cAx/41/4 (b,c)=(2k+1,2k+3) or (2k+3,2k+1)
    Ax2 x2+y2+(λx+μy)p2k+14(z,u)+g2k+32(z,u) (4;1,3,1,2)14(b,c,1,2) cAx/41/4 (b,c,λ,μ)=(2k+5,2k+3,1,0) or (2k+3,2k+5,0,1)
    Ax3 x2+y2+gk(z,u) (2;0,1,1,1)12(b,c,1,1) cAx/21/2 (b,c)=(k,k+1) or (k+1,k)
    Ax4 x2+y2+(λx+μy)pk2(z,u)+gk+1(z,u) (2;0,1,1,1)12(b,c,1,1) cAx/21/2 (b,c,λ,μ)=(k+2,k+1,1,0) or (k+1,k+2,0,1)

     | Show Table
    DownLoad: CSV
    Table 4.  Divisorial contractions to cD points with discrepancy one.
    No. defining equations weight typea(X,E) condition
    D1 x2+y2u+λyzk+gl(z,u) (b,b1,1,2) cD1 b=min{k1,l2}
    D2 x2+y2u+λyzk+g2l(z,u) (b,b,1,1) cD1 b=min{k,l}
    D3 {x2+ut+λyzk+g2b+2(z,u)y2+p2b(x,z,u)+t (b+1,b,1,1,2b+1) cD1 kb+2
    D4 x2+y2u+yhk(z,u)+g2b+1(x,z,u) (b+1,b,1,1) cD1 kb+1
    D5 {x2+yt+g2b(z,u)yu+pb(z,u)+t (b,b1,1,1,b+1) cD1 zbp(z,u)

     | Show Table
    DownLoad: CSV
    Table 5.  Divisorial contractions to cD points with discrepancies greater than one.
    No. defining equations weight typea(X,E) condition
    D6 x2+y2u+zk+g2b+1(x,y,z,u) (b+1,b,a,1) cDa ak=2b+1
    D7 {x2+yt+g2b+2(y,z,u)yu+zk+pb+1(z,u)+t (b+1,b,a,1,b+2) cDa ak=b+1
    D8 {x2+ut+λzb+14+gb+1(y,z,u)y2+μzb14+pb1(x,z,u)+t (b+12,b12,4,1,b) cD4 b+14N,λ=1,μ=0, or b14N,μ=1,λ=0.
    D9 {x2+ut+zb+12+gb+1(y,z,u)y2+pb1(x,z,u)+t (b+12,b12,2,1,b) cD2
    D10 x2+y2u+zb+g2b(y,z,u) (b,b,2,1) cD2
    D11 x2+y2u+yp3(z,u)+u3+g6(z,u) (3,3,1,2) cD42 z3p(z,u)
    D12 x2+y2u+z3+yu2+g6(y,z,u) (3,4,2,1) cD43

     | Show Table
    DownLoad: CSV
    Table 6.  Divisorial contractions to cD/r points with discrepancy one.
    No. defining equations (r;ai)weight typea(X,E) condition
    D13 x2+y3+gk(y,z,u) (3;0,2,1,1)13(3,2,4,1) cD/31/3 k=2 and zu2 or z3g, or k=3 and z2ug
    D14 x2+y3+z3+g4(y,z,u) (3;0,2,1,1)13(6,5,4,1) cD/31/3
    D15 x2+yzu+g2(y,z,u) (2;1,1,1,0)12(3,1,1,2) cD/21/2
    D16 x2+yzu+g3(y,z,u) (2;1,1,1,0)12(3,b,c,d) cD/21/2 (b,c,d)=(3,1,2)(1,1,4)
    D17 {x2+yt+g3(z,u)zu+y3+t (2;1,1,1,0,1)12(3,1,1,2,5) cD/21/2
    D18 x2+y2u+λyzk+gl(z,u) (2;1,1,1,0)12(b,b2,1,4) cD/21/2 b=min{k2,l21}
    D19 x2+y2u+λyzk+gl(z,u) (2;1,1,1,0)12(b,b,1,2) cD/21/2 b=min{k,l}
    D20 {x2+ut+λyzk+gb+2(z,u)y2+pb(x,z,u)+t (2;1,1,1,0,0)12(b+2,b,1,2,2b+2) cD/21/2 kb+4
    D21 x2+y2u+yhk(z,u)+gb+1(x,z,u) (2;1,1,1,0)12(b+2,b,1,2) cD/21/2 kb+2
    D22 {x2+yt+g2b(z,u)yu+zb+t (2;1,1,1,0,1)12(b,b2,1,2,b+2) cD/21/2

     | Show Table
    DownLoad: CSV
    Table 7.  Divisorial contractions to cD/r points with large disprepancies.
    No. defining equations (r;ai)weight typea(X,E) condition
    D23 x2+y2u+zm+gb+1(x,y,z,u) (2;1,1,1,0)12(b+2,b,a,2) cD/2a/2 ma=2b+2,aandbareodd
    D24 {x2+yt+gb+2(z,u)yu+zm+pb2+1(z,u)+t (2;1,1,1,0,1)12(b+2,b,a,2,b+4) cD/2a/2 ma=b+2ab(mod2)
    D25 x2+y2u+z4b+g4b(y,z,u) (2;1,1,1,0)(2b,2b,1,1) cD/21
    D26 x2+yzu+y4+zb+uc (2;1,1,1,0)(2,1,2,1) cD/21 b,c4 b is even
    D27 {x2+ut+y4+z4yz+u2+t (2;1,1,1,0,0)(2,1,1,1,3) cD/21
    D28 {x2+ut+g2b+2(y,z,u)y2+p2b(x,z,u)+t (2;1,1,1,0,0)(b+1,b,1,1,2b+1) cD/21 Either b is odd, or b is even and xzb1 or z2bp
    D29 {x2+ut+g2b+2(y,z,u)y2+p2b(x,z,u)+t (2;1,1,1,0,0)(b+1,b,2,1,2b+1) cD/22 xzb12 or zbp

     | Show Table
    DownLoad: CSV
    Table 8.  Divisorial contractions to cE points with discrepancy one.
    No. defining equations weight typea(X,E) condition
    E1 x2+y3+g4(y,z,u) (2,2,1,1) cE61 2y2g(y,z,u)=0
    E2 x2+xp2(z,u)+y3+g5(y,z,u) (3,2,1,1) cE6,71
    E3 x2+y3+gq6(y,z,u) (3,2,2,1) cE1
    E4 x2+y3+y2p2(z,u)+gq8(y,z,u) (4,3,2,1) cE1
    E5 x2+xp4(y,z,u)+y3+g9(y,z,u) (5,3,2,1) cE1
    E6 x2+y3+y2p3(z,u)+g10(y,z,u) (5,4,2,1) cE7,81
    E7 x2+y3+g12(y,z,u) (6,4,3,1) cE1
    E8 x2+y3+y2p4(z,u)+g14(y,z,u) (7,5,3,1) cE7,81
    E9 x2+xp7(y,z,u)+y3+g15(y,z,u) (8,5,3,1) cE7,81
    E10 x2+y3+g18(y,z,u) (9,6,4,1) cE7,81
    E11 x2+y3+y2p6(z,u)+g20(y,z,u) (10,7,4,1) cE81
    E12 x2+y3+g24(y,z,u) (12,8,5,1) cE81
    E13 x2+y3+g30(y,z,u) (15,10,6,1) cE81

     | Show Table
    DownLoad: CSV
    Table 9.  Divisorial contractions to cE points with discrepancy one, continued.
    No. defining equations weight typea(X,E) condition
    E14 {x2+y3+tz+g6(y,z,u)p4(x,y,z,u)+t (3,2,1,1,5) cE6,71 p(x,y,z,u) is irreducible
    E15 x2+xp2(z,u)+y3+g6(x,y,z,u) (4,2,1,1) cE61
    E16 {x2+y3+tp2(z,u)+g6(y,z,u)q3(y,z,u)+t (3,2,1,1,4) cE71 q(y,z,u) is irreducible
    E17 x2+y3+yz3+g6(y,z,u) (3,3,1,1) cE71 y2u2g
    E18 {x2+yt+g10(y,z,u)y2+p6(y,z,u)+t (5,3,2,1,7) cE7,81 y2+p(y,z,u) is irreducible

     | Show Table
    DownLoad: CSV
    Table 10.  Divisorial contractions to cE points with large disprepancies.
    No. defining equations weight typea(X,E) condition
    E19 x2+(y+p2(z,u))3+yu3+g6(z,u) (3,3,2,1) cE62 zp(z,u)
    E20 {x2+yt+g10(y,z,u)y2+p6(z,u)+t (5,3,2,2,7) cE72 gcd(p6,g10)=1
    E21 x2+y3+u7+g14(z,u) (7,5,3,2) cE7,82 yz3,z5 or z4ug(z,u)

     | Show Table
    DownLoad: CSV
    Table 11.  Divisorial contractions to cE/2 points.
    No. defining equations (r;ai)weight typea(X,E)
    E22 x2+y3+g3(y,z,u) (2;1,0,1,1)12(3,2,3,1) cE/21/2
    E23 x2+y3+g5(y,z,u) (2;1,0,1,1)12(5,4,3,1) cE/21/2
    E24 x2+xp52(y,z,u)+y3+g6(y,z,u) (2;1,0,1,1)12(7,4,3,1) cE/21/2
    E25 x2+y3+g9(y,z,u) (2;1,0,1,1)12(9,6,5,1) cE/21/2
    E26 x2+y3+z4+u8+g8(y,z,u) (2;1,0,1,1)(4,3,2,1) cE/21

     | Show Table
    DownLoad: CSV

    For the reader's convenience, we put these tables in Section 4. In those tables, we use the following notation: For a non-negative integer m, the notation gm represents a function gOW such that w(g)=m. The notation pm represents a function pOW which is homogeneous of weight m with respect to the weight w.

    The reference of each of the cases in Table 2, Table 3, Table 5, ..., Table 7, Table 9, ..., Table 11 is as follows:

    ● Case A1 is [13, Theorem 1.2 (ⅰ)]. Case A2 is [15, Theorem 2.6].

    ● Case Ax1–Ax4 are [9 Theorems 7.4, 7.9, 8.4 and 8.8 respectively.

    ● Cases D6 and D7 are [13, Theorem 1.2 (ⅱ)]. Case D8–D11 is [15] Theorem 2.1–2.4. Case D12 is [15, Theorem 2.7].

    ● Case D13 is [9, Theorems 9.9, 9.14, 9.20]. Case D14 is [9, Theorem 9.25].

    ● Case D16 is [10, Proposition 4.4]. Case D17 is [10, Proposition 4.7, 4.12]. Case D18 is [10, Proposition 4.9]. Case D18 is [10, Proposition 5.4]. Case D19 is [10, Propositions 5.8, 5.13, 5.22, 5.28, and 5.35]. Case D20 is [10, Propositions 5.18 and 5.25]. Case D21 is [10, Propositions 5.16 and 5.32]. Case D22 is [10, Propositions 5.9 and 5.36].

    ● Cases D23 and D24 is [13, Theorem 1.2(ⅱ)] and [11, Theorem 1.1 (ⅲ)]. Case D25–D28 is [11 Theorem 1.1 (ⅰ), (ⅰ'), (ⅱ'), (ⅲ), (ⅱ) respectively. Case D29 is [14, Theorem 2].

    ● Case E19–E21 is [15] Theorems 2.5, 2.9 and 2.10 respectively.

    ● Case E22 is [9, Theorems 10.11, 10.17, 10.22, 10.28, 10.33 and 10.41]. Case E23 is [9, Theorems 10.33 and 10.47]. Case E24 is [9, Theorems 10.54 and 10.61]. Case E25 is [9], Theorem 10.67]. Case E26 is [11, Theorem 1.2].

    Divisorial contractions to cD points of discrepancy one (Case D1–D5 in Table 4) and divisorial contractions to cE points of discrepancy one (Case E1–E18 in Table 8) was completely classified by Hayakawa in his two unpublished papers "Divisorial contractions to cD points" and "Divisorial contractions to cE points". We will briefly introduce how to derive this classification. For more detail, please contact the author or Professor Takayuki Hayakawa in Kanazawa University.

    Let (oX) be a germ of three-dimensional Gorenstein terminal singular point with type cD or cE.

    Step1: Construct a divisorial contraction X1X which contracts an exceptional divisor of discrepancy one to o. We refer to [2, Section 4, Section 6] for the explicit construction. X1X can be viewed as a weighted blow-up with respect to an explicit embedding and a explicit weight.

    Step2: Find all exceptional divisors E over X such that a(E,X)=1 and CenterXE=o. We know that exc(X1X) is an exceptional divisor of discrepancy one. Assume that Eexc(X1X). Then an easy computation on discrepancies shows that a(E,X1)<1. In particular, CenterX1E is a non-Gorenstein point. Since X1X is an explicit weighted blow-up, all non-Gorenstein points on X1 can be explicitly computed, and all exceptional divisors of discrepancy less than one can be explicitly write down. Say S is the set of exceptional divisors over X1 with discrepancy less than one. One can compute a(E,X) for ES. If a(E,X)>1, then we remove E form S. After that, S is a set consisting exceptional divisors over X of discrepancy one.

    Step3: For any exceptional divisor ES, the valuation of E on X can be calculated. One can construct a weighted blow-up YEX with respect to this valuation. If YE do not have terminal singularities, then we remove E form S. Now, YEX for all ES, together with X1X, are all divisorial contractions to o with discrepancy one.

    Definition 2.12. Let YX be a divisorial contraction which contracts a divisor E to a point P. We say that YX is a w-morphism if a(X,E)=1rP, where rP is the Cartier index of KX near P.

    Definition 2.13. The depth of a terminal singularity PX, dep(PX), is the minimal length of the sequence

    XmXm1X1X0=X,

    such that Xm is Gorenstein and XiXi1 is a w-morphism for all 1im.

    The generalized depth of a terminal singularity PX, gdep(PX), is the minimal length of the sequence

    XnXn1X1X0=X,

    such that Xn is smooth and XiXi1 is a w-morphism for all 1in. The variety Xn is called a feasible resolution of PX.

    The Gorenstein depth of a terminal singularity PX, depGor(PX), is defined by gdep(PX)dep(PX).

    For a terminal threefold we can define

    dep(X)=Pdep(PX),
    gdep(X)=Pgdep(PX)

    and

    depGor(X)=PdepGor(PX).

    Remark 2.14. In the above definition, the existence of a sequence

    XmXm1X1X0=X,

    such that Xm is Gorenstein follows from [10, Theorem 1.2]. The existence of a sequence

    XnXn1X1X0=X,

    such that Xn is smooth follows from [2, Theorem 2].

    Definition 2.15. Assume that YX is a w-morphism such that gdep(Y)=gdep(X)1. Then we say that YX is a strict w-morphism.

    Lemma 2.16. Assume that YX is a divisorial contraction which is a weighted blow-up with the weight w(x1,...,xn)=1r(a1,...,an) with respect to an embedding XAn(x1,...,xn)/G where G is a cyclic group of index r. Assume that E is an exceptional divisor over X and vE(x1,...,xn)=1r(b1,...,bn). Then CenterYEUi non-trivially if and only if biaibjaj for all 1jn. Here, U1, ..., Un denotes the canonical affine chart of the weighted blow-up on Y.

    Proof. Let y1, ..., yn be the local coordinates of Ui. Then we have the following change of coordinates formula:

    xi=yairi,xj=yajriyj if ij.

    One can see that

    vE(yi)=biai,vE(yj)=birajbirai if ij.

    We know that CenterEY intersects Ui non-trivially if and only if birajbirai0 for all ji, or, equivalently, bjajbiai for all ji.

    Corollary 2.17. Assume that YX and Y1X are two different w-morphisms over the same point. Let E and F be the exceptional divisors of YX and Y1X, respectively. Then there exists uOX such that vE(u)<vF(u).

    Proof. Let XAn(x1,...,xn)/G be the embedding so that Y1X can be obtained by the weighted blow-up with respect to the embedding. We may assume that (xn=0) defines a Du Val section. Then

    vE(xn)=a(E,X)=a(F,X)=vF(xn).

    It follows that an=bn=1 where (a1,...,an) and (b1,...,bn) are integers in Lemma 2.16. Now, since a(F,X)=a(E,X), one has that CenterY1E is a non-Gorenstein point (if Y1 is generically Gorenstein along CenterY1E, then an easy computation shows that a(E,X)>a(F,X)). It follows that CenterY1EUn is empty since an=1 implies that Un is Gorenstein. Thus, by Lemma 2.16 we know that there exists j so that bjaj<1. Hence, vE(u)>vF(u) if u=xj.

    We have the following factorization of steps of three-dimensional terminal MMP by Chen and Hacon [3].

    Theorem 2.18 ([3, Theorem 3.3]). Assume that either XX is a flip over V, or XV is a divisorial contraction to a curve such that the exceptional locus contains a non-Gorenstein singular point. Then there exists a diagram

    such that Y1X is a w-morphism, YkX is a divisorial contraction, Y1Y2 is a flip or a flop, and YiYi+1 is a flip for i>1. If XV is divisorial, then YkX is a divisorial contraction to a curve and XV is a divisorial contraction to a point.

    Remark 2.19. Notation as in the above theorem. From the construction of the diagram, we can state the following:

    (1) Let CY1 be a flipping/flopping curve of Y1Y2. Then CX is a flipping curve of XX. Here, we use the notation introduced in Section 2.1, so CX is the image of CY1 on X.

    (2) Assume that the exceptional locus of XV contains a non-Gorenstein point P which is not a cA/r or a cAx/r point. Then Y1X can be chosen to be any w-morphism over P. This statement follows from the proof of [3, Theorem 3.1].

    We have the following properties of the depth [3, Propositions 2.15, 3.8 and 3.9]:

    Lemma 2.20. Let X be a terminal threefold.

    1. If YX is a divisorial contraction to a point, then dep(Y)dep(X)1.

    2. If YX is a divisorial contraction to a curve, then dep(Y)dep(X).

    3. If XX is a flip, then dep(X)>dep(X).

    We have the following negativity lemma for flips.

    Lemma 2.21. Assume that XX is a (KX+D)-flip. Then, for all exceptional divisors E, one has that a(E,X,D)a(E,X,DX). The inequality is strict if CenterXE is contained in the flipping locus.

    Proof. It is a special case of [16, Lemma 3.38].

    What we really need is the following corollary of the negativity lemma.

    Corollary 2.22. Assume that XX is a (KX+D)-flip and CX is an irreducible curve which is not a flipping curve. Then (KX+D).C(KX+DX).CX. The inequality is strict if C intersects the flipping locus non-trivially.

    Proof. Let XϕWϕX be a common resolution such that C is not contained in the indeterminacy locus of ϕ. Then Lemma 2.21 implies that F=ϕ(KX+D)ϕ(KX+DX) is an effective divisor and is supported on exactly those exceptional divisors whose centers on X are contained in the flipping locus. Hence,

    (KX+D).C(KX+DX).CX=(ϕ(KX+D)ϕ(KX+DX)).CW=F.CW0.

    The last inequality is strict if and only if CW intersects F non-trivially, or, equivalently, C intersects the flipping locus non-trivially.

    Let YX be a divisorial contraction between Q-factorial terminal threefolds that contracts a divisor E to a point. We construct the diagram

    as follows: Let Z1Y be a w-morphism and let H|KX| be a Du Val section. According to [3, Lemma 2.7 (ⅱ)], we have a(E,X,H)=0. We run the (KZ1+HZ1+ϵEZ1)-MMP over X for some ϵ>0 such that (Z1,HZ1+ϵEZ1) is klt. Notice that a general curve inside EZ1 intersects the pair negatively, and a general curve in F intersects the pair positively where F=exc(Z1Y). Thus, after finitely many (KZ1+HZ1+ϵEZ1)-flips Z1...Zk, the MMP ends with a divisorial contraction ZkY1 which contracts EZk, and Y1X is a divisorial contraction which contracts FY1.

    Lemma 3.1. Keeping the above notation, assume that KZ1 is anti-nef over X and EZ1 is not covered by KZ1-trivial curves. Then ZiZi+1 is a KZi-flip or flop for all i and ZkY1 is a KZk-divisorial contraction. In particular, Y1, Z2, ..., Zk are all terminal.

    Proof. Assume first that k=1. If Z1Y1 is a KZ1-negative contraction, then we are done. Otherwise, Z1Y1 is a KZ1-trivial contraction. In this case, EZ1 is covered by KX1-trivial curves, which contradicts our assumption.

    Now assume that k>1. We know that Z1Z2 is a KZ1-flip or flop. Also, notice that a general curve on EZ1 is KZ1-negative. Hence a general curve on EZ2 is KZ2-negative by Corollary 2.22. Now the relative effective cone NE(Z2/X) is a two-dimensional cone. One of the boundaries of NE(Z2/X) corresponds to the flipped/flopped curve of Z1Z2 and is KZ2-non-negative. Since there is a KZ2-negative curve, we know that the other boundary of NE(Z2/X) is KZ2-negative. Therefore, if k=2, then Z2Y1 is a KZ2-divisorial contraction, and for k>2, Z2Z3 is a KZ2-flip. One can prove the statement by repeating this argument k2 more times.

    We are going to find the sufficient conditions for the assumptions of Lemma 3.1. Our final results are Lemmas 3.4 and 3.6.

    Let XAn(x1,...,xn)/G=W/G be the embedding such that YX is a weighted blow-up with respect to the weight w and this embedding. First, we show that after replacing W by a larger affine space, if necessary, we may assume that Z1YX can be viewed as a sequence of weighted blow-ups with respect to the embedding XW/G.

    Let V be a suitable open set which contains P=CenterYF such that Z1Y can be viewed as a weighted blow-up with respect to an embedding VAn1(y1,...,yn1)/G. For j=1, ..., n1, we define DjV to be the Weil divisor corresponding to yj=0. Then Dj,X=ϕDj is a Weil divisor on X. Since X is Q-factorial, Dj,X corresponds to a G-semi-invariant function sjOW. We can consider the embedding

    XW/G(xn+jsj=0)j=1,...,n1An(x1,...,xn+n1)/G=ˉW/G.

    Then, YX is also a weighted blow-up with respect the weight ˉw which is defined by ˉw(xj)=w(xj) if j<n, and ˉw(xn+j)=w(sj). The embedding XˉW is exactly what we need.

    Now, let WW be the first weighted blow-up. We may assume P=CenterYF is the origin of UiW. Let y1, ..., yn be the local coordinate system of Ui that is mentioned in Corollary 2.2. We know that E|Ui=(yi=0). Let f4, ..., fn be the defining equation of XW/G. Then f4, ..., fn define Y|Ui where fi=(ϕ|1Ui)(fi). Since Y has terminal singularities, the weighted embedding dimension of Y|Ui near P is less than 4. For 5jn, we may write fj=ξjyj+fj(y1,...y4) for some ξj which does not vanish on P. One can always assume that HY=(y3=0) and so i3.

    Let fj=fj|yi=y3=0. Then f4, ..., fn defines HE near P. If fj is irreducible as a G-semi-invariant function, then we let ηj=fj. Otherwise, let ηj be a G-semi-invariant irreducible factor of fj.

    Lemma 3.2. Assume that YX can be viewed as a four-dimensional weighted blow-up. Then η4=...=ηn=0 defines an irreducible component of HYE.

    Proof. Since YX can be viewed as a four-dimensional weighted blow-up, we know that i4 and fj=yj+fj|y3=yi=0 for all j>4. Hence, we have ηj=fj. One can see that the projection

    (η4=...=ηn=0)|HYEP(a1,...,an)P(a1,...,a4)(η4=0)|HYE

    is an isomorphism. Since η4 is an irreducible function, it defines an irreducible curve.

    Notice that ηj is a polynomial in y1, ..., yn. There exists ηjOW such that ηj=(ϕ|1Ui)(ηj). We assume that YX is a weighted blow-up with the weight 1r(a1,...,an) and Z1Y is a weighted blow-up with the weight 1r(a1,...,an).

    Lemma 3.3. Let Γ=(η4=...=ηn=0) and assume that Γ is an irreducible and reduced curve. Then

    KZ1.ΓZ1=a23vE(η4)...vE(ηn)rn3ma1...an+aivF(η4)...vF(ηn)rn4a1...an.

    Here, m is the integer in Lemma 2.1 corresponding to the weighted blow-up YX.

    Proof. Since ΓE and KZ1+HZ1 is numerically trivial over X, we only need to show that

    HZ1.ΓZ1=a23vE(η4)...vE(ηn)rn3ma1...anaivF(η4)...vF(ηn)rn4a1...an. (3.1)

    We know that HZ1.ΓZ1=H.ΓvF(HY)F.ΓZ1. We need to show that the first term of (3.1) equals H.Γ and the second term of (3.1) equals vF(HY)F.ΓZ1.

    We have an embedding YWPW(a1,...,an). Let Dj be the divisor on W which corresponds to ηj. Then Γ=D4..Dn.E.H is a weighted complete intersection, so ΓZ1=D4,Z1..Dn,Z1.EZ1.HZ1. To compute H.Γ, we view Γ as a curve inside P(a1,...,an) which is defined by H=D4=...=Dn=0. It follows that

    H.Γ=a23vE(η4)...vE(ηn)rn3ma1...an.

    To compute F.ΓZ1, one writes

    F.ΓZ1=F.(ψD4vF(η4)F)..(ψDnvF(ηn)F).(ψEvF(E)F).(ϕHYvF(HY)F)=(1)n1vF(η4)...vF(ηn)vF(E)vF(HY)Fn.

    Since Z1Y is a w-morphism, the integer λ in Section 2.2 is 1. Hence, we know that Fn=(1)n1rn1a1...an. Now, vF(E)=air and vF(HY)=a(Y,F)=1r, so

    vF(HY)F.ΓZ1=aivF(η4)...vF(ηn)rn4a1...an.

    Lemma 3.4. Notation and assumption as in Lemma 3.3. Assume that:

    (ⅰ) For all 4jn, there exists an integer δj so that xkjδj appears in ηj as a monomial for some positive integer kj. Moreover, the integers δ4, ..., δn are all distinct.

    (ⅱ) If ji, 3, δ4, ..., δn, then a3ajaj.

    Then KZ1.ΓZ10.

    Proof. Fix j4. From the construction and our assumption we know that that i, δ4, ..., δn are all distinct. One can see that rvE(ηj)=kjaδj and rvF(ηj)kjaδj. Thus, we have a relation

    rvE(ηj)aδjrvF(ηj)aδj.

    One can always assume that if j>4, ji, then δj=j. By interchanging the order of y1, ..., y4, we may assume that δ4=4. Now, if i<4, then we may assume that i=1. If i>4, then we may assume that δi=1. We can write

    a23vE(η4)...vE(ηn)rn3ma1...an=1maia3a2rvE(η4)aδ4...rvE(ηn)aδn

    and

    aivF(η4)...vF(ηn)rn4a1...an=1r1a2rvF(η4)aδ4...rvF(ηn)aδn.

    Since mai=r, a3a21a2 and rvE(ηj)aδjrvF(ηj)aδj, we know that

    a23vE(η4)...vE(ηn)rn3ma1...anaivF(η4)...vF(ηn)rn4a1...an,

    So KZ1.ΓZ10.

    Remark 3.5.

    (1) From the construction we know that if j5, ji, then one can choose δj=j.

    (2) If YX can be viewed as a four-dimensional weighted blow-up, then condition (ⅰ) of Lemma 3.4 always holds. Indeed, in this case one has i4, so η4 is a two-variable irreducible function, hence there exists δ44 such that yk4δ4η4 for some positive integer k4. One also has δj=j for all j>4. Thus, condition (ⅰ) of Lemma 3.4 holds.

    (3) If for ji, 3, δ4, ..., δn one has that ajai, then condition (ⅱ) of Lemma 3.4 holds. Indeed, by Lemma 2.1 we know that UiAn/τ,τ where τ corresponds to the vector v=1ai(a1,...,ai1,r,ai+1,an). Let ˉv be the vector corresponding to the cyclic action near PUi. Then vmˉv (mod Zn) and r=mai. Since Z1Y is a w-morphism, and since HY is defined by y3=0, we know that a3=1 and ˉva31r(a1,...,an) (mod Zn). One can see that a3ajaj (mod r). This implies that a3ajaj since

    ajaimai=r.

    Lemma 3.6. Assume that KZ1 is anti-nef over X and there exists uOX such that vE(u)<a(E,X)a(F,X)vF(u). Then ZiZi+1 is a KZi-flip or flop for all 1ik1 and ZkY1 is a terminal divisorial contraction.

    In particular, if there exists ji such that a3aj>aj, then the conclusion of this lemma holds.

    Proof. We only need to show that EZ1 is not covered by KZ1-trivial curves. Then the conclusion follows from Lemma 3.1.

    Assume that EZ1 is covered by KZ1-trivial curves. Since KZ1 is anti-nef, those KZ1-trivial curves are contained in the boundary of the relative effective cone NE(Z1/X). Hence, k=1 and Z1Y1 is a KZ1-trivial divisorial contraction. Notice that if CYE is a curve which does not contain P, then KZ1.CZ1=KY.CY<0, hence the curve CZ1 is not contracted by Z1Y1. Thus, Z1Y1 is a divisorial contraction to the curve CY1. Notice that, in this case, a(E,Y1)=0.

    By computing the discrepancy, one can see that the pull-back of FY1 on Z1 is FZ1+a(E,X)a(F,X)EZ1. It follows that for all uOX, one has that

    vE(u)a(E,X)a(F,X)vF(u).

    Hence, if there exists u such that vE(u)<a(E,X)a(F,X)vF(u), then EZ1 is not covered by KZ1-trivial curves, so ZkY1 is a terminal divisorial contraction.

    Now, by Lemma 3.7, we know that

    a(E,X)a(F,X)=ra3r+a3ai.

    Consider u=xj. For ji we know that vE(xj)=ajr and

    vF(xj)=vF(yjyajri)=ajr+ajairr=raj+ajairr.

    The inequality vE(u)a(E,X)a(F,X)vF(u) becomes

    ajrra3r+a3airaj+ajairr=1ra3(raj+ajai)r+a3ai,

    or, equivalently,

    aj(r+a3ai)a3(raj+ajai).

    This is equivalent to

    aja3aj.

    Hence, the condition a3aj>aj implies that vE(u)<a(E,X)a(F,X)vF(u).

    Lemma 3.7. One has that

    a(E,X)=a3r,a(F,X)=r+a3airr.

    Proof. Since a(E,X,H)=0, we know that a(E,X)=vE(H)=a3r. Then

    a(F,X)=1r+a3rair=r+a3airr.

    Lemma 3.8. Note that the assumption in Lemma 3.4 depends only on the first weighted blow-up YX. In other words, we can check whether the assumption holds or not by simply considering the embedding which defines the weighted blow-up YX instead of considering the (possibly) larger embedding which defines both YX and Z1Y. Likewise, to apply Lemma 3.6, we can simply look at the embedding that defines YX, if condition a3aj>aj already holds under this embedding.

    Notation 3.9.

    (1) We say that the condition (Ξ) holds if conditions (ⅰ) and (ⅱ) in Lemma 3.4 hold for all possible choices of Γ. We say that the condition (Ξ) holds if conditions (2) and (3) in Remark 3.5 hold for all possible choices of Γ. As explained in Remark 3.5, we know that the condition (Ξ) implies the condition (Ξ).

    (2) We say that the condition (Ξ) (resp. (Ξ)) holds if the condition (Ξ) (reps. (Ξ)) holds and the inequality in Lemma 3.4 is strict for all possible choices of Γ. Using the notation in Lemma 3.4, it is equivalent to say that either there exists ji, 3, δ4, ..., δn such that a3aj>aj, or there exists j4 such that

    rvE(ηj)aδj>rvF(ηj)aδj.

    (3) We say that the condition (Θu) holds for some function u if vE(u)<a(E,X)a(F,X)vF(u). We say that the condition (Θj) holds for some index j if a3aj>aj. In either case, Lemma 3.6 can be applied.

    Notation 3.10. We say that a divisorial contraction YX is linked to another divisorial contraction Y1X if the diagram

    exists, where Z1Y is a strict w-morphism over a non-Gorenstein point, ZkY1 is a divisorial contraction, and ZiZi+1 is a flip or a flop for all 1ik1. We use the notation YXY1 if YX is linked to Y1X.

    Furthermore, if all ZiZi+1 are all flips, or k=1, then we say that Y is negatively linked to Y1, and use the notation YXY1.

    Remark 3.11. At this point, it is not clear why Z1Y should be a divisorial contraction to a non-Gorenstein point. In fact, from the classification of divisorial contractions between terminal threefolds (cf. Tables in Section 4), one can see that if there are two different divisorial contractions YX and Y1X, then Y or Y1 always contain a non-Gorenstein point. It is natural to construct the diagram starting with the most singular point, which is always a non-Gorenstein point.

    Remark 3.12.

    (1) If (Ξ) or (Ξ) holds and (Θu) or (Θj) holds for some function u or index j, then by Lemmas 3.4 and 3.6 one has that YXY1.

    (2) Assume that (Ξ) or (Ξ) holds and (Θu) or (Θj) holds for some function u or index j. Then one has that YXY1.

    Lemma 3.13. Assume that

    X(x1(x1+p(x2,...,x4))+g(x2,...,x4)=0)A4/G,

    such that

    (1) vE(g)=a1r+vE(p)=2a1r1.

    (2) i=1, a2+a4=a1 and a3=1.

    Then YXY1.

    Proof. We know that a1>aj for j=2, ..., 4, so (Ξ) holds. Consider the embedding

    X(x1x5+g(x2,...,x4)=x5x1p(x2,...,x4)=0)A5(x1,...,x5)/G.

    Then YX can be viewed as a weighted blow-up with the weight 1r(a1,...,a5) with respect to this embedding, where a5=rvE(p). The origin of U1 is a cyclic quotient point of type 1a1(r,a2,...,a5). The only w-morphism is the weighted blow-up that corresponds to the weight w(y2,...y4)=1a1(a2,...,a4). One can see that a5=rvE(g)>rvE(p)=a5, hence (Θ5) holds. Moreover, one can see that η5=y5p(y2,0,y4), so rvF(η5)=rvE(η5)=rvE(p(x2,0,x4)), hence

    rvE(η5)a5>rvF(η5)a5.

    Thus, YXY1.

    Lemma 3.14. Assume that YX and Y1X are two divisorial contractions such that YXY1. Let E and F be the exceptional divisors of YX and Y1X, respectively. Assume that there exists uOX such that vF(u)=1r and vF(u)>0 where r is the Cartier index of CenterXE and u is the strict transform of u on Y. Then a(F,X)<a(E,X) if a(E,X)>1.

    Proof. Notice that we have

    vF(u)=vF(˜u)+vE(u)vF(E).

    Since vE(u)1r and vF(˜u)>0, we know that vF(E)<1. Thus

    a(F,X)=1r+a(E,X)vF(E)1r+a(E,X)r1r=a(E,X)+1a(E,X)r

    where r is the Cartier index of CenterYF. Hence, a(F,X)<a(E,X) if a(E,X)>1.

    The aim of this section is to prove the following proposition:

    Proposition 4.1. Let X be a terminal threefold and YX, YX be two different divisorial contractions to points over X. Then there exists Y1, ..., Yk, Y1, ..., Yk such that

    Y=Y1XY2X...XYk=YkX...XY1=Y.

    Proof. We need a case-by-case discussion according to the type of the singularity on X. Please see Propositions 4.2, 4.3, 4.5–4.10 and 4.15–4.18.

    We keep the notation in Section 3.

    In this subsection, we assume that X has cA/r singularities. Divisorial contractions over X are listed in Table 2.

    Proposition 4.2.

    (1) If YX is of type A1 with a>1, then YXY1 for some Y1X which is of type A1 with the discrepancy less than a.

    (2) If YX is of type A1 with a=1 and b>r, then YXY1 where Y1 is an A1 type weighted blow-up with the weight 1r(br,c+r,1,r). One also has Y1XY if we begin with Y1X and interchange the role of x and y. Moreover, YXY1 if and only if η4=y.

    (3) If YX is of type A2, then YXY1 where Y1X is a divisorial contraction of type A1.

    Proof. Assume first that YX is of type A1. We are going to prove (1) and (2). If both b and c are less than r, then a=k=1. In this case, there is exactly one divisorial contraction of type A1, so there is nothing to prove. Thus, we may assume that one of b or c, say b>r.

    The origin of the chart UxY is a cyclic quotient point. On this chart, one can choose (y1,...,y4)=(x,u,z,y) with i=1 and δ4=4. One can see that (Ξ) holds. Now the two action in Lemma 2.1 is given by

    τ=1b(r,c,a,r),τ=1b(β,β(b+c)r,baβr,bβ).

    Since Ux is terminal, there exists a vector τ"=1b(bδ,ϵ,1,δ) such that τaτ" (mod Z4) and τλτ" (mod Z4) for some integer λ. There is exactly one w-morphism over the origin of Ux which extracts the exceptional divisor F so that vF corresponds to the vector τ". One can also see that

    ϵb=vF(y)=vF(g)rkb

    where g is the strict transform of g on Ux, since if zrpuqg, then ap+qak and

    vF(zrpuq)=1b(rp+δq)=1ab(rap+δaq)rkaab=rkb

    for δar because δar(mod b) and b>r. Thus,

    a4=c<rkaaϵ=a3a4,

    hence (Θ4) holds, and there exists Y1X such that YXY1.

    We need to check whether (Ξ) holds or not. We have that f4=η4=y4+g. One always has that

    rvF(η4)a4=bϵbϵ=1.

    Now, g=0 if and only if

    rvE(η4)a4=rcrc=1,

    and δa=r if and only if

    a3a2=aδ=r=a2.

    Thus, (Ξ) holds if and only if g0 or a do not divide r.

    One can compute the discrepancy of Y1X using Lemma 3.7. We know that ai=bδ and a3=a, so

    a(F,X)=r+a3airr=rδa+barbar=a(E,X).

    If a=1, then r=δ, so a(F,X)=a(E,X)=1r. One can verify that Y1X is the weighted blow-up with the weight 1r(br,c+r,1,r). In this case, YXY1 if and only if g0. Hence, YXY1 if and only if η4=y. This proves (2).

    Now assume that a>1. We already know that δar. If δa>r, then a(F,X)<a(E,X) and YXY1, so (1) holds. Hence, one only needs to show that δar. If δa=r, then

    b=aβ+λr=a(β+λδ)

    where λ=baβr. One can see that bβ=(a1)β+λaδ. On the other hand, since τλτ" (mod Zn), we know that bβλδ (mod b). Hence, b divides

    bβλδ=(a1)(β+λδ).

    This is impossible since (a1)(β+λδ) is a positive integer and is less than b.

    Finally, assume that YX of type A2. In this case, one needs to look at the chart UxY, and we choose (y1,...,y4)=(x,z,y,u) with i=1 and δ4=4. One can see that (Ξ) holds. The origin of the chart Ux is a cAx/4 point of the form

    (x2y2+z3+u2+g(x,y,z,u)=0)A4(x,y,z,u)/14(1,1,2,3).

    From [9, Theorem 7.9], we know that there are exactly two w-morphisms over this point which are weighted blow-ups with the weights w±(x±y,xy,z,u)=14(5,1,2,3). For both these two w-morphisms, one has that a2=2, a3=3 and a2=2, so (Θ2) holds and (Ξ) holds since a3a2>a2. Thus, there exists Y1X so that YXY1. One can compute that the discrepancy of Y1X is one, so Y1X is of type A1. This proves (3).

    In this subsection, we assume that X has cAx/r singularities with r=2 or 4. Divisorial contractions over X are listed in Table 3.

    Proposition 4.3. (1) Assume that YX is of type Ax1 or Ax3. Then YX is the only divisorial contraction over X.

    (2) Assume that YX is of type Ax2 or Ax4. Then there are exactly two divisorial contractions over X. Let Y1X be another divisorial contraction. Then Y1X has the same type of YX, and one has that YXY1XY.

    Proof. The number of divisorial contractions follows from [9, Section 7, 8]. So we can assume that YX is of type Ax2 or Ax4, and Lemma 3.13 implies that YXY1.

    In this subsection, we assume that X has cD singularities. At first, we consider w-morphisms over X, which are listed in Table 4. Notice that for types D1, D2 or D5 in Table 4 there is at most one divisorial contraction over X which is of the given type. This is because the equations of type D1, D2 and D5 come from the normal form of cD-type singularities, which are unique, and the blowing-up weights are determined by the defining equations.

    Lemma 4.4. Assume that there exists two different divisorial contractions with discrepancy one over X. Then, one of the following holds:

    (1) One of the divisorial contractions is of type D1.

    (2) The two morphisms are of type D2 and D5, respectively.

    (3) Both of the divisorial contractions are of type D3.

    (4) Both of the divisorial contractions are of type D4.

    Proof. Assume that YX and Y1X are the two given divisorial contractions. It is enough to prove the following statements:

    (ⅰ) If YX is of type D4, then Y1X is of type D1 or D4.

    (ⅱ) If YX is of type D2, then Y1X is of type D1 or D5.

    (ⅲ) If YX is of type D5, then Y1X is not of type D3.

    Let E and F be the exceptional divisors of YX and Y1X, respectively. Then a(F,Y)<1, since otherwise a(F,X)>1. Thus, P=CenterFY is a non-Gorenstein point.

    First, assume that YX is of type D4. In this case, P may be the origin of Ux or the origin of Uy, and they are both cyclic quotient points. Exceptional divisors over P with discrepancy less than one are described in [4, Proposition 3.1]. The origin of Ux is a 1b+1(b,1,1) point. If P is this point, then, since z=0 defines a Du Val section, we have that vF(z)=a(F,X)=1. One can verify that vF(u)=vF(z)=1 and vF(x)=vF(y)=b. Now, vF(x)=b only when xpb(z,u)g(x,z,u) for some homogeneous polynomial p(z,u) of degree b. One can check that vF(x+p(z,u))=b+1. In this case, Y1X is also of type D4 after a change of coordinates xxp(z,u). If P is the origin of Uy, then it is a 1b(1,1,1) point. One can verify that vF(u)>1. This implies that Y1X is of type D1.

    Now, assume YX is of type D2. Then P is the origin of UyY, which is a cA/b point. Exceptional divisors of discrepancy less than one over P are described in [4, Proposition 3.4]. One can verify that if λ0 and k=b, then vF(u)=1. In this case, Y1X is of type D5. Otherwise, vF(u)=2, and so Y1X is of type D1.

    Finally, assume that YX is of type D5. One can see that Y1X cannot have type D3 since zbp(z,u). This finishes the proof.

    Proposition 4.5. Assume that there exist two different divisorial contractions with discrepancy one over X, say YX and Y1X.

    (1) If YX and Y1X are both of type D4, then YXY1XY.

    (2) If YX is of type D3 and Y1X is of type D1, then YXY1. If Y1X is of type D3, then there exists another divisorial contraction Y2X which is of type D1 so that YXY2XY1.

    (3) If YX is of type D2 and Y1X is of type D5, then YXY1.

    (4) If YX is of type D1 and Y1X is not of type D3, then YXY1.

    Proof. Assume first that YX and Y1X are both of type D4. Notice that, in this case, xpb(z,u)g(x,z,u). Thus, YXY1 by Lemma 3.13.

    Now assume that YX is of type D3. Consider the chart UtY which is defined by

    (x2+u+λyzktb+k2b2+g(z,u,t)=y2p(x,z,u)+t=0)
    A5(x,y,z,u,t)/12b+1(b+1,b,1,1,1).

    Notice that, using the notation in Section 3, we know that

    f4=η4=x2+u+g(0,u,0),f5=y2+p(x,0,u).

    η5 can be y±μub if p(x,0,u)=μ2u2b for some μC, and otherwise η5=f5. One can see that η4=η5=0 defines an irreducible and reduced curve. There is only one w-morphism over the origin of Ut which is defined by weighted blowing up the weight w(x,y,z,u,t)=12b+1(b+1,b,1,2b+2,2b). Now, in this case we choose (y1,...,y5)=(x,y,z,u,t) with i=5, δ4=4, δ5=2. One can see that (Ξ) holds and (Θ4) holds. Also, one has that rvE(η4)a4=2b+2 while rvF(η4)a4=1. Thus, (Ξ) holds and so there exists Y2X such that YXY2. One can compute that Y2X is of type D1. If Y1X is of type D1, then Y2=Y1 since there is at most one divisorial contraction with type D1. This proves statement (2).

    Now assume that YX is of type D2 and Y1X is of type D5. In this case, we consider the embedding corresponding to Y1X. Under this embedding, YX is given by the weighted blow-up with the weight (b,b,1,1,b) and the chart UyY is given by

    Uy=(x2t+g(y,z,u)=yu+zb+t=0)A5(x,y,z,u,t)/1b(0,1,1,1,0).

    We take (y1,...,y5)=(y,u,z,x,t) with δ4=4 and δ5=5. Then (Ξ) holds. The origin of Uy is a cA/b point and the weight w(y1,...,y5)=1b(b1,1,1,b,2b) defines a w-morphism over Uy. One can see that (Θ5) holds. Moreover, since a5=2b>b=a5, we know that (Ξ) holds. Thus, YXY1.

    Finally, assume that YX is of type D1 and Y1X is not of type D3. Let b and b1 be the integers in Table 4 corresponding to YX and Y1X, respectively. First, we claim that bb1. Indeed, if Y1X is of type D5, then zb1h(z,u), which implies that b1b+1. If Y1X is of type D2 or D4, then the inequality bb1 follows from Corollary 2.17. Now, the origin of the chart UuY is defined by

    (x2+y2+λyzkukb1+g(z,u)=0)A4(x,y,z,u)/12(b,b1,1,1),

    which is a cAx/2 point. We can take (y1,...,y4)=(u,y,z,x) with i=1 and δ4=4. w-morphisms over this point are fully described in [9, Section 8]. Since bb1, we know that 2kb1>b and the multiplicity of g(z,u) is greater than or equal to 2b. Hence, if F is the exceptional divisor of a w-morphism over Y, then

    vF(y)b>b1=vE(y).

    Thus, (Ξ) and (Θ2) holds and one has YXY1.

    Now we study divisorial contractions of discrepancy greater than one. All such divisorial contractions are listed in Table 5.

    Proposition 4.6. Assume that YX is a divisorial contraction with discrepancy a>1. Then there exists a divisorial contraction Y1X such that YXY1, and a(F,X)<a where F=exc(Y1X).

    Proof. First, notice that (Θu) holds in cases D6–D10 or D12, and (Θz) holds in case D11. This is because vE(u) or vE(z)=1 in those cases and a(E,X)a(F,X)=a>1.

    Now we list all cases in Table 5, write down the chart on Y we are looking at, and write down the variables y1, ..., yn. One can easily see that (Ξ) holds in all cases.

    (1) Assume that YX is of type D6. Consider the chart UxY and take (y1,...,y4)=(x,y,z,u) with δ4=4.

    (2) Assume that YX is of type D7. Consider the chart UtY and take (y1,...,y5)=(y,u,z,x,t), δ4=4 and δ5=1 or 2.

    (3) Assume that YX is of type D8 or D9. We consider the chart UtY and take (y1,...,y5)=(x,y,z,u,t) with δ4=4 and δ5=2.

    (4) Assume that YX is of type D10. We consider the chart UyY and take (y1,...,y4)=(y,u,z,x) with δ4=4.

    (5) Assume that YX is of type D11. We consider the chart UyY and (y1,...,y4)=(y,z,u+λy,x) for some λC with δ4=4.

    (6) Assume that YX is of type D12. We consider the chart UyY and (y1,...,y4)=(y,z,x+λy,u) for some λC with δ4=4.

    Now we know that there exists Y1X so that YXY1. Then Y1X is of one of types in Table 4 or Table 5. One can see that vF(z)=1 if Y1X is of types D1–D5, D7–D9 or D11 and vF(u)=1 if Y1X is of type D6, D10 or D12. Since CenterYF is the origin of the chart Ux, Uy or Ut, one can apply Lemma 3.14 to say that a(F,X)<a. This finishes the proof.

    In this subsection, we assume that X has cD/r singularities with r=2 or 3. We first study w-morphisms over X.

    Proposition 4.7. Assume that X has cD/3 singularities.

    (1) If YX is of type D14 or if YX is of type D13 and both zu2 and z2ug(y,z,u), then there is only one w-morphism over X.

    (2) If YX is of type D13 and zu2 or z2ug(y,z,u), then there are two or three w-morphisms over X. Say Y1X, ..., YkX are other w-morphisms with k=1 or 2. Then YXYiXY for all 1ik.

    Proof. The statement about the number of w-morphisms follows from [9, Section 9]. Now we may assume that YX is of type D13 and zu2 or z2ug(y,z,u). The chart UzY is defined by

    (x2+y3+g(y,z,u)=0)A4(x,y,z,u)/14(3,2,1,1)

    with u2 or zug(y,z,u). We can take (y1,...,y4)=(z,x,u+λz,y) for some λC with δ4=4. Now the w-morphism over Uz is a weighted blow-up with the weight w(y1,...y4)=14(3,5,1,2). One can see that (Θ2) and (Ξ) hold. Hence, we can get a divisorial contraction Y1X such that YXY1. One can compute that Y1X is also a w-morphism.

    If there are three w-morphisms over X, then the defining equation of X is of the form x2+y3+zu(z+u) as in [9, Section 9.A], so g(y,z,u)=u(z+u). One can make a change of coordinates uuz and again consider the weighted blow-up with the same weight 14(3,2,1,5). In this way, we can get a divisorial contraction Y2X which is different to Y1, and we also have that YXY2. This finishes the proof.

    Proposition 4.8. Assume that X has cD/2 singularities and YX is of type D15, D16 or D17.

    (1) If YX is of type D15, then there is only one w-morphism over X.

    (2) If YX is of type D17, then there exists exactly two w-morphisms over X. The other one, Y1X, is of type D16, and one has that YXY1.

    (3) If YX is of type D16 and there are no w-morphisms over X with type D17, then there are exactly three w-morphisms over X. They are all of type D16 and are negatively linked to each other.

    Proof. The statement about the number of w-morphisms follows from [10, Section 4]. Assume that YX is of type D17. Consider the chart UtY with (y1,...,y5)=(y,u,y+z,x,t) with δ4=4 and δ5=1. One can see that (Ξ) holds. Now the origin of Ut is a cyclic quotient point. Let F be the exceptional divisor of the w-morphism over Ut. Then one has that vF(y1,...y5)=15(6,2,1,3,3). One can see that (Θ1) holds.

    Assume that YX is of type D16 and there are no w-morphisms of type D17 over X. By [10, Section 4], we know that neither y4 nor z4g(y,z,u). Assume first that (b,c,d)=(1,1,4). Consider the chart UuY which has a cAx/4 singular point at the origin. We choose (y1,...,y4)=(y,u,y+z,x) with δ4=4. One can see that (Ξ) holds. Let w be the weight on Uu so that w(y1,...,y4)=14(5,2,1,3). Then the weighted blow-up with weight w gives a w-morphism. It follows that (Θ1) holds and also (Ξ) holds since a1=5>3=a1. Hence, there exists a w-morphism Y1X so that YXY1. If we interchange the roles of y and z, we can get another w-morphism Y2X with YXY2.

    Now assume that (b,c,d)=(3,1,2). Consider the chart UyY which is defined by

    (x2+zu+g(y,z,u)=0)A4(x,y,z,u)/13(0,1,1,2).

    Taking (y1,...,y4)=(y,u,y+z,x), then (Ξ) holds. Let w be the weight w(y1,...,y4)=13(1,5,1,3). Then the weighted blow-up with the weight w gives a w-morphism over Uy and (Θ2) and (Ξ) holds. If we take w to be another weight w(x,y,z,u)=13(3,1,4,2), then we get another w-morphism over Uy and (Θz) holds. Thus, we can get two different w-morphisms over X and Y is negatively linked to both of them.

    Proposition 4.9. Assume that X has cD/2 singularities, YX is of type D18–D22, and assume that there are two w-morphisms YX and Y1X.

    (1) Assume that YX is of type D18 and Y1X is not of type D20. Then YXY1.

    (2) Assume that both YX and Y1X are not of type D18. Then, one of the following holds:

    (2–1) YX is of type D19 and Y1X is of type D22. One has that YXY1.

    (2–2) Both YX and Y1X are of type D21 and YXY1XY.

    (2–3) Both YX and Y1X are of type D20 and there exists another w-morphism Y2X which is of type D18 so that YXY2XY1.

    Proof. The computation is similar to the proof of Proposition 4.5 after replacing the types D1–D5 by D18–D22, so we will omit the proof. Notice that an analog result of Lemma 4.4 can be proved by a similar computation, or can be directly followed by [10, Section 5].

    Now we consider non-w-morphisms over X. Notice that there is no divisorial contraction with discrepancy greater than 13 over cD/3 points. Divisorial contractions of discrepancy greater than 12 over cD/2 points are listed in Table 7.

    Proposition 4.10. Assume that r=2 and YX is a divisorial contraction with the discrepancy a2>1. Then there exists a divisorial contraction Y1X such that YXY1, and a(F,X)<a2 where F=exc(Y1X).

    Proof. First, assume that YX is of type D23. Consider the chart UxY which is defined by

    (x+y2u+zm+g(x,y,z,u)=0)A4(x,y,z,u)/1b+2(2,b,a,2).

    We take (y1,...,y4)=(x,y,z,u) with δ4=2 or 4. One can see that (Ξ) holds. Since a4=2 and a3, we know that (Θ4) holds. Thus, there exists Y1X such that YXY1. The origin of Ux is a cyclic quotient point. Let F be the exceptional divisor of the w-morphism over this point. Then vF(x)mb+2. It follows that

    a(F,X)=1b+2+a2vF(x)2+ma2b+4=1<a2.

    Assume that YX is of type D24. Consider the chart UtY which is defined by

    (x2+y+g(z,u,t)=yu+zm+p(z,u)+t=0)A4(x,y,z,u,t)/1b+4(c1,...,c5)

    where (c1,...,c5)=(b+2,b,a,2,2) if a, b are odd, and (c1,...,c5)=(b+3,b+2,a2,1,1) if a, b are even. We take (y1,...,y5)=(y,u,z,x,t) with δ4=4 and δ5=1 or 2. Then (Ξ) holds. Now the origin of Ut is a cyclic quotient point. Let F be the exceptional divisor over this point. Then vF(y1,...,y5)=1b+4(a1,...,a5) with a2+a4=b+4. It follows that a(a2+a4)>b+4=a2+a4, hence (Θj) holds for j=2 or 4. Thus, there exists Y1X so that YXY1. One has that

    a(F,X)=1b+4+a2vF(x)2+ma2b+812.

    Hence, Y1X is a w-morphism.

    Assume that YX is of type D25. The chart UyY is given by

    (x2+yu+z4b+g(y,z,u)=0)A4(x,y,z,u)/14b(0,2b1,1,2b+1).

    We take (y1,...,y4)=(y,u,z,x) with δ4=4. One can see that (Ξ) holds. The origin of Uy is a cA/4b point and there is only one w-morphism over this point. Let F be the exceptional divisor of the w-morphism. Then vF(y1,...,y4)=14b(2b1,2b+1,1,4b). Hence, (Θ2) holds. One can also compute that a(F,X)=12, hence there exists a w-morphism Y1X such that YXY1.

    Assume that YX is of type D26. The chart UzY is a cA/4 point given by

    (x2+yu+y4+z2b4+uczc4=0)A4(x,y,z,u)/14(0,1,1,3).

    We take (y1,...,y4)=(y,u,y+z,x) with δ4=4. One can see that (Ξ) holds. Now let w be the weight such that w(y1,...,y4)=14(1,3,1,4) if b=4 and w(y1,...,y4)=14(1,7,1,4) if b6. Hence, (Θ2) holds, and there exists Y1X such that YXY1. One can compute that a(F,X)=12.

    Assume that YX is of type D27. Consider the chart UtY which is defined by

    (x2+u+y4+z4=yz+u2+t=0)A5(x,y,z,u,t)/16(5,1,1,4,2).

    We take (y1,...,y5)=(u,y,y+z,x,t) with δ4=4 and δ5=1. In this case, (Ξ) holds. Now the origin of Ut is a cyclic quotient point. Let F be the exceptional divisor over this point which corresponds to a w-morphism. Then vF(y1,...,y5)=16(5,1,1,4,2). One can see that (Θ1) holds. Thus, there exists Y1X which extracts F so that YXY1. One can compute that a(F,X)=12.

    Finally, assume that YX is of type D28 or D29. The chart UtY is defined by

    (x2+u+g(y,z,u,t)=y2+p(x,z,u)+t=0)A4(x,y,z,u,t)/14b+2(1,1,a2b1,2,2),

    where a=2 in case D28 and a=4 in case D29. We take (y1,...,y5)=(x,y,z,u,t) with δ4=4 and δ5=2. Then (Ξ) holds. The origin of Ut is a cyclic quotient point. Let F be the exceptional divisor corresponding to the w-morphism over this point. Then vF(y1,...,y5)=14b+2(a1,...,a5) with a1+a2=4b+2, a2(2b+1a)1 (mod 4b+2) and a3=1. From the defining equation one can see that a4>1, hence (Θ4) holds. Thus, there exists a divisorial contraction Y1X which extracts F so that YXY1. We only need to show that a(F,X)<a2.

    Assume that YX is of type D28. In this case, a2 is the integer such that a2(2b1)1 (mod 4b+2). If b is odd, then a2=b since

    b(2b1)=2b2b=(4b+2)b12+1.

    One can see that a52b. If b is even, then a2=3b+1 since

    (3b+1)(2b1)=6b2b1=(4b+2)(32b1)+1.

    Hence, a1=b+1. Now, since xzb1 or z2bp(x,z,u), we also have that a52b. In either case we have

    a(F,X)=14b+2+a54b+22b+14b+2=12<1=a2.

    Finally, assume that YX is of type D29. We want to show that a5<4b+2. Then

    a(F,X)=14b+2+2a54b+214b+2+8b+24b+2<2=a2

    and we can finish the proof. If zbp(x,z,u), then a5b. If a2<2b+1, then a54b. Assume that zbp(x,z,u) and a22b+1. Then a12b+1 and xzb12p(x,z,u). Hence,

    a52b+1+b12<4b+2.

    In this subsection, we assume that X has cE singularities. First, we study w-morphisms over cE points. All w-morphisms over cE type points are listed in Tables 8 and 9.

    We assume that there exist two different w-morphisms over X, say YX and Y1X. Let F=exc(Y1X). Let P=CenterYF. One always has that a(F,Y)<1, so P is a non-Gorenstein point.

    Lemma 4.11. Assume that both YX and Y1X are of type E1–E13. Then YX is not of type E1 or E6.

    Proof. Assume that YX is of type E1. Then the only non-Gorenstein point on Y is the origin of

    Uy=(x2+y2+g(z,u)=0)A4(x,y,z,u)/12(0,1,1,1).

    This is a cAx/2 point. The exceptional divisor G of discrepancy less than one over this point is given by the weighted blow-up with the weight w(x,y,z,u)=12(2,3,1,1). One can compute that a(G,X)=2, hence there is only one w-morphism over X. Thus, YX is not of type E1.

    Assume that YX is of type E6. If X has cE8 singularities, then there is only one non-Gorenstein point on Y, namely the origin of Uy. If X has cE7 singularities, then the origin of Uz is also a non-Gorenstein point. Assume first that P is the origin of Uz. Then P is a cyclic quotient point of index two and there is only one exceptional divisor over P with discrepancy less than one. Hence, F should correspond to this exceptional divisor. One can compute that vF(x,y,z,u)=(3,3,1,1), so Y1X should be of type E17. Nevertheless, in this case one can see that vF(σ)vE(σ) for all σOX. This contradicts Corollary 2.17. Hence, P can not be the origin of Uz.

    We want to show that P is also not the origin of Uy. The chart Uy is defined by

    (x2+y(y+p(z,u))+g(y,z,u)=0)A4(x,y,z,u)/14(1,3,2,1).

    The origin of Uy is a cAx/4 point. Since (u=0) defines a Du Val section, we know that vF(u)=vF(y)+vF(u)=1. Hence, both vF(y) and vF(u)<1. This means that vF(y)3. Assume that YX and Y1X correspond to the same embedding XA4. Then, since vF(y)3, we know that Y1X is of type E1–E5. However, in those cases one always has that vF(σ)vE(σ) for all σOX. This contradicts Corollary 2.17. Thus, Y1X corresponds to a different embedding.

    Let ZY be a w-morphism over the origin of Uy. From the classification we know that ZY is a weighted blow-up with the weight w(x,y,z,u)=14(5,k,2,1) for k=3 or 7. One can compute that non-Gorenstein points on Z over Uy are cyclic quotient points. Let ˉZZ be an economic resolution over those cyclic quotient points. Then F appears on ˉZ since a(F,Y)<1. Moreover, ˉZX can be viewed as a sequence of weighted blow-ups with respect to the embedding XA4(x,y,z,u). We write (x1,...,x4)=(x,y,z,u) and let XA4(x1,...,x4) be the embedding corresponding to the weighted blow-up Y1X. One can always assume that x4=x4=u since vE(u)=1. We write xj=xj+qj. Since YX and Y1X correspond to different embeddings, there exists j<4 such that qj0 and vF(xj)>vF(xj)=vF(q). Since ˉZX can be viewed as a sequence of weighted blow-ups with respect to the embedding XA4(x,y,z,u), we know that the defining equation of ˉZ is of the form xj+qj+ˉh such that vF(xj)=vF(ˉh). Hence, there is exactly one j such that qj0, and the defining equation of X is of the form ξ(xj+qj)+h. One can see that either xj=z, or xj=y and qj=p.

    Now, if xj=z, then x1=x1=x and x2=x2=y. One can see that vF(x2)=vF(y)3. So, Y1X is of type E1–E5. In those cases, vF(xj)2, so vF(xj)=vF(qj)=1 and vF(xj)=2. Hence, Y1X is of type E3–E5 and vF(y)=vF(x2)2. Also, since vF(qj)=1, qj=λu for some λC. Therefore vF(z)=vF(xjqj)=1. But, then

    vF(z)vE(z)=12vF(y)vE(y).

    By Lemma 2.16, CenterYF can not be the origin of Uy. This leads to a contradiction.

    Finally, we assume that xj=y and qj=p. Notice that p=λ1zu+λ2u3, hence vF(p)2. If vF(z)=vF(x3)=1, then Y1X is of type E1 or E2, and so vF(xj)=2. However, we know that vF(p)2. This contradicts the assumption that vF(xj)>vF(qj)=vF(p). Hence, vF(z)2 and so vF(p)3. Since

    vF(y)=vF(xj)=vF(qj)=vF(p)3

    and vF(y)3 by the previous discussion, we know that vF(y)=3. Recall that we write

    Uy=(x2+y(y+p(z,u))+g(y,z,u)=0)A4(x,y,z,u)/14(1,3,2,1).

    Since vF(y)=3, vF(E)=vF(y)=34. This means that a(F,Y)=14, so F corresponds to a w-morphism over Uy. Nevertheless, as we mentioned before, w-morphisms over Uy can be obtained by a weighted blow-up with respect to the above embedding, hence Y1X and YX correspond to the same four-dimensional embedding, leading to a contradiction.

    Lemma 4.12. Assume that both YX and Y1X are of type E1–E13. If P is the origin of UxY, then YX is of type E2, E5 or E9, and Y1X has the same type. One has that YXY1XY.

    Proof. This assumption implies that the origin of Ux is contained in Y, so YX is of type E2, E5 or E9 and P is a cyclic quotient point. If YX is a weight blow-up with the weight (b,c,d,1), then b=c+d and

    Ux=(x+p(z,u)+g(x,y,z,u)=0)A4(x,y,z,u)/1b(1,c,d,1).

    Since a(F,Y)<1, F is the valuation described by [4, Proposition 3.1]. Hence, vF(y,z,u)=1b(c,d,a) with c+d=b and ab=a(F,Y). Since a(F,X)=a(F,Y)+vF(x)=1, we know that vF(x)=1ab. One can compute that

    vF(x,y,z,u)=(ba,c(acc)b,daddb,1).

    Since c<b and acc(mod b), we know that accb0, so vF(y)c=vE(y). Likewise, we know that vF(z)d=vE(z). One also has that vF(x)<vE(x) and vF(u)=vE(u).

    On the other hand, Corollary 2.17 says that there exists σOX such that vF(σ)>vE(σ). This can only happen when

    vF(x)=vF(p(z,u))<vF(x+p(z,u))=vF(g(x,y,z,u))

    and in this case one can choose σ=x+p(z,u)OX. Now, Y1X can be obtained by a weighted blow-up with respect to the embedding

    X(σ2σp(z,u)+y3+g(y,z,u)=0)A4(σ,y,z,u)

    and with the weight w1(σ,y,z,u)=(b1,c1,d1,1) where c1=c(acc)b and d1=daddb. Since c1c, d1d and b1=vF(v)>vE(v), by Lemma 2.16 we know that CenterY1E is the origin of U1,σ.

    Now, if we interchange Y and Y1, then the above argument yields that cc1 and dd1. Hence, c=c1 and d=d1 and so YX and Y1X are of the same type. One has that a=1 and c=c, d=d. Thus, F is the exceptional divisor of the w-morphism over P. Now we know that YXY1 by Lemma 3.13 and also Y1XY by the symmetry.

    Lemma 4.13. Assume that both YX and Y1X are of type E1–E13. Then P is not the origin of UyY.

    Proof. By Lemma 4.11, we know that YX is not of type E1 or E6, hence YX is of type E4, E8 or E11 and the origin of Uy is a cyclic quotient point. We assume that YX is a weighted blow-up with the weight (b,c,d,1). Following the same computation as in the proof of Lemma 4.12, we may write Y1X as a weighted blow-up with respect to the embedding

    X(x2+(σp(z,u))2σ+g(σ,z,u)=0)A4(x,σ,z,u)

    and with the weight (b1,c1,d1,1), such that CenterY1E is the origin of U1,σY1. Nevertheless, in this case one always has that b1<b since b>c. The symmetry between Y and Y1 yields that b>b1>b, which is impossible.

    Lemma 4.14. Assume that both YX and Y1X are of type E1–E13. Let XA4(x1,...,x4) be the embedding that corresponds to YX in Table 8. Then P is the origin of Ui for some i4.

    Proof. Assume that P is not the origin of Ui for all i4. Then Y has a non-Gorenstein point on UiUj for some ij. In this case, YX is of type E7 or E10–E13. For simplicity we assume that i=1 and j=2. If YX is a weighted blow-up with the weight (a1,...,a4), then we have the following observation:

    (1) a1=dk1 and a2=dk2 for some integers k1, k2 and d. We may assume that k2=2 and k1 is odd.

    (2) xk21 and xk12 appear in f where f is the defining equation of X. Moreover vE(xk21)=vE(xk12)=vE(f).

    (3) P is a cyclic quotient point of index d. On U1, the local coordinate system is given by (x1,x3,x4), where xl is the strict transform of xl on U1.

    Since F is a valuation of discrepancy less than 1 over P, we know that vF(x1,x3,x4)=1d(a1,a3,a4) with al<d for l=1, 3, and 4. One can compute that

    vF(x1,...,x4)=(k1a1,k2a1,1d(a3+a3a1),1d(a4+a4a1)),

    and Y1X can be obtained by the weighted blow-up with respect to the same embedding XA4(x1,...,x4) and with the weight vF. Nevertheless, one can easily see that vF(xl)vE(xl) for all 1l4. This contradicts Corollary 2.17.

    Proposition 4.15. Assume that both YX and Y1X are both of type E1–E13. Then YX is of type E2, E5 or E9, and Y1X has the same type. One has that YXY1XY.

    Proof. Let

    X(f(x,y,z,u)=0)A4(x,y,z,u)

    be the embedding corresponding to YX, and

    X(f1(x1,y1,z1,u1)=0)A4(x1,y1,z1,u1)

    be the embedding corresponding to Y1X. If CenterYF is the origin of UxY or CenterY1E is the origin of Ux1Y1, then the statement follows from Lemma 4.12. We do not consider these cases here. Then, Lemmas 16 and 17 imply that CenterYF=UzY and CenterY1E=Uz1Y1.

    We may assume that vE(f)vF(f1). Since vE(u)=vF(u1)=1, one can always assume that u=u1 and (u=0) defines a Du Val section. Lemma 2.16 implies that vE(z)>vE(z1) and vF(z)<vF(z1), hence zz1. We may write z1=z+h. If vE(h)vE(z), then we may replace z by z+h, which will lead to a contradiction. Hence, vE(h)<vE(z). Thus, h=λuk for some k<vE(z). Since (u=0) defines a Du Val section, we know that z4, yz3 or z5f. It follows that u4k, yu3k or u5k appear in either f or f1. This means that vE(f)4k, vE(y)+3k or 5k for some k<vE(z). One can easily check that for all the cases in Table 8 this inequality never holds. Thus, we get a contradiction.

    Proposition 4.16. Assume that YX is of type E14–E18. Then there exists Y1X which is of type E3 or E6 such that YXY1. Moreover, if YX is of type E15 or E17, then YXY1.

    Proof. Assume that YX is of type E14. Consider the chart

    Ut=(x2+y3+z+g(y,z,u,t)=p(x,y,z,u)+t=0)A5(x,y,z,u,t)/15(3,2,1,1,4).

    We choose (y1,...,y5)=(x,y,u,z,t) with δ4=1 and δ5=2 or 4, or δ4=2 and δ5=1 or 4. Then (Ξ) holds. Now, let F be the exceptional divisor that corresponds to the w-morphism over the origin of Ut. Then

    vF(y1,...,y5)=15(3,2,1,6,4).

    One can see that (Θ4) holds. Thus, there exists a divisorial contraction Y1X so that YXY1 which extracts F. One can compute that vF(x,y,z,u)=(3,2,2,1), so Y1X is of type E3.

    Assume that YX is of type E15. Consider the chart

    Ux=(x2+p(z,u)+y3+g(x,y,z,u)=0)A4(x,y,z,u)/14(3,2,1,1).

    We choose (y1,...,y4)=(x,z,u,y) with δ4=4. Then (Ξ) holds. Now the origin of Ux is a cAx/4 point. After a suitable change of coordinates, we may assume that u2p(z,u). Then the w-morphism over this point can be given by a weighted blow-up with the weight vF(y1,...,y4)=14(3,5,1,2). One can see that (Θ2) and (Ξ) hold. Hence, there exists a divisorial contraction Y1X such that YXY1. One also has vF(x,y,z,u)=(3,2,2,1), so Y1X is of type E3.

    Assume that YX is of type E16. Consider the chart

    Ut=(x2+y3+p(z,u)+g(y,z,u,t)=q(y,z,u)+t=0)A5(x,y,z,u,t)/14(3,2,1,1,3).

    We choose (y1,...,y5)=(y,z,u,x,t) with δ4=4 and δ5=1 or 2. Then (Ξ) holds. Now the origin of Ut is a cAx/4 point. After a suitable change of coordinates, we may assume that u2p(z,u). Then the weight vF(y1,...,y5)=14(2,5,1,3,3) defines a w-morphism over Ut. One can see that (Θ2) holds. Hence, there exists a divisorial contraction Y1X such that YXY1. One can compute that vF(x,y,z,u)=(3,2,2,1), so again Y1X is of type E3.

    Assume that YX is of type E17. Consider the chart

    Uy=(x2+y3+z3+g(y,z,u)=0)A4(x,y,z,u)/13(0,2,1,1).

    We can choose (y1,...,y4)=(y,z,u,x) with δ4=4. Then (Ξ) holds. The origin of Uy is a cD/3 point. Notice that y2u2g(y,z,u), so the w-morphism over Uy is given by the weighted blow-up with the weight vF(y1,...,y4)=13(2,4,1,3). One can see that (Θ2) and (Ξ) hold. One has that vF(x,y,z,u)=(3,2,1,1), so there exists Y1X which is of type E3 such that YXY1.

    Finally, assume that YX is of type E18. Consider the chart

    Ut=(x2+y+g(y,z,u,t)=y2+p(y,z,u)+t=0)A5(x,y,z,u,t)/17(5,3,2,1,6).

    We choose (y1,...,y5)=(y,z,u,x,t) with δ4=4 and δ5=1. Then (Ξ) holds. The origin of Ut is a cyclic quotient point. Let F be the exceptional divisor corresponding to the w-morphism over this point. Then

    vF(y1,...,y5)=17(10,2,1,5,6)

    (notice that the irreducibility of y2+p(y,z,u) implies that p(0,z,u)0, so vF(t)=67). One can see that (Θ1) holds. Thus, there exists a divisorial contraction Y1X which extracts F so that YXY1. One can compute that vF(x,y,z,u)=(5,4,2,1), and so Y1X is of type E6.

    Now we study divisorial contractions over cE points with discrepancy greater than one. Those divisorial contractions are given in Table 10.

    Proposition 4.17. Assume that YX is a divisorial contraction with discrepancy a>1. Then there exists a w-morphism Y1X such that YXY1.

    Proof. Assume first that YX is of type E19. The chart UyY is defined by

    (x2+(y+p(z,u))3+u3+g(y,z,u)=0)A4(x,y,z,u)/13(0,1,1,2).

    One can choose (y1,...,y4)=(x,y+p,z,u) with δ4=1. Then (Ξ) holds. The origin of Uy is a cD/3 point. The w-morphism over this point is given by the weighted blow-up with the weight

    w(y1,...,y4)=(3,4,1,2) or (6,4,1,5).

    One can see that (Θ4) holds. Thus, there exists Y1X such that YXY1. A direct computation shows that Y1X is a w-morphism.

    Assume that YX is of type E20. The chart UtY is defined by

    (x2+y+g(y,z,u,t)=y2+p(z,u)+t=0)A5(x,y,z,u,t)/17(5,3,2,2,6).

    We take (y1,...,y5)=(y,z,u,x,t) with δ4=4 and δ5=1. Then (Ξ) holds. The w-morphism over the origin of Ut is given by weighted blowing-up the weight w(y1,...,y5)=17(5,1,1,6,3). One can see that (Θ1) holds. Hence, there exists a divisorial contraction Y1X such that YXY1. One can compute that Y1X is a w-morphism.

    Finally assume that YX is of type E21. The chart UyY is defined by

    (x2+y+u7+g(y,z,u)=0)A4(x,y,z,u)/15(2,4,3,2).

    One can choose (y1,...,y4)=(y,z,u,x). Then (Ξ) holds. The w-morphism over Uy is given by the weighted blow-up with the weight w(y1,...,y4)=15(2,4,1,1). One can see that (Θ2) holds. Thus, there exists a divisorial contraction Y1X such that YXY1. One can compute that Y1X is a w-morphism.

    Finally, we need to study divisorial contractions over cE/2 points. All such divisorial contractions are listed in Table 11.

    Proposition 4.18. Let YX be a divisorial contraction.

    (1) Assume that there are two w-morphisms over X. Then:

    (1–1) If YX is of type E22 and there exists another w-morphism Y1X, then Y1X is of type E22 or E23 and YXY1.

    (1–2) If YX is of type E23, then there are exactly two w-morphisms. The other one, Y1X, is of type E22. Interchanging Y and Y1, we are back to Case (1–1).

    (1–3) If YX is of type E24, then there are exactly two w-morphisms. They are both of type E24 and are negatively linked to each other.

    (1–4) YX is not of type E25.

    (2) Assume that YX is of type E26. Then there is a w-morphism Y1X which is of type E22 such that YXY1.

    Proof. The statement about the number of w-morphisms follows from [9, Section 10]. First, assume that there exists two w-morphisms over X and YX is of type E22. Let F=exc(Y1X). The only non-Gorenstein point on Y is the origin of Uz, which is a cD/3 point defined by

    (x2+y3+g(y,z,u)=0)A4(x,y,z,u)/13(0,2,1,1).

    One can see that

    12=a(F,X)=a(F,Y)+12vF(z),

    hence a(F,Y)=vF(z)=13. Thus, F corresponds to a w-morphism over Uz. From Table 6 we know that

    vF(x,y,z,u+λz)=13(b,c,1,4)

    for some λC, where (b,c)=(3,2) or (6,5). Now one can choose (y1,...,y4)=(y,u+λz,u+ξz,x) with δ4=4, where ξC is a number so that u+ξz defines a Du Val section on X and ξλ. Thus, (Ξ) and (Θ2) hold and YXY1.

    Now assume that YX is of type E24. The chart UxY is defined by

    (x+p(y,z,u)+y3+g(x,y,z,u)=0)A4(x,y,z,u)/17(5,4,3,1).

    The origin is a cyclic quotient point and there is only one w-morphism over this point. Let F be the exceptional divisor corresponding to this w-morphism. By Lemma 3.13 we know that there exists a divisorial contraction Y1X which extracts F such that YXY1. One can compute that a(F,X)=12. Hence, Y1X is also a w-morphism.

    Finally, assume that YX is of type E26. Consider the chart UyY which is defined by

    (x2+y+z4+u8+g(y,z,u)=0)A4(x,y,z,u)/16(1,2,5,1).

    One can take (y1,...,y4)=(y,z,u,x) with δ4=4. One can see that (Ξ) holds. Now, let F be the exceptional divisor corresponding to the w-morphism over Ut. Then vF(y1,...,y4)=16(2,5,1,1). Hence, (Θ2) and (Ξ) hold. Thus, there exists a divisorial contraction Y1X which extracts F so that YXY1. One can compute that vF(x,y,z,u)=12(3,2,3,1), so Y1X is of type E22.

    We want to understand the change of singularities after running the minimal model program. The final result is the following proposition.

    Proposition 5.1.

    (1) Assume that YX is a divisorial contraction between terminal and Q-factorial threefolds.

    (1–1) If YX is a divisorial contraction to a point, then

    gdep(X)gdep(Y)+1anddep(X)dep(Y)+1.

    If YX is a divisorial contraction to a curve, then

    gdep(X)gdep(Y)anddep(X)dep(Y).

    (1-2) depGor(X)depGor(Y) and the inequality is strict if the non-isomorphic locus on X contains a Gorenstein singular point.

    (2) Assume that XX is a flip between terminal and Q-factorial threefolds.

    (2–1)

    gdep(X)>gdep(X)anddep(X)>dep(X).

    (2–2) depGor(X)depGor(X) and the inequality is strict if the non-isomorphic locus on X contains a Gorenstein singular point.

    Corollary 5.2. Assume that

    X0X1...Xk

    is a process of the minimal model program. Then:

    (1) ρ(X0/Xk)gdep(Xk)gdep(X0) and the equality holds if and only if XiXi+1 is a strict w-morphism for all i.

    (2) depGor(Xk)depGor(X0).

    In particular, if X is a terminal Q-factorial threefold and WX is a resolution of singularities, then ρ(W/X)gdep(X) and the equality holds if and only if W is a feasible resolution of X.

    Proof. Statement (2) easily follows from the inequalities in Proposition 5.1. Assume that the sequence contains m flips, then ρ(X0/Xk)=km. On the other hand, we know that

    gdep(Xi+1){gdep(Xi)1if XiXi+1 is a flipgdep(Xi)+1otherwise.

    It follows that gdep(Xk)gdep(X0)+k2m, hence gdep(Xk)gdep(X0)ρ(X0/Xk). Now gdep(Xk)gdep(X0)=ρ(X0/Xk) if and only if m=0 and gdep(Xi+1)=gdep(Xi)+1 for all i, which is equivalent to XiXi+1 being a strict w-morphism for all i.

    Now assume that X is a terminal Q-factorial threefold and WX is a resolution of singularities. We can run KW-MMP over X and the minimal model is X itself. Since gdep(W)=0, one has that ρ(W/X)gdep(X) and the equality holds if and only if W is a feasible resolution of X.

    The inequalities for the depth part are exactly Lemma 2.20. We only need to prove the inequalities for the generalized depth and the Gorenstein depth.

    Convention 5.3. Let S be a set consisting of birational maps between Q-factorial terminal threefolds. We say that ()S holds if, for all ZV inside S, one has that:

    (1) If ZV is a divisorial contraction to a point, then

    depGor(V)depGor(Z)depGor(V)(dep(Z)dep(V)+1).

    (2) If ZV is a divisorial contraction to a smooth curve, then

    depGor(V)depGor(Z)depGor(V)(dep(Z)dep(V)).

    (3) If ZV is a flip, then

    depGor(V)depGor(Z)depGor(V)(dep(Z)dep(V)1).

    (4) If ZV is a flop, then depGor(V)=depGor(Z).

    Moreover, if there exists a Gorenstein singular point PV such that P is not contained in the isomorphic locus of ZV, then depGor(V)>depGor(Z) unless VZ is a flop.

    We say that ()(1)S holds if statement (1) is true, but statements (2) and (3) are unknown.

    If VZ is a flip or a divisorial contraction, we denote the condition ()VZ=()S where S is the set containing only one element VZ.

    It is easy to see that if YX is a divisorial contraction, then ()YX holds if and only if the inequalities in Proposition 5.1 (1) hold. Likewise, if XX is a flip, then ()XX holds if and only if the inequalities in Proposition 5.1 (2) hold.

    Remark 5.4. If ZV is a flop, then the singularities on Z and V are the same by [17, Theorem 2.18]. Hence, statement (4) is always true.

    Convention 5.5. Given nZ0, we denote

    Sn={ϕ:ZV|ϕ is a flip, a flop or a divisorial contraction between  between Q-factorial terminal threefolds, gdep(Z)n}.

    Lemma 5.6. Assume that YX and Y1X are two divisorial contractions between terminal threefolds, such that YXY1. If ()Sgdep(Y)1 holds, then gdep(Y1)gdep(Y). Moreover, gdep(Y1)=gdep(Y) if and only if Y1XY.

    Proof. We have a diagram

    such that Z1Y is a strict w-morphism and ZiZi+1 is a flip or a flop for all 1ik1. Since gdep(Z1)=gdep(Y)1 and ()Sgdep(Y)1 holds, we know that gdep(Z2)gdep(Y)1. Repeating this argument k2 times, one can say that gdep(Zk)gdep(Y)1. Again, since ()Sgdep(Y)1 holds, we know that gdep(Y1)gdep(Zk)+1=gdep(Y).

    Now, gdep(Y1)=gdep(Y) if and only if all the inequalities above are equalities. This is equivalent to ZkY1 being a strict w-morphism and ZiZi+1 being a flop for all i=1, ..., k1, or k=1. In other words, we also have Y1XY.

    Corollary 5.7. Assume that YX is a strict w-morphism over PX and Y1X is another divisorial contraction over P. If ()Sgdep(Y) holds, then there exist divisorial contractions Y1X, ..., YkX such that

    Y1X...XYkXY.

    Proof. By Proposition 4.1, we know that there exists Y1, ..., Yl, Y1, ..., Yl such that

    Y1X...XYl=YlX...XY1XY.

    One can apply Lemma 5.6 to the sequence YXY1X...XYl and conclude that gdep(Yi)gdep(Y) for all i. Since YiXSgdep(Y) for all i=1, ..., l, one has gdep(Yi)gdep(X)1=gdep(Y). Thus, gdep(Yi)=gdep(Y) and Lemma 5.6 says that one has

    YlX...XY1XY.

    Now we can take k=l+l1 and let Yi=Yli+l for l<ik.

    Corollary 5.8. Assume that ()Sn2 holds. Assume that PX is a cA/r point or a cAx/r point such that gdep(X)=n. Then every w-morphism over P is a strict w-morphism. In particular, one can always assume that the morphism Y1X in Theorem 2.18 is a strict w-morphism.

    Proof. Propositions 4.2 and 4.3 say that if YX and Y1X are two different w-morphisms over P, then there exists Y2X, ..., YkX such that

    Y1XY2X...XYk=Y.

    We can assume that YX is a strict w-morphism, so gdep(Y)=n1. Lemma 5.6 implies that Y1 is also a strict w-morphism. Hence, every w-morphism over P is a strict w-morphism.

    Now assume that X is in the diagram in Theorem 2.18 and P is a non-Gorenstein point in the exceptional set of XW. If P is a cA/r or a cAx/r point, then we already know that every w-morphism over P is a strict w-morphism. Otherwise, by Remark 2.19 (2) we know that any w-morphism Y1X over P induces a diagram in Theorem 2.18. Hence, we can choose Y1X to be a strict w-morphism.

    Convention 5.9. Let DV be the set of symbols

    DV={Ai,Dj,Ek}iN,jN4,k=6,7,8.

    One can define an ordering on DV by

    Ai<Ai<Dj<Dj<Ek<Ekfor alli<i,j<j,k<k.

    Given DV, define

    T={X is a terminal GE(PX) for all Q-factorial threefold  non-Gorenstein point PX},
    T,n={XTgdep(X)n}

    and

    Tn=DVT,n.

    Here, GE(PX) denotes the type of the general elephant near P. That is, the type of a general Du Val section H|KX| near an analytic neighborhood of PX.

    Convention 5.10. Let T be a set of terminal threefolds. We say that the condition (Π)T holds if for all XT and for all strict w-morphisms YX over non-Gorenstein points of X, one has that dep(Y)=dep(X)1.

    Remark 5.11.

    (1) Assume that YX is a w-morphism over a non-Gorenstein point P. Then the general elephant of Y over X is better than the general elephant of X near P. This is because if H|KX| near P, then HY|KY| and HYH is a partial resolution by [3, Lemma 2.7].

    (2) One has that (Π)TA1 always holds since if GE(PX)=A1 for some non-Gorenstein point P, then P is a cyclic quotient point of index two (cf. [7, (6.4)]). In this case, there is only one w-morphism YX over P and Y is smooth over X. Hence,

    gdep(PX)=dep(PX)=1

    and gdep(Y)=dep(Y)=0 over X. Thus, YX is a strict w-morphism and also dep(Y)=dep(X)1.

    Lemma 5.12. Assume that YX is a strict w-morphism over a non-Gorenstein point P. If ()Sgdep(Y) holds and (Π)T holds for all <GE(PX), then dep(Y)=dep(X)1.

    Proof. By the definition we know that dep(Y)dep(X)1. Assume that dep(Y)>dep(X)1. Then there exists Y1X such that dep(Y1)=dep(X)1<dep(Y). Corollary 5.7 says that there exist divisorial contractions Y2X, ..., YkX such that Y1XY2X...XYkXY. From Remark 5.11 (1) we know that YiT for some <GE(PX) for all i, hence, if

    is the induced diagram of YiXYi+1, then dep(Zi,1)=dep(Yi)1. By Lemma 2.20 we know that

    dep(Yi+1)dep(Zi,ki)+1dep(Zi,1)+1=dep(Yi)

    for all i. Hence, dep(Y1)dep(Y). This leads to a contradiction as we assume that dep(Y)>dep(Y1).

    Lemma 5.13. Fix an integer n and assume that ()Sn1 holds. Then (Π)Tn holds.

    Proof. We need to show that for all XTn and for all strict w-morphism YX over a non-Gorenstein point, one has that dep(Y)=dep(X)1. By Remark 5.11 (2) we know that (Π)TA1 holds, hence (Π)TAm,n and (Π)TD4,n hold for all mN by Lemma 5.12 and by induction on m. Then, one can prove that (Π)TDm,n and (Π)TE6,n hold by again applying Lemma 5.12 and by induction on m. Statements (Π)TE7,n and (Π)TE8,n can be proved in the same way.

    Lemma 5.14. Fix an integer n and assume that ()Sn1 holds. Assume that we have a diagram

    such that gdep(X)=n, Y1X is a strict w-morphism, YiYi+1 is a flip or a flop for i=1, ..., k1, and YkX is a divisorial contraction. Then,

    (1) depGor(X)=depGor(Y1)depGor(X) and gdep(X)gdep(X).

    (2) ()YiYi+1 holds for i=1, ..., k1.

    (3) ()YkX holds.

    Moreover, if the non-isomorphic locus of XX on X contains a Gorenstein singular point, then depGor(X)<depGor(X).

    Proof. Since ()Sn1 holds, we know that (Π)Tn holds by Lemma 5.13. Hence, depGor(Y1)=depGor(X). We know that gdep(Y)=n1. Since ()Sn1 holds, one can prove that gdep(Yi)gdep(Y1)n1 for all i and hence ()YiYi+1 and ()YkX hold. Thus,

    gdep(X)gdep(Yk)1gdep(Y1)1gdep(X)

    and

    depGor(X)=depGor(Y1)...depGor(Yk)depGor(X).

    Now, if the non-isomorphic locus of XX on X contains a Gorenstein singular point, then either the non-isomorphic locus of YiYi+1 on Yi+1 contains a Gorenstein singular point or the non-isomorphic locus of YkX on X contains a Gorenstein singular point. Hence, at least one of the above inequalities is strict. Thus, one has depGor(X)<depGor(X).

    Lemma 5.15. Assume that ()Sn1 holds. Then ()(1)Sn holds.

    Proof. Let YX be a divisorial contraction to a point which belongs to Sn. We know that gdep(Y)=n. Assume first that YX is a strict w-morphism over a point PX. Notice that dep(Y)dep(X)1 by Lemma 2.20. If P is a non-Gorenstein point, then depGor(Y)=depGor(X) by Lemma 5.13, hence ()YX holds. If P is a Gorenstein point, then

    depGor(X)=gdep(X)dep(X)=gdep(Y)+1dep(X)=depGor(Y)+dep(Y)dep(X)+1.

    Moreover, since dep(PX)=0, dep(Y)dep(X)0, hence

    depGor(X)>depGor(Y)=depGor(X)(dep(Y)dep(X)+1).

    Thus, ()YX holds.

    In general, by Corollary 5.7 there exist Y1X, ..., YkX such that YkX is a strict w-morphism and one has YXY1X...XYk. By induction on k we may assume that ()YiX holds for all i (notice that gdep(Yi)gdep(Y)=n for all i by Lemma 5.6). Now we have a diagram

    By Lemma 5.14 we know that depGor(Y)depGor(Y1) and ()ZiZi+1, ()ZkY1 hold. Since ()Y1X holds, we know that depGor(X)depGor(Y1)depGor(Y) and

    depGor(Y)+(dep(Y)dep(X)+1)=depGor(Z1)+(dep(Z1)dep(X)+2)gdep(Z1)dep(X)+2gdep(Zk)dep(X)+2gdep(Y1)dep(X)+1depGor(Y1)+(dep(Y1)dep(X)+1)depGor(X).

    Moreover, if YX is a divisorial contraction to a Gorenstein point, then Y1X is also a divisorial contraction to a Gorenstein point. Hence, depGor(Y1)<depGor(X) and we also have depGor(Y)<depGor(X).

    Proof of Proposition 5.1. We need to say that ()Sn holds for all n and we will prove this by induction on n. If n=0, then S0 consists only smooth blow-downs and smooth flops. One can see that ()S0 holds. In general, assume that ()Sn1 holds. By Lemma 5.15 we know that ()(1)Sn holds. Hence, it is enough to show that, given a flip XX or a divisorial contraction to a curve XV such that gdep(X)=n, ()XX or ()XV holds.

    If XV is a smooth blow-down, then depGor(X)=depGor(V), and so there is nothing to prove. In general, we have a diagram as in Theorem 2.18:

    By Lemma 5.14 we know that depGor(X)depGor(X) and ()YiYi+1, ()YkX hold. One has that

    depGor(X)+(dep(X)dep(X))=gdep(X)dep(X)=gdep(Y1)dep(X)+1gdep(Yk)dep(X)+1gdep(X)dep(X)=depGor(X).

    Moreover, if XX is a flip, then either one of YiYi+1 is a flip or YkX is a divisorial contraction to a curve by [3, Remark 3.4]. This implies that either gdep(Y1)>gdep(Yk) or gdep(Yk)gdep(X). If XV is a divisorial contraction to a curve, then YkX is a divisorial contraction to a curve, hence one always has that gdep(Yk)gdep(X). In conclusion, we have

    depGor(X)depGor(X)(dep(X)dep(X)1).

    If XX is a flip, then one can see that ()XX holds. Now assume that XV is a divisorial contraction to a curve. Then XV is a divisorial contraction to a point. We know that ()XV holds since ()(1)Sn holds and gdep(X)n by Lemma 5.14. One can see that depGor(X)depGor(X)depGor(V) and

    depGor(X)depGor(X)(dep(X)dep(X)1)depGor(V)(dep(X)dep(V)+1)(dep(X)dep(X)1)=depGor(V)(dep(X)dep(V)).

    Thus, ()XV holds.

    In this section, we describe the difference between two different feasible resolutions of a terminal Q-factorial threefold. The final result is the diagrams in Theorem 1.1.

    Lemma 6.1. Assume that X is a terminal threefold and YX, Y1X are two different strict w-morphisms over PX such that YXY1. Let

    be the corresponding diagram. Then one of the following holds:

    (1) YXY1XY and k=1.

    (2) k=2 and Z1Z2 is a smooth flop.

    (3) k=2, Z1Z2 is a singular flop, PX is of type cA/r and both YX and Y1X are of type A1 in Table 2.

    Proof. Since YX and Y1X are both strict w-morphisms, we know that gdep(Y)=gdep(Y1). Hence, ZiZi+1 can not be a flip. Thus, if YXY1, then .

    Now assume that . According to the result in Section 4, one has that:

    (ⅰ) If is of type , then by Proposition 4.2 we know that and are both of type A1. Since , we know that is irreducible. Hence, there is only one -trivial curve on and so .

    (ⅱ) is not of type since one always has that if and are two different -morphisms over by Proposition 4.3.

    (ⅲ) is not of type by Proposition 4.5.

    (ⅳ) is not of type by Proposition 4.7.

    (ⅴ) If is of type , then by Propositions 4.8 and 4.9 we know that is of type D17 in Table 6.

    We know that is a -morphism over the origin of

    and we take . One can see that there are two curves contained in , namely and . We know that () holds. When computing the intersection number for one can take and , hence () holds in this case. This implies that the proper transform of on is a -trivial curve and is not contained in . Thus, is a flip along for some , and hence . This is a contradiction. Thus, this case will not happen.

    (ⅵ) If is of type , then by Proposition 4.15 and Proposition 4.16 we know that is of type E14, E16, or E18 in Table 9. In those cases, are five-dimensional weighted blow-ups. Let be a curve contained in such that the proper transform of on is a possibly -trivial curve. From the proof of Proposition 4.16, one can see that:

    (ⅵ–ⅰ) is of type E14. In this case, is a -morphism over the origin of

    We choose with and or , or and or . We also know that () holds. If both and , then () holds, which implies that . This contradicts our assumption. Hence, . Now, .

    (ⅵ–ⅱ) is of type E16. In this case is a -morphism over the origin of

    As in the proof of Proposition 4.16, we choose with and or . We know that () holds. Hence, if , then () holds, and so . This leads to a contradiction. Thus, and one has .

    (ⅵ–ⅲ) is of type E18. In this case, is a -morphism over the origin of

    We choose with and . We know that () holds. If , then we can choose . Then, since () holds, we know that () holds, and so . This contradicts our assumption. Hence, . In this case, divides and has singularities, so . One can see that .

    Now the origin of is a cyclic quotient point and the -morphism over this point is a weighted blow-up with weights , and , respectively. An easy computation shows that does not pass through any singular point of , hence is a smooth flop. Since there is only one -trivial curve, we know that .

    (ⅶ) If is of type , then by Proposition 4.18 we know that is of type E22 or E23 in Table 11. Moreover, if is of type E23, then is of type E22. Thus, interchanging and if necessary, we can always assume that is of type E22. As in the proof of Proposition 4.18, we know that is a -morphism over the origin of , which is a point defined by

    The only possible -trivial curve is the lifting of the curve on . Moreover, is defined by a weighted blow-up with the weight

    for some , where or . If , then one can see that holds by considering the function , hence , which contradicts our assumption. Hence, . In this case, an easy computation shows that does not pass through any singular point of , so is a smooth flop. Since there is only one -trivial curve, we know that .

    We need to construct a factorization of the flop in Lemma 6.1 (3). Assume that

    is a singularity with and , are two different strict -morphisms with the factorization

    such that is a flop.

    Remark 6.2. One always has that since, if , then there is only one -morphism over by classification Table 2.

    Lemma 6.3. is a flop over

    where . Moreover, and .

    Proof. We may assume that is a weighted blow-up with the weight with . The chart is defined by

    and the chart is defined by

    for some . As in the proof of Proposition 4.2, we know that is a weighted blow-up over the origin of with the weight . The flopping curve of is the strict transform of the curve such that and . One can see that intersects at the origin of , which is defined by

    It is easy to see that is contained in , and on we know that is defined by .

    We have the following change of coordinates formula:

    Also,

    One can see that

    It follows that

    If we choose an isomorphism

    then , where

    by noticing that

    Now we can choose .

    Finally, we know that is a weighted blow-up with the weight and is a -morphism over the origin of . Hence, the local picture of can be obtained by a similar computation, but interchanging the role of and . One then has that .

    In this subsection, we assume that

    with for some . Let and such that is a -factorial terminal flop. Let be the weight and . Then . Let be the chart

    and be the chart

    with the relations , and . Similarly, let be the chart

    and be the chart

    with the relations , and .

    Lemma 6.4. Let be a strict -morphism. Then:

    (1) The chart is -factorial.

    (2) If then is smooth. Otherwise, contains exactly one non--factorial point which is defined by where . One then has that .

    (3) All other singular points on are cyclic quotient points.

    Proof. From Table 2, we know that is a weighted blow-up with the weight with . Statement (3) follows from direct computations. One can compute that

    with . By [18, 2.2.7] we know that is -factorial if and only if is irreducible as a -invariant function. If is not -factorial, then for some non-unit -invariant functions and . Write

    Let

    Notice that if , then if because . Hence, we have the relation . Now let

    Then and can be viewed as -invariant functions on such that for , . This means that is not irreducible. Nevertheless, the chart is defined by and has -factorial singularities. This leads to a contradiction. Thus, is -factorial.

    The chart is defined by . When , we know that is a unit along , hence is smooth. If , then is a non-unit. In this case, the origin of is a non--factorial singularity.

    From the construction in [19], for any strict -morphism we have the following diagram

    where is a -factorization and is a composition of flips for , . By Lemma 6.4, we know that if , then . Otherwise, is a flop over the singularity .

    First we discuss the factorization of . If , then is covered by four affine charts , , and . The origin of and are cyclic quotient points and all other singular points are contained in . When , the chart has a non--factorial point and there are two charts and over this point. We fix the following notation for the latter discussion.

    with , , , and .

    with , , , and .

    with , , , and .

    If , define . In this case this chart is smooth. When we define

    with , , , and .

    with , , , and .

    On the other hand, we will see later that it is enough to assume that is a divisorial contraction over the origin of . This means that is covered by five affine charts , , , and where the latter four charts correspond to the weighted blow-up . Notice that and are the same divisor since and are isomorphic in codimension one. One can see that is a weighted blow-up with the weight . Again, we use the following notation:

    with , , , and .

    with , , , and .

    with , , , and .

    with , , , and .

    with , , , and .

    Lemma 6.5. Assume that is a divisorial contraction over the origin of . Then:

    (1) .

    (2) All the singular points on the non-isomorphic loci of on both and are cyclic quotient points.

    (3) The flip is of type IA in the convention of [20, Theorem 2.2].

    Proof. It is easy to see that for , and if . Since is smooth, the only singular point contained in the non-isomorphic locus of on is the origin of . This point is a cyclic quotient point of index , so it has generalized depth .

    On the other hand, singular points on the non-isomorphic locus of on are origins of and . They are cyclic quotient points of indices and , respectively. One can then see that

    Now we know that that the flipping curve contains only one singular point which is a cyclic quotient point. Also, the general elephant of the flip is of -type since it comes from the factorization of a flop over a point. Thus, the flip is of type IA by the classification from [20], Theorem 2.2].

    Lemma 6.6. Assume that is a flip of type IA. Assume that

    is the factorization in Theorem 2.18. Then:

    (1) If is a flop, then it is a Gorenstein flop.

    (2) If is a flip, or is a flop and is a flip, then the flip is of type IA.

    Proof. Let be the flipping curve. Since is a type IA flip, there is exactly one non-Gorenstein point which is contained in , and this point is a point. From the construction, we know that is a -morphism over this point. Also, there exists a Du Val section such that . We know that by [3, Lemma 2.7 (2)]. Hence, all non-Gorenstein point of is contained in . Now assume that contains a non-Gorenstein point. Then intersects non-trivially. Since , we know that , hence is a -negative curve and so is a flip. Thus, if is a flop, then it is a Gorenstein flop.

    If is a flip, then passes through a non-Gorenstein point since by [21, Theorem 0]. Since passes through exactly one non-Gorenstein point, passes through exactly one non-Gorenstein point and this point is contained in . Since is a -morphism over a point, an easy computation shows that contains only singularities. Also, we know that is a Du Val section which does not contain . Thus, is also a flip of type IA by the classification [20, Theorem 2.2].

    Now assume that is a flop and is a flip. Then the flipping curve of is contained in since all -negative curves over are contained in . Since is a flop, we know that contains only singularities [17, Theorem 2.18] and is also a Du Val section. If , then is a flip of type IA by the classification from [20, Theorem 2.2].

    Thus, we only need to prove that . Assume that . Then, since does not intersect (otherwise is a -negative curve and then is not a flop), we know that does not intersect the flopping curve . Let be a curve which intersects non-trivially. Then for some since the both curves are contracted by . Hence, for all divisors such that , we know that . On the other hand, is a -anti-flip since is not contained in . By Corollary 2.22 we know that , hence . Now we know that and is not numerically equivalent to a multiple of , hence we may write for some . Since

    and , we know that . Hence, is not contained in the boundary of the relative effective cone . Thus, cannot be the flipping curve of . This leads to a contradiction.

    Lemma 6.7. Assume that is a three-dimensional terminal -factorial flip which satisfies conditions (1)–(3) of Lemma 6.5. Then, the factorization in Theorem 2.18 for is one of the following diagrams:

    (1)

    where is a flip which also satisfies conditions (1)–(3) of Lemma 6.5 and is a strict -morphism.

    (2)

    where is a smooth flop and is a flip which also satisfies conditions (1)–(3) of Lemma 6.5 and is a strict -morphism.

    (3)

    where is a smooth flop and where is a smooth curve contained in the smooth locus of .

    Proof. We have the factorization

    such that is a flip or a flop and is a flip for all . One has that

    If , then and is a flop. By [3, Remark 3.4] we know that is a divisorial contraction to a curve. Since singular points on the non-isomorphic locus of are all cyclic quotient points and there is no divisorial contraction to a curve which passes through a cyclic quotient point [8, Theorem 5], we know that is a divisorial contraction to a curve contained in the smooth locus. Now, we also have that , hence is smooth over and so is also a smooth curve.

    Now assume that . Then , hence either or and is a flop. Also, is a -morphism. Since singular points on the non-isomorphic locus of are all cyclic quotient points, singular points on the exceptional divisor of are all cyclic quotient points. Since flops do not change singularities [17, Theorem 2.4], we know that the singular points on the non-isomorphic locus of are all cyclic quotient points for , ..., . If is a flop, then by Lemma 6.6 we know that it is a Gorenstein flop. Since cyclic quotient points are not Gorenstein, we know that the flop is in fact a smooth flop. Now assume that is a flip for or . Then again, by Lemma 6.6 we know that it is a flip of type IA. Hence, conditions (1)–(3) of Lemma 6.5 are satisfied for this flip.

    Corollary 6.8. Assume that is a three-dimensional terminal -factorial flip which satisfies conditions (1)–(3) of Lemma 6.5. Then we have a factorization

    where and are feasible resolutions, where is a smooth curve, and is a sequence of smooth flops.

    Proof. For convenience we denote the diagram by .

    We know that the factorization of is of the form (1)–(3) in Lemma 6.7. Notice that, since the only singular points on the non-isomorphic locus of are cyclic quotient points, the feasible resolutions and are uniquely determined. Hence, is also a feasible resolution of . If the factorization of is of type (3) in Lemma 6.7, then the non-isomorphic locus of contains no singular points. This means that induces a smooth blow-up on , and induces a smooth flop . Thus, exists.

    Now, if the factorization of is of type (1), then is a feasible resolution of . Since , we may utilize induction on and assume that exists, and then can be induced by . If the factorization of is of type (2), then again by induction we may assume that

    exists. Since is a strict -morphism, we know that . Also, since is a smooth flop, it induces a smooth flop . If we let and let be the composition , then we get the diagram .

    Definition 6.9. Let be a birational map between smooth threefolds. We say that is of type if it is a composition of smooth flops. We say that is of type if there exists the diagram

    such that:

    (1) for some smooth curve for , .

    (2) is a composition of smooth flops for , .

    (3) has the factorization

    such that is a birational map of type for some .

    Example 6.10. In the following diagrams, dashmaps stand for smooth flops and all other maps are blowing-down smooth curves.

    (1)

    (2)

    (3)

    In diagram (1), is of type . In both diagram (2) and (3), is of type .

    Definition 6.11. Let be a birational map between smooth threefolds. We say that has an -type factorization if there exists birational maps between smooth threefolds

    such that is of type for some .

    Proposition 6.12. Assume that is a -factorial terminal threefold and , are two different feasible resolutions. Then the birational map has an -type factorization.

    Proof. First, notice that if , then has either a cyclic quotient point of index 2, or a point defined by for or by [22, Corollary 3.4]. In those cases, there is exactly one feasible resolution (which is obtained by blowing-up the singular point). Hence, one may assume that . Let (resp. ) be the strict -morphism which is the first factor of (resp. ). If , then and are two different feasible resolutions of . In this case, the statement can be proved by induction on . Thus, we may assume that .

    Since both and are strict -morphisms, by Corollary 5.7 there exists a sequence of strict -morphisms , , ..., such that for , ..., . For each , let be a feasible resolution. Then is also a feasible resolution of and it is enough to prove that our statement holds for and , for all , ..., . Thus, we may assume that .

    We have the diagram

    Lemma 6.1 says that there are three possibilities. If , then and are two different feasible resolutions of . Since , again by induction on we know that has an -type factorization. Assume that and is a smooth flop. Then it induces a smooth flop where is a feasible resolution of for , . Also, we know that and are two feasible resolutions of , and and are two feasible resolutions of . Again, by induction on , we know that both and have -type factorizations, hence does as well.

    Finally, assume that we are in the case of Lemma 6.1 (3), namely that has a singularity and is a singular flop. By Lemma 6.3 we know that is a flop over

    We have the factorization of the flop

    We use the notation at the beginning of this subsection. Assume that . Then we can choose to be the weighted blow-up with the weight . In this case, is a divisorial contraction over the origin of for , since both and . When , let be the weighted blow-up with the weight . Then, is a divisorial contraction over the origin of , but is a divisorial contraction over the origin of . Since and by Remark 6.2, we know that for some unit . If we let , then we can write the defining equation of as where . One has that , and under this notation one also has that is a divisorial contraction over the origin of . In conclusion, Corollary 6.8 holds for both and .

    Let be a feasible resolution. If , then one can see that is of type . Assume that is a flop. Notice that is a -morphism since where . One has that

    where the second equality follows from Corollary 5.8. Now we know that is a flop over with . By induction on , we may assume that has an -type factorization where is a feasible resolution corresponding to the diagram in Corollary 6.8. One can see that can be connected by a diagram of the form for some . Finally, we know that and (resp. and ) are feasible resolutions of (resp. ). Again, by induction on , we may assume that and have -type factorizations. Hence, has an -type factorization.

    Remark 6.13. Assume that is a singularity. Then there exists feasible resolutions and such that the birational map is connected by for .

    In this section, we prove our main theorems. First, we recall some definitions which are defined in the introduction section.

    Definition 7.1. Let be a projective variety. We say that a resolution of singularities is a P-minimal resolution if for any smooth model one has that .

    Definition 7.2. Let be a birational map between smooth varieties. We say that this birational map is a P-desingularization of a flop if there exists a flop such that and are P-minimal resolutions.

    Proposition 7.3. Assume that is a threefold. Then, is a P-minimal resolution if and only if is a feasible resolution of a terminalization of . In particular, if is a terminal and -factorial threefold, then P-minimal resolutions of coincide with feasible resolutions.

    Proof. Let be a resolution of singularities and let be the -MMP over . Then is a terminalization of . We know that by Corollary 5.2. Assume first that is -minimal. Let be a feasible resolution of , then . Since is also a smooth resolution of , the inequality is an equality. Therefore, , which implies that is a sequence of strict -morphisms by Corollary 5.2, or, equivalently, is a feasible resolution.

    Conversely, assume that is not -minimal, but it is a feasible resolution of some which is a terminalization of . There exists a -minimal resolution such that . From the above argument, there exists a terminalization of such that is a feasible resolution. Hence, . However, since terminalizations are connected by flops [23, Theorem 1] and flops do not change singularities by [17, Theorem 2.4], we know that . This means that , so can not be a feasible resolution of . This is a contradiction.

    Proof of Theorem 1.1. Let be a threefold and , be two P-minimal resolutions. By Proposition 7.3, we know that (resp. ) is a feasible resolution of a terminalization (resp. ). If , then and are connected by flops [23, Theorem 1], hence is connected by P-desingularizations of terminal -factorial flops.

    Now assume that . Then and are two different feasible resolutions of . The first two paragraphs in the proof of Proposition 6.12 and Lemma 6.7 imply that can be also connected by P-desingularizations of terminal -factorial flops. Moreover, Proposition 6.12 says that those P-desingularizations of flops can be factorized into compositions of diagrams of the form . This finishes the proof.

    Remark 7.4. Assume that is a terminal -factorial threefold and , are two different P-minimal resolutions. We know that and can be connected by P-desingularizations of flops. Let be a P-desingularization of a flop which appears in the factorization of . Then, from the construction we know that .

    Now we compare an arbitrary resolution of singularities to a P-minimal resolution.

    Definition 7.5. Let be a birational map where is a smooth threefold and is a terminal threefold. We say that the birational map has a bfw-factorization if can be factorized into a composition of smooth blow-downs, P-desingularizations of flops, and strict -morphisms.

    Remark 7.6. If is a strict -morphism and is a smooth blow-down or a P-desingularization of a flop, then on the indeterminacy locus of is disjoint to the indeterminacy locus of since the former one lies on the smooth locus of and the latter one is a singular point. Hence, there exists where is a smooth blow-down or a P-desingularization of a flop, and is a strict -morphism. In other words, has a bfw-factorization if and only if there exists a birational map which is a composition of smooth blow-downs and P-desingularization of flops, where is a feasible resolution of .

    Proposition 7.7. Assume that a birational map has a bfw-factorization where is a smooth threefold and is a terminal threefold.

    (1) If is a flop, then there is a birational map which has a bfw-factorization.

    (2) If is a strict -morphism, then there exists a birational map which also has a bfw-factorization.

    (3) If is a flip or a divisorial contraction, then the induced birational map has a bfw-factorization.

    Proof. Assume first that is a flop. By Remark 7.6 we know that there exists a bfw-map , where is a feasible resolution of . Let be a feasible resolution of . Then is a composition of strict -morphisms and the induced birational map is a P-desingularization of the flop . It follows that the composition

    is a bfw-map. This proves (1).

    We will prove (2) and (3) by induction on . If , then is smooth. In this case, there is no strict -morphism or flip . Assume that is a divisorial contraction. If is smooth, then it is a smooth blow-down by [22, Theorem 3.3, Corollary 3.4], and if is singular, then should be a strict -morphism since in this case is terminal -factorial and is a P-minimal resolution of . Now we may assume that and statements (2) and (3) hold for threefolds with generalized depth less then .

    Let

    be a sequence of birational maps so that is a smooth blow-down, a P-desingularization of a flop, or a strict -morphism for all . By Remark 7.6 we can assume that is a strict -morphism. Now, given a strict -morphism , if , then there is nothing to prove. Otherwise, by Corollary 5.7 there exists a sequence of strict -morphisms , ..., such that

    For each , one has the factorization

    such that is a composition of flops, is a strict -morphism, and

    By the induction hypothesis, we know that if there exists a bfw-map , then there exists a bfw-map . Now one can prove statement (2) by induction on .

    Assume that is a flip. Then we have a factorization

    as in Theorem 2.18. Since is a strict -morphism by Corollary 5.8, there exists a bfw-map . Since

    the induction hypothesis implies that there exists a bfw-map .

    Finally, assume that is a divisorial contraction. If it is a smooth blow-down or a strict -morphism, then there is nothing to prove. Otherwise, there exists a diagram

    such that is a strict -morphism and is a flip or a flop for all . One has that

    hence there exists a bfw-map . If is a divisorial contraction to a curve, then is a divisorial contraction to a curve as in Theorem 2.18. In this case, we also have , so there exists a bfw-map . If is a divisorial contraction to a point, then the discrepancy of is less than the discrepancy of unless is a -morphism. Also, when is a -morphism, we know that by Lemma 5.6. Thus, we can prove statement (3) by induction on the generalized depth and the discrepancy of over .

    One can easily see the following corollary:

    Corollary 7.8. Assume that is a smooth threefold and is a birational map which is a composition of steps of MMP. Then, this birational map can be factorized into a composition of smooth blow-downs, P-desingularizations of flops, and strict -morphisms.

    Proof of Theorem 1.2. By Corollary 7.8 and Remark 7.6 we know that there exists a feasible resolution such that is a composition of smooth blow-downs and P-desingularizations of flops, where is a minimal model of over . By Proposition 7.3 we know that is also a P-minimal resolution of , hence the birational map is connected by P-desingularizations of flops. Thus, the composition is connected by smooth blow-downs and P-desingularizations of flops.

    Proof of Corollary 1.3. Let

    be a sequence of smooth blow-downs and P-desingularization of flops as in Theorem 1.2. We only need to show that if is a P-desingularization of a flop , then for all , ..., .

    By [24, Lemma 2.12] we know that for all . Since and have the same analytic singularities [17, Theorem 2.4], there exists a feasible resolution such that for all . Now, and are two different P-minimal resolutions of , so they can be connected by P-desingularizations of flops with smaller generalized depth by Remark 7.4. By induction on the generalized depth, one can see that . Hence, for all , ..., .

    This section is dedicated to exploring minimal resolutions for singularities in higher dimensions and the potential applications of our main theorems.

    In three dimensions, -minimal resolutions appear to be a viable generalization of minimal resolutions for surfaces. However, in higher dimensions, -minimal resolutions are not good enough. For example, let be a smooth flip (eg. a standard flip [1, Section 11.3]). Then, and are both -minimal resolutions of the underlying space, but is better than . It is reasonable to assume that is a minimal resolution, while is not. Inspired by Corollary 1.3, we define a new kind of minimal resolution:

    Definition 8.1. Let be a projective variety over complex numbers. We say that a resolution of singularities is a -minimal resolution if for any smooth model one has that for all .

    As stated in Corollary 1.3, B-minimal resolutions coincide with P-minimal resolutions in dimension three. Our main theorems say that B-minimal resolutions of threefolds satisfy certain nice properties. It is logical to anticipate that B-minimal resolutions of higher-dimensional varieties share similar properties.

    Conjecture. For any projective variety over the complex numbers, one has that:

    (1) B-minimal resolutions of exist.

    (2) Two different B-minimal resolutions are connected by desingularizations of -factorial terminal flops.

    (3) If is a B-minimal resolution and is an arbitrary resolution of singularities, then can be connected by smooth blow-downs, smooth flips, and desingularizations of -factorial terminal flops.

    Let

    where if is connected by P-desingularizations of -factorial terminal flops. For , we say that if there exist and so that and is a smooth blow-down. Then, Theorems 1.1 and 1.2 imply the following.

    Corollary 8.2. Given a threefold , let

    Then has a unique minimal element.

    In other words, if we consider the resolution of singularities inside , then there is a unique minimal resolution, which behaves similarly to the minimal resolution of a surface.

    As a consequence, inside the space the following strong factorization theorem holds.

    Theorem 8.3 (Strong factorization theorem for ). Assume that and are smooth threefolds which are birational to each other. Then there exists a smooth threefold such that inside one has for , .

    Proof. Let be a common resolution. Then for , . Since the minimal element of is itself, one has that for , .

    One can characterize a surface singularity by the information of exceptional curves on the minimal resolution. One may ask, does a similar phenomenon happen for higher-dimensional singularities? Since for higher-dimensional singularities there is no unique minimal resolution, what we really want to study is the following object.

    Definition 8.4. Let be a projective threefold over the complex numbers. We say that a divisorial valuation over is an almost essential valuation if for any P-minimal resolution one has that is an irreducible component of the exceptional locus of .

    This name comes from the "essential valuation" in the theories of arc spaces.

    Definition 8.5. Let be a variety. We say that a divisorial valuation over is an essential valuation if for any resolution of singularities one has that is an irreducible component of the exceptional locus of .

    From the definition, one can see that essential valuations are almost essential, but an almost essential valuation may not be essential.

    Example 8.6. Let for some . There is exactly one -morphism over the singular point, which is obtained by blowing-up the origin. There is only one singular point on , which is defined by . Blowing-up the singular point more times, we get a resolution of singularities . From the construction we know that is a unique feasible resolution of . Since is terminal and -factorial, is the unique P-minimal resolution of . Hence, almost essential valuations of are those divisorial valuations which appear on . One can compute that such that . On the other hand, by [25, Lemma 15] we know that essential valuations of are and . Hence, , ..., are almost essential valuations which are not essential.

    Notice that the set of essential valuations does not really characterize the singularity since it is independent of . The set of almost essential valuations carries more information of the singularity.

    Let be a smooth variety. The bounded derived category of coherent sheaves of , denoted by , is an interesting subject of investigation. One possible method to study is to construct a semi-orthogonal decomposition of (refer to [26] for more information). Orlov [27] proved that a smooth blow-down yields a semi-orthogonal decomposition. In particular, if is a smooth surface, then the -MMP is a series of smooth blow-downs, thereby resulting in a semi-orthogonal decomposition of .

    Now assume that is a smooth threefold and let

    be the process of -minimal model program. According to Corollary 7.8 and Remark 7.6, can be factored into a composition of smooth blow-downs and P-desingularizations of flops, where is a P-minimal resolution of . If every P-desingularizations of flops that appears in the factorization is a smooth flop, then the sequence induces a semi-orthogonal decomposition of since smooth flops are derived equivalent [28].

    Example 8.7. Let be the flip which is a quotient of an Atiyah flop by an -action [16, Example 2.7]. Then is smooth and has a singular point. Let be the smooth resolution obtained by blowing-up the singular point. Then is a sequence of MMP.

    The factorization of the flip is exactly diagram (3) in Lemma 6.7, namely the diagram

    where is a smooth flop and is a blow-down of a smooth curve. We know that there exists an equivalence of category and a semi-orthogonal decomposition . Hence, is a semi-orthogonal decomposition.

    In general, a P-desingularization of a flop may not be derived equivalent since and may not be isomorphic in codimension one. Nevertheless, due to the symmetry between and , one might expect that a semi-orthogonal decomposition on will result in a semi-orthogonal decomposition on . It still hopeful that our approach will be effective for all smooth threefolds.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The author thanks Jungkai Alfred Chen for his helpful comments. The author thanks the referees for carefully reading the paper and for their useful comments. The author is supported by KIAS individual Grant MG088901.

    The authors declare there is no conflicts of interest.



    [1] G. Bao, P. Li, Maxwell's Equations in Periodic Structures, Series on Applied Mathematical Sciences, Science Press, Beijing/Springer, Singapore, 2022.
    [2] W. Cai, Computational Methods for Electromagnetic Phenomena Electrostatics in Solvation, Scattering, and Electron Transport, Cambridge University Press, 2013.
    [3] L. Demkowicz, Computing With hp-Adaptive Finite Elements, I: One and Two-Dimensional Elliptic and Maxwell Problems, CRC Press, Taylor and Francis, 2006.
    [4] L. Demkowicz, J. Kurtz, D. Pardo, M. Paszenski, W. Rachowicz, A. Zdunek, Computing with hp-Adaptive Finite Elements, II: Frontiers: Three Dimensional Elliptic and Maxwell Problems with Applications, CRC Press, Taylor and Francis, 2007.
    [5] J. M. Jin, The Finite Element Method in Electromagnetics, 3rd edition, Piscataway, NJ, USA, IEEE Press, 2014.
    [6] P. Monk, Finite Element Methods for Maxwell's Equations, Oxford University Press, 2003.
    [7] A. Taflove, S. C. Haguess, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd edition, Artech House, Norwood, 2005.
    [8] A. Taflove, A. Oskooi, S. G. Johnson, Advances in FDTD Computational Electrodynamics: Photonics and Nanotechnology, Artech House, Boston, 2013.
    [9] T. J. Cui, D. R. Smith, R. Liu, Metamaterials: Theory, Design, and Applications, Springer, 2010.
    [10] N. Engheta, R. W. Ziolkowski, Electromagnetic Metamaterials: Physics and Engineering Explorations, Wiley-IEEE Press, 2006.
    [11] L. Solymar, E. Shamonina, Waves in Metamaterials, Oxford University Press, 2009.
    [12] Q. Bao, K. Loh, Graphene photonics, plasmonics, and broadband optoelectronic devices, ACS Nano, 6 (2012), 3677–3694. https://doi.org/10.1021/nn300989g doi: 10.1021/nn300989g
    [13] A. K. Geim, K. S. Novoselov, The rise of graphene, Nat. Mater., 6 (2007), 183–191. https://doi.org/10.1038/nmat1849
    [14] V. G. Veselago, The electrodynamics of substances with simultaneously negative values of and , Sov. Phys. Usp., 47 (1968), 509–514.
    [15] D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat-Nasser, S. Schultz, Composite medium with simultaneously negative permeability and permittivity, Phys. Rev. Lett., 84 (2000), 4184–4187. https://doi.org/10.1103/PhysRevLett.84.4184 doi: 10.1103/PhysRevLett.84.4184
    [16] A. Shelby, D. R. Smith, S. Schultz, Experimental verification of a negative index of refraction, Science, 292 (2001), 489–491. https://doi.org/10.1126/science.1058847 doi: 10.1126/science.1058847
    [17] J. B. Pendry, Negative refraction makes a perfect lens, Phys. Rev. Lett., 85 (2000), 3966–3969. https://doi.org/10.1103/PhysRevLett.85.3966 doi: 10.1103/PhysRevLett.85.3966
    [18] A. J. Holden, Towards some real applications for negative materials, Photonics Nanostruct. Fundam. Appl., 3 (2005), 96–99. https://doi.org/10.1016/j.photonics.2005.09.014 doi: 10.1016/j.photonics.2005.09.014
    [19] Y. Hao, R. Mittra, FDTD Modeling of Metamaterials: Theory and Applications, Artech House Publishers, 2008.
    [20] D. H. Werner, D. H. Kwon, Transformation Electromagnetics and Metamaterials: Fundamental Principles and Applications, Springer, 2013.
    [21] V. A. Bokil, Y. Cheng, Y. Jiang, F. Li, Energy stable discontinuous Galerkin methods for Maxwell's equations in nonlinear optical media, J. Comput. Phys., 350 (2017), 420–452. https://doi.org/10.1016/j.jcp.2017.08.009 doi: 10.1016/j.jcp.2017.08.009
    [22] N. Schmitt, C. Scheid, S. Lanteri, A. Moreau, J. Viquerat, A DGTD method for the numerical modeling of the interaction of light with nanometer scale metallic structures taking into account non-local dispersion effects, J. Comput. Phys., 316 (2016), 396–415. https://doi.org/10.1016/j.jcp.2016.04.020 doi: 10.1016/j.jcp.2016.04.020
    [23] J. Lin, An adaptive boundary element method for the transmission problem with hyperbolic metamaterials, J. Comput. Phys., 444 (2021), 110573. https://doi.org/10.1016/j.jcp.2021.110573 doi: 10.1016/j.jcp.2021.110573
    [24] B. Donderici, F. L. Teixeira, Mixed finite-element time-domain method for transient Maxwell equations in doubly dispersive media, IEEE Trans. Microwave Theory Tech., 56 (2008), 113–120. https://doi.org/10.1109/TMTT.2007.912217 doi: 10.1109/TMTT.2007.912217
    [25] J. Li, J. S. Hesthaven, Analysis and application of the nodal discontinuous Galerkin method for wave propagation in metamaterials, J. Comput. Phys., 258 (2014), 915–930. https://doi.org/10.1016/j.jcp.2013.11.018 doi: 10.1016/j.jcp.2013.11.018
    [26] C. Scheid, S. Lanteri, Convergence of a Discontinuous Galerkin scheme for the mixed time domain Maxwell's equations in dispersive media, IMA J. Numer. Anal., 33 (2013), 432–459. https://doi.org/10.1093/imanum/drs008 doi: 10.1093/imanum/drs008
    [27] Z. Xie, J. Wang, B. Wang, C. Chen, Solving Maxwell's equation in meta-materials by a CG-DG method, Commun. Comput. Phys., 19 (2016), 1242–1264. https://doi.org/10.4208/cicp.scpde14.35s doi: 10.4208/cicp.scpde14.35s
    [28] W. Li, D. Liang, The spatial fourth-order compact splitting FDTD scheme with modified energy-conserved identity for two-dimensional Lorentz model, J. Comput. Appl. Math., 367 (2020), 112428. https://doi.org/10.1016/j.cam.2019.112428 doi: 10.1016/j.cam.2019.112428
    [29] W. Li, D. Liang, Symmetric energy-conserved S-FDTD scheme for two-dimensional Maxwell's equations in negative index metamaterials, J. Sci. Comput., 69 (2016), 696–735. https://doi.org/10.1007/s10915-016-0214-9 doi: 10.1007/s10915-016-0214-9
    [30] X. Bai, H. Rui, New energy analysis of Yee scheme for metamaterial Maxwell's equations on non-uniform rectangular meshes, Adv. Appl. Math. Mech., 13 (2021), 1355–1383. https://doi.org/10.4208/aamm.OA-2020-0208 doi: 10.4208/aamm.OA-2020-0208
    [31] S. Nicaise, Stability and asymptotic properties of dissipative evolution equations coupled with ordinary differential equations, Math. Control Relat. Fields, 13 (2023), 265–302. https://doi.org/10.3934/mcrf.2021057 doi: 10.3934/mcrf.2021057
    [32] P. Fernandes, M. Raffetto, Well-posedness and finite element approximability of time-harmonic electromagnetic boundary value problems involving bianisotropic materials and metamaterials, Math. Model. Methods Appl. Sci., 19 (2009), 2299–2335. https://doi.org/10.1142/S0218202509004121 doi: 10.1142/S0218202509004121
    [33] P. Fernandes, M. Ottonello, M. Raffetto, Regularity of time-harmonic electromagnetic fields in the interior of bianisotropic materials and metamaterials, IMA J. Appl. Math., 79 (2014), 54–93. https://doi.org/10.1093/imamat/hxs039 doi: 10.1093/imamat/hxs039
    [34] P. Cocquet, P. Mazet, V. Mouysset, On the existence and uniqueness of a solution for some frequency-dependent partial differential equations coming from the modeling of metamaterials, SIAM J. Math. Anal., 44 (2012), 3806–3833. https://doi.org/10.1137/100810071 doi: 10.1137/100810071
    [35] E. T. Chung, P. Ciarlet Jr., A staggered discontinuous Galerkin method for wave propagation in media with dielectrics and metamaterials, J. Comput. Appl. Math., 239 (2013), 189–207. https://doi.org/10.1016/j.cam.2012.09.033 doi: 10.1016/j.cam.2012.09.033
    [36] M. Cassier, P. Joly, M. Kachanovska, Mathematical models for dispersive electromagnetic waves: An overview, Comput. Math. Appl., 74 (2017), 2792–2830. https://doi.org/10.1016/j.camwa.2017.07.025 doi: 10.1016/j.camwa.2017.07.025
    [37] J. Li, A literature survey of mathematical study of metamaterials, Int. J. Numer. Anal. Model., 13 (2016), 230–243.
    [38] U. Leonhardt, Optical conformal mapping, Science, 312 (2006), 1777–1780. https://doi.org/10.1126/science.1126493
    [39] J. B. Pendry, D. Schurig, D. R. Smith, Controlling electromagnetic fields, Science, 312 (2006), 1780–1782. https://doi.org/10.1126/science.1125907 doi: 10.1126/science.1125907
    [40] J. Li, Y. Huang, Time-Domain Finite Element Methods for Maxwell's Equations in Metamaterials, Springer Series in Computational Mathematics (vol.43), Springer, 2013.
    [41] H. Ammari, H. Kang, H. Lee, M. Lim, S. Yu, Enhancement of near cloaking for the full Maxwell equations, SIAM J. Appl. Math., 73 (2013), 2055–2076. https://doi.org/10.1137/120903610 doi: 10.1137/120903610
    [42] G. Bao, H. Liu, J. Zou, Nearly cloaking the full Maxwell equations: Cloaking active contents with general conducting layers, J. Math. Pures Appl., 101 (2014), 716–733. https://doi.org/10.1016/j.matpur.2013.10.010 doi: 10.1016/j.matpur.2013.10.010
    [43] H. Ammari, J. Garnier, V. Jugnon, H. Kang, H. Lee, M. Lim, Enhancement of near-cloaking. Part III: Numerical simulations, statistical stability, and related questions, Comtemp. Math., 577 (2012), 1–24.
    [44] H. Ammari, G. Ciraolo, H. Kang, H. Lee, G. W. Milton, Spectral theory of a Neumann-Poincaré-type operator and analysis of cloaking due to anomalous localized resonance, Arch. Ration. Mech. Anal., 208 (2013), 667–692. https://doi.org/10.1007/s00205-012-0605-5 doi: 10.1007/s00205-012-0605-5
    [45] R. V. Kohn, J. Lu, B. Schweizer, M. I. Weinstein, A variational perspective on cloaking by anomalous localized resonance, Commun. Math. Phys., 328 (2014), 1–27. https://doi.org/10.1007/s00220-014-1943-y doi: 10.1007/s00220-014-1943-y
    [46] R. V. Kohn, D. Onofrei, M. S. Vogelius, M. I. Weinstein, Cloaking via change of variables for the Helmholtz equation, Commun. Pure Appl. Math., 63 (2010), 973–1016. https://doi.org/10.1002/cpa.20326 doi: 10.1002/cpa.20326
    [47] A. Greenleaf, Y. Kurylev, M. Lassas, G. Uhlmann, Cloaking devices, electromagnetics wormholes and transformation optics, SIAM Rev., 51 (2009), 3–33. https://doi.org/10.1137/080716827 doi: 10.1137/080716827
    [48] F. Guevara Vasquez, G. W. Milton, D. Onofrei, Broadband exterior cloaking, Opt. Express, 17 (2009), 14800–14805. https://doi.org/10.1364/OE.17.014800
    [49] M. Lassas, M. Salo, L. Tzou, Inverse problems and invisibility cloaking for FEM models and resistor networks, Math. Mod. Meth. Appl. Sci., 25 (2015), 309–342. https://doi.org/10.1142/S0218202515500116 doi: 10.1142/S0218202515500116
    [50] S. C. Brenner, J. Gedicke, L. Y. Sung, An adaptive finite element method for two-dimensional transverse magnetic time harmonic Maxwell's equations with general material properties and general boundary conditions, J. Sci. Comput., 68 (2016), 848–863. https://doi.org/10.1007/s10915-015-0161-x doi: 10.1007/s10915-015-0161-x
    [51] J. J. Lee, A mixed method for time-transient acoustic wave propagation in metamaterials, J. Sci. Comput., 84 (2020). https://doi.org/10.1007/s10915-020-01275-0
    [52] S. Nicaise, J. Venel, A posteriori error estimates for a finite element approximation of transmission problems with sign changing coefficients, J. Comput. Appl. Math., 235 (2011), 4272–4282. https://doi.org/10.1016/j.cam.2011.03.028 doi: 10.1016/j.cam.2011.03.028
    [53] Z. Yang, L. L. Wang, Accurate simulation of ideal circular and elliptic cylindrical invisibility cloaks, Commun. Comput. Phys., 17 (2015), 822–849. https://doi.org/10.4208/cicp.280514.131014a doi: 10.4208/cicp.280514.131014a
    [54] Z. Yang, L. L. Wang, Z. Rong, B. Wang, B. Zhang, Seamless integration of global Dirichlet-to-Neumann boundary condition and spectral elements for transformation electromagnetics, Comput. Methods Appl. Mech. Engrg., 301 (2016), 137–163. https://doi.org/10.1016/j.cma.2015.12.020 doi: 10.1016/j.cma.2015.12.020
    [55] B. Wang, Z. Yang, L. L. Wang, S. Jiang, On time-domain NRBC for Maxwell's equations and its application in accurate simulation of electromagnetic invisibility cloaks, J. Sci. Comput., 86 (2021). https://doi.org/10.1007/s10915-020-01354-2
    [56] U. Leonhardt, T. Tyc, Broadband invisibility by non-Euclidean cloaking, Science, 323 (2009), 110–112. https://doi.org/10.1126/science.1166332 doi: 10.1126/science.1166332
    [57] R. Liu, C. Ji, J. J. Mock, J. Y. Chin, T. J. Cui, D. R. Smith, Broadband ground-plane cloak, Science, 323 (2009), 366–369. https://doi.org/10.1126/science.1166949
    [58] J. Li, Two new finite element schemes and their analysis for modeling of wave propagation in graphene, Results Appl. Math., 9 (2021), 100136. https://doi.org/10.1016/j.rinam.2020.100136 doi: 10.1016/j.rinam.2020.100136
    [59] Y. Wu, J. Li, Total reflection and cloaking by zero index metamaterials loaded with rectangular dielectric defects, Appl. Phys. Lett., 102 (2013), 183105. https://doi.org/10.1063/1.4804201 doi: 10.1063/1.4804201
    [60] J. Li, Well-posedness study for a time-domain spherical cloaking model, Comput. Math. Appl., 68 (2014), 1871–1881. https://doi.org/10.1016/j.camwa.2014.10.007 doi: 10.1016/j.camwa.2014.10.007
    [61] J. Li, Y. Huang, W. Yang, Well-posedness study and finite element simulation of time-domain cylindrical and elliptical cloaks, Math. Comput., 84 (2015), 543–562. https://doi.org/10.1090/s0025-5718-2014-02911-6 doi: 10.1090/s0025-5718-2014-02911-6
    [62] W. Yang, J. Li, Y. Huang, Mathematical analysis and finite element time domain simulation of arbitrary star-shaped electromagnetic cloaks, SIAM J. Numer. Anal., 56 (2018), 136–159. https://doi.org/10.1137/16M1093835 doi: 10.1137/16M1093835
    [63] J. Li, J. B. Pendry, Hiding under the carpet: a new strategy for cloaking, Phys. Rev. Lett., 101 (2008), 2039014. https://doi.org/10.1103/PhysRevLett.101.203901 doi: 10.1103/PhysRevLett.101.203901
    [64] J. Li, Y. Huang, W. Yang, A. Wood, Mathematical analysis and time-domain finite element simulation of carpet cloak, SIAM J. Appl. Math., 74 (2014), 1136–1151. https://doi.org/10.1137/140959250 doi: 10.1137/140959250
    [65] J. Li, C. W. Shu, W. Yang, Development and analysis of two new finite element schemes for a time-domain carpet cloak model, Adv. Comput. Math., 48 (2022), 24. https://doi.org/10.1007/s10444-022-09948-0 doi: 10.1007/s10444-022-09948-0
    [66] J. Li, Z. Liang, J. Zhu, X. Zhang, Anisotropic metamaterials for transformation acoustics and imaging, in Acoustic Metamaterials: Negative Refraction, Imaging, Sensing and Cloaking (eds. R. V. Craster and S. Guenneau), Springer Series in Materials Science, 166 (2013), 169–195. https://doi.org/10.1007/978-94-007-4813-2_7
    [67] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, et al., Electric field effect in atomically thin carbon films, Science, 306 (2004), 666–669. https://doi.org/10.1126/science.1102896 doi: 10.1126/science.1102896
    [68] F. Bonaccorso, Z. Sun, T. Hasan, A. C. Ferrari, Graphene photonics and optoelectronics, Nat. Photonics, 4 (2010), 611–622. https://doi.org/10.1038/nphoton.2010.186 doi: 10.1038/nphoton.2010.186
    [69] F. H. L. Koppens, D. E. Chang, F. J. Garcia de Abajo, Graphene plasmonics: A platform for strong light-matter interactions, Nano Lett., 11 (2011), 3370–3377. https://doi.org/10.1021/nl201771h doi: 10.1021/nl201771h
    [70] S. K. Tiwari, S. Sahoo, N. Wang, A. Huczko, Graphene research and their outputs: Status and prospect, J. Sci.: Adv. Mater. Devices, 5 (2020), 10–29. https://doi.org/10.1016/j.jsamd.2020.01.006 doi: 10.1016/j.jsamd.2020.01.006
    [71] G. Bal, P. Cazeaux, D. Massatt, S. Quinn, Mathematical models of topologically protected transport in twisted bilayer graphene, Multiscale Model. Simul., 21 (2023), 1081–1121. https://doi.org/10.1137/22M1505542 doi: 10.1137/22M1505542
    [72] Y. Hong, D. P. Nicholls, On the consistent choice of effective permittivity and conductivity for modeling graphene, JOSA A, 38 (2021), 1511–1520. https://doi.org/10.1364/JOSAA.430088 doi: 10.1364/JOSAA.430088
    [73] J. P. Lee-Thorp, M. I. Weinstein, Y. Zhu, Elliptic operators with honeycomb symmetry: Dirac points, edge states and applications to photonic graphene, Arch. Ration. Mech. Anal., 232 (2019), 1–63. https://doi.org/10.1007/s00205-018-1315-4 doi: 10.1007/s00205-018-1315-4
    [74] M. Maier, D. Margetis, M. Luskin, Dipole excitation of surface plasmon on a conducting sheet: finite element approximation and validation, J. Comput. Phys., 339 (2017), 126–145. https://doi.org/10.1016/j.jcp.2017.03.014 doi: 10.1016/j.jcp.2017.03.014
    [75] M. Maier, D. Margetis, M. Luskin, Generation of surface plasmon-polaritons by edge effects, Commun. Math. Sci., 16 (2018), 77–95.
    [76] J. H. Song, M. Maier, M. Luskin, Adaptive finite element simulations of waveguide configurations involving parallel 2D material sheets, Comput. Methods Appl. Mech. Eng., 351 (2019), 20–34. https://doi.org/10.1016/j.cma.2019.03.039 doi: 10.1016/j.cma.2019.03.039
    [77] J. Wilson, F. Santosa, P. A. Martin, Temporally manipulated plasmons on graphene, SIAM J. Appl. Math., 79 (2019), 1051–1074. https://doi.org/10.1137/18M1226889 doi: 10.1137/18M1226889
    [78] W. Yang, J. Li, Y. Huang, Time-domain finite element method and analysis for modeling of surface plasmon polaritons in graphene devices, Comput. Methods Appl. Mech. Eng., 372 (2020), 113349. https://doi.org/10.1016/j.cma.2020.113349 doi: 10.1016/j.cma.2020.113349
    [79] J. Li, L. Zhu, T. Arbogast, A new time-domain finite element method for simulation surface plasmon polaritons on graphene sheets, Comput. Math. Appl., 142 (2023), 268–383. https://doi.org/10.1016/j.camwa.2023.05.003 doi: 10.1016/j.camwa.2023.05.003
    [80] Y. Gong, N. Liu, Advanced numerical methods for graphene simulation with equivalent boundary conditions: a review, Photonics, 10 (2023), 712. https://doi.org/10.3390/photonics10070712 doi: 10.3390/photonics10070712
    [81] P. Li, L. J. Jiang, H. Bağci, Discontinuous Galerkin time-domain modeling of graphene nanoribbon incorporating the spatial dispersion effects, IEEE Trans. Antennas Propag., 66 (2018), 3590–3598. https://doi.org/10.1109/TAP.2018.2826567 doi: 10.1109/TAP.2018.2826567
    [82] L. Yang, J. Tian, K. Z. Rajab, Y. Hao, FDTD modeling of nonlinear phenomena in wave transmission through graphene, IEEE Antennas Wirel. Propag. Lett., 17 (2018), 126–129. https://doi.org/10.1109/LAWP.2017.2777530 doi: 10.1109/LAWP.2017.2777530
    [83] B. Alpert, L. Greengard, T. Hagstrom, Rapid evaluation of nonreflecting boundary kernels for time-domain wave propagation, SIAM J. Numer. Anal., 37 (2000), 1138–1164. https://doi.org/10.1137/S0036142998336916 doi: 10.1137/S0036142998336916
    [84] B. Engquist, A. Majda, Absorbing boundary conditions for the numerical simulation of waves, Math. Comput., 31 (1977), 629–651. https://doi.org/10.1073/pnas.74.5.1765 doi: 10.1073/pnas.74.5.1765
    [85] T. Hagstrom, T. Warburton, D. Givoli, Radiation boundary conditions for time-dependent waves based on complete plane wave expansions, J. Comput. Appl. Math., 234 (2010), 1988–1995. https://doi.org/10.1016/j.cam.2009.08.050 doi: 10.1016/j.cam.2009.08.050
    [86] J. P. Bérenger, A perfectly matched layer for the absorption of electromagnetic waves, J. Comput. Phys., 114 (1994), 185–200. https://doi.org/10.1006/jcph.1994.1159 doi: 10.1006/jcph.1994.1159
    [87] R. W. Ziolkowski, Maxwellian material based absorbing boundary conditions, Comput. Methods Appl. Mech. Eng., 169 (1999), 237–262. https://doi.org/10.1016/S0045-7825(98)00156-X doi: 10.1016/S0045-7825(98)00156-X
    [88] Y. Huang, J. Li, Z. Fang, Mathematical analysis of Ziolkowski's PML model with application for wave propagation in metamaterials, J. Comp. Appl. Math., 366 (2020), 112434. https://doi.org/10.1016/j.cam.2019.112434 doi: 10.1016/j.cam.2019.112434
    [89] J. Li, L. Zhu, Analysis and application of two novel finite element methods for solving Ziolkowski's PML model in the integro-differential form, SIAM J. Numer. Anal., 61 (2023), 2209–2236. https://doi.org/10.1137/22M1506936 doi: 10.1137/22M1506936
    [90] G. C. Cohen, P. Monk, Mur-Nédélec finite element schemes for Maxwell's equations, Comput. Methods Appl. Mech. Eng., 169 (1999), 197–217. https://doi.org/10.1016/S0045-7825(98)00154-6 doi: 10.1016/S0045-7825(98)00154-6
    [91] M. Chen, Y. Huang, J. Li, Development and analysis of a new finite element method for the Cohen-Monk PML model, Numer. Math., 147 (2021), 127–155. https://doi.org/10.1007/s00211-020-01166-4 doi: 10.1007/s00211-020-01166-4
    [92] J. L. Lions, J. Métral, O. Vacus, Well-posed absorbing layer for hyperbolic problems, Numer. Math., 92 (2002), 535–562. https://doi.org/10.1007/s002110100263 doi: 10.1007/s002110100263
    [93] G. Bao, H. Wu, Convergence analysis of the perfectly matched layer problems for time-harmonic Maxwell's equations, SIAM J. Numer. Anal., 43 (2005), 2121–2143. https://doi.org/10.1137/040604315 doi: 10.1137/040604315
    [94] J. H. Bramble, J. E. Pasciak, Analysis of a finite PML approximation for the three dimensional time-harmonic Maxwell and acoustic scattering problems, Math. Comput., 76 (2007), 597–614. https://doi.org/10.1090/S0025-5718-06-01930-2 doi: 10.1090/S0025-5718-06-01930-2
    [95] J. Chen, Z. Chen, An adaptive perfectly matched layer technique for 3-D time-harmonic electromagnetic scattering problems, Math. Comput., 77 (2007), 673–698. https://doi.org/10.1090/S0025-5718-07-02055-8 doi: 10.1090/S0025-5718-07-02055-8
    [96] Z. Chen, W. Zheng, PML method for electromagnetic scattering problem in a two-layer medium, SIAM J. Numer. Anal., 55 (2017), 2050–2084. https://doi.org/10.1137/16M1091757 doi: 10.1137/16M1091757
    [97] T. Hohage, F. Schmidt, L. Zschiedrich, Solving time-harmonic scattering problems based on the pole condition II: convergence of the PML method, SIAM J. Math. Anal., 35 (2003), 547–560. https://doi.org/10.1137/S0036141002406485 doi: 10.1137/S0036141002406485
    [98] M. Lassas, E. Somersalo, On the existence and convergence of the solution of PML equations, Computing, 60 (1998), 229–241. https://doi.org/10.1007/BF02684334 doi: 10.1007/BF02684334
    [99] J. H. Bramble, J. E. Pasciak, Analysis of a finite element PML approximation for the three dimensional time-harmonic Maxwell problem, Math. Comput., 77 (2008), 1–10. https://doi.org/10.1090/S0025-5718-07-02037-6 doi: 10.1090/S0025-5718-07-02037-6
    [100] Z. Chen, W. Zheng, Convergence of the uniaxial perfectly matched layer method for time-harmonic scattering problems in two-layered media, SIAM J. Numer. Anal., 48 (2010), 2158–2185. https://doi.org/10.1137/090750603 doi: 10.1137/090750603
    [101] Z. Chen, T. Cui, L. Zhang, An adaptive anisotropic perfectly matched layer method for 3-D time harmonic electromagnetic scattering problems, Numer. Math., 125 (2013), 639–677. https://doi.org/10.1007/s00211-013-0550-8 doi: 10.1007/s00211-013-0550-8
    [102] Y. Lin, K. Zhang, J. Zou, Studies on some perfectly matched layers for one-dimensional time-dependent systems, Adv. Comput. Math., 30 (2009), 1–35. https://doi.org/10.1007/s10444-007-9055-2 doi: 10.1007/s10444-007-9055-2
    [103] Y. Gao, P. Li, Electromagnetic scattering for time-domain Maxwell's equations in an unbounded structure, Math. Models Methods Appl. Sci., 27 (2017), 1843–1870. https://doi.org/10.1142/S0218202517500336 doi: 10.1142/S0218202517500336
    [104] P. Li, L. Wang, A. Wood, Analysis of transient electromagnetic scattering from a three-dimensional open cavity, SIAM J. Appl. Math., 75 (2015), 1675–1699. https://doi.org/10.1137/140989637 doi: 10.1137/140989637
    [105] C. Wei, J. Yang, B. Zhang, Convergence analysis of the PML method for time-domain electromagnetic scattering problems, SIAM J. Numer. Anal., 58 (2020), 1918–1940. https://doi.org/10.1137/19M126517X doi: 10.1137/19M126517X
    [106] C. Wei, J. Yang, B. Zhang, Convergence of the uniaxial PML method for time-domain electromagnetic scattering problems, ESAIM Math. Model. Numer. Anal., 55 (2021), 2421–2443. https://doi.org/10.1051/m2an/2021064 doi: 10.1051/m2an/2021064
    [107] F. L. Teixeira, W. C. Chew, Advances in the theory of perfectly matched layers, in Fast and Efficient Algorithms in Computational Electromagnetics, Artech House, Boston, 7 (2001), 283–346.
    [108] W. C. Chew, W. H. Weedon, A 3D perfectly matched medium from modified Maxwell's equations with stretched coordinates, Microwave Opt. Technol. Lett., 7 (1994), 599–604. https://doi.org/10.1002/mop.4650071304 doi: 10.1002/mop.4650071304
    [109] F. Collino, P. Monk, Optimizing the perfectly matched layer, Comput. Methods Appl. Mech. Eng., 164 (1998), 157–171. https://doi.org/10.1016/S0045-7825(98)00052-8 doi: 10.1016/S0045-7825(98)00052-8
    [110] T. Lu, P. Zhang, W. Cai, Discontinuous Galerkin methods for dispersive and lossy Maxwell's equations and PML boundary conditions, J. Comput. Phys., 200 (2004), 549–580.
    [111] J. L. Hong, L. H. Ji, L. H. Kong, Energy-dissipation splitting finite-difference time-domain method for Maxwell equations with perfectly matched layers, J. Comput. Phys., 269 (2014), 201–214. https://doi.org/10.1016/j.jcp.2014.03.025 doi: 10.1016/j.jcp.2014.03.025
    [112] J. S. Hesthaven, On the analysis and construction of perfectly matched layers for the linearized Euler equations, J. Comput. Phys., 142 (1998), 129–147. https://doi.org/10.1006/jcph.1998.5938 doi: 10.1006/jcph.1998.5938
    [113] F. Q. Hu, On absorbing boundary conditions for linearized euler equations by a perfectly matched layer, J. Comput. Phys., 129 (1996), 201–219. https://doi.org/10.1006/jcph.1996.0244 doi: 10.1006/jcph.1996.0244
    [114] F. Pled, C. Desceliers, Review and recent developments on the perfectly matched layer (PML) method for the numerical modeling and simulation of elastic wave propagation in unbounded domains, Arch. Comput. Methods Eng., 29 (2022), 471–518. https://doi.org/10.1007/s11831-021-09581-y doi: 10.1007/s11831-021-09581-y
    [115] D. Appelö, T. Hagstrom, G. Kreiss, Perfectly matched layers for hyperbolic systems: general formulation, well-posedness, and stability, SIAM J. Appl. Math., 67 (2006), 1–23. https://doi.org/10.1137/050639107 doi: 10.1137/050639107
    [116] D. Appelö, T. Hagstrom, A general perfectly matched layer model for hyperbolic-parabolic systems, SIAM J. Sci. Comput., 31 (2009), 3301–23. https://doi.org/10.1137/080713951 doi: 10.1137/080713951
    [117] D. H. Baffet, M. J. Grote, S. Imperiale, M. Kachanovska, Energy decay and stability of a perfectly matched layer for the wave equation, J. Sci. Comput., 81 (2019), 2237–2270. https://doi.org/10.1007/s10915-019-01089-9 doi: 10.1007/s10915-019-01089-9
    [118] Z. Yang, L. L. Wang, Y. Gao, A truly exact perfect absorbing layer for time-harmonic acoustic wave scattering problems, SIAM J. Sci. Comput., 43 (2021), A1027–A1061. https://doi.org/10.1137/19M1294071 doi: 10.1137/19M1294071
    [119] W. C. Chew, Q. H. Liu, Perfectly matched layers for elastodynamics: a new absorbing boundary condition, J. Comput. Acoust., 4 (1996), 341–349. https://doi.org/10.1142/S0218396X96000118 doi: 10.1142/S0218396X96000118
    [120] F. Collino, C. Tsogka, Application of the PML absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media, Geophysics, 88 (2001), 43–73.
    [121] K. Duru, L. Rannabauer, A. A. Gabriel, G. Kreiss, M. Bader, A stable discontinuous Galerkin method for the perfectly matched layer for elastodynamics in first order form, Numer. Math., 146 (2020), 729–782. https://doi.org/10.1007/s00211-020-01160-w doi: 10.1007/s00211-020-01160-w
    [122] S. Abarbanel, D. Gottlieb, A mathematical analysis of the PML method, J. Comput. Phys., 134 (1997), 357–363. https://doi.org/10.1006/jcph.1997.5717 doi: 10.1006/jcph.1997.5717
    [123] E. Bécache, P. Joly, On the analysis of Bérenger's perfectly matched layers for maxwell's equations, ESAIM: M2AN, 36 (2002), 87–119. https://doi.org/10.1051/m2an:2002004 doi: 10.1051/m2an:2002004
    [124] L. Zhao, A. C. Cangellaris, A general approach for the development of unsplit-field time-domain implementations of perfectly matched layers for FDTD grid truncation, IEEE Microwave Guid Wave Lett., 6 (1996), 209–211. https://doi.org/10.1109/75.491508 doi: 10.1109/75.491508
    [125] Y. Huang, M. Chen, J. Li, Development and analysis of both finite element and fourth-order in space finite difference methods for an equivalent Berenger's PML model, J. Comput. Phys., 405 (2020), 109154. https://doi.org/10.1016/j.jcp.2019.109154 doi: 10.1016/j.jcp.2019.109154
    [126] Y. Huang, J. Li, X. Liu, Developing and analyzing an explicit unconditionally stable finite element scheme for an equivalent Bérenger's PML model, ESAIM: M2AN, 57 (2023), 621–644. https://doi.org/10.1051/m2an/2022086 doi: 10.1051/m2an/2022086
    [127] D. Correia, J. M. Jin, 3D-FDTD-PML analysis of left-handed metamaterials, Microwave Opt. Technol. Lett., 40 (2004), 201–205. https://doi.org/10.1002/mop.11328 doi: 10.1002/mop.11328
    [128] S. A. Cummer, Perfectly matched layer behavior in negative refractive index materials, IEEE Antennas Wirel. Propag. Lett., 3 (2004), 172–175. https://doi.org/10.1109/LAWP.2004.833710 doi: 10.1109/LAWP.2004.833710
    [129] X. T. Dong, X. S. Rao, Y. B. Gan, B. Guo, W. Y. Yin, Perfectly matched layer-absorbing boundary condition for left-handed materials, IEEE Microwave Wireless Compon. Lett., 14 (2004), 301–303. https://doi.org/10.1109/LMWC.2004.827104 doi: 10.1109/LMWC.2004.827104
    [130] P. R. Loh, A. F. Oskooi, M. Ibanescu, M. Skorobogatiy, S. G. Johnson, Fundamental relation between phase and group velocity, and application to the failure of perfectly matched layers in backward-wave structures, Phys. Rev. E, 79 (2009), 065601. https://doi.org/10.1103/PhysRevE.79.065601 doi: 10.1103/PhysRevE.79.065601
    [131] Y. Shi, Y. Li, C. H. Liang, Perfectly matched layer absorbing boundary condition for truncating the boundary of the left-handed medium, Microwave Opt. Technol. Lett., 48 (2006), 57–63. https://doi.org/10.1002/mop.21260 doi: 10.1002/mop.21260
    [132] E. Bécache, P. Joly, M. Kachanovska, V. Vinoles, Perfectly matched layers in negative index metamaterials and plasmas, ESAIM: Proc. Surv., 50 (2015), 113–132. https://doi.org/10.1051/proc/201550006 doi: 10.1051/proc/201550006
    [133] E. Bécache, P. Joly, M. Kachanovska, V. Vinoles, On the analysis of perfectly matched layers for a class of dispersive media and application to negative index metamaterials, Math. Comput., 87 (2018), 2775–2810. https://doi.org/10.1090/mcom/3307 doi: 10.1090/mcom/3307
    [134] J. Li, L. Zhu, Analysis and FDTD simulation of a perfectly matched layer for the Drude metamaterial, Ann. Appl. Math., 38 (2022), 1–23.
    [135] Y. Huang, J. Li, X. Yi, H. Zhao, Analysis and application of a time-domain finite element method for the Drude metamaterial perfectly matched layer model, J. Comput. Appl. Math., 438 (2024), 115575. https://doi.org/10.1016/j.cam.2023.115575 doi: 10.1016/j.cam.2023.115575
    [136] E. Bécache, M. Kachanovska, Stable perfectly matched layers for a class of anisotropic dispersive models, Part I: necessary and sufficient conditions of stability, ESAIM: M2AN, 51 (2017), 2399–2434. https://doi.org/10.1051/m2an/2017019 doi: 10.1051/m2an/2017019
    [137] C. L. Holloway, E. F. Kuester, J. A. Gordon, J. O'Hara, J. Booth, D.R. Smith, An overview of the theory and applications of metasurfaces: The two-dimensional equivalents of metamaterials, IEEE Antennas Propag. Mag., 54 (2012), 10–35. https://doi.org/10.1109/MAP.2012.6230714 doi: 10.1109/MAP.2012.6230714
    [138] K. Achouri, C. Caloz, Design, concepts, and applications of electromagnetic metasurfaces, Nanophotonics, 7 (2018), 1095–1116. https://doi.org/10.1515/nanoph-2017-0119 doi: 10.1515/nanoph-2017-0119
    [139] S. B. Glybovski, S. A. Tretyakov, P. A. Belov, Y. S. Kivshar, C. R. Simovski, Metasurfaces: From microwaves to visible, Phys. Rep., 634 (2016), 1–72. https://doi.org/10.1016/j.physrep.2016.04.004 doi: 10.1016/j.physrep.2016.04.004
    [140] K. A. Lurie, An Introduction to the Mathematical Theory of Dynamic Materials, Springer, 2017.
    [141] P. A. Huidobro, E. Galiffi, S. Guenneau, R. V. Craster, J. B. Pendry, Fresnel drag in space-time-modulated metamaterials, Proc. Natl. Acad. Sci. U.S.A., 116 (2019), 24943–24948. https://doi.org/10.1073/pnas.1915027116 doi: 10.1073/pnas.1915027116
    [142] C. Caloz, Z. Deck-Léger, Spacetime metamaterials–Part II: theory and applications, IEEE Trans. Antennas Propag., 68 (2020), 1583–1598. https://doi.org/10.1109/TAP.2019.2944216 doi: 10.1109/TAP.2019.2944216
    [143] E. Galiffi, R. Tirole, S. Yin, H. Li, S. Vezzoli, P. A. Huidobro, et al., Photonics of time-varying media, Adv. Photonics, 4 (2022), 014002–014002. https://doi.org/10.1117/1.AP.4.1.014002 doi: 10.1117/1.AP.4.1.014002
    [144] N. Engheta, Four-dimensional optics using time-varying metamaterials, Science, 379 (2023), 1190–1191. https://doi.org/10.1126/science.adf1094 doi: 10.1126/science.adf1094
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1772) PDF downloads(148) Cited by(0)

Figures and Tables

Figures(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog