One-range Addition Theorems for Noninteger n Slater Functions using Complete Orthonormal Sets of Exponential Type Orbitals in Standard Convention

By the use of complete orthonormal sets of ${\psi ^{(\alpha^{\ast})}}$ -exponential type orbitals ( ${\psi ^{(\alpha^{\ast})}}$ -ETOs) with integer (for α ^* = α) and noninteger self-frictional quantum number α ^*(for α ^* ≠ α) in standard convention introduced by the author, the one-range addition theorems for ${\chi }$ -noninteger n Slater type orbitals ${(\chi}$ -NISTOs) are established. These orbitals are defined as follows $$\begin{array}{ll}\psi _{nlm}^{(\alpha^*)} (\zeta ,\vec {r}) = \frac{(2\zeta )^{3/2}}{\Gamma (p_l ^* + 1)} \left[{\frac{\Gamma (q_l ^* + )}{(2n)^{\alpha ^*}(n - l - 1)!}} \right]^{1/2}e^{-\frac{x}{2}}x^{l}_1 F_1 ({-[ {n - l - 1}]; p_l ^* + 1; x})S_{lm} (\theta ,\varphi )\\ \chi _{n^*lm} (\zeta ,\vec {r}) = (2\zeta )^{3/2}\left[ {\Gamma(2n^* + 1)}\right]^{{-1}/2}x^{n^*-1}e^{-\frac{x}{2}}S_{lm}(\theta ,\varphi ),\end{array}$$ where ${x=2\zeta r, 0<\zeta <\infty , p_l ^{\ast}=2l+2-\alpha ^{\ast}, q_l ^{\ast}=n+l+1-\alpha ^{\ast}, -\infty <\alpha ^{\ast} <3 , -\infty <\alpha \leq 2,_1 F_1 }$ is the confluent hypergeometric function and ${S_{lm} (\theta ,\varphi )}$ are the complex or real spherical harmonics. The origin of the ${\psi ^{(\alpha ^{\ast})} }$ -ETOs, therefore, of the one-range addition theorems obtained in this work for ${\chi}$ -NISTOs is the self-frictional potential of the field produced by the particle itself. The obtained formulas can be useful especially in the electronic structure calculations of atoms, molecules and solids when Hartree–Fock–Roothan approximation is employed.     

Three-loop on-shell charge renormalization without integration:\Lambda _{QED}^{\overline {MS} } to four loops

Using algebraic methods, we find the three-loop relation between the bare and physical couplings of one-flavourD-dimensional QED, in terms of Γ functions and a singleF ₃₂ series, whose expansion nearD=4 is obtained, by wreath-product transformations, to the order required for five-loop calculations. Taking the limitD→4, we find that the \(\overline {MS} \) coupling \(\bar \alpha (\mu )\) satisfies the boundary condition $$\begin{gathered} \frac{{\bar \alpha (m)}}{\pi } = \frac{\alpha }{\pi } + \frac{{15}}{{16}}\frac{{\alpha ^3 }}{{\pi ^3 }} + \left\{ {\frac{{11}}{{96}}\zeta (3) - \frac{1}{3}\pi ^2 \log 2} \right. \hfill \\ \left. { + \frac{{23}}{{72}}\pi ^2 - \frac{{4867}}{{5184}}} \right\}\frac{{\alpha ^4 }}{{\pi ^4 }} + \mathcal{O}(\alpha ^5 ), \hfill \\ \end{gathered} $$ wherem is the physical lepton mass and α is the physical fine structure constant. Combining this new result for the finite part of three-loop on-shell charge renormalization with the recently revised four-loop term in the \(\overline {MS} \) β-function, we obtain $$\begin{gathered} \Lambda _{QED}^{\overline {MS} } \approx \frac{{me^{3\pi /2\alpha } }}{{(3\pi /\alpha )^{9/8} }}\left( {1 - \frac{{175}}{{64}}\frac{\alpha }{\pi } + \left\{ { - \frac{{63}}{{64}}\zeta (3)} \right.} \right. \hfill \\ \left. { + \frac{1}{2}\pi ^2 \log 2 - \frac{{23}}{{48}}\pi ^2 + \frac{{492473}}{{73728}}} \right\}\left. {\frac{{\alpha ^2 }}{{\pi ^2 }}} \right), \hfill \\ \end{gathered} $$ at the four-loop level of one-flavour QED.     

