$\set{(x,y) \in \Z\times\N \where (y = 0 \wedge x\ge 0) \vee y > 0}$.
We interpret $<$ so that a point $u<{}$ another point $v$ iff $u$ is
vertically lower than $v$, or it's at the same vertical level but to the
left of $v$.
We haven't yet fully specified our model. Consider that we can express
$x~is~even$ in the language of arithmetic ($\exists y. x = y+y$).
Similarly for $x~is~odd$. And it will be a truth of arithmetic that
every number has one of these properties. So if we want to preserve that
sentence, $i$ must also be come out either "even" or "odd." If it's
even, then we'll need to include some further elements in the domain,
because nothing so far is such that double it is $i$. We'll need to add
$\frac{i}2$, and then add every other integer point on the same
horizontal line ($\dots, \frac{i}2-1, \dots, \frac{i}2+1, \dots$) to
take care of successors and predecessors. If on the other hand $i$ is
odd, then $i+1$ will be even, and we'll need to add $\frac{i+1}2$, and
every other integer point on *that* horizontal line. Similar reasoning
shows that we have to add points on the horizontal lines $\frac{i}3$ and
$\frac{2i}3$ (or else on lines slightly offset from those), and so on.
In the end, we'll need to include all the integers on every horizontal
line where $y$ is a positive rational. That is, now we'll have all these
points in the complex plane:
> $\set{(x,y) \in \Z\times\pmb{\Q} \where (y = 0 \wedge x\ge 0) \vee y > 0}$.
Have we fully characterized multiplication? We've said what twice $i$
is, and half $i$, and so on. But what is $i*i$? We don't want to use the
ordinary operation of complex multiplication because in true arithmetic
$\forall z. z*z*z*z=1 \cond z=1$. So to characterize multiplication
where non-standard elements are multiplied by other non-standard
elements, as when we multiply $i$ by $i$, we have to do something more
complicated. Currently all our non-standard elements are in the part of
the complex plane that lies above the axis $y=0$. Let's vertically
"squash" this region of the complex plane so that it occupies instead
the just the open vertical band above $y=0$ and below $y=1$. We can do
this by mapping every $y$-value through a monotonically increasing
continuous function whose domain is the first interval and whose range
is the second interval. One such function is graphed
[here](https://www.wolframalpha.com/input/?i=y%3D%28arctan+%28ln+x%29%29%2Fpi%2B1%2F2%2Cx+from+0+to+6).
We have to adjust the interpretation of $+$ to compensate, and also the
partial interpretation of $*$ that we've already identified. Our first
non-standard element $i$ might now appear at $(0,\frac12)$ in the
complex plane; but I will continue to refer to it as $i$. Our problem
was that we hadn't specified the interpretation of operations like
$i*i$. We can now do that, by adding *another copy* of the newly
squashed region, on top of the first copy. This new region will lie
above $y=1$ and below $y=2$ (none of the points on the horizontal lines
where $y=1$ or $y=2$ will themselves be included). The result of
squaring our $i$ element from the original non-standard block will be
the element at the corresponding location in the second block. The
result of multiplying the $i$ element from the first block by the $i$
element from the second block will then have to be an $i$ element in yet
a third block, above the previous two, and so on. In the end, we'll have
all these points in the complex plane:
> $\set{(x,y) \in \Z\times\R \where (y=0 \wedge x\ge 0) \vee (y>0 \wedge y \nin \N \wedge \operatorname{squash}^{-1}(y \bmod 1)\in \Q)}$,
where $\operatorname{squash}$ is whatever we use for our vertical
"squashing" function. Formulated differently, that's the same as:
> $\set{(x,0) \where x\in\N} \union \set{(x,n + \operatorname{squash}(y)) \where x\in\Z \wedge n\in\N \wedge y\in \Q^+}$.
Let's give this model a name, call it $\script N''$. The size of this
model's domain, like those of $\script N$ and $\script N'$, is still
countable. But $\script N''$ won't be isomorphic to the other models.
Although $\script N''$ is not isomorphic to $\script N$, it does make
all the same sentences true (it's elementarily equivalent). That is, it
satisfies all the sentences of true arithmetic. When we have a model
that satisfies a theory but isn't isomorphic to that theory's intended
model, we call it a **non-standard model**. So $\script N''$ is a
non-standard model of true arithmetic. Furthermore, up to isomorphism
this is the *only countable* non-standard model. What that means is that
any countable models not isomorphic the intended model have to be
isomorphic to this one.
Warning: I am sure that the model I've described has the right general shape to be the countable non-standard model of arithmetic. (Specifically, it has the right "order type," and it will satisfy many obvious truths of arithmetic.) But I haven't fully specified the interpretation of $+$ and $*$, and there are known difficult issues about doing so. So I can't warrant that what I've described is adequate in all respects. However, it should be good enough for you to have a working idea of what this model is like. See comments at the end of the page.
There seems to be a clear structural difference between our first two models and this new model $\script N''$. In the first two, a chain of predecessors always "bottoms out" at $0$. Whereas in $\script N''$, $i$ has a predecessor who has a predecessor who $\dots$, but this sequence never "bottoms out". When we learn more about order relations, we'll see that the official way to say this is that in $\script N''$, the domain is not **well-ordered** by $<$, but in the other models it was. And as it turns out, whether or not a relation is a well-order is one of those things you can only express in second-order logic. So it's not surprising that our first-order theories cannot capture this difference between the models.
This is connected to a question Sam asked in class. He asked us to consider the set $\set{\neg\forall x\theta x,\theta(0),\theta(S0),\theta(SS0),\dots}$. (Sets of this sort are called $\omega$-inconsistent.) Observe that for any finite subset of this, you could find *some* $\theta$ where that subset is satisfied on the intended model of arithmetic.
But there's no *single* $\theta$ for which all the finite subsets would be true, on the intended model. Still, every finite subset would be satisfi*able*, for *some* model, which interpreted symbols differently than the intended model does.
And so compactness then tells us that *the whole set* must be satisfiable. But intuitively it might seem like it shouldn't be. Doesn't the first sentence contradict what all the rest of the sentences together say? Well, *only* if we've managed to somehow introduce the requirement that the only objects there be are $0, S0, SS0, \dots$. Let's call objects of this sort, that are reachable by applying the successor function a finite number of times to $0$, "finitely generated." If we haven't excluded the possibility of there being additional objects that *weren't* finitely generated, then a model could satisfy all the sentences $\theta(0), \theta(S0), \theta(SS0), \dots$ while giving $\theta$ an interpretation that *doesn't* apply to the additional objects.
And as we were saying, in first-order logic it's *not* possible to require that the only objects there be were finitely generated. So there are no axioms we could state in a first-order language that would rule out the truth of every such set of sentences. As mentioned above, though, on the *intended* model $\script N$ of arithmetic, no such set is true --- there is no $\theta$ where $\script N$ satisfies the whole such set.
The preceding remarks still leave it open that there might be some *other* difference that the theory of true arithmetic could catch: some sentence that one of the models would make true and the other make false. But I report that there isn't. These models satisfy all the same first-order sentences in the language of arithmetic.
Let's press on and extend our model still further. Let $\script N'''$ include all the same elements from $\script N''$, and *also* all the integer points on the horizontal line where $y=\operatorname{squash}(\pi)$, too, and so on for every real. Now we'll have all *these* points in the complex plane:
> $\set{(x,0) \where x\in\N} \union \set{(x,n + \operatorname{squash}(y)) \where x\in\Z \wedge n\in\N \wedge y\in \R^+}$.
which can be simplified to:
> $\set{(x,0) \where x\in\N} \union \set{(x,y) \in \Z\times\pmb{\R^+} \where y\nin\N}$.
If we keep in mind that the cardinality of the reals exceeds the cardinality of $\Q$, it should be clear that *this* model's domain is *larger* than that of $\script N$ and the other models so far considered. There is no bijection between their domains; so they cannot be isomorphic. Yet this model too will also satisfy true arithmetic.
So we had $\script N, \script N', \script N'', \script N'''$. The first two were obviously isomorphic, and the fourth has a bigger domain so is *obviously* not isomorphic to any of the others. The third was *non*-obviously not isomorphic to the first two. $\script N''$ is a (the) countable non-standard model of true arithmetic, and $\script N'''$ is one uncountable non-standard model.
One of the important metalogical results introduced this week, the Loewenheim-Skolem theorem (in its "upward" form) will tell us that we should *expect* to find non-standard models like $\script N'''$. When a theory has any infinite model, it will also have models with larger domains. First-order theories are "blind" to the differences between different infinite cardinalities.
To read more, have a look at [Wikipedia](http://en.wikipedia.org/wiki/Peano_axioms).
Warning: This page contains a mix of
* things whose proof is here given or sketched
* things that have been proved but are here just being reported, and
* mere speculation that I'm hoping might be pedagogically useful.
Let me be explicit about which is which. The text beginning "Here's how we can persuade ourselves..." sketches a proof that, because of compactness, there must be models of true arithmetic that have at least one additional element in the domain, coming after all the standard numbers. It is provable but not here proven that this model will be non-isomorphic to the intended model, and hence it will be non-standard. We go on to sketch some of the additional characteristics that such a model would have. We just take it for granted that there can be such a model with a countable domain. It is provable that there can be, and that it has the "order type" of my model $\script N''$. I do not fully specify the $+$ and $*$ functions for this model, and it's mere speculation that the parts of them that I do sketch are sustainable. It has been proven that up to isomorphism, there is only one countable non-standard models of arithmetic; and that its $+$ and $*$ functions are not "arithmetically definable", that is, expressible in the language of arithmetic. It is known, because of the upward Loewenheim-Skolem theorem, that there are also uncountable models of arithmetic (for every cardinality); whether any of them is a "filling in" of the countable model in the way my $\script N'''$ is a "filling in" of $\script N''$ is merely my speculation.
I hope that some of this might help you fix your ideas; but I hope also you won't take anything stated here for a known fact that isn't.