A "surnatural numbers" as a largest model of the natural numbers One characteristic of the surreal numbers is that they are a monster model of the first-order theory of real numbers, according to Joel David Hamkins in this post. Thus they are real-closed, and every other real-closed field embeds into them with very nice properties. So we may ask if one can similarly create a monster model of the first-order theory of natural numbers as a kind of "surnaturals."
One natural question to ask is if the non-negative elements in Conway's ring of "omnific integers" would fit the bill, given that they are often promoted as the "surreal version of the integers" or something like that. However, it is rather easy to see that this is not true: the omnific integers have the interesting property that their field of fractions is the entire surreal number field, thus there are two omnific integers whose quotient is the square root of 2, whereas no such integers exist in any model of the naturals.
Thus in some sense, the omnific integers are not quite the most direct correspondant of the natural numbers within the surreal numbers. So one question is if the monster model can be built constructively, similarly to the surreals. I would suspect that such a construction exists, given that Joel David Hamkins was able to explicitly construct a monster model of all groups (!) in the above post, which would seem to suggest a monster model of the integers also exists, and thus the naturals as those non-negative integers.
How can one build such a model?
EDIT: I initially asked if the Grothendieck ring of the ordinals with commutative addition/multiplication could be a monster model for the integers, but it is apparently too small, since for $\omega$ we should have that $\omega$ is either even or odd. This means there should be some element $x$ such that $x + x = \omega$, or $x + x = \omega + 1$, so either $\omega/2$ or $(\omega + 1)/2$ should be in our set. Also, there will need to be some element in the monster model which is divisible by every standard finite number, and if that were $\omega$, we'd thus need to have $\omega \cdot q$ for every rational q. Thanks to Noah Schweber and Emil Jeřábek for pointing this out, and also to Emil for clearing up some confusion I had about whether the non-unique factorization of omnific integers necessarily implies they are not a model of the naturals (apparently it does not, but there are other reasons why, such as there are two omnifics whose ratio is sqrt(2), which is provably not true in any model of PA).
 A: I asked (and also answered) a more general version of this question a while ago. To summarize the answer, some results of Kanovei and Shelah have the following corollary:
Fact. In $\mathsf{ZFC}$ there is a uniform procedure for building 'set-saturated,' class-sized elementary extensions of arbitrary structures. That is to say there are formulas $S(M,L,x)$ and $F(M,L,f,x)$ in the language of set theory such that in any model $V \models \mathsf{ZFC}$ if $L \in V$ is a language and $M \in V$ is an $L$-structure, then the following hold (where $M^\ast = \{x \in V : V \models S(M,L,x)\}$):

*

*$M \subseteq M^\ast$,

*if $\varphi \in V$ is an $L$-formula with free variables $x_0,\dots,x_n$ and $\bar{a} \in M^\ast$ is an $n$-tuple, then $V \models F(M,L,\exists x_n\varphi,\bar{a})$ if and only if $V \models (\exists x \in M^\ast) F(M,L,\varphi,\bar{a}x)$ (where we are using some fixed coding of tuples in $\mathsf{ZFC}$),

*furthermore, if $\bar{c} \in M$ is an $(n+1)$-tuple, then $V \models F(M,L,\varphi,\bar{c})$ if and only if $V \models “M \models \varphi(\bar{c})”$ (in particular, if $\varphi$ is a sentence, then $V \models F(M,L,\varphi,\varnothing)$ if and only if $V \models “M \models \varphi”$),

*$F$ is compatible with Boolean combinations (i.e., $V\models F(M,L,\varphi\wedge \psi,\bar{a})$ if and only if $V\models F(M,L,\varphi,\bar{a})\wedge F(M,L,\psi,\bar{a})$, etc.), and

*if $A \subseteq M^\ast$ is a set and $p(x)$ is a finitely satisfiable set of $L_A$-formulas with free variable $x$, then there is $b \in M^\ast$ such that for any $\varphi(x,\bar{a}) \in p(x)$, $V \models F(M,L,\varphi,b\bar{a})$.

So to state it informally, $S(M,L,x)$ defines the universe of a class-sized elementary extension of $M$ and $F(M,L,f,x)$ is its truth predicate.
Applying this to the naturals tells us that there is a formula that defines a proper class monster model of $\mathrm{Th}(\mathbb{N})$ in any model of $\mathsf{ZFC}$.
One thing to note, though, is that without global choice (which makes my original question trivial), it's unclear whether there's always a definable isomorphism between different set-saturated class-sized models of a given theory. I believe this is related to an unanswered MathOverflow question of Hamkins. That said, if $M$ and $N$ are $L$-structures and $M \equiv N$, then there will be an isomorphism between $M^\ast$ and $N^\ast$ that is definable with certain parameters.
Another thing to note is that some constructions that model theorists commonly use with the monster model are unclear in the context of these class monster models. There isn't necessarily a good way to talk about arbitrary global types, for instance. You do, however, get a good homogeneity property: There is a subgroup $G$ of $\mathrm{Aut}(M^\ast)$ that can be represented as a class in a definable way which has the property that if $\bar{a}$ and $\bar{b}$ are set-sized tuples that realize the same type, then there is a $\sigma \in G$ such that $\sigma \bar{a} = \bar{b}$.
A: I think it should be said that the first order theory of the semi-ring $\mathbb{N}$ of non-negative integers is much more difficult to work with than that of the real ordered field.
The existence of certain proper elementary extensions of $\mathbb{N}$ usually relies on non-constructive methods, so it is unlikely that one could give a nice description of such an object within the class $\mathbf{No}$ of surreal numbers.
An easier task would be to try to identify models of Peano Arithmetic (PA) within $\mathbf{No}$ (the class of non-negative omnific integers only satisfies the fragment of PA called Open Induction, which only contains weak instances of the induction axioms). This seems already difficult enough, but for this one at least there is some literature. Specifically about integer parts of real-closed fields which are models of (the full integers version of) PA. I don't know this literature, but I know Salma Kuhlmann and Paola D'Aquino have look into related questions.