Light quark masses in quantum chromodynamics and chiral symmetry breaking

We derive model independent lower bounds for the sums of effective quark masses \(\bar m_u + \bar m_d \) and \(\bar m_u + \bar m_s \) . The bounds follow from the combination of the spectral representation properties of the hadronic axial currents two-point functions and their behavior in the deep euclidean region (known from a perturbative QCD calculation to two loops and the leading non-perturbative contribution). The bounds incorporate PCAC in the Nambu-Goldstone version. If we define the invariant masses \(\hat m\) by $$\bar m_i = \hat m_i \left( {{{\frac{1}{2}\log Q^2 } \mathord{\left/ {\vphantom {{\frac{1}{2}\log Q^2 } {\Lambda ^2 }}} \right. \kern-\nulldelimiterspace} {\Lambda ^2 }}} \right)^{{{\gamma _1 } \mathord{\left/ {\vphantom {{\gamma _1 } {\beta _1 }}} \right. \kern-\nulldelimiterspace} {\beta _1 }}} $$ and <F ²> is the vacuum expectation value of $$F^2 = \Sigma _a F_{(a)}^{\mu v} F_{\mu v(a)} $$ , we find, e.g., $$\hat m_u + \hat m_d \geqq \sqrt {\frac{{2\pi }}{3} \cdot \frac{{8f_\pi m_\pi ^2 }}{{3\left\langle {\alpha _s F^2 } \right\rangle ^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} }}} $$ ; with the value <α _u F ²?0.04GeV⁴, recently suggested by various analysis, this gives $$\hat m_u + \hat m_d \geqq 35MeV$$ . The corresponding bounds on \(\bar m_u + \bar m_s \) are obtained replacingm _π ² f _π bym _K ² f _K . The PCAC relation can be inverted, and we get upper bounds on the spontaneous masses, \(\hat \mu \) : $$\hat \mu \leqq 170MeV$$ where \(\hat \mu \) is defined by $$\left\langle {\bar \psi \psi } \right\rangle \left( {Q^2 } \right) = \left( {{{\frac{1}{2}\log Q^2 } \mathord{\left/ {\vphantom {{\frac{1}{2}\log Q^2 } {\Lambda ^2 }}} \right. \kern-\nulldelimiterspace} {\Lambda ^2 }}} \right)^d \hat \mu ^3 ,d = {{12} \mathord{\left/ {\vphantom {{12} {\left( {33 - 2n_f } \right)}}} \right. \kern-\nulldelimiterspace} {\left( {33 - 2n_f } \right)}}$$ .     

Neutral strange particle exclusive production in charged current high-energy antineutrino interactions

The cross section of the quasi-elastic reactions \(\bar v_\mu p \to \mu ^ + \Lambda (\Sigma ^0 )\) in the energy range 5–100 GeV is determined from Fermilab 15′ bubble chamber antineutrino data. TheQ ² analysis of quasi-elastic Λ events yieldsM _A=1.0±0.3 GeV/c² for the axial mass value. With zero µΛ K ⁰ events observed, the 90% confidence level upper limit \(\sigma (\bar v_\mu p \to \mu ^ + \Lambda {\rm K}^0 )< 2.0 \cdot 10^{ - 40} cm^2 \) is obtained. At the same time, we found that the cross section of reaction \(\bar v_\mu p \to \mu ^ + \Lambda {\rm K}^0 + m\pi ^0 \) is equal to \(\left( {3.9\begin{array}{*{20}c} { + 1.6} \\ { - 1.3} \\ \end{array} } \right) \cdot 10^{ - 40} cm^2 \) .     

The NLS Equation in Dimension One with Spatially Concentrated Nonlinearities: the Pointlike Limit

