Imports
/- Copyright (c) 2026 Gregory J. Loges. All rights reserved. Released under Apache 2.0 license as described in the file LICENSE. Authors: Gregory J. Loges -/ module public import Physlib.QuantumMechanics.Operators.Multiplication public import Physlib.QuantumMechanics.HilbertSpaces.SpaceD.PolyBddSchwartzSubmodule public import Physlib.SpaceAndTime.Space.Norm.Regularized

Position operators

i. Overview

In this module we introduce several position operators for quantum mechanics on Space d.

ii. Key results

Definitions:

    positionCLM : (components of) the position vector operator acting on Schwartz maps 𝓢(Space d, ℂ) by multiplication by xᵢ.

    radiusRegPowCLM : operator acting on Schwartz maps by multiplication by (‖x‖² + ε²)^(s/2), a smooth regularization of ‖x‖ˢ.

    positionOperator : a self-adjoint multiplication operator acting on SpaceDHilbertSpace d.

    readiusRegPowOperator : a self-adjoint multiplication operator acting on SpaceDHilbertSpace d.

Notation:

    𝐱 for positionCLM

    𝐫₀ for radiusRegPowCLM

    𝐫 for radiusPowLM

iii. Table of contents

    A. Schwartz operators

      A.1. Position vector

      A.2. Radius powers (regularized)

      A.3. Radius powers

        A.3.1. As limit of regularized operators

    B. Unbounded operators

      B.1. Position vector

      B.2. Radius powers (regularized)

      B.3. Radius powers

        B.3.1. As limit of regularized operators

iv. References

@[expose] public section

A. Schwartz operators

A.1. Position vector

Component i of the position operator is the continuous linear map from 𝓢(Space d, ℂ) to itself which maps ψ to xᵢψ.

