arXiv:1509.01519v3 [math.AC] 4 May 2017

THE SUPPORT OF LOCAL COHOMOLOGY MODULES
MORDECHAI KATZMAN AND WENLIANG ZHANG
Abstract. We describe the support of F -finite F -modules over polynomial rings R of prime characteristic. Our description yields an algorithm to compute the support of such modules; the complexity of our
algorithm is also analyzed. To the best of our knowledge, this is the
first algorithm to avoid extensive use of Gröbner bases and hence of
substantial practical value. We also use the idea behind this algorithm
to prove that the support of HIj (S) is Zariski closed for each ideal I of
S where R is noetherian commutative ring of prime characteristic with
finitely many isolated singular points and S = R/gR (g ∈ R).

1. Introduction
Local cohomology is a powerful tool introduced by Alexander Grothendieck
in the 1960’s ([Har67]) and it has since yielded many geometric and algebraic insights. From an algebraic point of view, given an ideal I in
a commutative ring R, local cohomology modules HiI (−) (i ≥ 0) arise
as right-derived functors of the torsion functor on R-modules given by
ΓI (M ) = {a ∈ M | I k a = 0 for some k ≥ 0}. A central question in the
theory of local cohomology is to determine for which values of i does the
local cohomology module HiI (M ) vanish. This question is both useful and
difficult even in the case where R is a regular local ring and M = R, and this
case has been studied intensely since the introduction of local cohomology
(e.g., cf. [Har68], [PS73] and [Ogu73]).
The aim of this paper is to describe the support of local cohomology
modules in prime characteristic. Specifically, we first study the support of
F -finite F -modules over polynomial rings R and show a computationally
feasible method for computing these without the need to compute generating roots. To the best of our knowledge, this is the first computationally feasible algorithm for calculating the support of these modules in
prime characteristic. We then apply this to the calculation of supports
of local
cohomology
modules
and of iterated local cohomology modules



in
i2
i1
thus, for example, giving an effective method
HI1 HI2 . . . HIn (R) . . .
for determining the vanishing of Lyubeznik numbers.
2010 Mathematics Subject Classification. 13D45, 13A35.
Key words and phrases. local cohomology, prime characteristic.
M.K. gratefully acknowledges support from EPSRC grant EP/J005436/1. W.Z. is
partially supported by NSF grants DMS #1405602/#1606414.
1

2

MORDECHAI KATZMAN AND WENLIANG ZHANG

Our methods are interesting both from theoretical and practical points
of view. A careful analysis of the algorithms resulting from these methods
(see Section 4 below) shows that
(a) the degrees of the polynomials appearing in the calculations have a
low upper bound, and, furthermore,
(b) when the method is applied to the calculation of supports of local
cohomology modules, if the input is given by polynomials with integer coefficients, then the calculation of supports modulo different
primes p involves polynomials whose degrees can be bounded from
above by a constant times p, that constant being independent of p.
In [Lyu97]) Gennady Lyubeznik described an algorithm for computing
the support of F -finite F -modules. That algorithm requires the calculation
for roots of these modules, and this relies on the repeated calculation of
Grobner bases; these are often too complex to be computed in practice.
Our algorithm consists of an iterative procedure (as described in section 3)
which produces a quotient of a finite-rank free module with the same support
as the given F -finite F -module. To find the support itself one needs to find
a presentation for this submodule as a cokernel of a matrix: the support is
then defined by the ideal of maximal minors of that matrix.
Crucially, the iterative procedure above does not require the calculation
of Gröbner bases, and consists essentially of matrix multiplications together
with the listing of terms of polynomials whose degrees are bounded by a
constant (independent of p) times p. The final step of the algorithm, finding
a presentation of a finitely generated module, requires the calculation of one
module of syzygies, hence the calculation of one Gröber basis.
It is this that makes our algorithm a practical tool for computing supports
of F -finite F -modules.1
The reason why we are able to compute and analyze in characteristic p the
support of F -finite F -modules is the existence of the eth iterated Frobenius
e
endomorphism f e : R → R, taking a ∈ R to ap (e ≥ 0). The usefulness of
these lies in the fact that given an R-module M , we may endow it with a new
R-module structure via f e : let F∗e M denote the additive Abelian group M
denoting its elements {F∗e m | m ∈ M }, and endow F∗e M with the R-module
e
structure is given by aF∗e m = F∗e ap m for all a ∈ R and m ∈ M .
This also allows us to define the eth Frobenius functors from the category
of R-modules to itself given by FRe (M ) = F∗e R ⊗R M and viewing this as
a R-module via the identification of F∗e R with R: the resulting R-module
structure on FRe (M ) satisfies a(F∗e b ⊗ m) = F∗e ab ⊗ m and F∗e ap b ⊗ m =
F∗e b ⊗ am for all a, b ∈ R and m ∈ M .
We will be interested in this construction mainly for regular rings and
henceforth in this paper R will denote a regular ring of characteristic p > 0.
1The various algorithms in this paper have been incorporated in the “FSing” package

of Macaulay 2[GS].

THE SUPPORT OF LOCAL COHOMOLOGY MODULES

3

Recall that an F -finite FR -module M is an R-module obtained as a direct
limit of a direct limit system of the form
U

F 2 (U )

F 1 (U )

M−
→ FR1 (M ) −−R−−→ FR2 (M ) −−R−−→ . . .
where M is a finitely generated module and U is an R-linear map (cf. [Lyu97]).
The main interest in F -finite F -modules follows from the fact that local cohomology modules are F -finite F -modules, as we now explain.
The jth local cohomology module of M with support on an ideal I ⊂ R is
defined as
(1)

