# Calculating the error in <x^2> from the error in <x> (Molecular Dynamics)

Hi, does anyone know of an easy way to calculate the error in <x^2> from the error in <x>? I am running a molecular dynamics simulation and trying to work out the error in the fluctuation of kinetic energy <dEk> = <3/2NT^2> - <3/2NT>^2 from the error in <T>.

BvU
Homework Helper
You say you have the error in <T> ? Is'n that the same as ##<\sigma^2> = <(T-<T>)^2> = <T^2> - (<T>)^2 \ \ ## ?

You say you have the error in <T> ? Is'n that the same as ##<\sigma^2> = <(T-<T>)^2> = <T^2> - (<T>)^2 \ \ ## ?
Yes, I have the standard error in the mean of T but I also need the error in the mean of T^2

BvU
Homework Helper
So you don't just want the spread <dEk> but also the error in this spread ? Doesn't that depend on the size of the system ?

 adding: can't you use something like $${\sigma_{T^2}\over T^2 } = 2 {\sigma_T \over T} \ \ \ \rm ?$$

So you don't just want the spread <dEk> but also the error in this spread ? Doesn't that depend on the size of the system ?

 adding: can't you use something like $${\sigma_{T^2}\over T^2 } = 2 {\sigma_T \over T} \ \ \ \rm ?$$
I think so, this is what I'm having difficulty with. There must be error in <dEk> because there is error in <T>?

 adding: can't you use something like $${\sigma_{T^2}\over T^2 } = 2 {\sigma_T \over T} \ \ \ \rm ?$$[/QUOTE]
I know of this formula but I wasn't sure if it would work with the error being in the average of T...

BvU
Homework Helper
I suppose you are doing something very sophisticated ?
The Maxwell Boltzmann distribution has well-defined characteristics with 'exact' ##\sigma##. It's only when you generate samples, that the error in such ##\sigma## (determined from the sample) comes into the picture ?

I suppose you are doing something very sophisticated ?
The Maxwell Boltzmann distribution has well-defined characteristics with 'exact' ##\sigma##. It's only when you generate samples, that the error in such ##\sigma## (determined from the sample) comes into the picture ?
Well I'm using a molecular dynamics program to determine a set of T values over time, so I get some fluctuation about the average T and I've worked out the error in the average of T. The problem is that I now need the error in <dEk> which I think needs the error in T but I'm not completely sure..

I think you want to use the propagation of errors formula.

https://en.m.wikipedia.org/wiki/Propagation_of_uncertainty
I've looked at this carefully but I am still unsure and can't seem to find any information of whether these formulas work for errors in the mean of quantities. For example, I'm not sure that

d<T^2>/<T^2> = 2d<T>/<T>

Dale
Mentor
2021 Award
Hmm, I can't answer for certain, but when I have been in similar situations I would just generate 10000 normally distributed numbers and see if the formula works.

It probably depends on the mean, so you might try mean 100, variance 9 or something similar.

Kieran
Hmm, I can't answer for certain, but when I have been in similar situations I would just generate 10000 normally distributed numbers and see if the formula works.

It probably depends on the mean, so you might try mean 100, variance 9 or something similar.

Dale
Mentor
2021 Award
You are welcome. Sometimes a few simulations are nearly as good as a proof.

BvU
Homework Helper
Well I'm using a molecular dynamics program to determine a set of T values over time, so I get some fluctuation about the average T and I've worked out the error in the average of T. The problem is that I now need the error in <dEk> which I think needs the error in T but I'm not completely sure..
Looks a bit like sample mean sigma ##\sigma_m## is ##\displaystyle \sigma_{\rm\ population}\over \sqrt N## and since you only have the sample it's ##\displaystyle \sigma_m## is ##\displaystyle\sigma_{\rm\ sample}\over \sqrt {n-1}##.

Fill us in on the result of Dale's suggestion ! My bet is on the simplest guess: double the relative error in T to get the relative error in T2

FactChecker
I think using (x+ε)2 = x2+2xε+ε2 would give the answer that the error would be 2$\bar{x}$ε + ε2