def positionCLM : 𝓢(Space d, ) →L[] 𝓢(Space d, ) := SchwartzMap.smulLeftCLM (Complex.ofRealCLM ∘L coordCLM i)
@[inherit_doc positionCLM] notation "𝐱" => positionCLM@[inherit_doc positionCLM] notation "𝐱[" d' "]" => positionCLM (d := d')lemma positionCLM_apply_fun (ψ : 𝓢(Space d, )) : 𝐱 i ψ = (fun x : Space d x i) ψ := d:i:Fin dψ:𝓢(Space d, )((𝐱 i) ψ) = (fun x => x.val i) ψ d:i:Fin dψ:𝓢(Space d, )x✝:Space d((𝐱 i) ψ) x✝ = ((fun x => x.val i) ψ) x✝ d:i:Fin dψ:𝓢(Space d, )x✝:Space d((smulLeftCLM (Complex.ofRealCLM ∘SL coordCLM i)) ψ) x✝ = ((fun x => x.val i) ψ) x✝ erw [smulLeftCLM_apply_apply (g := Complex.ofRealCLM (coordCLM i)) (d:i:Fin dψ:𝓢(Space d, )x✝:Space dHasTemperateGrowth (Complex.ofRealCLM (coordCLM i)) All goals completed! 🐙)d:i:Fin dψ:𝓢(Space d, )x✝:Space d(Complex.ofRealCLM (coordCLM i)) x✝ ψ x✝ = ((fun x => x.val i) ψ) x✝ All goals completed! 🐙@[simp] lemma positionCLM_apply (ψ : 𝓢(Space d, )) (x : Space d) : 𝐱 i ψ x = x i * ψ x := d:i:Fin dψ:𝓢(Space d, )x:Space d((𝐱 i) ψ) x = (x.val i) * ψ x All goals completed! 🐙

A.2. Radius powers (regularized)

The radius operator to power s, regularized by ε ≠ 0, is the continuous linear map from 𝓢(Space d, ℂ) to itself which maps ψ to (‖x‖² + ε²)^(s/2) • ψ.

def radiusRegPowCLM {d : } (ε : ˣ) (s : ) : 𝓢(Space d, ) →L[] 𝓢(Space d, ) := SchwartzMap.smulLeftCLM (Complex.ofReal normRegularizedPow d ε s)
@[inherit_doc radiusRegPowCLM] notation "𝐫₀" => radiusRegPowCLM@[inherit_doc radiusRegPowCLM] notation "𝐫₀[" d' "]" => radiusRegPowCLM (d := d')lemma radiusRegPowCLM_apply_fun {d : } (ε : ˣ) (s : ) (ψ : 𝓢(Space d, )) : 𝐫₀ ε s ψ = fun x (x ^ 2 + ε ^ 2) ^ (s / 2) ψ x := d:ε:ˣs:ψ:𝓢(Space d, )((𝐫₀ ε s) ψ) = fun x => (x ^ 2 + ε ^ 2) ^ (s / 2) ψ x d:ε:ˣs:ψ:𝓢(Space d, )x:Space d((𝐫₀ ε s) ψ) x = (x ^ 2 + ε ^ 2) ^ (s / 2) ψ x d:ε:ˣs:ψ:𝓢(Space d, )x:Space d((smulLeftCLM (Complex.ofReal normRegularizedPow d (↑ε) s)) ψ) x = ((x ^ 2 + ε ^ 2) ^ (s / 2)) * ψ x d:ε:ˣs:ψ:𝓢(Space d, )x:Space dHasTemperateGrowth (Complex.ofReal normRegularizedPow d (↑ε) s) exact HasTemperateGrowth.comp (d:ε:ˣs:ψ:𝓢(Space d, )x:Space dHasTemperateGrowth Complex.ofReal All goals completed! 🐙) (normRegularizedPow_hasTemperateGrowth d ε s)All goals completed! 🐙@[simp] lemma radiusRegPowCLM_comp_eq {d : } (ε : ˣ) (s t : ) : 𝐫₀[d] ε s ∘L 𝐫₀ ε t = 𝐫₀ ε (s+t) := d:ε:ˣs:t:𝐫₀ ε s ∘SL 𝐫₀ ε t = 𝐫₀ ε (s + t) d:ε:ˣs:t:ψ:𝓢(Space d, )x:Space d((𝐫₀ ε s ∘SL 𝐫₀ ε t) ψ) x = ((𝐫₀ ε (s + t)) ψ) x All goals completed! 🐙@[simp] lemma radiusRegPowCLM_zero {d : } (ε : ˣ) : 𝐫₀ ε 0 = ContinuousLinearMap.id 𝓢(Space d, ) := d:ε:ˣ𝐫₀ ε 0 = ContinuousLinearMap.id 𝓢(Space d, ) d:ε:ˣx✝¹:𝓢(Space d, )x✝:Space d((𝐫₀ ε 0) x✝¹) x✝ = ((ContinuousLinearMap.id 𝓢(Space d, )) x✝¹) x✝ All goals completed! 🐙lemma positionSqCLM_eq {d : } (ε : ˣ) : i, 𝐱 i ∘L 𝐱 i = 𝐫₀ ε 2 - ε.1 ^ 2 ContinuousLinearMap.id 𝓢(Space d, ) := d:ε:ˣ i, 𝐱 i ∘SL 𝐱 i = 𝐫₀ ε 2 - ε ^ 2 ContinuousLinearMap.id 𝓢(Space d, ) d:ε:ˣx✝¹:𝓢(Space d, )x✝:Space d((∑ i, 𝐱 i ∘SL 𝐱 i) x✝¹) x✝ = ((𝐫₀ ε 2 - ε ^ 2 ContinuousLinearMap.id 𝓢(Space d, )) x✝¹) x✝ All goals completed! 🐙

A.3. Radius powers

The radius operator to power s is the linear map from 𝓢(Space d, ℂ) to Space d → ℂ that maps ψ to x ↦ ‖x‖ˢψ(x) (which is 'nearly' Schwartz for general s).

d✝:μ:Measure (Space d✝)i:Fin d✝d:s:x✝¹:x✝:𝓢(Space d, )(fun x => x ^ s) (x✝¹ x✝) = (fun x => x ^ s) (RingHom.id ) x✝¹ x✝; All goals completed! 🐙
@[inherit_doc radiusPowLM] notation "𝐫" => radiusPowLM@[inherit_doc radiusPowLM] notation "𝐫[" d' "]" => radiusPowLM (d := d')lemma radiusPowLM_apply_fun {d : } (s : ) (ψ : 𝓢(Space d, )) : 𝐫 s ψ = fun x x ^ s ψ x := rflAll goals completed! 🐙

x ↦ ‖x‖ˢψ(x) is smooth away from x = 0.

