For separation of variables, put and set . The boundary conditions are , . Multiplying by and using integration by parts gives
A nonzero solution therefore has : the zero case forces constant, and then forces . Writing , the eigenfunctions and their temporal factors are
These sine functions are orthogonal on , with squared norm . Expanding the initial condition and multiplying by gives, using the permitted interchange of the sum and integral,
Here and . Thus the Fourier series solution is
The exponential decay permits differentiation of the series for , and each term satisfies both boundary conditions. At it is the mixed-end Fourier expansion of .
The gained heat is the increase in internal energy. Since the paper specifies no cross-sectional area, first give the answer per unit cross-sectional area; multiply by for a rod of cross-sectional area . The initial temperature has integral , whereas . At this gives
Numerically . This is positive and less than the eventual gain , as expected from relaxation towards the uniform steady state.
For a nonzero integer , direct integration gives
whereas . The product-to-sum formulas
therefore establish Fourier orthogonality: for positive integers ,
Here is the Kronecker delta. In particular the constant function is orthogonal to every nonconstant sine and cosine in this system. Integrating the given Fourier series, then multiplying it by or and integrating, isolates its Fourier coefficients. Using the permitted interchange of the sum and integral gives
If is an even function, then is an odd function, while is an even function. Thus
The Fourier series consequently reduces to a Fourier cosine series.
Pair adjacent terms in the finite sum:
There are pairs, each equal to . The result is .
For , use the finite geometric series with ratio :
The second equality follows by factoring in numerator and denominator. Taking real parts proves the finite trigonometric sum identity
When , every cosine is , so the sum is . This also gives the continuous limiting value at the excluded arguments.
Let . Its gradient cannot vanish on the ellipsoid, since the origin is not on the surface. Thus the Lagrange multiplier equations capture every constrained stationary point:
The temperature is a sum of squares, so . Equality requires and . The surface constraint then gives the four coldest points:
These are all such points, because a sum of nonnegative terms can vanish only when each term vanishes.
Group consecutive terms in pairs. The th pair, for , is
There are such pairs in this finite sum, so the answer is .
Use the unitary discrete Fourier transform
The transformed operators obey . The hopping terms contribute to the number coefficient. In the pairing sum, the odd sine part cancels between and . Thus
where
The factor one half in the full momentum sum avoids counting an anomalous pair twice. For , for every allowed momentum. Set
Real even angles ensure consistency between opposite momenta, and preserves the canonical commutation relations. The anomalous coefficient becomes . Normal ordering then gives
with
The field-polarized planar ferromagnetic spin wave has energy gap . The momenta run over the whole Brillouin zone once, including the self-inverse momenta and, for even , ; these use single-mode squeezing rather than introducing an extra independent partner.
The zero-field case has a zero-mode obstruction to bosonic Bogoliubov diagonalization. Its nonzero-momentum modes have , but , so no finite squeezing angle exists for the uniform mode. With , its Hamiltonian is , not an identically zero oscillator. For example, its matrix element between the vacuum and the two-boson number state is , whereas would vanish. Consequently the requested oscillator diagonalization is literal for nonzero field; at zero field treat the collective orientation separately or introduce a field regulator before taking a limit.
For nonzero field, the spin-wave ground state is the vacuum of all . The inverse Bogoliubov transformation gives and vanishing one-point oscillator averages. Translational invariance therefore gives the depletion
Rotating back to the physical axes yields, to the retained spin-wave order,
The zero transverse components also follow from rotation by about the physical axis in a symmetry-preserving ground state. The nonzero depletion describes zero-point quantum fluctuations: the original boson vacuum is not the ground state because the anomalous terms create pairs. It is nonnegative and must satisfy for the expansion to be reliable. At strong field it tends to zero, with in the infinite-chain limit.
For fixed , the momentum sum has a finite thermodynamic limit obtained by . At zero field the uniform mode already has no oscillator vacuum, as explained in part d. Even after removing that mode, the nonzero modes behave as
Thus assigning a finite ordered moment by simply setting in the infinite momentum sum is invalid. The divergence indicates failure of this expansion about a fixed planar orientation; it is not a negative physical magnetization. At finite size an exact zero-field state chosen to have definite total has zero planar one-point magnetization by the continuous rotation symmetry. Establishing its full correlations requires a treatment beyond the divergent ordered-state approximation.
With no pairing, the Hamiltonian is diagonal in the fermionic occupation numbers. The Pauli exclusion principle allows each momentum mode to be either empty or occupied once. Without a fixed fermion number constraint, occupying a mode lowers the energy precisely when its single-particle energy is negative. Therefore, with the Fock vacuum, a ground state is
Choose a fixed ordering of the fermionic creation operators; changing it alters only the overall phase. Any zero-energy modes may be occupied or empty without changing the ground-state energy.
To identify the occupied modes, complete the square:
Here is the dispersion parameter, not inverse temperature. If , every mode has positive energy and the unique ground state is the Fock vacuum, with . If , there are still no negative-energy modes, but have zero energy if they lie on the allowed momentum grid. Each such mode doubles the ground-state degeneracy.
For , define
Both are positive, and the occupied Fermi sea consists of the two intervals . With , the requested finite-size ground-state energy is explicitly
The factor two counts opposite momenta. Endpoints are excluded from this sum because they have zero energy; when an endpoint is allowed, its optional occupation adds ground-state degeneracy. A narrow negative-energy interval can contain no allowed finite-size momenta, in which case the sum is empty. The zero-momentum mode is always empty since .
Since every quasiparticle energy is positive, the ground state is the vacuum of all the annihilation operators, including . Its zero-momentum mode is therefore empty. For each positive momentum, the inverse Bogoliubov transformation gives
Both annihilate the normalized pair state . For example, applying produces the two cancelling terms . Thus
The factors have even fermion parity and commute for disjoint momentum pairs. Each normalized factor is a superposition of zero and two fermions, whose pair-occupation probability is . Therefore
There is no zero-mode contribution. For fixed , the sum converges since .
A normalized simple excited state is for . Its excitation energy is , and its total momentum is : both terms in add momentum . Acting on the pair gives
The affected pair now contains exactly one fermion instead of mean number ; all other pairs are unchanged. Consequently
This is the fermion number of a Bogoliubov excitation. Its change in mean fermion number lies between zero and one for , with the value one at . The Bogoliubov quasiparticle coherently mixes adding a fermion and removing the opposite-momentum fermion; it is not an eigen-excitation of the original fermion number operator. Indeed,
so fermion number is not conserved when , while fermion parity is conserved. The generally noninteger change is an expectation value in a coherent paired background, rather than a fractional eigenvalue of the fermion number operator. In the weak-pairing limit at fixed , it tends to one, as expected for an ordinary added fermion.
Evaluate the mixed radix representation from its most significant end, using Horner's method. There is no need to compute the positional weights separately.
def digsToInt(xs: Array[Int], bs: Array[Int]): Int = {
  require(xs.length == bs.length)
  val N = xs.length
  var i = N
  var acc = 0
  while (i > 0) {
    i -= 1
    acc = xs(i) + bs(i) * acc
  }
  acc
}
The digit conditions are additional preconditions. The loop invariant at the loop test is
An empty sum is zero and an empty product is one, so initialization at is correct. For the next digit, writing gives
which preserves the loop invariant. The loop variant is . At termination , giving
The code executes exactly multiplications, meeting the bound, and uses auxiliary space. When it returns zero. For the time representation in the source, it computes .
A zero-address instruction has no explicit operand addresses, so its operands and result location must be implicit in the instruction set architecture. A stack machine supplies these naturally: operations remove their inputs from the top of an operand stack and push the result. A fixed-processor register design is logically possible, so “almost certainly” is appropriately weaker than a proof that only a stack machine could work.
Write the stack from bottom to top, with its top on the right. Let stand for any unchanged prefix. Suitable effects are
Thus add pops two operands and pushes one sum, reducing stack occupancy by one. The special const c instruction pushes its immediate constant and increases occupancy by one. Only const needs an explicit value in its encoding.
Generate each reciprocal factorial from the preceding one, instead of recomputing a factorial for every summand:
def eb(n: Int): Double = {
  require(n >= 0)
  var i = 0
  var term = 1.0
  var sum = 1.0
  while (i < n) {
    i += 1
    term /= i.toDouble
    sum += term
  }
  sum
}
The loop invariant is , , and . Initially and . The division produces , which is then added to the previous partial sum. The loop variant is . Consequently the returned value is , using exactly divisions and additions, no multiplications, and time complexity. The slightly weaker precondition also permits the useful zero-term boundary case. For Double, small positive terms eventually cease to change sum because of floating-point rounding.
Let be the set of permutations fixing object . If a specified collection of objects is fixed, the remaining objects can be arranged in ways. There are such collections. Applying the inclusion-exclusion principle to exclude all fixed objects gives
For , the terms and cancel. Thus the number of derangements is , interpreting the remaining sum as empty when . The full inclusion-exclusion formula also gives .
By the preceding convolution, has a Poisson distribution with parameter . Directly from its probability mass function, shifting the index in the absolutely convergent sums gives and . Since , . This is also the identity .