e

HjI (M ) = lim
ExtjR (R/I [p ] , M )
→
e

where maps in the direct limit system are induced by the surjections R/I [p
e
R/I [p ] . If we apply this with M = R, we obtain

e+1 ]

e

HjI (R) = lim
ExtjR (R/I [p ] , R)
→
e

∼
ExtjR (FRe (R/I), FRe R)
= lim
→
e


j
e
∼
Ext
(R/I,
R)
F
= lim
R
R
→
e

e
where we use the facts that FRe (R) ∼
= R, FRe (R/I) ∼
= R/I [p ] , and that,
e
since R is regular, the Frobenius functor FR (−) is exact and thus commutes
with the computation of cohomology. This shows that HjI (R) are F -finite
F -modules, and we may apply our F -finite F -module machinery to them.
Finally, in section 7 we turn our attention to hypersurfaces and describe
the support of their local cohomology modules, which turn out to be closed.2
Given a fixed g ∈ R, one can ask for the locus of primes P ⊆ R for which
g
→ HiI (RP ) is injective and the locus
the multiplication by g map HiI (RP ) −
of primes for which this is surjective. We show that these two loci are
Zariski closed by describing explicitly the defining ideals of these loci, and
we use these to describe the defining ideal of the (Zariski closed) support for
HiI (R/gR). We also extend the Zariski-closedness of HiI (R/gR) to the case
when R has finitely many isolated singular points.
The methods used for the various calculations in this paper are described
in section 2.

2. Prime characteristic tools
Definition 2.1. Let e ≥ 0. Let T be a commutative ring of prime characteristic p.
2The fact that the support is closed was simultaneously and independently also discov-

ered by Mel Hochster and Luis Núñez-Betancourt in [HNB] using a different method.

→

4

MORDECHAI KATZMAN AND WENLIANG ZHANG
e

(a) Given any matrix (or vector) A with entries in T , we define A[p ]
to be the matrix obtained from A by raising its entries to the pe th
power.
e
(b) Given any submodule K ⊆ T α , we define K [p ] to be the R-submodule
e
of T α generated by {v [p ] | v ∈ K}.
Henceforth in this section, T will denote a regular ring with the property
that F∗e T are intersection flat T -modules for all e ≥ 0, i.e., for any family of
T -modules {Mλ }λ∈Λ ,
\
\
F∗e T ⊗T
Mλ =
F∗e T ⊗T Mλ .
λ∈Λ

λ∈Λ

F∗e T

are free T -modules (e. g. , polyThese include rings T for which
nomial rings and power series rings with F -finite coefficient rings,) and
also all complete regular rings (cf. [Kat08, Proposition 5.3]). These rings
have that property that for any collection of submodules {Lλ }λ∈Λ of T α ,
[pe ]
T
T
[pe ]
= λ∈Λ Lλ : indeed, the regularity of T implies that for
λ∈Λ Lλ
e
any submodule L ⊆ T α , L[p ] can be identified with FTe (L) and and the
intersection-flatness of F∗e T implies
\
\
\
\
FTe (
Lλ ) = F∗e T ⊗T
Lλ =
F∗e T ⊗T Lλ =
FTe (Lλ ).
λ∈Λ

λ∈Λ

λ∈Λ

λ∈Λ

The theorem below extends the Ie (−) operation defined on ideals in
[Kat08, Section 5] and in [BMS08, Definition 2.2] (where it is denoted
e
(−)[1/p ] ) to submodules of free R-modules.

Theorem 2.2. Let e ≥ 1. Given a submodule K ⊆ T α there exists a
e
minimal submodule L ⊆ T α for which K ⊆ L[p ] . We denote this minimal
submodule Ie (K).
Proof. Let L be the intersection of all submodules M ⊆ T α for which K ⊆
e
e
M [p ] . The intersection-flatness of T implies that K ⊆ L[p ] and clearly, L
is minimal with this property.

When F∗e T is T -free, this is a straightforward generalization of the calculation of Ie for ideals. To do so, fix a free basis B for F∗e T and note that
P
[pe ]
every element v ∈ T α can be expressed uniquely in the form v = b∈B ub b
where ub ∈ T α for all b ∈ B.
Proposition 2.3. Let e ≥ 1.
(a) For any submodules V1 , . . . , Vℓ ⊆ Rn , Ie (V1 + · · · + Vℓ ) = Ie (V1 ) +
· · · + Ie (Vℓ ).
(b) Let B be a free basis for F∗e T . Let v ∈ Rα and let
X [pe]
v=
ub b
b∈B

THE SUPPORT OF LOCAL COHOMOLOGY MODULES

5

be the unique expression for v where ub ∈ T α for all b ∈ B. Then
Ie (T v) is the submodule W of T α generated by {ub | b ∈ B}.
Proof. The proof of this proposition is a straightforward modification of the
proofs of propositions 5.2 and 5.6 in [Kat08] and Lemma 2.4 in [BMS08].
Clearly, Ie (V1 +· · ·+Vℓ ) ⊇ Ie (Vi ) for all 1 ≤ i ≤ ℓ, hence Ie (V1 +· · ·+Vℓ ) ⊇
Ie (V1 ) + · · · + Ie (Vℓ ). On the other hand
e

e

e