d:s:n:ℕ∞ψ:𝓢(Space d, )x:Space dhx:x 0h: (x : Space d), x ^ s = inner x x ^ (s / 2)ContDiffAt (↑n) (fun x => x ^ s) x d:s:n:ℕ∞ψ:𝓢(Space d, )x:Space dhx:x 0h: (x : Space d), x ^ s = inner x x ^ (s / 2)ContDiffAt (↑n) (fun x => inner x x ^ (s / 2)) x exact ContDiffAt.rpow_const_of_ne (d:s:n:ℕ∞ψ:𝓢(Space d, )x:Space dhx:x 0h: (x : Space d), x ^ s = inner x x ^ (s / 2)ContDiffAt (↑n) (fun x => inner x x) x All goals completed! 🐙) (inner_self_ne_zero.mpr hx)

x ↦ ‖x‖ˢψ(x) is strongly measurable.

d:s:ψ:𝓢(Space d, )StronglyMeasurable fun x => x ^ s ψ x exact StronglyMeasurable.smul (f := fun x : Space d => x ^ s) (d:s:ψ:𝓢(Space d, )StronglyMeasurable fun x => x ^ s All goals completed! 🐙) ψ.continuous.stronglyMeasurable

x ↦ ‖x‖ˢψ(x) is square-integrable provided s is not too negative.

d:s:ψ:𝓢(Space d, )a:hψ✝:ψ PolyBddSchwartzMap d ah:0 < d + 2 * (a + s)hd:d > 0this:NeZero dlintegral_lt_top_aux: {S : Set (Space d)} {C p : }, IntegrableOn (fun x => x ^ p) S volume {G : Space d ENNReal}, ∫⁻ (x : Space d) in S, G x ∫⁻ (x : Space d) in S, C ^ 2‖ₑ * x ^ p‖ₑ ∫⁻ (x : Space d) in S, G x < C:hC_pos:0 < ChC: (x : Space d), x ^ (s.toNat + d) * ψ x Cx:Space dhx:1 xhx':0 < x:ψ x C * x ^ (-(s.toNat + d))s s.toNat exact (Int.le_ceil s).trans (d:s:ψ:𝓢(Space d, )a:hψ✝:ψ PolyBddSchwartzMap d ah:0 < d + 2 * (a + s)hd:d > 0this:NeZero dlintegral_lt_top_aux: {S : Set (Space d)} {C p : }, IntegrableOn (fun x => x ^ p) S volume {G : Space d ENNReal}, ∫⁻ (x : Space d) in S, G x ∫⁻ (x : Space d) in S, C ^ 2‖ₑ * x ^ p‖ₑ ∫⁻ (x : Space d) in S, G x < C:hC_pos:0 < ChC: (x : Space d), x ^ (s.toNat + d) * ψ x Cx:Space dhx:1 xhx':0 < x:ψ x C * x ^ (-(s.toNat + d))s s.toNat All goals completed! 🐙)
A.3.1. As limit of regularized operators

Neighborhoods of "0" in the non-zero reals, i.e. those sets containing (-ε,0) ∪ (0,ε) ⊆ ℝˣ for some ε > 0.

abbrev nhdsZeroUnits : Filter ˣ := comap (Units.coeHom ) (nhds 0)
instance : NeBot nhdsZeroUnits := d:μ:Measure (Space d)i:Fin dnhdsZeroUnits.NeBot d:μ:Measure (Space d)i:Fin dt:Set ht:t nhds 0 a, (Units.coeHom ) a t d:μ:Measure (Space d)i:Fin dt:Set ht:t nhds 0ε:hε_pos:ε > 0:Metric.ball 0 ε t a, (Units.coeHom ) a t use Units.mk0 (ε / 2) (d:μ:Measure (Space d)i:Fin dt:Set ht:t nhds 0ε:hε_pos:ε > 0:Metric.ball 0 ε tε / 2 0 All goals completed! 🐙) exact (d:μ:Measure (Space d)i:Fin dt:Set ht:t nhds 0ε:hε_pos:ε > 0:Metric.ball 0 ε t(Units.coeHom ) (Units.mk0 (ε / 2) ) Metric.ball 0 ε All goals completed! 🐙)

𝐫[ε,s] ψ converges pointwise to 𝐫[s] ψ as ε → 0 except perhaps at x = 0.