In the present paper, we study the following scaled nonlinear Schrödinger equation (NLS) in one space dimension: $$ i\frac{\rm d}{{\rm d}t}\psi^{\varepsilon}(t)=-\Delta\psi^{\varepsilon}(t) +\frac{1}{\varepsilon}V\left(\frac{x}{\varepsilon} \right)|\psi^{\varepsilon}(t)|^{2\mu}\psi^{\varepsilon}(t)\quad \varepsilon > 0\,\quad V\in L^1(\mathbb{R},(1+|x|){\rm d}x) \cap L^\infty(\mathbb{R}).$$ This equation represents a nonlinear Schrödinger equation with a spatially concentrated nonlinearity. We show that in the limit \({\varepsilon\to 0}\) the weak (integral) dynamics converges in \({H^1(\mathbb{R})}\) to the weak dynamics of the NLS with point-concentrated nonlinearity: $$ i\frac{{\rm d}}{{\rm d}t} \psi(t) =H_{\alpha} \psi(t) .$$ where H _α is the Laplacian with the nonlinear boundary condition at the origin \({\psi'(t,0+)-\psi'(t,0-)=\alpha|\psi(t,0)|^{2\mu}\psi(t,0)}\) and \({\alpha=\int_{\mathbb{R}}V{\rm d}x}\) . The convergence occurs for every \({\mu\in \mathbb{R}^+}\) if V ≥  0 and for every  \({\mu\in (0,1)}\) otherwise. The same result holds true for a nonlinearity with an arbitrary number N of concentration points.     

Universality of the Local Regime for the Block Band Matrices with a Finite Number of Blocks

We consider the block band matrices, i.e. the Hermitian matrices $H_N$ , $N=|\Lambda |W$ with elements $H_{jk,\alpha \beta }$ , where $j,k \in \Lambda =[1,m]^d\cap \mathbb {Z}^d$ (they parameterize the lattice sites) and $\alpha , \beta = 1,\ldots , W$ (they parameterize the orbitals on each site). The entries $H_{jk,\alpha \beta }$ are random Gaussian variables with mean zero such that $\langle H_{j_1k_1,\alpha _1\beta _1}H_{j_2k_2,\alpha _2\beta _2}\rangle =\delta _{j_1k_2}\delta _{j_2k_1} \delta _{\alpha _1\beta _2}\delta _{\beta _1\alpha _2} J_{j_1k_1},$ where $J=1/W+\alpha \Delta /W$ , $\alpha < 1/4d$ . This matrices are the special case of Wegner’s $W$ -orbital models. Assuming that the number of sites $|\Lambda |$ is finite, we prove universality of the local eigenvalue statistics of $H_N$ for the energies $|\lambda _0|< \sqrt{2}$ .     

Time decay of solutions to the cauchy problem for time-dependent Schrödinger-Hartree equations

We consider the time-dependent Schrödinger-Hartree equation (1) $$iu_t + \Delta u = \left( {\frac{1}{r}*|u|^2 } \right)u + \lambda \frac{u}{r},(t, x) \in \mathbb{R} \times \mathbb{R}^3 ,$$ (2) $$u(0,x) = \phi (x) \in \Sigma ^{2,2} ,x \in \mathbb{R}^3 ,$$ where λ≧0 and \(\Sigma ^{2,2} = \{ g \in L^2 ;\parallel g\parallel _{\Sigma ^{2,2} }^2 = \sum\limits_{|a| \leqq 2} {\parallel D^a g\parallel _2^2 + \sum\limits_{|\beta | \leqq 2} {\parallel x^\beta g\parallel _2^2< \infty } } \} \) . We show that there exists a unique global solutionu of (1) and (2) such that $$u \in C(\mathbb{R};H^{1,2} ) \cap L^\infty (\mathbb{R};H^{2,2} ) \cap L_{loc}^\infty (\mathbb{R};\Sigma ^{2,2} )$$ with $$u \in L^\infty (\mathbb{R};L^2 ).$$ Furthermore, we show thatu has the following estimates: $$\parallel u(t)\parallel _{2,2} \leqq C,a.c. t \in \mathbb{R},$$ and $$\parallel u(t)\parallel _\infty \leqq C(1 + |t|)^{ - 1/2} ,a.e. t \in \mathbb{R}.$$      

A Strong Central Limit Theorem for a Class of Random Surfaces