(Ie (V1 ) + · · · + Ie (Vℓ ))[p ] = Ie (V1 )[p ] + · · · + Ie (Vℓ )[p ] ⊇ V1 + · · · + Vℓ
and the minimality of Ie (V1 + · · · + Vℓ ) implies that Ie (V1 + · · · + Vℓ ) ⊆
Ie (V1 ) + · · · + Ie (Vℓ ) and (a) follows.
e
Clearly v ∈ W [p ] , and so Ie (T v) ⊆ W . On the other hand, let W be a
P
e
[pe ]
submodule of T α such that v ∈ W [p ] . Write v = si=1 ri wi for ri ∈ T and
P
pe
wi ∈ W for all 1 ≤ i ≤ s, and for each such i write ri = b∈B rbi
b where
rbi ∈ T for all b ∈ B. Now
!
s
X [pe]
X X
e
e
[p ]
p
b
ub b = v =
wi
rbi
b∈B

b∈B

i=1

[pe ]

and
are direct sums, we comparePcoefficients and obtain ub =
P sincee these
[pe ]
p
s
for all b ∈ B and so ub = ( si=1 rbi wi ) for all b ∈ B hence
i=1 rbi wi
ub ∈ W for all b ∈ B.

The behavior of the Ie operation under localization and completion will
be crucial for obtaining the results of this paper. To investigate this we need
the following generalization of [LS01, Lemma 6.6].
Lemma 2.4. Let T be a completion of T at a prime ideal P . Let α ≥ 0 and
e
e
let W be a submodule of T α . For all e ≥ 0, W [p ] ∩ T = (W ∩ T )[p ] .
Proof. If T is local with maximal ideal P , the result follows from a straightforward modification of the proof of [LS01, Lemma 6.6].
We now reduce the general case to the previous case which implies that
e
e
W [p ] ∩ TP = (W ∩ TP )[p ] . Intersecting with T now gives
e

e

e

e

W [p ] ∩ T = (W ∩ TP )[p ] ∩ T = (W ∩ TP ∩ T )[p ] = (W ∩ T )[p ] .

Lemma 2.5 (cf. [Mur13]). Let T be a localization of T or a completion at
a prime ideal.
For all e ≥ 1, and all submodules V ⊆ T α , Ie (V ⊗T T) exists and equals
Ie (V ) ⊗T T.
e

Proof. Let L ⊆ T α be a submodule, such that L[p ] ⊇ V ⊗T T. We clearly
e
e
have L[p ] ∩ T α = (L ∩ T α )[p ] when T is a localization of T and when T is
a completion of T this follows from the previous Lemma. We deduce that
(L ∩ T α ) ⊇ Ie (V ⊗T T ∩ T α ) and hence L ⊇ (L ∩ T α ) ⊗T T ⊇ Ie (V ⊗T T ∩
T α ) ⊗T T.

6

MORDECHAI KATZMAN AND WENLIANG ZHANG

But since Ie (V ⊗T T ∩ T α ) ⊗T T satisfies
e

e

(Ie (V ⊗T T∩α ) ⊗T T)[p ] = Ie (V ⊗T T∩T α )[p ] ⊗T T ⊇ (V ⊗T T∩T α )⊗T T ⊇ V ⊗T T
we deduce that Ie (V ⊗T T ∩ T α ) ⊗T T is the smallest submodule K ⊆ T α
e
for which K [p ] ⊇ V ⊗T T. We conclude that Ie (V ⊗T T) equals Ie (V ⊗T T ∩
T α ) ⊗T T.
We always have
Ie (V ⊗T T) = Ie (V ⊗T T ∩ T α ) ⊗T T ⊇ Ie (V ) ⊗T T.
On the other hand
e

e

(Ie (V ) ⊗T T)[p ] = Ie (V )[p ] ⊗T T ⊇ V ⊗T T
hence Ie (V ⊗T T) ⊆ Ie (V ) ⊗T T and thus Ie (V ⊗T T) = Ie (V ) ⊗T T.

3. Calculation of supports of FR -finite FR -modules
We begin by recalling the following result from [Lyu97, Proposition 2.3].
Remark 3.1 (Vanishing of γt ). Let M be an FR -finite FR -module with a
generating homomorphsim γ : M → FR (M ). Let γt denote the composition
M → FR (M ) → · · · → FRt (M ). We may assume that M has a presentaA

→ Rβ → M → 0 and write the generating homomorphism as
tion: Rα −
U
Coker(A) −
→ Coker(A[p] ) where U is a β × β matrix with entries in R. Then
t
γt is the composition of Coker(A) → · · · → Coker(A[p ] ). Note that M = 0
if and only if there is a t such that γt = 0. We have
γt = 0 ⇔ Im(U [p

t−1 ]

⇔ It (Im U [p

t

◦ · · · ◦ U [p] ◦ U ) ⊆ Im(A[p ] )

t−1 ]

◦ · · · ◦ U [p] ◦ U ) ⊆ Im A