d:s:ψ:𝓢(Space d, )x:Space dhx:x 0hpow:x ^ s = (x ^ 2 + 0 ^ 2) ^ (s / 2)Tendsto (fun ε => ((𝐫₀ ε s) ψ) x) nhdsZeroUnits (nhds ((𝐫 s) ψ x)) d:s:ψ:𝓢(Space d, )x:Space dhx:x 0hpow:x ^ s = (x ^ 2 + 0 ^ 2) ^ (s / 2)Tendsto (fun ε => ((x ^ 2 + ε ^ 2) ^ (s / 2)) * ψ x) nhdsZeroUnits (nhds (((x ^ 2 + 0 ^ 2) ^ (s / 2)) * ψ x)) d:s:ψ:𝓢(Space d, )x:Space dhx:x 0hpow:x ^ s = (x ^ 2 + 0 ^ 2) ^ (s / 2)Tendsto (fun ε => (x ^ 2 + ε ^ 2) ^ (s / 2)) nhdsZeroUnits (nhds ((x ^ 2 + 0 ^ 2) ^ (s / 2))) refine Tendsto.rpow_const ?_ (Or.inl <| d:s:ψ:𝓢(Space d, )x:Space dhx:x 0hpow:x ^ s = (x ^ 2 + 0 ^ 2) ^ (s / 2)x ^ 2 + 0 ^ 2 0 All goals completed! 🐙) All goals completed! 🐙

𝐫[ε,s] ψ converges pointwise to 𝐫[s] ψ as ε → 0 provided 𝐫[ε,s] ψ 0 is bounded.

d:s:ψ:𝓢(Space d, )hs:0 sthis:0 ^ s = (0 ^ 2) ^ (s / 2)Tendsto (fun i => ((i ^ 2) ^ (s / 2)) * ψ 0) nhdsZeroUnits (nhds (((0 ^ 2) ^ (s / 2)) * ψ 0)) d:s:ψ:𝓢(Space d, )hs:0 sthis:0 ^ s = (0 ^ 2) ^ (s / 2)Tendsto (fun i => (i ^ 2) ^ (s / 2)) nhdsZeroUnits (nhds ((0 ^ 2) ^ (s / 2))) exact Tendsto.rpow_const (Tendsto.pow tendsto_comap 2) (Or.inr <| d:s:ψ:𝓢(Space d, )hs:0 sthis:0 ^ s = (0 ^ 2) ^ (s / 2)0 s / 2 All goals completed! 🐙) d:s:ψ:𝓢(Space d, ):ψ 0 = 0Tendsto (fun i => ((𝐫₀ i s) ψ) 0) nhdsZeroUnits (nhds ((𝐫 s) ψ 0)) All goals completed! 🐙 d:s:ψ:𝓢(Space d, )h:0 s ψ 0 = 0x:Space dhx:NeZero xTendsto (fun i => ((𝐫₀ i s) ψ) x) nhdsZeroUnits (nhds ((𝐫 s) ψ x)) All goals completed! 🐙

a.e. version of radiusRegPow_tendsto_radiusPow