This paper is concerned with d = 2 dimensional lattice field models with action ${V(\nabla\phi(\cdot))}$ , where ${V : \mathbf{R}^d \rightarrow \mathbf{R}}$ is a uniformly convex function. The fluctuations of the variable ${\phi(0) - \phi(x)}$ are studied for large |x| via the generating function given by ${g(x, \mu) = \ln \langle e^{\mu(\phi(0) - \phi(x))}\rangle_{A}}$ . In two dimensions ${g'' (x, \mu) = \partial^2g(x, \mu)/\partial\mu^2}$ is proportional to ${\ln\vert x\vert}$ . The main result of this paper is a bound on ${g''' (x, \mu) = \partial^3 g(x, \mu)/\partial \mu^3}$ which is uniform in ${\vert x \vert}$ for a class of convex V. The proof uses integration by parts following Helffer–Sjöstrand and Witten, and relies on estimates of singular integral operators on weighted Hilbert spaces.     

Uniqueness of Self-similar Solutions to Smoluchowski’s Coagulation Equations for Kernels that are Close to Constant

We consider self-similar solutions to Smoluchowski’s coagulation equation for kernels \(K=K(x,y)\) that are homogeneous of degree zero and close to constant in the sense that $$\begin{aligned} -\varepsilon \le K(x,y)-2 \le \varepsilon \Big ( \Big (\frac{x}{y}\Big )^{\alpha } + \Big (\frac{y}{x}\Big )^{\alpha }\Big ) \end{aligned}$$ for \(\alpha \in [0,1)\) . We prove that self-similar solutions with given mass are unique if \(\varepsilon \) is sufficiently small which is the first such uniqueness result for kernels that are not solvable. Our proof relies on a contraction argument in a norm that measures the distance of solutions with respect to the weak topology of measures.     

Localization in the ground state of a disordered array of quantum rotators

We consider the zero-temperature behavior of a disordered array of quantum rotators given by the finite-volume Hamiltonian: $$H_\Lambda = - \mathop \Sigma \limits_{x \in \Lambda } \frac{{h(x)}}{2}\frac{{\partial ^2 }}{{\partial \varphi (x)^2 }} - J\mathop \Sigma \limits_{\left\langle {x,y} \right\rangle \in \Lambda } \cos (\varphi (x) - \varphi (y))$$ , wherex,y∈Z ^d , 〈,〉 denotes nearest neighbors inZ ^d ;J>0 andh={h(x)>0,x∈Z ^d } are independent identically distributed random variables with common distributiondμ(h), satisfying ∫h ^?δ dμ(h)<∞ for some δ>0. We prove that for anym>0 it is possible to chooseJ(m) sufficiently small such that, if 0<J<J(m), for almost every choice ofh and everyx∈Z ^d the ground state correlation function satisfies $$\left\langle {\cos (\varphi (x) - \varphi (y))} \right\rangle \leqq C_{x,h,J} e^{ - m\left| {x - y} \right|} $$ for ally∈Z ^d withC _x,h,J <∞.     

Liouville-Type Theorems for the Forced Euler Equations and the Navier–Stokes Equations

In this paper we study the Liouville-type properties for solutions to the steady incompressible Euler equations with forces in ${\mathbb {R}^N}$ . If we assume “single signedness condition” on the force, then we can show that a ${C^1 (\mathbb {R}^N)}$ solution (v, p) with ${|v|^2+ |p| \in L^{\frac{q}{2}}(\mathbb {R}^N),\,q \in (\frac{3N}{N-1}, \infty)}$ is trivial, v = 0. For the solution of the steady Navier–Stokes equations, satisfying ${v(x) \to 0}$ as ${|x| \to \infty}$ , the condition ${\int_{\mathbb {R}^3} |\Delta v|^{\frac{6}{5}} dx < \infty}$ , which is stronger than the important D-condition, ${\int_{\mathbb {R}^3} |\nabla v|^2 dx < \infty}$ , but both having the same scaling property, implies that v = 0. In the appendix we reprove Theorem 1.1 (Chae, Commun Math Phys 273:203–215, 2007), using the self-similar Euler equations directly.     

Multi-quark structure ofU(3.1)

We interpret the recently observedU(3.1) mesons with the \(\Lambda \bar p\) + pions decays as the bound state of \(\Lambda ,\bar p\) andX ₀(1480). TheX ₀(1480) is a mesonium with \(Q^2 \bar Q^2 \) structures observed in γγ reactions and \(\bar pn\) annihilations. With this interpretation, we can understand its decay modes. Furthermore, we predict the ratio of \(\sigma (\Lambda \bar p\pi ^ + \pi ^ - )/\sigma (\Lambda \bar p\pi ^ + \pi ^ + )\) to be ?3.1 for centrally produced events and that the width of \(U^ - (\Lambda \bar p\pi ^ + \pi ^ - )\) to be greater than that of \(U^ + (\Lambda \bar p\pi ^ + \pi ^ + )\) . Both predictions seem to be in reasonable accord with the available data. We call for the detection of the \(\Lambda \bar p\pi ^ - \pi ^ - \) mode to verify the present interpretation.     

On the freund-witten adelic formula for Veneziano amplitudes

On the basis of the analysis of the adele group (Tate's formula), a regularization for the divergent infinite product ofp-adic Г-functions $$\Gamma _p (\alpha ) = \frac{{1 - p^{\alpha - 1} }}{{[ - p^{ - \alpha } }}$$ is proposed, and the adelic formula is proved $$reg\coprod\limits_{p = 2}^\infty {\Gamma _p (\alpha )} = \frac{{\zeta (\alpha )}}{{\zeta (1 - \alpha )}}$$ whereζ(α) is the Riemannζ-function.     

On the rate of quantum ergodicity I: Upper bounds

One problem in quantum ergodicity is to estimate the rate of decay of the sums $$S_k (\lambda ;A) = \frac{1}{{N(\lambda )}}\sum\limits_{\sqrt {\lambda _j } \leqq \lambda } {\left| {(A\varphi _j ,\varphi _j ) - \bar \sigma _A } \right|^k } $$ on a compact Riemannian manifold (M, g) with ergodic geodesic flow. Here, {λ _j ,? _j } are the spectral data of the Δ of(M, g), A is a 0-th order ψDO, $\bar \sigma _A $ is the (Liouville) average of its principal symbol and $N(\lambda ) = \# \{ j:\sqrt {\lambda _j } \leqq \lambda \} $ . ThatS _k (λ;A)=o(1) is proved in [S, Z.1, CV.1]. Our purpose here is to show thatS _k (λ;A)=O((logλ) ^?k/2 ) on a manifold of (possibly variable) negative curvature. The main new ingredient is the central limit theorem for geodesic flows on such spaces ([R, Si]).     

Charm Production in Antiproton-Proton Annihilation

We study the production of charmed mesons (D) and baryons (?? _c ) in antiproton-proton ${(\bar{p}p)}$ annihilation close to their respective production thresholds. The elementary charm production process is described by either baryon/meson exchange or by quark/gluon dynamics. Effects of the interactions in the initial and final states are taken into account rigorously. The calculations are performed in close analogy to our earlier study on ${\bar{p}p \to \bar{\Lambda} \Lambda}$ and ${\bar{p} p \to \bar{K} K}$ by connecting the processes via SU(4) flavor symmetry. Our predictions for the ${\bar{\Lambda}_c \Lambda_c}$ production cross section are in the order of 1 to 7 mb, i.e. a factor of around 10?C70 smaller than the corresponding cross sections for ${\bar{\Lambda} \Lambda}$ However, they are 100 to 1000 times larger than predictions of other model calculations in the literature. On the other hand, the resulting cross sections for ${\bar{D} D}$ production are found to be in the order of 10^?2 ?C 10^?1 ??b and they turned out to be comparable to those obtained in other studies.     

Structure of Neutron-Rich Be Hyper Isotopes

The deformation change of ${^{9}_\Lambda}$ Be and the low-lying states of ${^{12}_{\Lambda}}$ Be are studied by using the antisymmetrized molecular dynamics for hypernuclei (HyperAMD). In ${^{9}_{\Lambda}}$ Be, the Λ hyperon in p orbit enhances nuclear quadrupole deformation, while the Λ hyperon in s orbit reduces it. In ${^{12}_{\Lambda}}$ Be, the ground state parity inverted in ¹¹Be is reverted in ${^{12}_{\Lambda}}$ Be by adding a Λ hyperon as an impurity (impurity effect).     

Three-loop relation of quark\overline {MS} and pole masses

We calculate, exactly, the next-to-leading correction to the relation between the \(\overline {MS} \) quark mass, \(\bar m\) , and the scheme-independent pole mass,M, and obtain $$\begin{gathered} \frac{M}{{\bar m(M)}} \approx 1 + \frac{4}{3}\frac{{\bar \alpha _s (M)}}{\pi } + \left[ {16.11 - 1.04\sum\limits_{i = 1}^{N_F - 1} {(1 - M_i /M)} } \right] \hfill \\ \cdot \left( {\frac{{\bar \alpha _s (M)}}{\pi }} \right)^2 + 0(\bar \alpha _s^3 (M)), \hfill \\ \end{gathered} $$ as an accurate approximation forN _F?1 light quarks of massesM _i <M. Combining this new result with known three-loop results for \(\overline {MS} \) coupling constant and mass renormalization, we relate the pole mass to the \(\overline {MS} \) mass, \(\bar m\) (μ), renormalized at arbitrary μ. The dominant next-to-leading correction comes from the finite part of on-shell two-loop mass renormalization, evaluated using integration by parts and checked by gauge invariance and infrared finiteness. Numerical results are given for charm and bottom \(\overline {MS} \) masses at μ=1 GeV. The next-to-leading corrections are comparable to the leading corrections.     

Error Analysis on Photonic Qubit Rotations Implemented by Wave Plates

Optical Poincare sphere rotations $e^{-i\theta\sigma_{x}/2}$ , $e^{-i\theta\sigma_{y}/2}$ and $e^{-i\theta\sigma_{z}/2}$ can be realized by wave-plate combinations. Errors due to combinations with non-ideal wave plates are discussed for three specific combinations (θ=π) by trace distance. The result shows that different settings of combinations affect trace distance: (i) trace distance for $e^{-i\pi\sigma_{x}/2}$ equals that for $e^{-i\pi\sigma_{z}/2}$ , but both of them are smaller than that for $e^{-i\pi\sigma_{y}/2}$ , when optics-axis random errors are considered; (ii) trace distance for $e^{-i\pi\sigma_{x}/2}$ also equals that for $e^{-i\pi\sigma_{z}/2}$ , but both of them are larger than that for $e^{-i\pi\sigma_{y}/2}$ , when phase-shift random errors are considered. The method outlined in this paper is general and is useful to analyze other combinations.     

Non-Isothermal Boundary in the Boltzmann Theory and Fourier Law

In the study of the heat transfer in the Boltzmann theory, the basic problem is to construct solutions to the following steady problem: $$v \cdot \nabla _{x}F =\frac{1}{{\rm K}_{\rm n}}Q(F,F),\qquad (x,v)\in \Omega \times \mathbf{R}^{3}, \quad \quad (0.1) $$ v · ? x F = 1 K n Q ( F , F ) , ( x , v ) ∈ Ω × R 3 , ( 0.1 ) $$F(x,v)|_{n(x)\cdot v<0} = \mu _{\theta}\int_{n(x) \cdot v^{\prime}>0}F(x,v^{\prime})(n(x)\cdot v^{\prime})dv^{\prime},\quad x \in\partial \Omega,\quad \quad (0.2) $$ F ( x , v ) | n ( x ) · v < 0 = μ θ ∫ n ( x ) · v ′ > 0 F ( x , v ′ ) ( n ( x ) · v ′ ) d v ′ , x ∈ ? Ω , ( 0.2 ) where Ω is a bounded domain in ${\mathbf{R}^{d}, 1 \leq d \leq 3}$ R d , 1 ≤ d ≤ 3 , K_n is the Knudsen number and ${\mu _{\theta}=\frac{1}{2\pi \theta ^{2}(x)} {\rm exp} [-\frac{|v|^{2}}{2\theta (x)}]}$ μ θ = 1 2 π θ 2 ( x ) exp [ - | v | 2 2 θ ( x ) ] is a Maxwellian with non-constant(non-isothermal) wall temperature θ(x). Based on new constructive coercivity estimates for both steady and dynamic cases, for ${|\theta -\theta_{0}|\leq \delta \ll 1}$ | θ - θ 0 | ≤ δ ? 1 and any fixed value of K_n, we construct a unique non-negative solution F _s to (0.1) and (0.2), continuous away from the grazing set and exponentially asymptotically stable. This solution is a genuine non-equilibrium stationary solution differing from a local equilibrium Maxwellian. As an application of our results we establish the expansion ${F_s=\mu_{\theta_0}+\delta F_{1}+O(\delta ^{2})}$ F s = μ θ 0 + δ F 1 + O ( δ 2 ) and we prove that, if the Fourier law holds, the temperature contribution associated to F ₁ must be linear, in the slab geometry.