bring back the "correlated nucleon tail"#421
Merged
sjgardiner merged 9 commits intomasterfrom Dec 20, 2024
Merged
Conversation
The correlated tail from previous work disappeared with the most recent modifications to the LFG model here. Forcing the computed removal energy to be nonnegative makes sense for momenta drawn from the bulk of the LFG distribution; negative removal energy ordinarily means the nucleon is not bound (i.e. outside the nucleus). However, for momenta drawn from the correlated tail, they will *always* be negative relative to the calculation from the LFG model: their associated KE is much larger because their *internal* (pairwise) potential energy offsets it. Obviously the LFG computation can't be aware of that. Here we take cases where the LFG calculation winds up negative and simply reassign them a (more) sensible removal energy. In principle we would like something specific to SRC correlated pairs. In the absence of that, we fall back on the mean binding energy per nucleon, which is a (bad) proxy for the one-nucleon separation energy.
… its strength at 20% Otherwise identical to AR23_20i
… the LocalFGM nuclear model to turn use of the new Maxwell-Boltzmann distribution on/off even when not forcing positive removal energies. This allows the default configuration in config/LocalFGM.xml to retain its previous behavior. Adjust the new N24_20i configuration to turn this option on. Include the new parameter in the AR23_20i configuration with the default value of false (no change in behavior expected). Remove default arguments to LocalFGM::MaxwellBoltzmannRemovalE() since this member function is currently only called in a single place with explicit arguments.
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
3 participants
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
The correlated tail from previous work disappeared with the most recent modifications to the LFG model here. Forcing the computed removal energy to be nonnegative makes sense for momenta drawn from the bulk of the LFG distribution; negative removal energy ordinarily means the nucleon is not bound (i.e. outside the nucleus). However, for momenta drawn from the correlated tail, they will always be negative relative to the calculation from the LFG model: their associated KE is much larger because their internal (pairwise) potential energy offsets it. Obviously the LFG computation can't be aware of that.
Here we take cases where the LFG calculation winds up negative and simply reassign them a (more) sensible removal energy. The removal energy is sampled from a Maxwell-Boltzmann distribution, as is done for the "correlated spectral function" in spectral function fits. Here we use rough parameters taken from estimates by Artur Ankowski given in private communication to J. Wolcott (Ankowski paper forthcoming). The values thrown by the MB distribution are capped such that the off-shell struck nucleon in the initial state has nonnegative mass, which otherwise causes NaNs when trying to boost into its rest frame.
Because the default behavior of the "AR23" CMC would otherwise be changing under the feet of folks who want to use it, this PR also includes some changes to the AR23 configuration so that it should preserve the previous behavior. Here we also include a "N24" (NOvA 2024) CMC that turns it on.