d:inst✝:NeZero ds:ψ:𝓢(Space d, ){x | ¬Tendsto (fun ε => ((𝐫₀ ε s) ψ) x) nhdsZeroUnits (nhds ((𝐫 s) ψ x))} {0} d:inst✝:NeZero ds:ψ:𝓢(Space d, )x:Space dhx:x {x | ¬Tendsto (fun ε => ((𝐫₀ ε s) ψ) x) nhdsZeroUnits (nhds ((𝐫 s) ψ x))}x {0} d:inst✝:NeZero ds:ψ:𝓢(Space d, )x:Space dhx:x {x | ¬Tendsto (fun ε => ((𝐫₀ ε s) ψ) x) nhdsZeroUnits (nhds ((𝐫 s) ψ x))}hx':x {0}False All goals completed! 🐙
d:inst✝:NeZero ds:ψ:𝓢(Space d, )φ:Space d t₁:Set (Space d) := {x | ¬Tendsto (fun ε => ((𝐫₀ ε s) ψ) x) nhdsZeroUnits (nhds (φ x))}t₂:Set (Space d) := {x | φ x (𝐫 s) ψ x}t₁ {0} = t₂ {0} d:inst✝:NeZero ds:ψ:𝓢(Space d, )φ:Space d t₁:Set (Space d) := {x | ¬Tendsto (fun ε => ((𝐫₀ ε s) ψ) x) nhdsZeroUnits (nhds (φ x))}t₂:Set (Space d) := {x | φ x (𝐫 s) ψ x}x:Space dx t₁ {0} x t₂ {0} d:inst✝:NeZero ds:ψ:𝓢(Space d, )φ:Space d t₁:Set (Space d) := {x | ¬Tendsto (fun ε => ((𝐫₀ ε s) ψ) x) nhdsZeroUnits (nhds (φ x))}t₂:Set (Space d) := {x | φ x (𝐫 s) ψ x}0 t₁ {0} 0 t₂ {0}d:inst✝:NeZero ds:ψ:𝓢(Space d, )φ:Space d t₁:Set (Space d) := {x | ¬Tendsto (fun ε => ((𝐫₀ ε s) ψ) x) nhdsZeroUnits (nhds (φ x))}t₂:Set (Space d) := {x | φ x (𝐫 s) ψ x}x:Space dhx:NeZero xx t₁ {0} x t₂ {0} d:inst✝:NeZero ds:ψ:𝓢(Space d, )φ:Space d t₁:Set (Space d) := {x | ¬Tendsto (fun ε => ((𝐫₀ ε s) ψ) x) nhdsZeroUnits (nhds (φ x))}t₂:Set (Space d) := {x | φ x (𝐫 s) ψ x}0 t₁ {0} 0 t₂ {0} All goals completed! 🐙 d:inst✝:NeZero ds:ψ:𝓢(Space d, )φ:Space d t₁:Set (Space d) := {x | ¬Tendsto (fun ε => ((𝐫₀ ε s) ψ) x) nhdsZeroUnits (nhds (φ x))}t₂:Set (Space d) := {x | φ x (𝐫 s) ψ x}x:Space dhx:NeZero xx t₁ {0} x t₂ {0} d:inst✝:NeZero ds:ψ:𝓢(Space d, )φ:Space d t₁:Set (Space d) := {x | ¬Tendsto (fun ε => ((𝐫₀ ε s) ψ) x) nhdsZeroUnits (nhds (φ x))}t₂:Set (Space d) := {x | φ x (𝐫 s) ψ x}x:Space dhx:NeZero xx t₁ x t₂ d:inst✝:NeZero ds:ψ:𝓢(Space d, )φ:Space d t₁:Set (Space d) := {x | ¬Tendsto (fun ε => ((𝐫₀ ε s) ψ) x) nhdsZeroUnits (nhds (φ x))}t₂:Set (Space d) := {x | φ x (𝐫 s) ψ x}x:Space dhx:NeZero xhLim:Tendsto (fun ε => ((𝐫₀ ε s) ψ) x) nhdsZeroUnits (nhds ((𝐫 s) ψ x))x t₁ x t₂ All goals completed! 🐙

B. Unbounded operators

B.1. Position vector

The operator on SpaceDHilbertSpace d acting by multiplication by fun x ↦ xᵢ.

def positionOperator : SpaceDHilbertSpace d μ →ₗ.[] SpaceDHilbertSpace d μ := 𝓜 μ (Complex.ofRealCLM ∘L Space.coordCLM i)
@[inherit_doc positionOperator] notation "𝓧" => positionOperatorlemma positionOperator_hasDenseDomain : (𝓧 μ i).HasDenseDomain := mulOperator_hasDenseDomain (d:μ:Measure (Space d)i:Fin dAEStronglyMeasurable (⇑(Complex.ofRealCLM ∘SL coordCLM i)) μ All goals completed! 🐙)lemma positionOperator_isSelfAdjoint [IsFiniteMeasureOnCompacts μ] : IsSelfAdjoint (𝓧 μ i) := mulOperator_isSelfAdjoint_ofReal (d:μ:Measure (Space d)i:Fin dinst✝:IsFiniteMeasureOnCompacts μAEStronglyMeasurable (⇑(Complex.ofRealCLM ∘SL coordCLM i)) μ All goals completed! 🐙) (d:μ:Measure (Space d)i:Fin dinst✝:IsFiniteMeasureOnCompacts μ(starRingEnd ) (Complex.ofRealCLM ∘SL coordCLM i) = (Complex.ofRealCLM ∘SL coordCLM i) d:μ:Measure (Space d)i:Fin dinst✝:IsFiniteMeasureOnCompacts μx✝:Space d((starRingEnd ) (Complex.ofRealCLM ∘SL coordCLM i)) x✝ = (Complex.ofRealCLM ∘SL coordCLM i) x✝; All goals completed! 🐙)lemma positionOperator_isUnbounded [IsFiniteMeasureOnCompacts μ] : (𝓧 μ i).IsUnbounded := LinearPMap.IsSelfAdjoint.isUnbounded (positionOperator_isSelfAdjoint μ i)

