Newer
Older
and so we can write:
ts(beta=beta).rank(u)
var("q", domain="real")
c_in_terms_of_q = c_lower_bound + q
c=\chern_1(u) = \sage{c_in_terms_of_q}
\qquad 0 \leq q \coloneqq \chern_1^{\beta}(u) \leq \chern_1^{\beta}(v)
\end{equation}
Furthermore, $\chern_1 \in \ZZ$ so we only need to consider
$q \in \frac{1}{n} \ZZ \cap [0, \chern_1^{\beta}(F)]$,
where $n$ is the denominator of $\beta$.
For the next subsections, we consider $q$ to be fixed with one of these values,
and we shall be varying $\chern_0(E) = r$ to see when certain inequalities fail.
This section studies the numerical conditions that $u$ must satisfy as per
lemma \ref{lem:num_test_prob1}
(or corrolary \ref{cor:num_test_prob2})
\subsubsection{Size of pseudo-wall: $\chern_2^P(u)>0$ }
\label{subsect-d-bound-radiuscond}
This condition refers to condition
\ref{item:radiuscond:lem:num_test_prob1}
from lemma \ref{lem:num_test_prob1}
(or corrolary \ref{cor:num_test_prob2}).
In the case where we are tackling problem \ref{problem:problem-statement-2}
(with $\beta = \beta_{-}$), this condition, when expressed as a bound on $d$,
amounts to:
\begin{align}
\label{eqn:radius-cond-betamin}
\chern_2^{\beta_{-}}(u) &\geq 0 \\
d &\geq \beta_{-}q + \frac{1}{2} \beta_{-}^2r
\end{align}
Semistability of the Semistabilizer:
\texorpdfstring{
$\Delta(u) \geq 0$
This condition refers to condition
\ref{item:bgmlvu:lem:num_test_prob1}
from lemma \ref{lem:num_test_prob1}
(or corrolary \ref{cor:num_test_prob2}).
# First Bogomolov-Gieseker form expression that must be non-negative:
bgmlv2 = Δ(u)
Expressing $\Delta(u)\geq 0$ in terms of $q$ as defined in eqn \ref{eqn-cintermsofm}
we get the following:
\begin{sagesilent}
bgmlv2_with_q = (
bgmlv2
.expand()
.subs(c == c_in_terms_of_q)
)
\end{sagesilent}
\begin{equation}
\sage{0 <= bgmlv2_with_q}
\end{equation}
\noindent
This can be rearranged to express a bound on $d$ as follows
(recall from condition \ref{item:rankpos:lem:num_test_prob1}
in lemma \ref{lem:num_test_prob1} or corrolary
\ref{cor:num_test_prob2} that $r>0$):
bgmlv2_d_ineq = (
(0 <= bgmlv2_with_q)/2/r # rescale assuming r > 0
+ d # Rearrange for d
).expand()
bgmlv2_d_upperbound = bgmlv2_d_ineq.rhs()
\label{eqn-bgmlv2_d_upperbound}
\sage{bgmlv2_d_ineq}
\end{equation}
\begin{sagesilent}
# Seperate out the terms of the lower bound for d
bgmlv2_d_upperbound_without_hyp = (
bgmlv2_d_upperbound
bgmlv2_d_upperbound_const_term = (
bgmlv2_d_upperbound_without_hyp
bgmlv2_d_upperbound_linear_term = (
bgmlv2_d_upperbound_without_hyp
- bgmlv2_d_upperbound_const_term
bgmlv2_d_upperbound_exp_term = (
bgmlv2_d_upperbound
- bgmlv2_d_upperbound_without_hyp
\end{sagesilent}
Viewing equation \ref{eqn-bgmlv2_d_upperbound} as a lower bound for $d$ in term
$\sage{bgmlv2_d_upperbound_const_term}$,
a linear term
$\sage{bgmlv2_d_upperbound_linear_term}$,
and a hyperbolic term
$\sage{bgmlv2_d_upperbound_exp_term}$.
Notice that in the context of problem \ref{problem:problem-statement-2}
($\beta = \beta_{-}$),
the constant and linear terms match up with the ones
for the bound found for $d$ in subsubsection \ref{subsect-d-bound-radiuscond}.
Semistability of the Quotient:
\texorpdfstring{
$\Delta(v-u) \geq 0$
Δ(v-u) ≥ 0
\label{subsect-d-bound-bgmlv3}
This condition refers to condition
\ref{item:bgmlvv-u:lem:num_test_prob1}
from lemma \ref{lem:num_test_prob1}
(or corrolary \ref{cor:num_test_prob2}).
Expressing $\Delta(v-u)\geq 0$ in term of $q$ and rearranging as a bound on
$d$ yields:
1165
1166
1167
1168
1169
1170
1171
1172
1173
1174
1175
1176
1177
1178
1179
1180
1181
1182
1183
1184
1185
1186
# Third Bogomolov-Gieseker form expression that must be non-negative:
bgmlv3 = Δ(v-u)
bgmlv3_with_q = (
bgmlv3
.expand()
.subs(c == c_in_terms_of_q)
)
var("r_alt",domain="real") # r_alt = r - R temporary substitution
bgmlv3_with_q_reparam = (
bgmlv3_with_q
.subs(r == r_alt + R)
/r_alt # This operation assumes r_alt > 0
).expand()
bgmlv3_d_ineq = (
((0 <= bgmlv3_with_q_reparam)/2 + d) # Rearrange for d
.subs(r_alt == r - R) # Resubstitute r back in
.expand()
)
# Check that this equation represents a bound for d
bgmlv3_d_upperbound = bgmlv3_d_ineq.rhs() # Keep hold of lower bound for d
# Seperate out the terms of the lower bound for d
bgmlv3_d_upperbound_without_hyp = (
bgmlv3_d_upperbound
.subs(1/(R-r) == 0)
)
bgmlv3_d_upperbound_const_term = (
bgmlv3_d_upperbound_without_hyp
.subs(r==0)
)
bgmlv3_d_upperbound_linear_term = (
bgmlv3_d_upperbound_without_hyp
- bgmlv3_d_upperbound_const_term
).expand()
bgmlv3_d_upperbound_exp_term = (
bgmlv3_d_upperbound
- bgmlv3_d_upperbound_without_hyp
).expand()
# Verify the simplified forms of the terms that will be mentioned in text
var("chb1v chb2v",domain="real") # symbol to represent ch_1^\beta(v)
var("psi phi", domain="real") # symbol to represent ch_1^\beta(v) and
# ch_2^\beta(v)
(
# keep hold of this alternative expression:
bgmlv3_d_upperbound_const_term_alt := (
)
)
.subs(phi == v.twist(beta).ch[2]) # subs real val of ch_1^\beta(v)
.expand()
)
# Keep hold of this alternative expression:
bgmlv3_d_upperbound_exp_term_alt :=
(
+ (C - q)^2/2
- D*R
)/(r-R)
)
.subs(phi == v.twist(beta).ch[2]) # subs real val of ch_1^\beta(v)
.expand()
)
assert bgmlv3_d_upperbound_exp_term == (
(
# Keep hold of this alternative expression:
bgmlv3_d_upperbound_exp_term_alt2 :=
(
(psi - q)^2/2/(r-R)
)
)
.subs(psi == v.twist(beta).ch[1]) # subs real val of ch_1^\beta(v)
.expand()
\def\psi{\chern_1^{\beta}(v)}
\def\phi{\chern_2^{\beta}(v)}
\begin{equation*}
\label{eqn-bgmlv3_d_upperbound}
d \leq
\sage{bgmlv3_d_upperbound_linear_term}
+ \sage{bgmlv3_d_upperbound_const_term_alt}
+ \sage{bgmlv3_d_upperbound_exp_term_alt2}
\qquad
\text{where }r>R
\end{equation*}
For $r=R$, $\Delta(v-u)\geq 0$ is always true, and for $r<R$ it gives a lower
bound on $d$, but it is weaker than the one given by the lower bound in
subsubsection \ref{subsect-d-bound-radiuscond}.
Viewing the right hand side of equation \ref{eqn-bgmlv3_d_upperbound}
as a function of $r$, the linear and constant terms almost match up with the
ones in the previous section, up to the
$\chern_2^{\beta}(v)$ term.
However, when specializing to problem \ref{problem:problem-statement-2} again
(with $\beta = \beta_{-}$), then we have $\chern^{\beta}_2(v) = 0$.
And so in this context, the linear and constant terms do match up with the
previous subsubsections.
\subsubsection{All Bounds on $d$ Together for Problem
\ref{problem:problem-statement-2}}
As mentioned in passing, when specializing to solutions $u$ of problem
\ref{problem:problem-statement-2}, the bounds on
$d=\chern^{\beta_{-}}_2(u)$ induced by conditions
\ref{item:bgmlvu:lem:num_test_prob2},
\ref{item:bgmlvv-u:lem:num_test_prob2}, and
\ref{item:radiuscond:lem:num_test_prob1}
from corrolary \ref{cor:num_test_prob2} have the same constant and linear
terms in $r$, but different hyperbolic terms.
These give bounds with the same assymptotes when we take $r\to\infty$
(for any fixed $q=\chern_1^{\beta_{-}}(u)$).
% redefine \beta (especially coming from rendered SageMath expressions)
% to be \beta_{-} for the rest of this subsubsection
\bgroup
\let\originalbeta\beta
\renewcommand\beta{{\originalbeta_{-}}}
\bgroup
% redefine \psi in sage expressions (placeholder for ch_1^\beta(F)
\def\psi{\chern_1^{\beta}(F)}
d &>&
\frac{1}{2}\beta^2 r
&+ \beta q,
\phantom{+}& % to keep terms aligned
&\qquad\text{when\:} r > 0
\label{eqn:radiuscond_d_bound_betamin}
\sage{bgmlv2_d_upperbound_linear_term}
&+ \sage{bgmlv2_d_upperbound_const_term}
+& \sage{bgmlv2_d_upperbound_exp_term},
&\qquad\text{when\:} r > 0
\label{eqn:bgmlv2_d_bound_betamin}
\sage{bgmlv3_d_upperbound_linear_term}
&+ \sage{bgmlv3_d_upperbound_const_term_alt.subs(phi == 0)}
% ^ ch_2^\beta(F)=0 for beta_{-}
+& \sage{bgmlv3_d_upperbound_exp_term_alt2},
\label{eqn:bgmlv3_d_bound_betamin}
positive_radius_condition = (
(
(0 > - u.twist(beta).ch[2])
+ d # rearrange for d
.subs(solve(q == u.twist(beta).ch[1], c)[0]) # express c in term of q
.expand()
)
def beta_min(chern):
ts = stability.Tilt()
return min(
map(
lambda soln: soln.rhs(),
solve(
(ts.degree(chern))
.expand()
.subs(ts.alpha == 0),
beta
)
)
)
v_example = Chern_Char(3,2,-2)
q_example = 7/3
def plot_d_bound(
v_example,
q_example,
ymax=5,
ymin=-2,
xmax=20,
aspect_ratio=None
):
1375
1376
1377
1378
1379
1380
1381
1382
1383
1384
1385
1386
1387
1388
1389
1390
1391
1392
1393
1394
1395
1396
1397
1398
1399
1400
1401
1402
1403
1404
1405
1406
# Equations to plot imminently representing the bounds on d:
eq2 = (
bgmlv2_d_upperbound
.subs(R == v_example.ch[0])
.subs(C == v_example.ch[1])
.subs(D == v_example.ch[2])
.subs(beta = beta_min(v_example))
.subs(q == q_example)
)
eq3 = (
bgmlv3_d_upperbound
.subs(R == v_example.ch[0])
.subs(C == v_example.ch[1])
.subs(D == v_example.ch[2])
.subs(beta = beta_min(v_example))
.subs(q == q_example)
)
eq4 = (
positive_radius_condition.rhs()
.subs(q == q_example)
.subs(beta = beta_min(v_example))
)
example_bounds_on_d_plot = (
plot(
eq3,
(r,v_example.ch[0],xmax),
color='green',
linestyle = "dashed",
legend_label=r"upper bound: $\Delta(v-u) \geq 0$",
)
+ plot(
eq2,
(r,0,xmax),
color='blue',
linestyle = "dashed",
legend_label=r"upper bound: $\Delta(u) \geq 0$"
)
+ plot(
eq4,
(r,0,xmax),
color='orange',
linestyle = "dotted",
legend_label=r"lower bound: $\mathrm{ch}_2^{\beta_{-}}(u)>0$"
)
)
example_bounds_on_d_plot.ymin(ymin)
example_bounds_on_d_plot.ymax(ymax)
example_bounds_on_d_plot.axes_labels(['$r$', '$d$'])
if aspect_ratio:
example_bounds_on_d_plot.set_aspect_ratio(aspect_ratio)
return example_bounds_on_d_plot
\end{sagesilent}
\begin{figure}
\centering
\centering
\sageplot[width=\linewidth]{plot_d_bound(v_example, 0, ymin=-0.5)}
\caption{$q = 0$ (all bounds other than green coincide on line)}
\label{fig:d_bounds_xmpl_min_q}
\end{subfigure}%
\centering
\sageplot[width=\linewidth]{plot_d_bound(v_example, 4, ymin=-3, ymax=3)}
\caption{$q = \chern^{\beta}(F)$ (all bounds other than blue coincide on line)}
\label{fig:d_bounds_xmpl_max_q}
\end{subfigure}
\caption{
Bounds on $d\coloneqq\chern_2(E)$ in terms of $r\coloneqq\chern_0(E)$ for fixed, extreme,
values of $q\coloneqq\chern_1^{\beta}(E)$.
Where $\chern(F) = (3,2,-2)$.
}
\label{fig:d_bounds_xmpl_extrm_q}
\end{figure}
Recalling that $q \coloneqq \chern^{\beta}_1(E) \in [0, \chern^{\beta}_1(F)]$,
it is worth noting that the extreme values of $q$ in this range lead to the
tightest bounds on $d$, as illustrated in figure
(\ref{fig:d_bounds_xmpl_extrm_q}).
In fact, in each case, one of the weak upper bounds coincides with one of the
weak lower bounds, (implying no possible destabilizers $E$ with
$\chern_0(E)=\vcentcolon r>R\coloneqq\chern_0(F)$ for these $q$-values).
This indeed happens in general since the right hand sides of
(eqn \ref{eqn:bgmlv2_d_bound_betamin}) and
(eqn \ref{eqn:radiuscond_d_bound_betamin}) match when $q=0$.
In the other case, $q=\chern^{\beta}_1(F)$, it is the right hand sides of
(eqn \ref{eqn:bgmlv3_d_bound_betamin}) and
(eqn \ref{eqn:radiuscond_d_bound_betamin}) which match.
The more generic case, when $0 < q\coloneqq\chern_1^{\beta}(E) < \chern_1^{\beta}(F)$
for the bounds on $d$ in terms of $r$ is illustrated in figure
(\ref{fig:d_bounds_xmpl_gnrc_q}).
The question of whether there are pseudo-destabilizers of arbitrarily large
rank, in the context of the graph, comes down to whether there are points
$(r,d) \in \ZZ \oplus \frac{1}{2} \ZZ$ (with large $r$)
% TODO have a proper definition for pseudo-destabilizers/walls
that fit above the yellow line (ensuring positive radius of wall) but below the
blue and green (ensuring $\Delta(E), \Delta(G) > 0$).
These lines have the same assymptote at $r \to \infty$
(eqns \ref{eqn:bgmlv2_d_bound_betamin},
\ref{eqn:bgmlv3_d_bound_betamin},
As mentioned in the introduction (sec \ref{sec:intro}), the finiteness of these
solutions is entirely determined by whether $\beta$ is rational or irrational.
Some of the details around the associated numerics are explored next.
\begin{figure}
\centering
\sageplot[
width=\linewidth
]{plot_d_bound(v_example, 2, ymax=4, ymin=-2, aspect_ratio=1)}
\caption{
Bounds on $d\coloneqq\chern_2(E)$ in terms of $r\coloneqq \chern_0(E)$ for a fixed
value $\chern_1^{\beta}(F)/2$ of $q\coloneqq\chern_1^{\beta}(E)$.
Where $\chern(F) = (3,2,-2)$.
}
\label{fig:d_bounds_xmpl_gnrc_q}
\end{figure}
\subsection{Bounds on Semistabilizer Rank \texorpdfstring{$r$}{r}}
Now, the inequalities from the above (TODO REF) will be used to find, for
each given $q=\chern^{\beta}_1(E)$, how large $r$ needs to be in order to leave
no possible solutions for $d$. At that point, there are no Chern characters
$(r,c,d)$ that satisfy all inequalities to give a pseudowall.
\subsubsection{All Semistabilizers Left of $V_v$ for Rational beta}
The strategy here is similar to what was shown in (sect
\ref{sec:twisted-chern}).
One specialization here is to use that $\ell\coloneqq c_1(H)$ generates $NS(X)$, so that
in fact, any Chern character can be written as
$\left(r,c\ell,\frac{e}{2}\ell^2\right)$ for $r,c,e\in\ZZ$.
var("a_v b_q n") # Define symbols introduce for values of beta and q
beta_value_expr = (beta == a_v/n)
q_value_expr = (q == b_q/n)
\end{sagesilent}
\renewcommand{\aa}{{a_v}}
\newcommand{\bb}{{b_q}}
Suppose $\beta = \frac{\aa}{n}$ for some coprime $n \in \NN,\aa \in \ZZ$.
Then fix a value of $q$:
\begin{equation}
\in
\frac{1}{n} \ZZ
\cap [0, \chern_1^{\beta}(F)]
\end{equation}
as noted at the beginning of this section (\ref{sec:refinement}).
Substituting the current values of $q$ and $\beta$ into the condition for the
radius of the pseudo-wall being positive
(eqn \ref{eqn:positive_rad_d_bound_betamin}) we get:
\begin{equation}
\label{eqn:positive_rad_condition_in_terms_of_q_beta}
\sage{positive_radius_condition.subs([q_value_expr,beta_value_expr]).factor()}
\in
\frac{1}{2n^2}\ZZ
\end{equation}
# placeholder for the specific values of k (start with 1):
var("kappa", domain="real")
assymptote_gap_condition1 = (kappa/(2*n^2) < bgmlv2_d_upperbound_exp_term)
assymptote_gap_condition2 = (kappa/(2*n^2) < bgmlv3_d_upperbound_exp_term_alt2)
)
assert r_upper_bound1.lhs() == r
r_upper_bound2 = (
* (r-R) * 2*n^2 / kappa + R
)
assert r_upper_bound2.lhs() == r
\end{sagesilent}
\begin{theorem}[Bound on $r$ \#1]
\label{thm:rmax_with_uniform_eps}
Let $v = (R,C\ell,D\ell^2)$ be a fixed Chern character. Then the ranks of the
pseudo-semistabilizers for $v$ with
$\chern_1^\beta = q$
are bounded above by the following expression.
\bgroup
\def\psi{\chern_1^{\beta}(F)}
\sage{r_upper_bound1.subs(kappa==1).rhs()}, \:\:
\sage{r_upper_bound2.subs(kappa==1).rhs()}
Taking the maximum of this expression over
$q \in [0, \chern_1^{\beta}(F)]\cap \frac{1}{n}\ZZ$
would give an upper bound for the ranks of pseudo-semistabilizers for $v$.
Both $d$ and the lower bound in
(eqn \ref{eqn:positive_rad_condition_in_terms_of_q_beta})
are elements of $\frac{1}{2n^2}\ZZ$.
So, if any of the two upper bounds on $d$ come to within
$\frac{1}{2n^2}$ of this lower bound, then there are no solutions for
Hence any corresponding $r$ cannot be a rank of a
pseudo-semistabilizer for $v$.
To avoid this, we must have,
considering equations
\ref{eqn:bgmlv2_d_bound_betamin},
\ref{eqn:bgmlv3_d_bound_betamin},
\ref{eqn:positive_rad_d_bound_betamin}.
\bgroup
\let\originalepsilon\epsilon
\renewcommand\epsilon{{\originalepsilon_{v}}}
\begin{sagesilent}
var("epsilon")
var("chbv") # symbol to represent \chern_1^{\beta}(v)
# Tightness conditions:
bounds_too_tight_condition1 = (
bgmlv2_d_upperbound_exp_term
< epsilon
)
bounds_too_tight_condition2 = (
bgmlv3_d_upperbound_exp_term_alt.subs(chbv==0)
< epsilon
)
\end{sagesilent}
\bgroup
\def\psi{\chern_1^{\beta}(F)}
\begin{equation}
\min\left(
\sage{bgmlv2_d_upperbound_exp_term},
\sage{bgmlv3_d_upperbound_exp_term_alt2}
\noindent
This is equivalent to:
\bgroup
\def\psi{\chern_1^{\beta}(F)}
\label{eqn:thm-bound-for-r-impossible-cond-for-r}
r \leq
r_upper_bound1.subs(kappa==1).rhs()
r_upper_bound2.subs(kappa==1).rhs()
\egroup % end scope where epsilon redefined
q_sol = solve(
r_upper_bound1.subs(kappa==1).rhs()
== r_upper_bound2.subs(kappa==1).rhs()
, q)[0].rhs()
r_upper_bound_all_q = (
r_upper_bound1.rhs()
.expand()
.subs(q==q_sol)
.subs(psi**2 == Delta/nu^2)
.subs(1/psi**2 == nu^2/Delta)
)
\end{sagesilent}
Luke Naylor
committed
\begin{corrolary}[Bound on $r$ \#2]
\label{cor:direct_rmax_with_uniform_eps}
Luke Naylor
committed
Let $v$ be a fixed Chern character and
Luke Naylor
committed
Then the ranks of the pseudo-semistabilizers for $v$
are bounded above by the following expression.
\bgroup
\let\originalDelta\Delta
\renewcommand\Delta{{\originalDelta(v)}}
\begin{equation*}
\sage{r_upper_bound_all_q.expand()}
\end{equation*}
Luke Naylor
committed
\egroup
\end{corrolary}
\begin{proof}
\bgroup
\def\psi{\chern_1^{\beta}(F)}
\let\originalDelta\Delta
The ranks of the pseudo-semistabilizers for $v$ are bounded above by the
maximum over $q\in [0, \chern_1^{\beta}(F)]$ of the expression in theorem
\ref{thm:rmax_with_uniform_eps}.
Noticing that the expression is a maximum of two quadratic functions in $q$:
\begin{equation*}
f_1(q)\coloneqq\sage{r_upper_bound1.subs(kappa==1).rhs()} \qquad
f_2(q)\coloneqq\sage{r_upper_bound2.subs(kappa==1).rhs()}
\end{equation*}
These have their minimums at $q=0$ and $q=\chern_1^{\beta}(F)$ respectively.
It suffices to find their intersection in
$q\in [0, \chern_1^{\beta}(F)]$, if it exists,
and evaluating on of the $f_i$ there.
The intersection exists, provided that
$f_1(\chern_1^{\beta}(F))>f_2(\chern_1^{\beta}(F))=R$,
or equivalently,
$R \leq n^2{\chern_1^{\beta}(F)}^2$.
Setting $f_1(q)=f_2(q)$ yields
$q=\sage{q_sol.expand()}$.
And evaluating $f_1$ at this $q$-value gives:
$\sage{r_upper_bound_all_q.expand().subs([nu==1,Delta==psi^2])}$.
Finally, noting that $\originalDelta(v)=\psi^2\ell^2$, we get the bound as
\begin{example}[$v=(3, 2\ell, -2)$ on $\PP^2$]
\label{exmpl:recurring-second}
Just like in example \ref{exmpl:recurring-first}, take
$\ell=c_1(\mathcal{O}(1))$ as the standard polarization on $\PP^2$, so
that $m=2$, $\beta=\sage{recurring.b}$,
giving $n=\sage{recurring.b.denominator()}$.
\begin{sagesilent}
recurring.n = recurring.b.denominator()
recurring.bgmlv = recurring.chern.Q_tilt()
corrolary_bound = (
r_upper_bound_all_q.expand()
.subs(Delta==recurring.bgmlv)
.subs(nu==1) ## \ell^2=1 on P^2
.subs(R==recurring.chern.ch[0])
.subs(n==recurring.n)
)
\end{sagesilent}
Using the above corrolary \ref{cor:direct_rmax_with_uniform_eps}, we get that
the ranks of tilt semistabilizers for $v$ are bounded above by
$\sage{corrolary_bound} \approx \sage{float(corrolary_bound)}$,
which is much closer to real maximum 25 than the original bound 144.
\end{example}
1757
1758
1759
1760
1761
1762
1763
1764
1765
1766
1767
1768
1769
1770
1771
1772
1773
1774
1775
1776
1777
\begin{example}[extravagant example: $v=(29, 13\ell, -3/2)$ on $\PP^2$]
\label{exmpl:extravagant-second}
Just like in example \ref{exmpl:extravagant-first}, take
$\ell=c_1(\mathcal{O}(1))$ as the standard polarization on $\PP^2$, so
that $m=2$, $\beta=\sage{extravagant.b}$,
giving $n=\sage{extravagant.b.denominator()}$.
\begin{sagesilent}
extravagant.n = extravagant.b.denominator()
extravagant.bgmlv = extravagant.chern.Q_tilt()
corrolary_bound = (
r_upper_bound_all_q.expand()
.subs(Delta==extravagant.bgmlv)
.subs(nu==1) ## \ell^2=1 on P^2
.subs(R==extravagant.chern.ch[0])
.subs(n==extravagant.n)
)
\end{sagesilent}
Using the above corrolary \ref{cor:direct_rmax_with_uniform_eps}, we get that
the ranks of tilt semistabilizers for $v$ are bounded above by
$\sage{corrolary_bound} \approx \sage{float(corrolary_bound)}$,
which is much closer to real maximum $\sage{extravagant.actual_rmax}$ than the
original bound 215296.
%% refinements using specific values of q and beta
These bound can be refined a bit more by considering restrictions from the
possible values that $r$ take.
Furthermore, the proof of theorem \ref{thm:rmax_with_uniform_eps} uses the fact
that, given an element of $\frac{1}{2n^2}\ZZ$, the closest non-equal element of
$\frac{1}{2}\ZZ$ is at least $\frac{1}{2n^2}$ away. However this a
conservative estimate, and a larger gap can sometimes be guaranteed if we know
this value of $\frac{1}{2n^2}\ZZ$ explicitly.
The expressions that will follow will be a bit more complicated and have more
parts which depend on the values of $q$ and $\beta$, even their numerators
$\aa,\bb$ specifically. The upcoming theorem (TODO ref) is less useful as a
`clean' formula for a bound on the ranks of the pseudo-semistabilizers, but has a
purpose in the context of writing a computer program to find
pseudo-semistabilizers. Such a program would iterate through possible values of
$q$, then iterate through values of $r$ within the bounds (dependent on $q$),
which would then determine $c$, and then find the corresponding possible values
for $d$.
Firstly, we only need to consider $r$-values for which $c\coloneqq\chern_1(E)$ is
\begin{equation}
c =
\sage{c_in_terms_of_q.subs([q_value_expr,beta_value_expr])}
\in \ZZ
\end{equation}
\noindent
That is, $r \equiv -\aa^{-1}\bb$ mod $n$ ($\aa$ is coprime to
$n$, and so invertible mod $n$).
positive_radius_condition
.rhs()
.subs([q_value_expr,beta_value_expr])
.factor()
.numerator()
\end{sagesilent}
\noindent
Let $\aa^{'}$ be an integer representative of $\aa^{-1}$ in $\ZZ/n\ZZ$.
Next, we seek to find a larger $\epsilon$ to use in place of $\epsilon_F$ in the
proof of theorem \ref{thm:rmax_with_uniform_eps}:
\begin{lemmadfn}[
Finding a better alternative to $\epsilon_v$:
$\epsilon_{v,q}$
]
\label{lemdfn:epsilon_q}
Suppose $d \in \frac{1}{2}\ZZ$ satisfies the condition in
eqn \ref{eqn:positive_rad_condition_in_terms_of_q_beta}.
That is:
\begin{equation*}
\sage{positive_radius_condition.subs([q_value_expr,beta_value_expr]).factor()}
\end{equation*}
\noindent
Then we have:
\begin{equation}
\label{eqn:epsilon_q_lemma_prop}
d - \frac{(\aa r + 2\bb)\aa}{2n^2}
\geq \epsilon_{v,q} \geq \epsilon_v > 0
\end{equation}
\noindent
Where $\epsilon_{v,q}$ is defined as follows:
\begin{equation*}
\frac{k_{q}}{2n^2}
\end{equation*}
with $k_{v,q}$ being the least $k\in\ZZ_{>0}$ satisfying $k \equiv -\aa\bb \mod n$
\end{lemmadfn}
\begin{proof}
- \frac{
(\aa r+2\bb)\aa
}{
2n^2
}
\quad \text{for some } x \in \ZZ
\span \span \span \span \span
\label{eqn:finding_better_eps_problem}
\aa^2 \aa^{-1}\bb - 2\aa\bb
\mod n
\label{eqn:better_eps_problem_k_mod_gcd2n2_a2mn}
\label{eqn:better_eps_problem_k_mod_n}
In our situation, we want to find the least $k$ satisfying
eqn \ref{eqn:finding_better_eps_problem}.
Since such a $k$ must also satisfy eqn \ref{eqn:better_eps_problem_k_mod_n},
we can pick the smallest $k_{q,1} \in \ZZ_{>0}$ which satisfies this new condition
(a computation only depending on $q$ and $\beta$, but not $r$).
We are then guaranteed that $k_{v,q}$ is less than any $k$ satisfying eqn
\ref{eqn:finding_better_eps_problem}, giving the first inequality in eqn
\ref{eqn:epsilon_q_lemma_prop}.
Furthermore, $k_{v,q}\geq 1$ gives the second part of the inequality:
$\epsilon_{v,q}\geq\epsilon_v$, with equality when $k_{v,q}=1$.
\begin{theorem}[Bound on $r$ \#3]
\label{thm:rmax_with_eps1}
Let $v$ be a fixed Chern character, with $\frac{a_v}{n}=\beta\coloneqq\beta(v)$
rational and expressed in lowest terms.
Then the ranks $r$ of the pseudo-semistabilizers $u$ for $v$ with
$\chern_1^\beta(u) = q = \frac{b_q}{n}$
are bounded above by the following expression:
var("delta", domain="real") # placeholder symbol to be replaced by k_{q,i}
\def\kappa{k_{v,q}}
\def\psi{\chern_1^{\beta}(F)}
\begin{align*}
\min
\left(
\sage{r_upper_bound1.rhs()}, \:\:
\sage{r_upper_bound2.rhs()}
\right)
\end{align*}
\egroup
Where $k_{v,q}$ is defined as in definition/lemma \ref{lemdfn:epsilon_q},
and $R = \chern_0(v)$
Furthermore, if $\aa \not= 0$ then
$r \equiv \aa^{-1}b_q (\mod n)$.
\begin{example}[$v=(3, 2\ell, -2)$ on $\PP^2$]
\label{exmpl:recurring-third}
Just like in examples \ref{exmpl:recurring-first} and
\ref{exmpl:recurring-second},
take $\ell=c_1(\mathcal{O}(1))$ as the standard polarization on $\PP^2$, so that
$\beta=\sage{recurring.b}$, giving $n=\sage{recurring.b.denominator()}$
and $\chern_1^{\sage{recurring.b}}(F) = \sage{recurring.twisted.ch[1]}$.
The (non-exclusive) upper bounds for $r\coloneqq\chern_0(u)$ of a tilt semistabilizer $u$ of $v$
in terms of the possible values for $q\coloneqq\chern_1^{\beta}(u)$ are as follows:
\begin{sagesilent}
import numpy as np
def bound_comparisons(example):
n = example.b.denominator()
a_v = example.b.numerator()
def theorem_bound(v_twisted, q_val, k):
return int(min(
n^2*q_val^2/k
,
v_twisted.ch[0]
+ n^2*(v_twisted.ch[1] - q_val)^2/k
))
def k(n, a_v, b_q):
n = int(n)
a_v = int(a_v)
b_q = int(b_q)
k = -a_v*b_q % n
return k if k > 0 else k + n
b_qs = list(range(example.twisted.ch[1]*n+1))
qs = list(map(lambda x: x/n,b_qs))
ks = list(map(lambda b_q: k(n, a_v, b_q), b_qs))
theorem2_bounds = [
theorem_bound(example.twisted, q_val, 1)
for q_val in qs
]
theorem3_bounds = [
theorem_bound(example.twisted, q_val, k)
for q_val, k in zip(qs,ks)
return qs, theorem2_bounds, theorem3_bounds
qs, theorem2_bounds, theorem3_bounds = bound_comparisons(recurring)
\directlua{ table_width = 3*4+1 }
\begin{tabular}{l\directlua{for i=0,table_width-1 do tex.sprint([[|c]]) end}}
$q=\chern_1^\beta(u)$
\directlua{for i=0,table_width-1 do
local cell = [[&$\noexpand\sage{qs[]] .. i .. "]}$"
tex.sprint(cell)