I will propose an answer my question. Please let me know if I have spoken accurately.
Within the hyperreal field (that is, the most traditional hyperreal field as described by some of the original texts on the subject by Robinson, Goldblatt, Keisler, and Henle), can we express ln(ϵ) and ln(ω) as a power series in terms of ϵ and/or ω? (That is, where ϵ is infinitesimal, ω is infinite, and ϵ = 1/ω.)
Yes:
ln(ω) = (1-ϵ) + ((1-ϵ)^2)/2 + ((1-ϵ)^3)/3 + ((1-ϵ)^4)/4 + ....
ln(ϵ) = -(1-ϵ) + -((1-ϵ)^2)/2 + -((1-ϵ)^3)/3 + -((1-ϵ)^4)/4 + ....
These results are indeed related by the log identity ln(ω) = -ln(ϵ), but can both be derived independently from the Taylor series for natural logarithm
ln(1 + x) = x - (x^2)/2 + (x^3)/4 - (x^4)/4 + ...
whose radius of convergence is |x| < 1 and x = 1.
PROOF:
First, find ln(ϵ).
Since the infinitesimal ϵ is in the radius of convergence, we can directly let x = -1 + ϵ:
ln(ϵ) = ln(1 + (-1+ϵ)) = (-1+ϵ) - ((-1+ϵ)^2)/2 + ((-1+ϵ)^3)/3 - ((-1+ϵ)^4)/4 + ...
We can rewrite this expression more neatly as:
ln(ϵ) = -(1-ϵ) - ((1-ϵ)^2)/2 - ((1-ϵ)^3)/3 - ((1-ϵ)^4)/4 - ...
After expanding the right side and gathering all the standard terms, we see an expression that clearly diverges to negative infinity as expected:
ln(ϵ) = -1 - 1/2 - 1/3 - 1/4 - .... + (all of the non-standard terms)
Second, find ln(ω).
To do this, I will derive a transformed Taylor series for ln(m) whose radius of convergence is any real number m > 1.
Assuming m > 1,
we have 1 > 1/m > 0,
which can be rewritten as 0 > -1 + 1/m > -1.
Therefore, (-1 + 1/m) is within the radius of convergence of the Taylor series used above.
Thus, we can rewrite ln(m) as:
ln(m) = ln((1/m)^-1) = -ln(1/m) = -ln(1 + (-1 + 1/m))
and can now use the given Taylor series using x = -1 + 1/m:
ln(m) = -ln( 1 + (-1 + 1/m)) = (1-1/m) + ((1-1/m)^2)/2 + ((1-1/m)^3)/3 + ((1-1/m)^4)/4 + ....
This is a convenient result in and of itself because radius of convergence of this Taylor series for ln(m) is any real m > 1.
By transfer principle, since ω > 1, we can use ln(m) to find:
ln(ω) = (1-1/ω) + ((1-1/ω)^2)/2 + ((1-1/ω)^3)/3 + ((1-1/ω)^4)/4 + ....
and since 1/ω = ϵ, we have:
ln(ω) = (1-ϵ) + ((1-ϵ)^2)/2 + ((1-ϵ)^3)/3 + ((1-ϵ)^4)/4 + ....
which is indeed the opposite sign of the result for ln(ϵ):
ln(ϵ) = -(1-ϵ) + -((1-ϵ)^2)/2 + -((1-ϵ)^3)/3 + -((1-ϵ)^4)/4 + ....
In both cases, after expanding the terms and gathering all the standard terms, we see expressions that clearly diverge to infinity as expected:
ln(ϵ) = -1 - 1/2 - 1/3 - 1/4 - .... + (all of the non-standard terms in terms of ϵ)
ln(ω) = 1 + 1/2 + 1/3 + 1/4 + .... + (all of the non-standard terms in terms of ϵ)