B.2. Radius powers (regularized)

The operator on SpaceDHilbertSpace d acting by multiplication by fun x ↦ (‖x‖² + ε²)^(s/2).

def radiusRegPowOperator (ε : ˣ) (s : ) : SpaceDHilbertSpace d μ →ₗ.[] SpaceDHilbertSpace d μ := 𝓜 μ (Complex.ofReal normRegularizedPow d ε s)
@[inherit_doc radiusRegPowOperator] notation "𝓡₀" => radiusRegPowOperatorlemma radiusRegPowOperator_hasDenseDomain (ε : ˣ) (s : ) : (𝓡₀ μ ε s).HasDenseDomain := mulOperator_hasDenseDomain (d:μ:Measure (Space d)ε:ˣs:AEStronglyMeasurable (Complex.ofReal normRegularizedPow d (↑ε) s) μ All goals completed! 🐙)lemma radiusRegPowOperator_isSelfAdjoint [IsFiniteMeasureOnCompacts μ] (ε : ˣ) (s : ) : IsSelfAdjoint (𝓡₀ μ ε s) := mulOperator_isSelfAdjoint_ofReal (d:μ:Measure (Space d)inst✝:IsFiniteMeasureOnCompacts με:ˣs:AEStronglyMeasurable (Complex.ofReal normRegularizedPow d (↑ε) s) μ All goals completed! 🐙) (d:μ:Measure (Space d)inst✝:IsFiniteMeasureOnCompacts με:ˣs:(starRingEnd ) Complex.ofReal normRegularizedPow d (↑ε) s = Complex.ofReal normRegularizedPow d (↑ε) s d:μ:Measure (Space d)inst✝:IsFiniteMeasureOnCompacts με:ˣs:x✝:Space d((starRingEnd ) Complex.ofReal normRegularizedPow d (↑ε) s) x✝ = (Complex.ofReal normRegularizedPow d (↑ε) s) x✝; All goals completed! 🐙)lemma radiusRegPowOperator_isUnbounded [IsFiniteMeasureOnCompacts μ] (ε : ˣ) (s : ) : (𝓡₀ μ ε s).IsUnbounded := LinearPMap.IsSelfAdjoint.isUnbounded (radiusRegPowOperator_isSelfAdjoint μ ε s)

B.3. Radius powers

The operator on SpaceDHilbertSpace d acting by multiplication by fun x ↦ ‖x‖ˢ.

