Elegant users forum

Chao, The "nonlinear" method in elegant is just Bassetti-Erskine, which I think is what you are using as well. The linear method is just the linear term of B-E. You can use longitudinal slicing in elegant for comparison with your code, since it does that as well. The particle_tune command uses NAFF ...

Chao, This plot compares results with and without space charge, using 1ps longitudinal slices. Seems that space charge has a significant impact on the vertical plane. Since the runs don't take much time, I increased the number of particles to 10 million, but it's excessive. comparison.png --Michael

Chao,

Here's what I get. No sign of the "noisy" behavior. I think this was a result of the above-mentioned bug in the older versions of the code.
--Michael

Chao,

I'm running some tests now. I wonder if you using the latest version of elegant. In version 2025.2 we fixed several bugs related to space charge. In addition, we added longitudinal slicing.

--Michael

Samira, I think the issue with the 3Hz ramp is that the quantum excitation rate is just not high enough to get the emittance and energy spread up to the equilibrium values in the time available. The beam goes through minimum emittance at t~=0.1s and then starts rising fairly slowly. Just as a check,...

Chao, The problem is that the beam is not matched to the lattice at the injection point. You can fix this with sddsmatchtwiss train.sdds rematched.sdds -xplane=filename=lattice.twi -yplane=filename=lattice.twi mv train.sdds old-train.sdds mv rematched.sdds train.sdds where "lattice.twi" is a file wi...

Chao,

I see that the emittance initially increases, then goes down. Can you check to see if particles are lost from the simulation?

--Michael

Esraa,

You can add a chirp (R65 nonzero) and compression (R56 nonzero) using EMATRIX elements. Please see attached files for an example

When optimizing both the chirp and R56, you'll want to place reasonable limits on the maximum allowable energy spread and R56.

--Michael

Samira,

It might be because your path does not include your present working directory. You can either add it to your path, or modify setUpCases as to change "setUpCase" to "./setUpCase".

--Michael

Samira,

Looking at your files, I don't see how the SREFFECTS element gets set up.

--Michael