⇔ I1 (U I1 (· · · I1 (U I1 (Im U ))) ⊆ Im A
⇔

I1 (U I1 (· · · I1 (U I1 (Im U ))) + Im A
=0
Im A

where we made repeated use of the facts that for any submodule M ⊆ Rβ
ℓ
we have Iℓ+1 (M ) = I1 (Iℓ (M )) and also Iℓ (U [p ] M ) = U Iℓ (M ).
Theorem 3.2.
(a) If
Ie (Im U [p

e−1 ]

e

◦ · · · ◦ U [p] ◦ U ) = Ie+1 (Im U [p ] ◦ · · · ◦ U [p] ◦ U )

then
(2)

Ie (Im U [p

e−1 ]

◦ · · · ◦ U [p] ◦ U ) = Ie+j (Im U [p

e+j−1 ]

for all j ≥ 0.
(b) There exists an integer e such that (2) holds.

◦ · · · ◦ U [p] ◦ U )

THE SUPPORT OF LOCAL COHOMOLOGY MODULES

7

e−1

Proof. Write Ve = Ie (Im U [p ] ◦ · · · ◦ U [p] ◦ U ). First we claim that if
Ve = Ve+1 then Ve = Ve+j for all j ≥ 0; we proceed by induction on j ≥ 0.
Using again the facts that for any submodule W ⊆ Rβ we have Iℓ+1 (M ) =
ℓ
I1 (Iℓ (M )) and that Iℓ (U [p ] M ) = U Iℓ (M ), we deduce that, if j ≥ 1, then
e+j−1
] ◦ · · · ◦ U [p] ◦ U )) = I (U V
Ve+j = I1 (Ie+(j−1) (Im U [p
1
e+j−1 ) and this, by
the induction hypothesis, equals I1 (U Ve ) = Ve+1 .
Next, we wish to show that for each prime ideal p there exists an integer
ep such that
(3)

Vep Rp = Vep +1 Rp .

and that for this ep , Vep +j Rp = Vep Rp for all j ≥ 0. After completing at p,
we have assume that our ring is a complete regular local ring, and we let E
cp . As this ring is complete
denote the injective hull of the residue field of R
and regular, there is a natural Frobenius map on E which we denote T ,
which can be extended to a Frobenius map on direct sums of E by letting
T act coordinate-wise; we denote these Frobenius maps also with T .
We now consider the Frobenius map Θ = U t T on E β ; in [KZ14, Lemma
e−1
3.6] it is shown that annE β Ie (Im U [p ] ◦ · · · ◦ U [p] ◦ U Rp )t ⊆ E β consists
of all elements killed by Θe . Now (cf. [HS77, Proposition 1.11] and [Lyu97,
e−1
Proposition 4.4]) show that there is an integer ep such that Iep (Im U [p ] ◦
ep
ep
· · · ◦ U [p] ◦ U Rp ) = Iep +1 (Im U [p ] ◦ · · · ◦ U [p] ◦ U Rp ) and Iep (Im U [p −1] ◦ · · · ◦
ep +j−1
U [p] ◦ U Rp ) = Iep +j (Im U [p  ] ◦ · · · ◦ U [p] ◦ U Rp ) forall j ≥ 0. Crucially,
e−1

Lemma 2.5 implies that Ie Im U [p ] ◦ · · · ◦ U [p] ◦ U Rp = Ve Rp for all
e ≥ ep and so (3) holds.
Consider the following subsets of Spec(R):
Pt = {p ∈ Spec(R) | Vt Rp = Vt+1 Rp } = Spec R \ Supp

Vt
.
Vt+1

These form an increasing sequence of open subsets of Spec R, and since for
each prime ideal p there is an integer tp such that
S

Vtp Rp = Vtp +1 Rp ,

we have t Pt = Spec(R). Now the quasicompactness of Spec R, guarantees
the existence of an integer e such that Pe = Spec(R); clearly that e satisfies
(2).

e−1

e

Corollary 3.3. If Ie (Im U [p ] ◦· · ·◦U [p] ◦U ) = Ie+1 (Im U [p ] ◦· · ·◦U [p] ◦U ),
then
 Im I (U [pe−1 ] ◦ · · · ◦ U [p] ◦ U ) + Im A 
e
= SuppR (M).
SuppR
Im A
4. Our algorithm and its complexity

Henceforth in this paper R will denote a polynomial ring over a field K
of prime characteristic p.

8

MORDECHAI KATZMAN AND WENLIANG ZHANG

Let M be an F -finite F -module with a generating morphism M → FR (M ).
Let A be an α×β matrix, which gives a presentation of M , i.e. Coker A ∼
= M.
U
[p]
Let U be a β × β matrix, for which the map Coker A −
→ Coker A is isomorphic to a generating morphism M → FR (M ). We compute the support
of M as follows.
Initialize L = Rβ .
Compute L′ = I1 (U L).
Check whether L′ ⊆ L. If this holds, let L = L′ and go to (2).
A
Compute a presentation for L+Im
Im A as the cokernel of a matrix W
with entries in R.
(5) Compute the ideal J of maximal minors of W .
(6) Output J: this is the defining ideal of the support of M.
(1)
(2)
(3)
(4)

In the rest of this section we discuss the complexity of this algorithm for
computing supports of F -finite F -modules.
We start with the following observation, relevant to the complexity of
steps (2) and (3) above.
Remark 4.1. Let δ be the largest degree of an entry in U and for any j ≥ 0 let
δj be the largest degree of a polynomial in a generator of Lj . The calculation
of I1 (−) as described in Proposition 2.3 implies that δj+1 ≤ (δj + δ)/p hence
δe ≤

1
1
δ
δ0
.
+ δ( + · · · + e ) ≤
pe
p
p
p−1

The iterative calculation I1 (U L) in step (2) involves a matrix multiplication and a I1 (−) operation, which, in view of Proposition 2.3, amounts to
collecting terms in polynomials. The complexity of this step depends on
(a) the size β of U , which is an input to the algorithm and does not
depend on p,
(b) the total number of terms occurring in each of the coordinates of a
set of generators of L.
In the worst case scenario, if the maximal degree of an entry in Up is Cp,
the total number of terms in (b) is bounded by




Cp + n − 1 β
= O pβ(n−1) .
n−1

In practice, the number of terms is much lower than this worst case.
Checking the inclusion in step (3) does not require computing Gröbner
bases: the discussion above shows that both L and L′ are generated by
vectors whose coordinates are polynomials whose degrees are bounded by
D := δ/(p − 1), i.e., L and L′ are given by generators in (R≤D )β where R≤D
denotes the K-vector space of polynomials of degrees at most D in R. Thus
L′ ⊆ L can be checked by checking whether each given generator of L′ is in
the sub-vector space of (R≤D )β spanned by the generators of L.

THE SUPPORT OF LOCAL COHOMOLOGY MODULES

9

Computing the presentation in step (4) of the algorithm involves computing the syzygies of the generators of L + Im A; this involves computing
a Gröbner bases for L + Im A.
In order to assess the practical advantage of our algorithm, we computed the support of 100 F -finite F -modules with randomly generated generating morphism C → FR1 (C) where C is a quotient of R2 and R =
Z/2Z[x1 , . . . , x5 ]. We denote t1 the time in seconds required by our algorithm to compute the support and t2 the time in seconds required to
compute a root using Grobner bases. The following is a plot of log t2 as a
function of log t1
log t2
+

+

6

+
+

+
+

4

+

++
+

+
++
+
+

+
+
+ +
+
+ ++ + +
+ + +
+
+
++
+
++
+ +++++
++
+
+
+
+
+
+++
++ +
++

+
+

++
++
+

+
+ +

+

2
+

+
++
+

−4

−2
++

+
+

+

+

++
+

+ ++
+
+

+
+

2

log t1
4

+

−2

+
++ +
+
+

+

+

+

+

−4

This suggests that for this characteristic and rank, t2 is approximately t21 .3
To further illustrate the effectiveness of our algorithm we compute the
following example.
Example 4.2. Consider three generic degree-2 polynomials in t: F1 (t) =
x0 + x1 t + x2 t2 , F2 (t) = y0 + y1 t + y2 t2 , F1 (t) = z0 + z1 t + z2 t2 . For any two
polynomials F (t), G(t) let Res(F, G) denote their Sylvester resultant, e. g. ,


x0 x1 x2 0
 0 x0 x1 x2 

Res(F1 , F2 ) = det 
 y0 y1 y2 0 
0 y0 y1 y2

Let I denote the ideal generated by Res(F1 , F2 ), Res(F1 , F3 ), Res(F2 , F3 ), Res(F1 +
F2 , F3 ) in the polynomial ring R over a field k whose variables are the x,y,
and zs above. In [Lyu95] it was asked whether H4I (R) = 0 and this was
settled in prime characteristic p > 2 (cf. [Kat97]) and in characteristic zero
3The Macaulay2 code used to produce this data and the data itself is available at [KZ].

10

MORDECHAI KATZMAN AND WENLIANG ZHANG

(cf. [Yan99, Theorem 3].) We used an implementation of our algorithm [KZ]
with Macaulay2 ([GS]) to settle the remaining case of characteristic 2: a
20-second run calculated the support of H4I (R) to be empty.
5. The complexity of computing the support of local
cohomology
We can use the fact that local cohomology modules HiJ (R) are F -finite
F -modules and apply the results in the previous sections to compute their
supports. To do so we would need to exhibit a generating morphism for
these; a standard choice of generating morphism is given by the R-linear
map Extj (R/J, R) → Extj (R/J [p] , R) induced by the natural surjection
R/J [p] → R/J (cf. [Lyu97, Proposition 1.11]). The calculation of this induced map would normally involve finding a free resolution F• for R/J,
extending the quotient map R/J [p] → R/J to a map of free resulutions
FR1 (F• ) → F• , applying Hom(−, R) to both resolutions, and finally computing the induced map of cohomologies. Each of these steps involves computing
multiple Gröbner bases, and, therefore, potentially unfeasible even in simple
instances. For example, even when J is generated by one f ∈ R, the map of
f p−1

Ext1 (R/f R, R) → Ext1 (R/f p R, R) above is isomorphic to R/f −−−→ R/f p
and even when f is, say, a random polynomial of degree 5 in 5 variables,
expanding f p−1 is not feasible beyond the first few primes (e.g., the M2
server Habanero crashes after p = 13.) Clearly, if the calculation of a generating morphism is infeasible, the algorithm of the preceding section cannot
be used to calculate the support of the given local cohomology module.
However, in those cases when a generating morphism can actually be
computed, the method above is guaranteed to produce generating morphisms which yield “bounded” generating morphism working over increasing prime characteristics p, as we describe below. Let S = Z[x1 , . . . , xn ] and
let J ⊆ S be an ideal. For any prime p let Jp denote the image of J in
Rp = Z/pZ[x1 , . . . , xn ]. An interesting and natural question arising in this
context is the description of the properties of the local cohomology module
HjJp (Rp ) as p ranges over all primes and we now turn our attention to these.
For different choices of prime p the matrices A and U above will be different and this could result in different values of δe which are unbounded
as p ranges over all primes. We now show that this is not the case. Let
Up denote the square matrix that induces the map Extj (Rp /Jp , Rp ) →
[p]
Extj (Rp /Jp , Rp ) and δp to denote the maximal degree of entries in Up .
We also denote L0,p = Rp and Li+1,p = I1 (Up Li,p ) and use δe,p to denote the
largest degree of a polynomial in a generator of Le,p .
A

A

A

1
2
s
Theorem 5.1. Let 0 → S bs −→
−→
S → S/J → 0 be a
→
S b1 −
· · · −−
free resolution of S/J. Let ∆ denote the maximal degree of any entry in

THE SUPPORT OF LOCAL COHOMOLOGY MODULES

11

A1 , . . . , As . Let p be a prime integer which is also a regular element on S/J.
Then δe,p ≤ 2j∆ for all integers e ≥ 1.
Proof. Since p is a regular element on S/J, tensoring the free resolution of
S/J with R/pR produces a free resolution of Rp /Jp . Hence the maximal
degree of entries in the maps of this free resolution of Rp /Jp is at most ∆.
b
b
Let θj denote the map Rpj → Rpj in the following commutative diagram
[p]
induced by Rp /Jp → Rp /Jp
0

/ Rbs
p
O

As,p

/ ···

/ Rbs
p

Aj,p

/ ···

/ Rb1
p

A1,p

[p]
As,p

/ ···

/ Rpbj

/ Rp
O

/ Rp /Jp
O

/0

/ Rp /Jp[p]

/0

=

θj

θs

0

/ Rpbj
O

[p]
Aj,p

/ ···

/ Rb1
p

[p]
A1,p

/ Rp

An easy induction on j shows that the maximal degree of entries in θj is at
[p]
most jp∆. The map Up : Extj (Rp /Jp , Rp ) → Extj (Rp /Jp , Rp ) is induced
by the transpose of θj and hence the maximal degree in an entry of Up is
also bounded by jp∆ (cf. Remark 4.1.) Now
jp∆
≤ 2j∆,
δe,p ≤
p−1
for all e ≥ 1.

Corollary 5.2. There is an integer N , independent of p, such that δe,p ≤ N
for all e and all p.
In particular, there is an integer N ′ , independent of p, such that min{e |
Le,p = Le+1,p } ≤ N ′ for all p, i.e. for each prime integer p, the number of
steps required to compute the stable value Le,p is bounded by N ′ .
Proof. The second statement follows immediately from the first since, once
the degree is bounded, the number of steps will be bounded by the number
of monomials with the bonded degree.
To prove the first statement, it suffices to note there are only finitely
many associated prime ideals of S/J in S and hence p is a regular element
on S/J for almost all p.

6. Iterated local cohomology modules
Let f1 , . . . , fm be a L
sequence of elements in R and let N be an R-module.
We will write Ki :=
1≤j1 <···<ji ≤m Nj1 ···ji to denote the i-th term of the
Koszul (co)complex K• (M ; f ) (where each Nj1 ···ji = N ), and we will use
H i (N ; f ) to denote the i-th Koszul (co)homology.
Proposition 6.1. Let M be an FR -finite FR -module with a generating hoϕ
momorphism M −
→ FR (M ) and let I = (f1 , . . . , fm ) be an ideal of R. Then
HIi (M) admits a generating homomorphism
H i (M ; f ) → FR (H i (M ; f )).

12

MORDECHAI KATZMAN AND WENLIANG ZHANG

Proof. Consider the following commutative diagram:
..
.O

..
.O
FR (φ1 )

FR (φi )

/ FR (K 1 )
/ ···
O FR (δ1 )

0

φ1

/ K1

0

..
.O
FR (φm )

/ ···
/ FR (K i )
O FR (δi )

/ FR (K m ) = FR (M )
O

φi
δ1

/ ···

/ Ki

/0

φm
δi

/ Km = M

/ ···

/

where the bottom row is the Koszul (co)complex of M on f and
φi :

M

1≤j1 <···<ji ≤m

⊕1≤j

p−1
<···<j ≤m ϕ◦(fj1 ···fji )

i
Mj1 ···ji −−−−1−−−−−
−−−−−−−−−−−→ FR (

M

Mj1 ···ji ).

1≤j1 <···<ji ≤m

It follows from [Lyu97, 1.10(c)] that the φi are generating morphisms of
Mfj1 ···fji . Therefore taking direct limit of each row of the diagram produces
the Čech complex Č(M; f ). Since taking direct limits preserves exactness,
lim(H i (M ; f ) → FR (H i (M ; f )) → · · · ) = HIi (M). Our conclusion follows.
−→

Combining what we have so far in this section, we now have an algorithm
to compute the support of HIi11 · · · HIiss (R). For example, the case s = 2 relevant to the calculation of Lyubeznik numbers in handled as follows. Start
with a generating morphism Exti2 (R/I2 , R) → FR (Exti2 (R/I2 , R)). Using
Proposition 6.1, we know that the Koszul cohomology H i1 (Exti2 (R/I2 , R); f )
(with I1 = (f )) is a generating homomorphism of HIi11 HIi22 (R). We may then
apply Corollary 3.3 to compute the support of HIi11 HIi22 (R).
7. The support of local cohomology of hypersurfaces
Throughout this section R denotes a regular ring of prime characteristic
p, I ⊆ R an ideal, and g ∈ R some fixed element.
Following [Lyu97, §2] we write


3φ
2φ
FR
FR
φ
2
i
i
i
1
i
HI (R) = lim ExtR (R/I, R) −
→ FR ExtR (R/I, R) −−→ FR ExtR (R/I, R) −−→ . . .
→

φ

where FRe (−) denotes the eth Frobenius functor, and φ : ExtiR (R/I, R) −
→
i
i
1
[p]
∼
FR ExtR (R/I, R) = ExtR (R/I , R) is the R-linear map induced by the
A

i
surjection R/I [p] → R/I. For all i ≥ 0 we fix a presentation Rαi −→
Rβi where Ai is a βi × αi matrix with entries in R. We can now find a

φ

βi × βi matrix Ui with entries in R for which the map φ : ExtiR (R/I, R) −
→
FR1 ExtiR (R/I, R) is isomorphic to the map Ui : Coker Ai → FR1 (Coker Ai ) =
[p]
Coker Ai given by multiplication by Ui .

THE SUPPORT OF LOCAL COHOMOLOGY MODULES

13

Theorem 7.1. For any i ≥ 0 consider the map g : HiI (RP ) → HiI (RP ) given
by multiplication by g. Let Ii denote the set of primes P ⊂ R for which the
map g is not injective and let Si denote the set of primes P ⊂ R for which
the map g is not surjective. For ℓ, e, j ≥ 0 write
[pe+j−1 ]

(ℓ)

Vej = Uℓ

[pe+j−2 ]

Uℓ

[pe ]

· · · Uℓ

.

Then
(i)

(a)

Ii

is closed and equal to Supp

(ker V0η :Rβ g)
(i)

ker V0η

(b) Si is closed and equal to Supp S

j≥0



for some η > 0,

Rβ
(i)

j

gRβ + Im A[p ] :Rβ V0j

(c) the support of HiI (R/gR) is closed and equal to Ii ∪ Si .

 , and

Proof. Fix some i ≥ 0 and write β, A and U for βi , Ai and Ui . The map
g : HiI (R) → HiI (R) can be described as a map of direct limit systems
(4)

Coker A

U

/ Coker A[p]U

g



Coker A

[p]

e−1 ]

U [p /

/ ...

Coker A[p

g
U



[pe ]
eU
]

/ ... .

g
U [p]

/ Coker A[p]

e−1 ]
U [p



/ Coker A[peU]

/ ...

[pe ]

/ ...

(i)

For any e, j ≥ 0 abbreviate Vej = Vej , and note that it is the matrix
e

e+j

corresponding to the composition map Coker A[p ] → Coker A[p ] in the
direct limits in (4). Any element in HiI (RP ) can be represented by an element
[pe ]
a ∈ Coker AP for some e ≥ 0, and this element represents the zero element
[pe+j ]

if and only if there exists a j ≥ e for which Vej a ∈ Im AP
only if
a ∈ (Im A[p

e+j ]

, i.e., if and

:Rβ Vej )P .
j

Consider the kernels Kj of the maps V0j : Coker A → Coker A[p ] ; these form
an ascending chain of submodules of Coker A and hence stabilize for all j
e
e+j
beyond some η ≥ 0. Note that the map Vej : Coker A[p ] → Coker A[p ]
is obtained by applying the exact functor FRe (−) to the map V0j , hence the
kernels of the maps Vej also stabilize for j ≥ η.
To prove (a) we now note that an element in HiI (RP ) represented a ∈
[pe ]
Coker AP is multiplied by g to zero
 a ∈ (ker Veη :Rβ g)P
 if and only if
(ker Veη :Rβ g)
= 0, i.e., if g is not
and so g is injective if and only iff
ker Veη
P

a zero divisor on Rβ / ker Veη P . But Rβ / ker Veη = FRe (Rβ / ker V0η ) and,
since R is regular, FRe (Rβ / ker V0η ) and Rβ / ker V0η have the same associated
primes, so we deduce that multiplication
by g is injective if and only if g is

not a zero divisor on Rβ / ker V0η P . We deduce that for a prime P ⊂ R,

14

MORDECHAI KATZMAN AND WENLIANG ZHANG

multiplication by g on HiI (RP ) is injective if and only if


(ker V0η :Rβ g)
=0
ker V0η
P
so Ii = Supp

(ker V0η :Rβ g)
.
ker V0η

To prove (b) we now note that an element in HiI (RP ) represented a ∈
[pe ]
Coker AP is in the image
 of g if and only if there exists a j ≥ 0 such that
[pe+j ]
β
hence g is surjective if for all e ≥ 0,
Vej a ∈ gR + Im AP

[
e+j
gRβ + Im A[p ] :Rβ Vej
= Rβ P .
P

j≥0

Furthermore,

[p]


e+1+j ]
e+j
=
gp Rβ + Im A[p
:Rβ Ve+1,j
gRβ + Im A[p ] :Rβ Vej


e+1+j ]
⊆
gRβ + Im A[p
:Rβ Ve+1,j
so for for all e ≥ 0,
[

gRβ + Im A[p

e+j ]

:Rβ Vej

j≥0

if and only if

[

j≥0

j

gRβ + Im A[p ] :Rβ V0j





P

P

= Rβ P

= Rβ P .

We conclude that g is not surjective if and only if P ∈ Supp Rβ /
To prove (c) consider the long exact sequence

S

j≥0




j
gRβ + Im A[p ] :Rβ V0j .

g

g

→ Hi+1
· · · → HiI (R) −
→ HiI (R) → HiI (R/gR) → Hi+1
I (R) → . . .
I (R) −
g

induced by the short exact sequence 0 → R −
→ R → R/gR
 → 0. Note that
g
i
i
i
HI (R/gR)P = 0 if and only if both HI (R) −
→ HI (R)
is surjective and
P


g
→ Hi+1
Hi+1
and the result follows.
I (R)
I (R) −
P

Question 7.2. Theorem 3.2 gives us an effective method for the calculation
of Ii . However, we do not know how to compute Si , hence we ask the
following: is there an effective method to bound the value of e for which
 

[
e
j
(i)
(i)
gRβ + Im A[p ] :Rβ V0j = gRβ + Im A[p ] :Rβ V0e ?
j≥0

It turns out that part of our Theorem 7.1 can be extended to the case of
isolated singular points.

THE SUPPORT OF LOCAL COHOMOLOGY MODULES

15

Corollary 7.3. Let R be a noetherian commutative ring of prime characteristic that has finitely many isolated singular points. Let g ∈ R be a
nonzerodivisor. Then Supp(HIj (R/gR)) is Zariski-closed for each integer j
and ideal I of R.
Proof. Let
T {m1 , . . . , mt } denotes the set of isolated singular points of R.
Set a = ti=1 mi . Let {f1 , . . . , fs } be a set of generators of a. It follows
from Theorem 7.1 that SuppRf (HjI (Rfk /gRfk )) is closed, i.e. it has finitely
k
many minimal associated primes. By the bijection between the set of associated primes of HjI (R/gR) that do not contain fk and Ass(HjI (Rfk /gRfk )),
it follows that the minimal associated primes of HjI (R/gR) are contained
in the union of {m1 , . . . , mt } and the set of minimal associated primes of

HjI (Rfk /gRfk ) which is a finite set.
The proof of Corollary 7.3 can also be used to prove the following result
which is of independent interest.
Proposition 7.4. Let R be either
(1) a noetherian commutative ring of prime characteristic, or
(2) of finite type over a field of characteristic 0.
Suppose that R has finitely many isolated singular points. Then HjI (R) has
only finitely many associated primes for each integer j and each ideal I of
R.
Proof. Let {m1 , . . . , mt } denotes the set of isolated singular points of R. Set
T
a = ti=1 mi . Let {f1 , . . . , fs } be a set of generators of a. It follows from
our assumptions on R that Rfk is either a noetherian regular ring of prime
characteristic or a regular ring of finite type over a field of characteristic 0
(cf. [Lyu93, Corollary 3.6]). Consequently, Ass(HjI (Rfk )) is finite for each
generator fk . Since there is a bijection between the set of associated primes
of HjI (R) that do not contain fk and Ass(HjI (Rfk )), it follows that
Ass(HjI (R))

⊆

s
[

Ass(HjI (Rfk ))

k=1

[

{m1 , . . . , mt }.


Acknowledgement
We would like to thank the anonymous referee whose suggestions and
criticisms improved this paper substantially.
References
[BMS08] M. Blickle, M. Mustaţă, and K. Smith: Discreteness and rationality of
F-thresholds, Michigan Math. J. 57 (2008), 43–61.
[GS]
D. R. Grayson and M. E. Stillman: Macaulay2, a software system for research in algebraic geometry.

16

[Har67]

MORDECHAI KATZMAN AND WENLIANG ZHANG

R. Hartshorne: Local cohomology, A seminar given by A. Grothendieck, Harvard University, Fall, vol. 1961, Springer-Verlag, Berlin, 1967. MR0224620 (37
#219)

[Har68]

R. Hartshorne: Cohomological dimension of algebraic varieties, Ann. of Math.
(2) 88 (1968), 403–450. 0232780 (38 #1103)
[HS77] R. Hartshorne and R. Speiser: Local cohomological dimension in characteristic p, Ann. of Math. (2) 105 (1977), no. 1, 45–79. MR0441962 (56 #353)
[HNB]
M. Hochster and L. Núñez-Betancourt: Support of local cohomology modules over hypersurfaces and rings with FFRT, preliminary preprint.
[Kat97] M. Katzman: The cohomological dimension of a resultant variety in prime
characteristic p, J. Algebra 191 (1997), no. 2, 510–517. 1448806 (98a:13023)
[Kat08] M. Katzman: Parameter-test-ideals of Cohen-Macaulay rings, Compos. Math.
144 (2008), no. 4, 933–948. MR2441251 (2009d:13030)
[KZ]
M. Katzman and W. Zhang:
A macaulay2 library for computing
supports
of
F-finite
F-modules,
Freely
available
from
http://www.katzman.staff.shef.ac.uk/Fsupport/.
[KZ14] M. Katzman and W. Zhang: Annihilators of Artinian modules compatible
with a Frobenius map, J. Symbolic Comput. 60 (2014), 29–46. 3131377
[Lyu93] G. Lyubeznik: Finiteness properties of local cohomology modules (an application of D-modules to commutative algebra), Invent. Math. 113 (1993), no. 1,
41–55. 1223223 (94e:13032)
[Lyu95] G. Lyubeznik: Minimal resultant systems, J. Algebra 177 (1995), no. 2, 612–
616. 1355218 (96i:13013)
[Lyu97] G. Lyubeznik: F -modules: applications to local cohomology and D-modules in
characteristic p > 0, J. Reine Angew. Math. 491 (1997), 65–130. MR1476089
(99c:13005)

[LS01]

G. Lyubeznik and K. E. Smith: On the commutation of the test ideal with
localization and completion, Trans. Amer. Math. Soc. 353 (2001), no. 8, 3149–
3180 (electronic). MR1828602 (2002f:13010)
[Mur13] S. Murru: On the upper semi-continuity of HSL numbers, arXiv:1302.1124.
[Ogu73] A. Ogus: Local cohomological dimension of algebraic varieties, Ann. of Math.
(2) 98 (1973), 327–365. 0506248 (58 #22059)
[PS73] C. Peskine and L. Szpiro: Dimension projective finie et cohomologie locale.
Applications à la démonstration de conjectures de M. Auslander, H. Bass et
A. Grothendieck, Inst. Hautes Études Sci. Publ. Math. (1973), no. 42, 47–119.
MR0374130 (51 #10330)

[Yan99] Z. Yan: Minimal resultant systems, J. Algebra 216 (1999), no. 1, 105–123.
1694582 (2000d:14022)

Department of Pure Mathematics, University of Sheffield, Hicks Building,
Sheffield S3 7RH, United Kingdom
E-mail address: M.Katzman@sheffield.ac.uk
Department of Mathematics, Statistics, and Computer Science, University
of Illinois at Chicago, 851 S. Morgan Street, Chicago, IL 60607-7045
E-mail address: wlzhang@uic.edu