def radiusPowOperator (s : ) : SpaceDHilbertSpace d μ →ₗ.[] SpaceDHilbertSpace d μ := 𝓜 μ (Complex.ofReal fun x x ^ s)
@[inherit_doc radiusPowOperator] notation "𝓡" => radiusPowOperatorlemma radiusPowOperator_hasDenseDomain (s : ) : (𝓡 μ s).HasDenseDomain := mulOperator_hasDenseDomain (Measurable.aestronglyMeasurable (d:μ:Measure (Space d)s:Measurable (Complex.ofReal fun x => x ^ s) All goals completed! 🐙))lemma radiusPowOperator_isSelfAdjoint [IsFiniteMeasureOnCompacts μ] (s : ) : IsSelfAdjoint (𝓡 μ s) := mulOperator_isSelfAdjoint_ofReal (Measurable.aestronglyMeasurable (d:μ:Measure (Space d)inst✝:IsFiniteMeasureOnCompacts μs:Measurable (Complex.ofReal fun x => x ^ s) All goals completed! 🐙)) (d:μ:Measure (Space d)inst✝:IsFiniteMeasureOnCompacts μs:((starRingEnd ) Complex.ofReal fun x => x ^ s) = Complex.ofReal fun x => x ^ s d:μ:Measure (Space d)inst✝:IsFiniteMeasureOnCompacts μs:x✝:Space d((starRingEnd ) Complex.ofReal fun x => x ^ s) x✝ = (Complex.ofReal fun x => x ^ s) x✝; All goals completed! 🐙)lemma radiusPowOperator_isUnbounded [IsFiniteMeasureOnCompacts μ] (s : ) : (𝓡 μ s).IsUnbounded := LinearPMap.IsSelfAdjoint.isUnbounded (radiusPowOperator_isSelfAdjoint μ s)private lemma add_floor_toNat_pos_aux (d : ) (s : ) : 0 < d + 2 * (1 - d / 2 - s.toNat + s) := d:s:0 < d + 2 * (1 - d / 2 - s.toNat + s) d:s:n: := 1 - d / 2 - s0 < d + 2 * (1 - d / 2 - s.toNat + s) d:s:n: := 1 - d / 2 - shn₁:1 - d / 2 - s < n + 10 < d + 2 * (1 - d / 2 - s.toNat + s) d:s:n: := 1 - d / 2 - shn₁:1 - d / 2 - s < n + 1hn₂:n n.toNat0 < d + 2 * (1 - d / 2 - s.toNat + s) All goals completed! 🐙lemma radiusPowLM_apply_polyBddSchwartz_memHS {d : } {s : } (ψ : PolyBddSchwartzSubmodule d 1 - d / 2 - s.toNat) : MemHS (𝐫[d] s ((polyBddSchwartzEquiv volume).symm ψ)) := let f := (polyBddSchwartzEquiv volume).symm ψ radiusPowLM_apply_memHS s f.1 1 - d / 2 - s.toNat f.2 (add_floor_toNat_pos_aux d s)lemma radiusPowOperator_domain_ge {d : } (s : ) : PolyBddSchwartzSubmodule d 1 - d / 2 - s.toNat (radiusPowOperator volume s).domain := d:s:PolyBddSchwartzSubmodule d (↑1 - d / 2 - s.toNat) volume (𝓡 volume s).domain d:s:ψ:(SpaceDHilbertSpace d):ψ PolyBddSchwartzSubmodule d (↑1 - d / 2 - s.toNat) volumeψ (𝓡 volume s).domain d:s:ψ:(SpaceDHilbertSpace d):ψ PolyBddSchwartzSubmodule d (↑1 - d / 2 - s.toNat) volumef:(PolyBddSchwartzMap d 1 - d / 2 - s.toNat) := (polyBddSchwartzEquiv volume).symm ψ, ψ (𝓡 volume s).domain d:s:ψ:(SpaceDHilbertSpace d):ψ PolyBddSchwartzSubmodule d (↑1 - d / 2 - s.toNat) volumef:(PolyBddSchwartzMap d 1 - d / 2 - s.toNat) := (polyBddSchwartzEquiv volume).symm ψ, MemHS ((ofReal fun x => x ^ s) ψ) d:s:ψ:(SpaceDHilbertSpace d):ψ PolyBddSchwartzSubmodule d (↑1 - d / 2 - s.toNat) volumef:(PolyBddSchwartzMap d 1 - d / 2 - s.toNat) := (polyBddSchwartzEquiv volume).symm ψ, (𝐫 s) f =ᵐ[volume] (ofReal fun x => x ^ s) ψd:s:ψ:(SpaceDHilbertSpace d):ψ PolyBddSchwartzSubmodule d (↑1 - d / 2 - s.toNat) volumef:(PolyBddSchwartzMap d 1 - d / 2 - s.toNat) := (polyBddSchwartzEquiv volume).symm ψ, MemHS ((𝐫 s) f) d:s:ψ:(SpaceDHilbertSpace d):ψ PolyBddSchwartzSubmodule d (↑1 - d / 2 - s.toNat) volumef:(PolyBddSchwartzMap d 1 - d / 2 - s.toNat) := (polyBddSchwartzEquiv volume).symm ψ, (𝐫 s) f =ᵐ[volume] (ofReal fun x => x ^ s) ψ d:s:ψ:(SpaceDHilbertSpace d):ψ PolyBddSchwartzSubmodule d (↑1 - d / 2 - s.toNat) volumef:(PolyBddSchwartzMap d 1 - d / 2 - s.toNat) := (polyBddSchwartzEquiv volume).symm ψ, (a : Space d), ((polyBddSchwartzEquiv volume) f) a = f a (𝐫 s) (↑f) a = ((ofReal fun x => x ^ s) ψ) a All goals completed! 🐙 d:s:ψ:(SpaceDHilbertSpace d):ψ PolyBddSchwartzSubmodule d (↑1 - d / 2 - s.toNat) volumef:(PolyBddSchwartzMap d 1 - d / 2 - s.toNat) := (polyBddSchwartzEquiv volume).symm ψ, MemHS ((𝐫 s) f) All goals completed! 🐙