Pony maths universe thing

[atomic7732]

I'm going to simulate a 2-dimensional universe based on 6 elementary particles, each named after the elements of harmony.

Forces

The force of attraction between two particles (bonded or not) can be expressed as the equation:

 F = k*((p1*-p2)/r^2)

Where:

• k is the pony universe constant (0.74386705412065310000 e+idk yet)
• p1 is a particle's charge
• p2 is another particle's charge
• r is the distance between the particles

The force of bonding repulsion can be expressed:

F = 1/(10x)

Where x is the difference in charge. (Yeah, this isn't going to work, must figure out how to get them to stableize around .1 unit away)

The properties of each particle can be expressed by several values: charge, stability, (will come up with more if needed).

Charge is used to determine force of attraction.
Stability is used to determine how a particle decays.

Particle properties
Loyalty particle (l)
Charge: -1
Base stability: 0
Size: 1 pu (particle unit)

Kindness particle (k)
Charge: 2
Base stability: 0
Size: 3.542 pu

Generosity particle (g)
Charge: -3
Base stability: 0

Laughter particle (la)
Charge: 4
Base stability: -0.41843553120

Magic particle (m)
Charge: -7
Base stability: -1.28353298071

Honesty particle (h)
Charge: 7
Base stability: -5.5132

Extra rules
When two particles occupy the same space (within .2 units) and are still being attracted to each other are judged to have bonded. (This is for determining rarity and commonality of compound particles)

Neutral compounds will not be able to physically bond to other particles.

Compounds will be calculated as one particle (for overall movement) as well as separate particles to simulate reactions in case a particle does indeed break away.

Particle stability averages with all particles bonded with, but only unstable particles decay.

Charge is never created nor destroyed.

A particle's stability while in a compound is determined by taking the average of the particles' stability in the compound and then multiplying by 1/10 of the overall charge.

To determine what a particle decays into, take stability value: (c/|s|)/2 = charge of average particle, then use largest smaller particles for remaining charge.

Particle-lifetime is determined using the following formula: s*5 (units of time)

A decay product is propelled a random direction at a force of 2 units per iteration?

Absolute charge is a measure of "mass" where the absolute values of the charges of all particles are added together of a compound or system.

Naming compounds and formulae
Overall neutral compounds containing only 2 types of particles should be named in the following way:
adj.x n.y

Where x is the particle that occurs more often.

Still-charged (and extremely reactive) compounds containing only 2 types of particles should append the overall charge to the end and the number of particle y there is. This is so that one may destinguish kl from k2l3, like so:
adj.x n.y (+/-z,w)
loyal kindness (+1,1) vs loyal kindness (+1,2)

Formulae: the larger absolute charge particle appears first, then all others of the same type (positive or negative), then list the largest opposite charge, and then others.

Predicted compounds
loyal kindness (0), kl2, 4, -, l + k + k, ?

Observed compounds
Sim 1: Not yet run
name, formula, abs charge, lifetime (before reacted/decay/otherwise destroyed; stability), created how, bond radii

Observed reactions

Reaction tendencies

[vh]

can they be ponies instead?
kol.

also nids anti-elements of harmony and explosions

[atomic7732]

Nah, there won't be anti-matter.

Also, not sure if I should do 2^r or r^2. Probably the latter... I'll start with that.